getfem-commits
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Getfem-commits] (no subject)


From: Markus Bürg
Subject: [Getfem-commits] (no subject)
Date: Mon, 28 Aug 2017 07:29:52 -0400 (EDT)

branch: mb-transInversion
commit 0346c9e8236ef967537ea9d19b255306a240ae60
Author: mb <address@hidden>
Date:   Mon Aug 28 13:29:47 2017 +0200

    Check complexity_ instead of degree.
---
 src/bgeot_geometric_trans.cc | 9 +--------
 1 file changed, 1 insertion(+), 8 deletions(-)

diff --git a/src/bgeot_geometric_trans.cc b/src/bgeot_geometric_trans.cc
index 18e4008..fb2655b 100644
--- a/src/bgeot_geometric_trans.cc
+++ b/src/bgeot_geometric_trans.cc
@@ -649,18 +649,11 @@ 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 = get_degree(*a) * get_degree(*b) == 1 ? true : false;
       complexity_ = a->complexity() * b->complexity();
+      is_lin = complexity_ == 1 ? true : false;
 
       size_type n1 = a->nb_points(), n2 = b->nb_points();
       trans.resize(n1 * n2);



reply via email to

[Prev in Thread] Current Thread [Next in Thread]