[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Getfem-commits] (no subject)
|
From: |
Yves Renard |
|
Subject: |
[Getfem-commits] (no subject) |
|
Date: |
Fri, 21 Jun 2019 09:24:46 -0400 (EDT) |
branch: master
commit f9c81ff11675b6aa87bcd9ce4d0b48f041d35b49
Author: Yves Renard <address@hidden>
Date: Fri Jun 21 14:42:53 2019 +0200
minor fix on extend_vector
---
src/getfem/getfem_mesh_fem.h | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/src/getfem/getfem_mesh_fem.h b/src/getfem/getfem_mesh_fem.h
index 49e4bc2..dbcfac2 100644
--- a/src/getfem/getfem_mesh_fem.h
+++ b/src/getfem/getfem_mesh_fem.h
@@ -257,8 +257,9 @@ namespace getfem {
template<typename VECT1, typename VECT2>
void extend_vector(const VECT1 &V1, const VECT2 &V2) const {
- if (is_reduced()) {
- size_type qqdim = gmm::vect_size(V1) / nb_dof();
+ size_type nbd = nb_dof();
+ if (is_reduced() && nbd) {
+ size_type qqdim = gmm::vect_size(V1) / nbd;
if (qqdim == 1)
gmm::mult(extension_matrix(), V1, const_cast<VECT2 &>(V2));
else