Skip to content
Snippets Groups Projects

More changes in use compile time matrix sizes in LIE

@@ -193,11 +193,12 @@ void HydroMechanicsLocalAssemblerFracture<ShapeFunctionDisplacement,
pressure_size);
// Projection of the body force vector at the element.
// GlobalDim-1 is the element's dimension.
Eigen::MatrixXd const global2local_rotation =
R.topLeftCorner(_element.getDimension(), GlobalDim);
R.template topLeftCorner<GlobalDim - 1, GlobalDim>();
auto const& gravity_vec =
(global2local_rotation * _process_data.specific_body_force).eval();
Eigen::Matrix<double, ShapeFunctionPressure::DIM, 1> const gravity_vec =
(global2local_rotation * _process_data.specific_body_force);
ParameterLib::SpatialPosition x_position;
x_position.setElementID(_element.getID());
Loading