[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Getfem-commits] [getfem-commits] branch master updated: Fix bug in seco
From: |
Konstantinos Poulios |
Subject: |
[Getfem-commits] [getfem-commits] branch master updated: Fix bug in second derivative of inverse operator |
Date: |
Sat, 03 Oct 2020 19:00:11 -0400 |
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 <logari81@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");
}
};
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Getfem-commits] [getfem-commits] branch master updated: Fix bug in second derivative of inverse operator,
Konstantinos Poulios <=