From 8e92243c9dbbbdaebcb734e4174f44597834d69f Mon Sep 17 00:00:00 2001
From: Christophe Geuzaine <cgeuzaine@ulg.ac.be>
Date: Sat, 5 May 2007 11:36:32 +0000
Subject: [PATCH] *** empty log message ***
---
Plugin/SphericalRaise.cpp | 11 ++++-------
1 file changed, 4 insertions(+), 7 deletions(-)
diff --git a/Plugin/SphericalRaise.cpp b/Plugin/SphericalRaise.cpp
index add290864e..1bed1d4f3f 100644
--- a/Plugin/SphericalRaise.cpp
+++ b/Plugin/SphericalRaise.cpp
@@ -1,4 +1,4 @@
-// $Id: SphericalRaise.cpp,v 1.25 2007-05-05 10:24:53 geuzaine Exp $
+// $Id: SphericalRaise.cpp,v 1.26 2007-05-05 11:36:32 geuzaine Exp $
//
// Copyright (C) 1997-2007 C. Geuzaine, J.-F. Remacle
//
@@ -105,7 +105,7 @@ static void sphericalRaiseList(Post_View * v, List_T * list, int nbElm,
if(!nbElm)
return;
- if(timeStep > v->NbTimeStep - 1){
+ if(timeStep < 0 || timeStep > v->NbTimeStep - 1){
Msg(GERROR, "Invalid TimeStep (%d) in View[%d]", timeStep, v->Index);
return;
}
@@ -117,7 +117,7 @@ static void sphericalRaiseList(Post_View * v, List_T * list, int nbElm,
// compute d=(x-Xc,y-Yc,z-Zc)
// norm d
// get nodal value val at xyz
- // compute (x,y,z)_new = (x,y,z)_old + raise*val*(dx,dy,dz)
+ // compute (x,y,z)_new = (x,y,z)_old + (offset+raise*val)*(dx,dy,dz)
nb = List_Nbr(list) / nbElm;
for(i = 0; i < List_Nbr(list); i += nb) {
@@ -130,10 +130,7 @@ static void sphericalRaiseList(Post_View * v, List_T * list, int nbElm,
d[1] = y[j] - center[1];
d[2] = z[j] - center[2];
norme(d);
- if(timeStep < 0)
- coef = raise + offset;
- else
- coef = raise * val[nbNod * timeStep + j] + offset;
+ coef = offset + raise * val[nbNod * timeStep + j];
x[j] += coef * d[0];
y[j] += coef * d[1];
z[j] += coef * d[2];
--
GitLab