[Getfem-commits] (no subject)

[Liang Jin Lim]

branch: thread_safe_fem_computation
commit a194c885d058f9663be0ff0383238da988624576
Author: lj <address@hidden>
Date:   Wed Aug 2 11:31:16 2017 +0200

    Hierarchical fem should recompute again grad and hessian, instead of taking 
the copy from the initial fem.
---
 src/getfem_fem.cc | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/src/getfem_fem.cc b/src/getfem_fem.cc
index 43b6440..d4aff1e 100644
--- a/src/getfem_fem.cc
+++ b/src/getfem_fem.cc
@@ -828,6 +828,8 @@ namespace getfem {
 
   thierach_femi::thierach_femi(ppolyfem fi1, ppolyfem fi2)
     : fem<base_poly>(*fi1) {
+    grad_computed_ = false;
+    hess_computed_ = false;
     GMM_ASSERT1(fi2->target_dim()==fi1->target_dim(), "dimensions mismatch.");
     GMM_ASSERT1(fi2->basic_structure(0) == fi1->basic_structure(0),
                 "Incompatible elements.");