...
 
Commits (2)
......@@ -8,6 +8,7 @@
#include <duneuro/eeg/source_model_interface.hh>
#include <duneuro/eeg/unfitted_partial_integration_source_model.hh>
#include <duneuro/eeg/unfitted_patch_based_venant_source_model.hh>
#include <duneuro/eeg/unfitted_truncated_spatial_venant_source_model.hh>
namespace duneuro
{
......@@ -31,6 +32,11 @@ namespace duneuro
Vector>>(solver.functionSpace().getGFS(), solver.subTriangulation(),
solver.elementSearch(), config.get<std::size_t>("compartment"), scaleToBBox,
config);
} else if (type == "truncated_spatial_venant") {
return Dune::Std::make_unique<UnfittedTruncatedSpatialVenantSourceModel<
typename Solver::Traits::FunctionSpace::GFS, typename Solver::Traits::SubTriangulation,
Vector>>(solver.functionSpace().getGFS(), solver.subTriangulation(),
solver.elementSearch(), config.get<std::size_t>("compartment"), config);
} else {
DUNE_THROW(duneuro::SourceModelException, "unknown source model of type \"" << type
<< "\"");
......@@ -56,6 +62,11 @@ namespace duneuro
Vector>>(solver.functionSpace().getGFS(), solver.subTriangulation(),
solver.elementSearch(), config.get<std::size_t>("compartment"), scaleToBBox,
config);
} else if (type == "truncated_spatial_venant") {
return Dune::Std::make_unique<UnfittedTruncatedSpatialVenantSourceModel<
typename Solver::Traits::FunctionSpace::GFS, typename Solver::Traits::SubTriangulation,
Vector>>(solver.functionSpace().getGFS(), solver.subTriangulation(),
solver.elementSearch(), config.get<std::size_t>("compartment"), config);
} else {
DUNE_THROW(duneuro::SourceModelException, "unknown source model of type \"" << type
<< "\"");
......
......@@ -9,6 +9,7 @@
#include <duneuro/eeg/udg_subtraction_source_model.hh>
#include <duneuro/eeg/unfitted_partial_integration_source_model.hh>
#include <duneuro/eeg/unfitted_patch_based_venant_source_model.hh>
#include <duneuro/eeg/unfitted_truncated_spatial_venant_source_model.hh>
namespace duneuro
{
......@@ -37,6 +38,11 @@ namespace duneuro
typename Solver::Traits::FunctionSpace, typename Solver::Traits::SubTriangulation,
Vector>>(solver.functionSpace(), solver.subTriangulation(), solver.elementSearch(),
config.get<std::size_t>("compartment"), config, solverConfig);
} else if (type == "truncated_spatial_venant") {
return Dune::Std::make_unique<UnfittedTruncatedSpatialVenantSourceModel<
typename Solver::Traits::FunctionSpace::GFS, typename Solver::Traits::SubTriangulation,
Vector>>(solver.functionSpace().getGFS(), solver.subTriangulation(),
solver.elementSearch(), config.get<std::size_t>("compartment"), config);
} else {
DUNE_THROW(duneuro::SourceModelException, "unknown source model of type \"" << type
<< "\"");
......@@ -62,6 +68,11 @@ namespace duneuro
Vector>>(solver.functionSpace().getGFS(), solver.subTriangulation(),
solver.elementSearch(), config.get<std::size_t>("compartment"), scaleToBBox,
config);
} else if (type == "truncated_spatial_venant") {
return Dune::Std::make_unique<UnfittedTruncatedSpatialVenantSourceModel<
typename Solver::Traits::FunctionSpace::GFS, typename Solver::Traits::SubTriangulation,
Vector>>(solver.functionSpace().getGFS(), solver.subTriangulation(),
solver.elementSearch(), config.get<std::size_t>("compartment"), config);
} else {
DUNE_THROW(duneuro::SourceModelException, "unknown source model of type \"" << type
<< "\"");
......
This diff is collapsed.