From 3a2ff94cf2cae08f9daecffc731cb7389cb79e90 Mon Sep 17 00:00:00 2001
From: Jean-Francois Remacle <jean-francois.remacle@uclouvain.be>
Date: Thu, 9 Aug 2001 22:05:51 +0000
Subject: [PATCH] *** empty log message ***
---
Plugin/LevelsetPlugin.cpp | 20 +++++++++++---------
1 file changed, 11 insertions(+), 9 deletions(-)
diff --git a/Plugin/LevelsetPlugin.cpp b/Plugin/LevelsetPlugin.cpp
index 15ae87dd31..5d50855246 100644
--- a/Plugin/LevelsetPlugin.cpp
+++ b/Plugin/LevelsetPlugin.cpp
@@ -1,4 +1,4 @@
-// $Id: LevelsetPlugin.cpp,v 1.14 2001-08-09 21:15:30 remacle Exp $
+// $Id: LevelsetPlugin.cpp,v 1.15 2001-08-09 22:05:51 remacle Exp $
#include "LevelsetPlugin.h"
#include "List.h"
@@ -40,7 +40,7 @@ Post_View *GMSH_LevelsetPlugin::execute (Post_View *v)
a cut of the actual view with a levelset.
*/
int k,i,nb,edtet[6][2] = {{0,1},{0,2},{0,3},{1,2},{1,3},{2,3}};
- double Xpi[6],Ypi[6],Zpi[6];
+ double Xpi[6],Ypi[6],Zpi[6],myValsi[6];
Post_View *View;
// for all scalar simplices
@@ -69,7 +69,7 @@ Post_View *GMSH_LevelsetPlugin::execute (Post_View *v)
nx++;
}
}
- /*
+
if(nx == 4)
{
double xx = Xp[3];
@@ -85,7 +85,7 @@ Post_View *GMSH_LevelsetPlugin::execute (Post_View *v)
Zp[2] = zz;
myVals[2] = vv;
}
- */
+
double v1[3] = {Xp[2]-Xp[0],Yp[2]-Yp[0],Zp[2]-Zp[0]};
double v2[3] = {Xp[1]-Xp[0],Yp[1]-Yp[0],Zp[1]-Zp[0]};
double gr[3];
@@ -95,16 +95,18 @@ Post_View *GMSH_LevelsetPlugin::execute (Post_View *v)
gradSimplex(X,Y,Z,VAL,gr);
prosca(gr,n,&xx);
- if(xx > 0){
+ if(n[0] + n[1] + n[2] > 0){
for(k=0;k<nx;k++){
Xpi[k] = Xp[k];
Ypi[k] = Yp[k];
Zpi[k] = Zp[k];
+ myValsi[k] = myVals[k];
}
for(k=0;k<nx;k++){
Xp[k] = Xpi[nx-k-1];
Yp[k] = Ypi[nx-k-1];
Zp[k] = Zpi[nx-k-1];
+ myVals[k] = myValsi[nx-k-1];
}
}
@@ -118,10 +120,10 @@ Post_View *GMSH_LevelsetPlugin::execute (Post_View *v)
}
if(nx == 4)
{
- for(k=1;k<4;k++)List_Add(View->ST, &Xp[k %4]);
- for(k=1;k<4;k++)List_Add(View->ST, &Yp[k % 4]);
- for(k=1;k<4;k++)List_Add(View->ST, &Zp[k % 4]);
- for(k=1;k<4;k++)List_Add(View->ST, &myVals[k %4]);
+ for(k=2;k<5;k++)List_Add(View->ST, &Xp[k %4]);
+ for(k=2;k<5;k++)List_Add(View->ST, &Yp[k % 4]);
+ for(k=2;k<5;k++)List_Add(View->ST, &Zp[k % 4]);
+ for(k=2;k<5;k++)List_Add(View->ST, &myVals[k %4]);
View->NbST++;
}
}
--
GitLab