diff --git a/dune/microstructure/matrix_operations.hh b/dune/microstructure/matrix_operations.hh index 98c1650242b85d062ab1325dac654065545728c8..35b8083331acfc3de7672425458d48e76f04b8ef 100644 --- a/dune/microstructure/matrix_operations.hh +++ b/dune/microstructure/matrix_operations.hh @@ -67,8 +67,15 @@ namespace MatrixOperations { return rankoneTensorproduct(U_plus_k_x_ecs, e_1); } - static MatrixRT crossSectionDirectionScaling(double w, MatrixRT M){ - return {M[0], M[1], w*M[2]}; +// static MatrixRT crossSectionDirectionScaling(double w, MatrixRT M){ +// return {M[0], M[1], w*M[2]}; +// } + + static MatrixRT crossSectionDirectionScaling(double w, MatrixRT M){ + return {{M[0][0], M[0][1], w*M[0][2]}, + {M[1][0], M[1][1], w*M[1][2]}, + {M[2][0], M[2][1], w*M[2][2]} + }; } static double trace (MatrixRT M){