[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);