From 9dcd522a4181b6e03709410b60680f5db45a89a4 Mon Sep 17 00:00:00 2001
From: Ling Wu <l.wu@ulg.ac.be>
Date: Thu, 19 Jan 2012 17:20:42 +0000
Subject: [PATCH] symmetric
---
Solver/STensor43.h | 9 ++++++---
1 file changed, 6 insertions(+), 3 deletions(-)
diff --git a/Solver/STensor43.h b/Solver/STensor43.h
index 906f7451f5..7ea38c1f24 100644
--- a/Solver/STensor43.h
+++ b/Solver/STensor43.h
@@ -31,10 +31,13 @@ class STensor43 {
for (int j = 0; j < 3; j++)
for (int k = 0; k < 3; k++)
for (int l = 0; l < 3; l++)
+ {
+ _val[getIndex(i, j, k, l)]= 0.;
if ((i==k)&&(j==l))
- _val[getIndex(i, j, k, l)]=v;
- else
- _val[getIndex(i, j, k, l)]=0.0;
+ _val[getIndex(i, j, k, l)]+=0.5*v;
+ if ((i==l)&&(j==k))
+ _val[getIndex(i, j, k, l)]+=0.5*v;
+ }
}
inline double &operator()(int i, int j,int k, int l)
{
--
GitLab