This is an automated email from the git hooks/post-receive script. logari81 pushed a commit to branch master in repository getfem.
The following commit(s) were added to refs/heads/master by this push: new 3507069 Fix bug in second derivative of inverse operator 3507069 is described below commit 3507069aa22f5d232bd3cf5b3e3d45c500e5c5ba Author: Konstantinos Poulios <logar...@gmail.com> AuthorDate: Sun Oct 4 01:00:02 2020 +0200 Fix bug in second derivative of inverse operator --- src/getfem_generic_assembly_functions_and_operators.cc | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/getfem_generic_assembly_functions_and_operators.cc b/src/getfem_generic_assembly_functions_and_operators.cc index 6081028..1c889ac 100644 --- a/src/getfem_generic_assembly_functions_and_operators.cc +++ b/src/getfem_generic_assembly_functions_and_operators.cc @@ -331,7 +331,7 @@ namespace getfem { gmm::copy(__mat_aux1().as_vector(), result.as_vector()); } - // Derivative : -M^{-1}{ik}M^{-1}{lj} (comes from H -> M^{-1}HM^{-1}) + // Derivative : -M^{-1}{ik}M^{-1}{lj} (comes from H -> -M^{-1}HM^{-1}) void derivative(const arg_list &args, size_type, base_tensor &result) const { // to be verified size_type N = args[0]->sizes()[0]; @@ -374,7 +374,7 @@ namespace getfem { for (size_type j = 0; j < N; ++j) for (size_type i = 0; i < N; ++i, ++it) *it = __mat_aux1()(i,k)*__mat_aux1()(l,m)*__mat_aux1()(n,j) - + __mat_aux1()(i,m)*__mat_aux1()(m,k)*__mat_aux1()(l,j); + + __mat_aux1()(i,m)*__mat_aux1()(n,k)*__mat_aux1()(l,j); GA_DEBUG_ASSERT(it == result.end(), "Internal error"); } };