[Getfem-commits] r5361 - /trunk/getfem/src/getfem_generic_assembly.cc

[logari81]

Author: logari81
Date: Mon Aug  1 11:41:11 2016
New Revision: 5361

URL: http://svn.gna.org/viewcvs/getfem?rev=5361&view=rev
Log:
small fix in ga_instruction_hess for 1D problems

Modified:
    trunk/getfem/src/getfem_generic_assembly.cc

Modified: trunk/getfem/src/getfem_generic_assembly.cc
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/src/getfem_generic_assembly.cc?rev=5361&r1=5360&r2=5361&view=diff
==============================================================================
--- trunk/getfem/src/getfem_generic_assembly.cc (original)
+++ trunk/getfem/src/getfem_generic_assembly.cc Mon Aug  1 11:41:11 2016
@@ -3193,7 +3193,8 @@
       size_type target_dim = Z.sizes()[1];
       size_type N = size_type(round(sqrt(scalar_type(Z.sizes()[2]))));
       size_type Qmult = qdim / target_dim;
-      GA_DEBUG_ASSERT((qdim == 1 && t.sizes()[0] == N && t.sizes()[1] == N) ||
+      GA_DEBUG_ASSERT((qdim == 1 && N == 1 && t.sizes()[0] == 1) ||
+                      (qdim == 1 && t.sizes()[0] == N && t.sizes()[1] == N) ||
                       (t.sizes()[1] == N && t.sizes()[2] == N
                        && t.sizes()[0] == qdim), "dimensions mismatch");
       GA_DEBUG_ASSERT(gmm::vect_size(coeff) == ndof*Qmult,