search

mikaem / Oasis

55 stars 48 forks source link

Computational cost Wale LES model #16

Open sjengq opened 6 years ago

sjengq commented 6 years ago

I notice a significant increase in computational time when enabling the Wale LES model (i.e. time steps take a factor 5 longer compared to no LES model). In addition I get the following notification when enabling the Wale model: WARNING: The number of integration points for each cell will be: 512 Consider using the option 'quadrature_degree' to reduce the number of points

Is this the anticipated behavior, or might I be doing something wrong?

Kind regards, Sjeng

mikaem commented 6 years ago

This seems to be a new problem with the model. Do you know which line triggers the warning? We should, as the message suggests, use 'quadrature_degree' to reduce the number of points.

sjengq commented 6 years ago

I use a slightly older version of OASIS, so the line numbers might differ slightly. The problem is triggered in Wale.py, line 38: nut_(). The warning is then mentioned in ffc/analysis.py(365)_check_quadrature_degree() num_points = ((degree + 1 + 1) // 2)**top_dim=512 (degree =14, top_dim=3).