From bcaa0d85e5221d256cb2b7388a4ba24bb12272f8 Mon Sep 17 00:00:00 2001
From: Laura Pfeiffer <laura.pfeiffer@kit.edu>
Date: Fri, 8 Dec 2023 09:36:51 +0100
Subject: [PATCH] change evaluation of potential
---
src/electrophysiology/solvers/LinearImplicitSolver.cpp | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/src/electrophysiology/solvers/LinearImplicitSolver.cpp b/src/electrophysiology/solvers/LinearImplicitSolver.cpp
index 09a8fd3c8..26dc1d0f7 100644
--- a/src/electrophysiology/solvers/LinearImplicitSolver.cpp
+++ b/src/electrophysiology/solvers/LinearImplicitSolver.cpp
@@ -474,14 +474,16 @@ void SemiImplicitSolverOnCells::SolvePDEOnCells(IElphyAssemble &A){
for (cell c = (*potential).cells(); c != (*potential).cells_end(); ++c) {
ScalarElement E((*potential),*c);
+ /*
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));
- }
+ (*vcw_c)[0](c(),0) = V / A; //
+ }*/
+ (*vcw_c)[0](c(),0) = E.Value(c.LocalCenter(),(*potential));
}
}
void SemiImplicitSolverOnCells::SolvePDE(IElphyAssemble &A){
--
GitLab