Commit a0d4a780 authored by Thomas KUMMER's avatar Thomas KUMMER
Browse files

euler lifev update

parent af9dd5dd
......@@ -314,9 +314,9 @@ public:
for (int p (0); p < n; ++p)
{
int i = f0_->blockMap().GID (p);
int j = f0_->blockMap().GID (p + n);
int k = f0_->blockMap().GID (p + 2 * n);
int i = f_->blockMap().GID (p);
int j = f_->blockMap().GID (p + n);
int k = f_->blockMap().GID (p + 2 * n);
F *= 0.;
F(0,0) = 1.; F(1,1) = 1.; F(2,2) = 1.;
......@@ -328,7 +328,7 @@ public:
// f0.normalize();
//
auto f = F * f0;
(*f_)[i] = f(0);
(*f_)[i] = (*f0_)[i];//f(0);
(*f_)[j] = f(1);
(*f_)[k] = f(2);
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment