const Tensor<2, dim, ADNumberType> C5 = A * a;
const Tensor<2, dim, ADNumberType> C6 = A / a;
- const ADNumberType det_A = determinant(A);
- const ADNumberType tr_A = trace(A);
- const Tensor<2, dim, ADNumberType> A_inv = invert(A);
- const Tensor<2, dim, ADNumberType> A_T = transpose(A);
- const Tensor<2, dim, ADNumberType> A_adj = adjugate(A);
- const Tensor<2, dim, ADNumberType> A_cof = cofactor(A);
- const Tensor<2, dim, ADNumberType> A_orth =
- project_onto_orthogonal_tensors(A);
- const ADNumberType A_l1_norm = l1_norm(A, 1e-8);
- const ADNumberType A_linf_norm = linfty_norm(A);
+ const ADNumberType det_A = determinant(A);
+ const ADNumberType tr_A = trace(A);
+ const Tensor<2, dim, ADNumberType> A_inv = invert(A);
+ const Tensor<2, dim, ADNumberType> A_T = transpose(A);
+ const Tensor<2, dim, ADNumberType> A_adj = adjugate(A);
+ const Tensor<2, dim, ADNumberType> A_cof = cofactor(A);
+ const ADNumberType A_l1_norm = l1_norm(A);
+ const ADNumberType A_linf_norm = linfty_norm(A);
const ADNumberType A_ddot_B = double_contract<0, 0, 1, 1>(A, B);
const Tensor<2, dim, ADNumberType> A_dot_B = contract<1, 0>(A, B);
DEAL_II_CONSTEXPR const auto dummy_8 = cofactor(a);
deallog << "Deteriminant before orthogonalization: " << determinant(a)
<< std::endl;
- const auto dummy_9 = project_onto_orthogonal_tensors(a, 1e-8);
+ const auto dummy_9 = project_onto_orthogonal_tensors(a, 0);
deallog << "Deteriminant after orthogonalization: " << determinant(dummy_9)
<< std::endl;
Assert(determinant(dummy_9) - 1. < 1e-8, ExcInternalError());