Commit 9b7a5789 authored by laura.pfeiffer's avatar laura.pfeiffer
Browse files

added Evaluation points

parent 69743b40
......@@ -57,7 +57,7 @@ ElphyPolynomialDegree = 1
StartTime = 0.0
EndTime = 0.1
#EndTime = 0.1
EndTime=0.0001
EndTime = 0.3
......
......@@ -7,16 +7,16 @@ constexpr std::array<double, 2> potentialValues{-84.584430, -85.298003};
std::string ElphyVentricleProblem::Evaluate(const Vector &solution) const {
std::string evaluation{};
auto values = getEvaluationAt(solution, {Point(82.53800201416016, -88.54129791259766, -55.96120071411133)});
auto values = getEvaluationAt(solution, {Point(74.837097, -81.447998, -44.074902)});
evaluation += "V(-33.4778, -27.0649, 44.7976) = " + std::to_string(values[0]) + "\n";
evaluation += "V(26.297, 10.486, 7.529) = " + std::to_string(values[1]) + "\n";
evaluation += "V(74.837097, -81.447998, -44.074902) = " + std::to_string(values[0]) + "\n";
//evaluation += "V(26.297, 10.486, 7.529) = " + std::to_string(values[1]) + "\n";
return evaluation;
}
std::vector<double> ElphyVentricleProblem::EvaluationResults(const Vector &solution) const {
auto values = getEvaluationAt(solution, {Point(82.53800201416016, -88.54129791259766, -55.96120071411133)});
auto values = getEvaluationAt(solution, {Point(74.837097, -81.447998, -44.074902)});
for (int i = 0; i < values.size(); ++i) {
values[i] -= potentialValues[i];
......@@ -26,10 +26,13 @@ std::vector<double> ElphyVentricleProblem::EvaluationResults(const Vector &solut
void ElphyVentricleProblem::InitializeEvaluationPoints() {
evaluationPoints.emplace_back(Point(82.53800201416016, -88.54129791259766, -55.96120071411133));// Apex
//evaluationPoints.emplace_back(Point(82.538, -88.5413, -55.9612));
evaluationPoints.emplace_back(Point(19.492599487304688, -74.51229858398438, -31.519500732421875));// Innenwand Ventricle angeregt
//evaluationPoints.emplace_back(Point(19.4926, -74.5123, -31.5195));
evaluationPoints.emplace_back(Point(74.837097, -81.447998, -44.074902));// Apex excited left ventricle
evaluationPoints.emplace_back(Point(74.327103, -86.796303, -49.017601));//Apex midwall left ventricle
evaluationPoints.emplace_back(Point(39.611698, -83.080399, -34.947701));// Innenwand between ventricles bottom
evaluationPoints.emplace_back(Point(18.007799, -63.959301, -19.1492));// Innenwand between ventricles middle
evaluationPoints.emplace_back(Point(-0.33377, -44.208099, -7.27446));// Innenwand between ventricles top
evaluationPoints.emplace_back(Point(64.5504, -27.874399, -25.565701));// Innenwand Mitte linkes Ventrikel
evaluationPoints.emplace_back(Point(-8.8439, -95.517899, -19.521999)); //Innenwand Mitte rechtes Ventrikel
}
......
......@@ -8,7 +8,7 @@
class ElphyVentricleProblem : public ElphyProblem {
std::string meshName;
public:
explicit ElphyVentricleProblem() : ElphyProblem("Ventricle", true), meshName("BiVentriclefine") {
explicit ElphyVentricleProblem() : ElphyProblem("Ventricle", true), meshName("BiVentricle") {
Tensor SV(0.2073, 0.0, 0.0,
0.0, 0.0921, 0.0,
0.0, 0.0, 0.023);
......
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