[Getfem-commits] (no subject)

[Markus Bürg]

branch: mb-transInversion
commit 27e02f9ba7fd21198e374f00e0625d1eef318ca0
Author: mb <address@hidden>
Date:   Wed Aug 23 13:36:04 2017 +0200

    Return linear flag correctly.
---
 src/bgeot_geometric_trans.cc | 9 ++++++++-
 1 file changed, 8 insertions(+), 1 deletion(-)

diff --git a/src/bgeot_geometric_trans.cc b/src/bgeot_geometric_trans.cc
index 83e5b0d..954b6b2 100644
--- a/src/bgeot_geometric_trans.cc
+++ b/src/bgeot_geometric_trans.cc
@@ -649,10 +649,17 @@ namespace bgeot {
   /* direct product transformation                                        */
   /* ******************************************************************** */
 
+  size_type get_degree(const poly_geometric_trans &poly_trans) {
+    auto degree = 1;
+    for (auto &&poly : poly_trans.trans) degree *= poly.real_degree();
+
+    return degree;
+  }
+
   struct cv_pr_t_ : public poly_geometric_trans {
     cv_pr_t_(const poly_geometric_trans *a, const poly_geometric_trans *b) {
       cvr = convex_ref_product(a->convex_ref(), b->convex_ref());
-      is_lin = false;
+      is_lin = get_degree(*a) * get_degree(*b) == 1 ? true : false;
       complexity_ = a->complexity() * b->complexity();
 
       size_type n1 = a->nb_points(), n2 = b->nb_points();