From 5d525152faf2119aebf6b99b94b0ba4bf3318774 Mon Sep 17 00:00:00 2001
From: Matti Pellika <matti.pellikka@tut.fi>
Date: Wed, 3 Mar 2010 06:44:53 +0000
Subject: [PATCH] Fixed orientation issue with overlapping chains.
---
Geo/ChainComplex.cpp | 25 +++++++++++++++----------
1 file changed, 15 insertions(+), 10 deletions(-)
diff --git a/Geo/ChainComplex.cpp b/Geo/ChainComplex.cpp
index da486aec68..5dedbf9cad 100644
--- a/Geo/ChainComplex.cpp
+++ b/Geo/ChainComplex.cpp
@@ -648,22 +648,27 @@ int Chain::writeChainMSH(const std::string &name)
void Chain::createPView()
{
std::vector<MElement*> elements;
- std::map<int, std::vector<double> > data;
-
+ std::map<int, std::vector<double> > data;
+ MElementFactory factory;
+
for(citer cit = _cells.begin(); cit != _cells.end(); cit++){
Cell* cell = (*cit).first;
int coeff = (*cit).second;
- MElement *e = cell->getImageMElement();
- cell->setDeleteImage(false);
-
- if(cell->getDim() > 0 && coeff < 0) e->revert(); // flip orientation
- for(int i = 0; i < abs(coeff); i++) elements.push_back(e);
-
+ MElement* e = cell->getImageMElement();
+ std::vector<MVertex*> v;
+ for(int i = 0; i < e->getNumVertices(); i++){
+ v.push_back(e->getVertex(i));
+ }
+ MElement* ne = factory.create(e->getTypeForMSH(), v);
+
+ if(cell->getDim() > 0 && coeff < 0) ne->revert(); // flip orientation
+ for(int i = 0; i < abs(coeff); i++) elements.push_back(ne);
+
std::vector<double> coeffs (1,abs(coeff));
- data[e->getNum()] = coeffs;
+ data[ne->getNum()] = coeffs;
}
-
+
int max[4];
for(int i = 0; i < 4; i++) max[i] = _model->getMaxElementaryNumber(i);
int entityNum = *std::max_element(max,max+4) + 1;
--
GitLab