From 16b5f8e94a1f8f1cee6bcddf89b33e54253bb2f2 Mon Sep 17 00:00:00 2001
From: Thomas Toulorge <thomas.toulorge@mines-paristech.fr>
Date: Tue, 12 Jun 2012 10:01:18 +0000
Subject: [PATCH] Fixed ellipse/ellipsoid tensor visualization
---
Post/PViewVertexArrays.cpp | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/Post/PViewVertexArrays.cpp b/Post/PViewVertexArrays.cpp
index 055814b03b..c28df12cfa 100644
--- a/Post/PViewVertexArrays.cpp
+++ b/Post/PViewVertexArrays.cpp
@@ -1077,7 +1077,7 @@ static void addTensorElement(PView *p, int iEnt, int iEle, int numNodes, int typ
for (int k = 0; k < 3; k++) {
vval[k][0] = xyz[i][k];
for (int j = 0; j < 3; j++) {
- vval[k][j+1] = V(j,k)*S(j);
+ vval[k][j+1] = V(k,j)*S(j);
}
}
double lmax = std::max(S(0), std::max(S(1), S(2)));
@@ -1098,9 +1098,9 @@ static void addTensorElement(PView *p, int iEnt, int iEle, int numNodes, int typ
}
tensor.eig(S, imS, V, rightV, false);
for (int j = 0; j < 3; j++) {
- vval[0][j+1] += V(j,0)*S(j)/numNodes;
- vval[1][j+1] += V(j,1)*S(j)/numNodes;
- vval[2][j+1] += V(j,2)*S(j)/numNodes;
+ vval[0][j+1] += V(0,j)*S(j)/numNodes;
+ vval[1][j+1] += V(1,j)*S(j)/numNodes;
+ vval[2][j+1] += V(2,j)*S(j)/numNodes;
}
vval[0][0] += xyz[i][0]/numNodes;
vval[1][0] += xyz[i][1]/numNodes;
--
GitLab