From c5bd25dfd3af51e7c1520f3c0470cca3e34039a1 Mon Sep 17 00:00:00 2001
From: Gauthier Becker <gauthierbecker@gmail.com>
Date: Tue, 24 Jan 2012 08:20:34 +0000
Subject: [PATCH] Better initialization for STensor43 sysmmetric identity
---
Solver/STensor43.h | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/Solver/STensor43.h b/Solver/STensor43.h
index 5e3ecfe7f0..c764fde277 100644
--- a/Solver/STensor43.h
+++ b/Solver/STensor43.h
@@ -37,7 +37,8 @@ class STensor43 {
else
_val[getIndex(i, j, k, l)]=0.0;
}
- inline void symmetricIdentity(const double v = 1.0)
+ // Symmetric identity tensor
+ STensor43(const double vik, const double vil)
{
for (int i = 0; i < 3; i++)
for (int j = 0; j < 3; j++)
@@ -46,9 +47,9 @@ class STensor43 {
{
_val[getIndex(i, j, k, l)]= 0.;
if ((i==k)&&(j==l))
- _val[getIndex(i, j, k, l)]+=0.5*v;
+ _val[getIndex(i, j, k, l)]+=0.5*vik;
if ((i==l)&&(j==k))
- _val[getIndex(i, j, k, l)]+=0.5*v;
+ _val[getIndex(i, j, k, l)]+=0.5*vil;
}
}
inline double &operator()(int i, int j,int k, int l)
--
GitLab