Currently when using a polynomial description of the material properties, we assume that we have a polynomial of order four. This is often not the case and we end up wasting a lot of flops. This PR makes the polynomial order a template parameter. The correct order is computed at runtime by doing a first read of the material parameters before calling MaterialProperties. The new template parameter p_order propagates everywhere because we need to template MaterialProperties and almost every class owns a reference to MaterialProperties.
Rombur
commented
4 months ago
Currently when using a polynomial description of the material properties, we assume that we have a polynomial of order four. This is often not the case and we end up wasting a lot of flops. This PR makes the polynomial order a template parameter. The correct order is computed at runtime by doing a first read of the material parameters before calling
MaterialProperties. The new template parameterp_orderpropagates everywhere because we need to templateMaterialPropertiesand almost every class owns a reference toMaterialProperties.