Skip to content
Snippets Groups Projects
Commit 76b11a3f authored by Oliver Sander's avatar Oliver Sander Committed by sander@PCPOOL.MI.FU-BERLIN.DE
Browse files

bugfix: local stress must be a std::vector

[[Imported from SVN: r4082]]
parent e8bc81d1
No related branches found
No related tags found
No related merge requests found
...@@ -299,10 +299,13 @@ getResultantForce(const BoundaryPatch<GridType>& boundary, ...@@ -299,10 +299,13 @@ getResultantForce(const BoundaryPatch<GridType>& boundary,
double pos = nIt->intersectionSelfLocal().corner(0); double pos = nIt->intersectionSelfLocal().corner(0);
Dune::array<RigidBodyMotion<3>,2> localSolution = {sol[indexSet.template subIndex<1>(*eIt,0)], std::vector<RigidBodyMotion<3> > localSolution(2);
sol[indexSet.template subIndex<1>(*eIt,1)]}; localSolution[0] = sol[indexSet.template subIndex<1>(*eIt,0)];
Dune::array<RigidBodyMotion<3>,2> localRefConf = {dynamic_cast<RodLocalStiffness<typename GridType::LeafGridView, double>* >(this->localStiffness_)->referenceConfiguration_[indexSet.template subIndex<1>(*eIt,0)], localSolution[1] = sol[indexSet.template subIndex<1>(*eIt,1)];
dynamic_cast<RodLocalStiffness<typename GridType::LeafGridView, double>* >(this->localStiffness_)->referenceConfiguration_[indexSet.template subIndex<1>(*eIt,1)]};
std::vector<RigidBodyMotion<3> > localRefConf(2);
localRefConf[0] = dynamic_cast<RodLocalStiffness<typename GridType::LeafGridView, double>* >(this->localStiffness_)->referenceConfiguration_[indexSet.template subIndex<1>(*eIt,0)];
localRefConf[1] = dynamic_cast<RodLocalStiffness<typename GridType::LeafGridView, double>* >(this->localStiffness_)->referenceConfiguration_[indexSet.template subIndex<1>(*eIt,1)];
FieldVector<double, blocksize> strain = dynamic_cast<RodLocalStiffness<typename GridType::LeafGridView, double>* >(this->localStiffness_)->getStrain(localSolution, *eIt, pos); FieldVector<double, blocksize> strain = dynamic_cast<RodLocalStiffness<typename GridType::LeafGridView, double>* >(this->localStiffness_)->getStrain(localSolution, *eIt, pos);
FieldVector<double, blocksize> referenceStrain = dynamic_cast<RodLocalStiffness<typename GridType::LeafGridView, double>* >(this->localStiffness_)->getStrain(localRefConf, *eIt, pos); FieldVector<double, blocksize> referenceStrain = dynamic_cast<RodLocalStiffness<typename GridType::LeafGridView, double>* >(this->localStiffness_)->getStrain(localRefConf, *eIt, pos);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment