Commit 69ed9498 authored by thomaskummer's avatar thomaskummer
Browse files

euler lifev update

parent 3e873dae
...@@ -1078,53 +1078,58 @@ EMSolver<Mesh, ElectroSolver>::computeDeformedFiberDirection (boost::shared_ptr< ...@@ -1078,53 +1078,58 @@ EMSolver<Mesh, ElectroSolver>::computeDeformedFiberDirection (boost::shared_ptr<
// dUdy = feSpacePtr->gradientRecovery(disp, 1); // dUdy = feSpacePtr->gradientRecovery(disp, 1);
// dUdz = feSpacePtr->gradientRecovery(disp, 2); // dUdz = feSpacePtr->gradientRecovery(disp, 2);
// int n = f_->epetraVector().MyLength() / 3; int n = f_->epetraVector().MyLength() / 3;
// int i (0); int j (0); int k (0); int i (0); int j (0); int k (0);
// MatrixSmall<3,3> F; VectorSmall<3> f0; MatrixSmall<3,3> F; VectorSmall<3> f0;
*f_ = *f0_; return; // *f_ = *f0_; return;
// for (int p (0); p < n; p++) for (int p (0); p < n; p++)
// { {
// i = f_->blockMap().GID (p); i = f_->blockMap().GID (p);
// j = f_->blockMap().GID (p + n); j = f_->blockMap().GID (p + n);
// k = f_->blockMap().GID (p + 2 * n); k = f_->blockMap().GID (p + 2 * n);
//
//// F(0,0) = 1.0 + dUdx[i]; // F(0,0) = 1.0 + dUdx[i];
//// F(0,1) = dUdy[i]; // F(0,1) = dUdy[i];
//// F(0,2) = dUdz[i]; // F(0,2) = dUdz[i];
//// F(1,0) = dUdx[j]; // F(1,0) = dUdx[j];
//// F(1,1) = 1.0 + dUdy[j]; // F(1,1) = 1.0 + dUdy[j];
//// F(1,2) = dUdz[j]; // F(1,2) = dUdz[j];
//// F(2,0) = dUdx[k]; // F(2,0) = dUdx[k];
//// F(2,1) = dUdy[k]; // F(2,1) = dUdy[k];
//// F(2,2) = 1.0 + dUdz[k]; // F(2,2) = 1.0 + dUdz[k];
//
//// this is right // this is right
//// F(0,0) = 1.0 + dUdx[i]; // F(0,0) = 1.0 + dUdx[i];
//// F(0,1) = dUdx[j]; // F(0,1) = dUdx[j];
//// F(0,2) = dUdx[k]; // F(0,2) = dUdx[k];
//// F(1,0) = dUdy[i]; // F(1,0) = dUdy[i];
//// F(1,1) = 1.0 + dUdy[j]; // F(1,1) = 1.0 + dUdy[j];
//// F(1,2) = dUdy[k]; // F(1,2) = dUdy[k];
//// F(2,0) = dUdz[i]; // F(2,0) = dUdz[i];
//// F(2,1) = dUdz[j]; // F(2,1) = dUdz[j];
//// F(2,2) = 1.0 + dUdz[k]; // F(2,2) = 1.0 + dUdz[k];
//
// F *= 0.; if ( f_->blockMap().LID ( f_->blockMap().GID (i) ) != -1 )
// F(0,0) = 1.; F(1,1) = 1.; F(2,2) = 1.; {
// //(*f_) ( f_->blockMap().GID (j) ) = value;
// f0(0) = (*f0_)[i];
// f0(1) = (*f0_)[j]; F *= 0.;
// f0(2) = (*f0_)[k]; F(0,0) = 1.; F(1,1) = 1.; F(2,2) = 1.;
////
//// f0.normalize(); f0(0) = (*f0_)[i];
//// f0(1) = (*f0_)[j];
// auto f = F * f0; f0(2) = (*f0_)[k];
// (*f_)[i] = f(0); //
// (*f_)[j] = f(1); // f0.normalize();
// (*f_)[k] = f(2); //
// } auto f = F * f0;
(*f_)[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