Skip to content
Snippets Groups Projects
Commit 9dcd522a authored by Ling Wu's avatar Ling Wu
Browse files

symmetric

parent 59424969
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
Showing
with 6 additions and 3 deletions
Solver/STensor43.h
+ 6
− 3
View file @ 9dcd522a
...@@ -31,10 +31,13 @@ class STensor43 { ...@@ -31,10 +31,13 @@ class STensor43 {
for (int j = 0; j < 3; j++) for (int j = 0; j < 3; j++)
for (int k = 0; k < 3; k++) for (int k = 0; k < 3; k++)
for (int l = 0; l < 3; l++) for (int l = 0; l < 3; l++)
{
_val[getIndex(i, j, k, l)]= 0.;
if ((i==k)&&(j==l)) if ((i==k)&&(j==l))
_val[getIndex(i, j, k, l)]=v; _val[getIndex(i, j, k, l)]+=0.5*v;
else if ((i==l)&&(j==k))
_val[getIndex(i, j, k, l)]=0.0; _val[getIndex(i, j, k, l)]+=0.5*v;
}
} }
inline double &operator()(int i, int j,int k, int l) inline double &operator()(int i, int j,int k, int l)
{ {
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment