Skip to content
Snippets Groups Projects
Commit df6a7112 authored by Christian Wieners's avatar Christian Wieners
Browse files

added Mittelvwert für V

parent e3e0039d
No related branches found
No related tags found
1 merge request!273fixed hopefully problem with computing displacement by depending mechassamble...
Pipeline #164614 failed
...@@ -474,7 +474,14 @@ void SemiImplicitSolverOnCells::SolvePDEOnCells(IElphyAssemble &A){ ...@@ -474,7 +474,14 @@ void SemiImplicitSolverOnCells::SolvePDEOnCells(IElphyAssemble &A){
for (cell c = (*potential).cells(); c != (*potential).cells_end(); ++c) { for (cell c = (*potential).cells(); c != (*potential).cells_end(); ++c) {
ScalarElement E((*potential),*c); ScalarElement E((*potential),*c);
(*vcw_c)[0](c(),0) = E.Value(c.LocalCenter(),(*potential)); double A = 0;
double V = 0;
for (int q = 0; q < E.nQ(); ++q) {
double w = E.QWeight(q);
A += w;
V += w * E.Value(q, *potential, 0);
(*vcw_c)[0](c(),0) = V / A; // E.Value(c.LocalCenter(),(*potential));
}
} }
} }
void SemiImplicitSolverOnCells::SolvePDE(IElphyAssemble &A){ void SemiImplicitSolverOnCells::SolvePDE(IElphyAssemble &A){
......
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