From df6a71120113c288c4952108253ad56429828822 Mon Sep 17 00:00:00 2001
From: wieners <christian.wieners@kit.edu>
Date: Mon, 27 Nov 2023 15:41:46 +0100
Subject: [PATCH] =?UTF-8?q?added=20Mittelvwert=20f=C3=BCr=20V?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
---
src/electrophysiology/solvers/LinearImplicitSolver.cpp | 9 ++++++++-
1 file changed, 8 insertions(+), 1 deletion(-)
diff --git a/src/electrophysiology/solvers/LinearImplicitSolver.cpp b/src/electrophysiology/solvers/LinearImplicitSolver.cpp
index b0b65bb2a..09a8fd3c8 100644
--- a/src/electrophysiology/solvers/LinearImplicitSolver.cpp
+++ b/src/electrophysiology/solvers/LinearImplicitSolver.cpp
@@ -474,7 +474,14 @@ void SemiImplicitSolverOnCells::SolvePDEOnCells(IElphyAssemble &A){
for (cell c = (*potential).cells(); c != (*potential).cells_end(); ++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){
--
GitLab