def inject_Example(module, exceptions, interfaces, CONFIG_H):
'''injects the user code into the module'''
# first the discretization
GridType = 'Dune::SGrid< 2, 2 >'
GridLayerType = 'Dune::Stuff::Grid::ChooseLayer::leaf'
RangeFieldType = 'double'
dimRange = '1'
polOrder = '1'
SpaceBackendType = 'Dune::GDT::ChooseSpaceBackend::pdelab'
LaBackendType = 'Dune::Stuff::LA::ChooseBackend::istl_sparse'
if 'istl_sparse' in LaBackendType:
MatrixType = 'Dune::Stuff::LA::IstlRowMajorSparseMatrix< ' + RangeFieldType + ' >'
VectorType = 'Dune::Stuff::LA::IstlDenseVector< ' + RangeFieldType + ' >'
elif 'eigen_sparse' in LaBackendType:
MatrixType = 'Dune::Stuff::LA::EigenRowMajorSparseMatrix< ' + RangeFieldType + ' >'
VectorType = 'Dune::Stuff::LA::EigenDenseVector< ' + RangeFieldType + ' >'
DiscretizationName = 'Dune::HDD::LinearElliptic::Discretizations::CG'
DiscretizationFullName = (DiscretizationName + '< '
+ GridType + ', ' + GridLayerType + ', '
+ RangeFieldType + ', ' + dimRange + ', ' + polOrder + ', '
+ SpaceBackendType + ', '
+ LaBackendType + ' >')
discretization = inject_StationaryDiscretizationImplementation(
module, exceptions, interfaces, CONFIG_H,
DiscretizationName,
Traits={'VectorType': VectorType,
'OperatorType': 'Dune::Pymor::Operators::LinearAffinelyDecomposedContainerBased< ' + MatrixType + ', ' + VectorType + ' >',
'FunctionalType': 'Dune::Pymor::Functionals::LinearAffinelyDecomposedVectorBased< ' + VectorType + ' >',
'ProductType': 'Dune::Pymor::Operators::LinearAffinelyDecomposedContainerBased< ' + MatrixType + ', ' + VectorType + ' > '},
template_parameters=[GridType, GridLayerType, RangeFieldType, dimRange, polOrder, SpaceBackendType, LaBackendType ])
# then add the example
LinearellipticExampleCG = module.add_class('LinearellipticExampleCG',
template_parameters=['Dune::SGrid< 2, 2 >'],
custom_name='LinearellipticExampleCG')
LinearellipticExampleCG.add_method('static_id',
retval('std::string'),
[], is_const=True, throw=[exceptions['Exception']])
LinearellipticExampleCG.add_method('write_config_file',
None, [], is_const=True, throw=[exceptions['Exception']])
LinearellipticExampleCG.add_constructor([], throw=[exceptions['Exception']])
LinearellipticExampleCG.add_method('initialize', None,
[param('const std::vector< std::string >', 'arguments')],
is_const=True, throw=[exceptions['Exception']])
LinearellipticExampleCG.add_method('discretization_and_return_ptr',
retval(DiscretizationFullName + ' *', caller_owns_return=True),
[], is_const=True, throw=[exceptions['Exception']],
custom_name='discretization')
def add_example(GridType, space_backend, la_backend, name):
# build all types needed for the discretization
dimRange = '1'
polOrder = '1'
grid_layer = 'Dune::Stuff::Grid::ChooseLayer::leaf'
MatrixType = 'Dune::Stuff::LA::'
VectorType = 'Dune::Stuff::LA::'
if 'eigen_sparse' in la_backend:
MatrixType += 'EigenRowMajorSparseMatrix'
VectorType += 'EigenDenseVector'
elif 'istl_sparse' in la_backend:
MatrixType += 'IstlRowMajorSparseMatrix'
VectorType += 'IstlDenseVector'
MatrixType += '< ' + RangeFieldType + ' >'
VectorType += '< ' + RangeFieldType + ' >'
OperatorType = 'Dune::Pymor::Operators::LinearAffinelyDecomposedContainerBased< ' + MatrixType + ', ' + VectorType + ' >'
ProductType = OperatorType
FunctionalType = 'Dune::Pymor::Functionals::LinearAffinelyDecomposedVectorBased< ' + VectorType + ' >'
DiscretizationName = 'Dune::HDD::LinearElliptic::Discretizations::CG'
DiscretizationType = (DiscretizationName + '< ' + GridType + ', ' +
grid_layer + ', ' + RangeFieldType + ', ' +
dimRange + ', ' + polOrder + ', ' +
space_backend + ', ' + la_backend + '>')
inject_StationaryDiscretizationImplementation(module,
exceptions,
interfaces,
CONFIG_H,
DiscretizationName,
Traits={
'VectorType':
VectorType,
'OperatorType':
OperatorType,
'FunctionalType':
FunctionalType,
'ProductType':
ProductType
},
template_parameters=[
GridType, grid_layer,
RangeFieldType,
dimRange, polOrder,
space_backend,
la_backend
])
# then create the example
Example = module.add_class(
'PbSpe10Example',
template_parameters=[GridType, space_backend, la_backend],
custom_name=name)
Example.add_constructor([
param('const std::string&', 'num_grid_elements'),
param('const bool', 'logging_enabled')
],
throw=exceptions)
Example.add_constructor(
[param('const std::string&', 'num_grid_elements')],
throw=exceptions)
Example.add_constructor([], throw=exceptions)
Example.add_method('pb_discretization_and_return_ptr',
retval(DiscretizationType + ' *',
caller_owns_return=True), [],
is_const=True,
throw=exceptions,
custom_name='discretization')
Example.add_method('visualize_grid',
None,
[param('const std::string&', 'filename_prefix')],
is_const=True,
throw=exceptions)
Example.add_method('visualize_problem',
None,
[param('const std::string&', 'filename_prefix')],
is_const=True,
throw=exceptions)
Example.add_method('visualize_problem',
None, [
param('const std::string&', 'filename_prefix'),
param('const Dune::Pymor::Parameter&', 'mu')
],
is_const=True,
throw=exceptions)
if space_backend == space_backend_fem:
Example.add_method(
'visualize_darcy_velocity',
None, [
param('const {}&'.format(VectorType), 'vector'),
param('const std::string&', 'filename'),
param('const std::string&', 'name'),
param('const Dune::Pymor::Parameter&', 'mu')
],
is_const=True,
throw=exceptions)
def inject_Example(module, exceptions, interfaces, CONFIG_H):
ssize_t = CONFIG_H['DUNE_STUFF_SSIZE_T']
'''injects the user code into the module'''
# first the discretization
GridType = 'Dune::ALUGrid< 2, 2, Dune::simplex, Dune::conforming >'
#GridType = 'Dune::SGrid< 2, 2 >'
RangeFieldType = 'double'
dimRange = '1'
polOrder = '1'
la_backend = 'Dune::Stuff::LA::ChooseBackend::eigen_sparse'
if la_backend == 'Dune::Stuff::LA::ChooseBackend::eigen_sparse':
MatrixType = 'Dune::Stuff::LA::EigenRowMajorSparseMatrix< ' + RangeFieldType + ' >'
VectorType = 'Dune::Stuff::LA::EigenDenseVector< ' + RangeFieldType + ' >'
elif la_backend == 'Dune::Stuff::LA::ChooseBackend::istl_sparse':
MatrixType = 'Dune::Stuff::LA::IstlRowMajorSparseMatrix< ' + RangeFieldType + ' >'
VectorType = 'Dune::Stuff::LA::IstlDenseVector< ' + RangeFieldType + ' >'
space_backend = 'Dune::GDT::ChooseSpaceBackend::fem'
OperatorType = 'Dune::Pymor::Operators::LinearAffinelyDecomposedContainerBased< ' + MatrixType + ', ' + VectorType + ' >'
ProductType = OperatorType
FunctionalType = 'Dune::Pymor::Functionals::LinearAffinelyDecomposedVectorBased< ' + VectorType + ' >'
DiscretizationName = 'Dune::HDD::LinearElliptic::Discretizations::BlockSWIPDG'
DiscretizationFullName = (DiscretizationName + '< '
+ GridType + ', '
+ RangeFieldType + ', '
+ dimRange + ', ' + polOrder + ', '
+ la_backend + '>')
LocalDiscretizationName = 'Dune::HDD::LinearElliptic::Discretizations::SWIPDG'
local_layer = 'Dune::Stuff::Grid::ChooseLayer::local'
LocalDiscretizationFullName = (LocalDiscretizationName
+ '< ' + GridType + ', '
+ local_layer + ', '
+ RangeFieldType + ', '
+ dimRange + ', ' + polOrder + ', '
+ space_backend + ', '
+ la_backend + ' >')
inject_StationaryDiscretizationImplementation(module, exceptions, interfaces, CONFIG_H,
LocalDiscretizationName,
Traits={'VectorType': VectorType,
'OperatorType': OperatorType,
'FunctionalType': FunctionalType,
'ProductType': ProductType},
template_parameters=[GridType, local_layer, RangeFieldType,
dimRange, polOrder, space_backend, la_backend])
oversampled_layer = 'Dune::Stuff::Grid::ChooseLayer::local_oversampled'
OversampledDiscretizationFullName = (LocalDiscretizationName
+ '< ' + GridType + ', '
+ oversampled_layer + ', '
+ RangeFieldType + ', '
+ dimRange + ', ' + polOrder + ', '
+ space_backend + ', '
+ la_backend + ' >')
inject_StationaryDiscretizationImplementation(module, exceptions, interfaces, CONFIG_H,
LocalDiscretizationName,
Traits={'VectorType': VectorType,
'OperatorType': OperatorType,
'FunctionalType': FunctionalType,
'ProductType': ProductType},
template_parameters=[GridType, oversampled_layer, RangeFieldType,
dimRange, polOrder, space_backend, la_backend])
# discretization = inject_StationaryDiscretizationImplementation(
inject_StationaryMultiscaleDiscretizationImplementation(
module, exceptions, interfaces, CONFIG_H,
DiscretizationName,
Traits={'VectorType': VectorType,
'OperatorType': OperatorType,
'FunctionalType': FunctionalType,
'ProductType': ProductType,
'LocalDiscretizationType': LocalDiscretizationFullName,
'OversampledDiscretizationType': OversampledDiscretizationFullName},
template_parameters=[GridType, RangeFieldType,
dimRange, polOrder, la_backend])
# then add the example
def add_example(name):
Example = module.add_class(name, template_parameters=[GridType], custom_name=name)
Example.add_constructor([param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels'),
param('const ' + ssize_t, 'debug_log_levels'),
param('const bool', 'enable_warnings'),
param('const bool', 'enable_colors'),
param('const std::string', 'info_color'),
param('const std::string', 'debug_color'),
param('const std::string', 'warn_color')],
throw=exceptions)
Example.add_constructor([param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels'),
param('const ' + ssize_t, 'debug_log_levels'),
param('const bool', 'enable_warnings'),
param('const bool', 'enable_colors'),
param('const std::string', 'info_color'),
param('const std::string', 'debug_color')],
throw=exceptions)
Example.add_constructor([param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels'),
param('const ' + ssize_t, 'debug_log_levels'),
param('const bool', 'enable_warnings'),
param('const bool', 'enable_colors'),
param('const std::string', 'info_color')],
throw=exceptions)
Example.add_constructor([param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels'),
param('const ' + ssize_t, 'debug_log_levels'),
param('const bool', 'enable_warnings'),
param('const bool', 'enable_colors')],
throw=exceptions)
Example.add_constructor([param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels'),
param('const ' + ssize_t, 'debug_log_levels'),
param('const bool', 'enable_warnings')],
throw=exceptions)
Example.add_constructor([param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels'),
param('const ' + ssize_t, 'debug_log_levels')],
throw=exceptions)
Example.add_constructor([param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels')],
throw=exceptions)
Example.add_constructor([param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference')],
throw=exceptions)
Example.add_constructor([param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products')],
throw=exceptions)
Example.add_constructor([param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers')],
throw=exceptions)
Example.add_constructor([param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements')],
throw=exceptions)
Example.add_constructor([param('const std::string', 'partitioning')],
throw=exceptions)
Example.add_constructor([], throw=exceptions)
Example.add_method('discretization_and_return_ptr',
retval(DiscretizationFullName + ' *', caller_owns_return=True),
[], is_const=True, throw=exceptions,
custom_name='discretization')
Example.add_method('visualize',
None,
[param('const std::string&', 'filename_prefix')],
is_const=True, throw=exceptions)
Example.add_method('project',
retval(VectorType),
[param('const std::string', 'expression')], is_const=True, throw=exceptions)
Example.add_method('compute_error',
retval(RangeFieldType),
[param('const ' + VectorType + '&', 'solution'),
param('const std::string', 'product_type'),
param('const Dune::Pymor::Parameter', 'mu'),
param('const Dune::Pymor::Parameter', 'mu_product')],
is_const=True, throw=exceptions)
Example.add_method('compute_error',
retval(RangeFieldType),
[param('const ' + VectorType + '&', 'solution'),
param('const std::string', 'product_type'),
param('const Dune::Pymor::Parameter', 'mu')],
is_const=True, throw=exceptions)
Example.add_method('compute_error',
retval(RangeFieldType),
[param('const ' + VectorType + '&', 'solution'),
param('const std::string', 'product_type')],
is_const=True, throw=exceptions)
Example.add_method('compute_error',
retval(RangeFieldType),
[param('const ' + VectorType + '&', 'solution')],
is_const=True, throw=exceptions)
Example.add_method('available_estimators',
retval('std::vector< std::string >'),
[],
is_const=True, throw=exceptions)
Example.add_method('estimate',
retval(RangeFieldType),
[param('const ' + VectorType + '&', 'vector'),
param('const std::string', 'type'),
param('const Dune::Pymor::Parameter', 'mu_hat'),
param('const Dune::Pymor::Parameter', 'mu_bar'),
param('const Dune::Pymor::Parameter', 'mu')],
is_const=True, throw=exceptions)
Example.add_method('estimate',
retval(RangeFieldType),
[param('const ' + VectorType + '&', 'vector'),
param('const std::string', 'type')],
is_const=True, throw=exceptions)
Example.add_method('available_local_estimators',
retval('std::vector< std::string >'),
[],
is_const=True, throw=exceptions)
Example.add_method('estimate_local',
retval('std::vector< ' + RangeFieldType + ' >'),
[param('const ' + VectorType + '&', 'vector'),
param('const std::string', 'type'),
param('const Dune::Pymor::Parameter', 'mu_hat'),
param('const Dune::Pymor::Parameter', 'mu_bar'),
param('const Dune::Pymor::Parameter', 'mu')],
is_const=True, throw=exceptions)
Example.add_method('estimate_local',
retval('std::vector< ' + RangeFieldType + ' >'),
[param('const ' + VectorType + '&', 'vector'),
param('const std::string', 'type')],
is_const=True, throw=exceptions)
Example.add_method('pb_project_global_to_oversampled',
retval(VectorType + '*', caller_owns_return=True),
[param('const ' + VectorType + '&', 'global_vector'),
param('const ' + ssize_t, 'subdomain')],
is_const=True, throw=exceptions,
custom_name='project_global_to_oversampled')
Example.add_method('pb_project_global_to_local',
retval(VectorType + '*', caller_owns_return=True),
[param('const ' + VectorType + '&', 'global_vector'),
param('const ' + ssize_t, 'subdomain')],
is_const=True, throw=exceptions,
custom_name='project_global_to_local')
Example.add_method('pb_project_oversampled_to_local',
retval(VectorType + '*', caller_owns_return=True),
[param('const ' + VectorType + '&', 'oversampled_vector'),
param('const ' + ssize_t, 'subdomain')],
is_const=True, throw=exceptions,
custom_name='project_oversampled_to_local')
Example.add_method('solve_for_local_correction',
retval(VectorType),
[param('std::vector< ' + VectorType + ' >&', 'local_vectors'),
param('const ' + ssize_t, 'subdomain'),
param('const Dune::Pymor::Parameter', 'mu')],
is_const=True, throw=exceptions)
Example.add_method('solve_for_local_correction',
retval(VectorType),
[param('std::vector< ' + VectorType + ' >&', 'local_vectors'),
param('const ' + CONFIG_H['DUNE_STUFF_SSIZE_T'], 'subdomain')],
is_const=True, throw=exceptions)
Example.add_method('solve_oversampled',
retval(VectorType),
[param('const ' + ssize_t, 'subdomain'),
param('const std::string&', 'boundary_value_type'),
param('const ' + VectorType + '&', 'boundary_value_type'),
param('const Dune::Pymor::Parameter', 'mu')],
is_const=True, throw=exceptions)
Example.add_method('solve_oversampled',
retval(VectorType),
[param('const ' + ssize_t, 'subdomain'),
param('const std::string&', 'boundary_value_type'),
param('const ' + VectorType + '&', 'boundary_value_type')],
is_const=True, throw=exceptions)
Example.add_method('visualize_on_coarse_grid',
None,
[param('const std::vector< double >&', 'vector'),
param('const std::string&', 'filename'),
param('const std::string&', 'name')],
is_const=True, throw=exceptions)
Example.add_method('alpha',
retval(RangeFieldType),
[param('const Dune::Pymor::Parameter&', 'mu_1'),
param('const Dune::Pymor::Parameter&', 'mu_2')],
is_const=True, throw=exceptions)
Example.add_method('gamma',
retval(RangeFieldType),
[param('const Dune::Pymor::Parameter&', 'mu_1'),
param('const Dune::Pymor::Parameter&', 'mu_2')],
is_const=True, throw=exceptions)
add_example('OS2015MultiscaleExample')
add_example('OS2015AcademicExample')
def add_example(GridType, space_backend, la_backend, name):
# build all types needed for the discretization
dimRange = '1'
polOrder = '1'
grid_layer = 'Dune::Stuff::Grid::ChooseLayer::leaf'
MatrixType = 'Dune::Stuff::LA::'
VectorType = 'Dune::Stuff::LA::'
if 'eigen_sparse' in la_backend:
MatrixType += 'EigenRowMajorSparseMatrix'
VectorType += 'EigenDenseVector'
elif 'istl_sparse' in la_backend:
MatrixType += 'IstlRowMajorSparseMatrix'
VectorType += 'IstlDenseVector'
MatrixType += '< ' + RangeFieldType + ' >'
VectorType += '< ' + RangeFieldType + ' >'
OperatorType = 'Dune::Pymor::Operators::LinearAffinelyDecomposedContainerBased< ' + MatrixType + ', ' + VectorType + ' >'
ProductType = OperatorType
FunctionalType = 'Dune::Pymor::Functionals::LinearAffinelyDecomposedVectorBased< ' + VectorType + ' >'
DiscretizationName = 'Dune::HDD::LinearElliptic::Discretizations::CG'
DiscretizationType = (DiscretizationName + '< '
+ GridType + ', ' + grid_layer + ', '
+ RangeFieldType + ', '
+ dimRange + ', ' + polOrder + ', '
+ space_backend + ', ' + la_backend + '>')
inject_StationaryDiscretizationImplementation(module, exceptions, interfaces, CONFIG_H,
DiscretizationName,
Traits={'VectorType': VectorType,
'OperatorType': OperatorType,
'FunctionalType': FunctionalType,
'ProductType': ProductType},
template_parameters=[GridType, grid_layer, RangeFieldType,
dimRange, polOrder, space_backend, la_backend])
# then create the example
Example = module.add_class('PbSpe10Example', template_parameters=[GridType, space_backend, la_backend], custom_name=name)
Example.add_constructor([param('const std::string&', 'num_grid_elements'),
param('const bool', 'logging_enabled')],
throw=exceptions)
Example.add_constructor([param('const std::string&', 'num_grid_elements')],
throw=exceptions)
Example.add_constructor([],
throw=exceptions)
Example.add_method('pb_discretization_and_return_ptr',
retval(DiscretizationType + ' *', caller_owns_return=True),
[], is_const=True, throw=exceptions,
custom_name='discretization')
Example.add_method('visualize_grid',
None,
[param('const std::string&', 'filename_prefix')],
is_const=True, throw=exceptions)
Example.add_method('visualize_problem',
None,
[param('const std::string&', 'filename_prefix')],
is_const=True, throw=exceptions)
Example.add_method('visualize_problem',
None,
[param('const std::string&', 'filename_prefix'),
param('const Dune::Pymor::Parameter&', 'mu')],
is_const=True, throw=exceptions)
if space_backend == space_backend_fem:
Example.add_method('visualize_darcy_velocity',
None,
[param('const {}&'.format(VectorType), 'vector'),
param('const std::string&', 'filename'),
param('const std::string&', 'name'),
param('const Dune::Pymor::Parameter&', 'mu')],
is_const=True, throw=exceptions)
def inject_Example(module):
'''injects the user code into the module'''
namespace = module.add_cpp_namespace('Example')
AnalyticalProblem = namespace.add_class('AnalyticalProblem')
AnalyticalProblem.add_constructor([])
AnalyticalProblem.add_constructor([param('const int', 'dd')])
if __name__ == '__main__':
# prepare the module
module, pybindgen_filename, config_h_filename = prepare_python_bindings(sys.argv[1:])
# add all of libdunepymor
module, exceptions, interfaces, CONFIG_H = inject_lib_dune_pymor(module, config_h_filename)
# add example user code
inject_Example(module)
# add the users discretization
MatrixType = 'Dune::Stuff::LA::CommonDenseMatrix< double >'
VectorType = 'Dune::Stuff::LA::CommonDenseVector< double >'
OperatorType = 'Dune::Pymor::Operators::LinearAffinelyDecomposedContainerBased< ' + MatrixType + ', ' + VectorType + ' >'
discretization = inject_StationaryDiscretizationImplementation(
module, exceptions, interfaces, CONFIG_H,
'Example::SimpleDiscretization',
Traits={'VectorType': VectorType,
'OperatorType': OperatorType,
'FunctionalType': 'Dune::Pymor::Functionals::LinearAffinelyDecomposedVectorBased< ' + VectorType + ' >',
'ProductType': OperatorType})
# and add the custom constructor to the discretization
discretization.add_constructor([param('const Example::AnalyticalProblem *', 'prob', transfer_ownership=True)])
# and finally write the pybindgen .cc file
finalize_python_bindings(module, pybindgen_filename)
def inject_Example(module, exceptions, interfaces, CONFIG_H):
ssize_t = CONFIG_H['DUNE_STUFF_SSIZE_T']
'''injects the user code into the module'''
# first the discretization
GridType = 'Dune::ALUGrid< 2, 2, Dune::simplex, Dune::conforming >'
#GridType = 'Dune::SGrid< 2, 2 >'
RangeFieldType = 'double'
dimRange = '1'
polOrder = '1'
la_backend = 'Dune::Stuff::LA::ChooseBackend::eigen_sparse'
if la_backend == 'Dune::Stuff::LA::ChooseBackend::eigen_sparse':
MatrixType = 'Dune::Stuff::LA::EigenRowMajorSparseMatrix< ' + RangeFieldType + ' >'
VectorType = 'Dune::Stuff::LA::EigenDenseVector< ' + RangeFieldType + ' >'
elif la_backend == 'Dune::Stuff::LA::ChooseBackend::istl_sparse':
MatrixType = 'Dune::Stuff::LA::IstlRowMajorSparseMatrix< ' + RangeFieldType + ' >'
VectorType = 'Dune::Stuff::LA::IstlDenseVector< ' + RangeFieldType + ' >'
space_backend = 'Dune::GDT::ChooseSpaceBackend::fem'
OperatorType = 'Dune::Pymor::Operators::LinearAffinelyDecomposedContainerBased< ' + MatrixType + ', ' + VectorType + ' >'
ProductType = OperatorType
FunctionalType = 'Dune::Pymor::Functionals::LinearAffinelyDecomposedVectorBased< ' + VectorType + ' >'
DiscretizationName = 'Dune::HDD::LinearElliptic::Discretizations::BlockSWIPDG'
DiscretizationFullName = (DiscretizationName + '< ' + GridType + ', ' +
RangeFieldType + ', ' + dimRange + ', ' +
polOrder + ', ' + la_backend + '>')
LocalDiscretizationName = 'Dune::HDD::LinearElliptic::Discretizations::SWIPDG'
local_layer = 'Dune::Stuff::Grid::ChooseLayer::local'
LocalDiscretizationFullName = (LocalDiscretizationName + '< ' + GridType +
', ' + local_layer + ', ' + RangeFieldType +
', ' + dimRange + ', ' + polOrder + ', ' +
space_backend + ', ' + la_backend + ' >')
inject_StationaryDiscretizationImplementation(
module,
exceptions,
interfaces,
CONFIG_H,
LocalDiscretizationName,
Traits={
'VectorType': VectorType,
'OperatorType': OperatorType,
'FunctionalType': FunctionalType,
'ProductType': ProductType
},
template_parameters=[
GridType, local_layer, RangeFieldType, dimRange, polOrder,
space_backend, la_backend
])
oversampled_layer = 'Dune::Stuff::Grid::ChooseLayer::local_oversampled'
OversampledDiscretizationFullName = (LocalDiscretizationName + '< ' +
GridType + ', ' + oversampled_layer +
', ' + RangeFieldType + ', ' +
dimRange + ', ' + polOrder + ', ' +
space_backend + ', ' + la_backend +
' >')
inject_StationaryDiscretizationImplementation(
module,
exceptions,
interfaces,
CONFIG_H,
LocalDiscretizationName,
Traits={
'VectorType': VectorType,
'OperatorType': OperatorType,
'FunctionalType': FunctionalType,
'ProductType': ProductType
},
template_parameters=[
GridType, oversampled_layer, RangeFieldType, dimRange, polOrder,
space_backend, la_backend
])
# discretization = inject_StationaryDiscretizationImplementation(
inject_StationaryMultiscaleDiscretizationImplementation(
module,
exceptions,
interfaces,
CONFIG_H,
DiscretizationName,
Traits={
'VectorType': VectorType,
'OperatorType': OperatorType,
'FunctionalType': FunctionalType,
'ProductType': ProductType,
'LocalDiscretizationType': LocalDiscretizationFullName,
'OversampledDiscretizationType': OversampledDiscretizationFullName
},
template_parameters=[
GridType, RangeFieldType, dimRange, polOrder, la_backend
])
# then add the example
def add_example(name):
Example = module.add_class(name,
template_parameters=[GridType],
custom_name=name)
Example.add_constructor([
param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels'),
param('const ' + ssize_t, 'debug_log_levels'),
param('const bool', 'enable_warnings'),
param('const bool', 'enable_colors'),
param('const std::string', 'info_color'),
param('const std::string', 'debug_color'),
param('const std::string', 'warn_color')
],
throw=exceptions)
Example.add_constructor([
param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels'),
param('const ' + ssize_t, 'debug_log_levels'),
param('const bool', 'enable_warnings'),
param('const bool', 'enable_colors'),
param('const std::string', 'info_color'),
param('const std::string', 'debug_color')
],
throw=exceptions)
Example.add_constructor([
param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels'),
param('const ' + ssize_t, 'debug_log_levels'),
param('const bool', 'enable_warnings'),
param('const bool', 'enable_colors'),
param('const std::string', 'info_color')
],
throw=exceptions)
Example.add_constructor([
param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels'),
param('const ' + ssize_t, 'debug_log_levels'),
param('const bool', 'enable_warnings'),
param('const bool', 'enable_colors')
],
throw=exceptions)
Example.add_constructor([
param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels'),
param('const ' + ssize_t, 'debug_log_levels'),
param('const bool', 'enable_warnings')
],
throw=exceptions)
Example.add_constructor([
param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels'),
param('const ' + ssize_t, 'debug_log_levels')
],
throw=exceptions)
Example.add_constructor([
param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference'),
param('const ' + ssize_t, 'info_log_levels')
],
throw=exceptions)
Example.add_constructor([
param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products'),
param('const bool', 'with_reference')
],
throw=exceptions)
Example.add_constructor([
param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers'),
param('const std::vector< std::string >', 'products')
],
throw=exceptions)
Example.add_constructor([
param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements'),
param('const ' + ssize_t, 'oversampling_layers')
],
throw=exceptions)
Example.add_constructor([
param('const std::string', 'partitioning'),
param('const ' + ssize_t, 'num_refinements')
],
throw=exceptions)
Example.add_constructor([param('const std::string', 'partitioning')],
throw=exceptions)
Example.add_constructor([], throw=exceptions)
Example.add_method('discretization_and_return_ptr',
retval(DiscretizationFullName + ' *',
caller_owns_return=True), [],
is_const=True,
throw=exceptions,
custom_name='discretization')
Example.add_method('visualize',
None,
[param('const std::string&', 'filename_prefix')],
is_const=True,
throw=exceptions)
Example.add_method('project',
retval(VectorType),
[param('const std::string', 'expression')],
is_const=True,
throw=exceptions)
Example.add_method(
'compute_error',
retval(RangeFieldType), [
param('const ' + VectorType + '&', 'solution'),
param('const std::string', 'product_type'),
param('const Dune::Pymor::Parameter', 'mu'),
param('const Dune::Pymor::Parameter', 'mu_product')
],
is_const=True,
throw=exceptions)
Example.add_method('compute_error',
retval(RangeFieldType), [
param('const ' + VectorType + '&', 'solution'),
param('const std::string', 'product_type'),
param('const Dune::Pymor::Parameter', 'mu')
],
is_const=True,
throw=exceptions)
Example.add_method('compute_error',
retval(RangeFieldType), [
param('const ' + VectorType + '&', 'solution'),
param('const std::string', 'product_type')
],
is_const=True,
throw=exceptions)
Example.add_method('compute_error',
retval(RangeFieldType),
[param('const ' + VectorType + '&', 'solution')],
is_const=True,
throw=exceptions)
Example.add_method('available_estimators',
retval('std::vector< std::string >'), [],
is_const=True,
throw=exceptions)
Example.add_method('estimate',
retval(RangeFieldType), [
param('const ' + VectorType + '&', 'vector'),
param('const std::string', 'type'),
param('const Dune::Pymor::Parameter', 'mu_hat'),
param('const Dune::Pymor::Parameter', 'mu_bar'),
param('const Dune::Pymor::Parameter', 'mu')
],
is_const=True,
throw=exceptions)
Example.add_method('estimate',
retval(RangeFieldType), [
param('const ' + VectorType + '&', 'vector'),
param('const std::string', 'type')
],
is_const=True,
throw=exceptions)
Example.add_method('available_local_estimators',
retval('std::vector< std::string >'), [],
is_const=True,
throw=exceptions)
Example.add_method('estimate_local',
retval('std::vector< ' + RangeFieldType + ' >'), [
param('const ' + VectorType + '&', 'vector'),
param('const std::string', 'type'),
param('const Dune::Pymor::Parameter', 'mu_hat'),
param('const Dune::Pymor::Parameter', 'mu_bar'),
param('const Dune::Pymor::Parameter', 'mu')
],
is_const=True,
throw=exceptions)
Example.add_method('estimate_local',
retval('std::vector< ' + RangeFieldType + ' >'), [
param('const ' + VectorType + '&', 'vector'),
param('const std::string', 'type')
],
is_const=True,
throw=exceptions)
Example.add_method(
'pb_project_global_to_oversampled',
retval(VectorType + '*', caller_owns_return=True), [
param('const ' + VectorType + '&', 'global_vector'),
param('const ' + ssize_t, 'subdomain')
],
is_const=True,
throw=exceptions,
custom_name='project_global_to_oversampled')
Example.add_method(
'pb_project_global_to_local',
retval(VectorType + '*', caller_owns_return=True), [
param('const ' + VectorType + '&', 'global_vector'),
param('const ' + ssize_t, 'subdomain')
],
is_const=True,
throw=exceptions,
custom_name='project_global_to_local')
Example.add_method(
'pb_project_oversampled_to_local',
retval(VectorType + '*', caller_owns_return=True), [
param('const ' + VectorType + '&', 'oversampled_vector'),
param('const ' + ssize_t, 'subdomain')
],
is_const=True,
throw=exceptions,
custom_name='project_oversampled_to_local')
Example.add_method(
'solve_for_local_correction',
retval(VectorType), [
param('std::vector< ' + VectorType + ' >&', 'local_vectors'),
param('const ' + ssize_t, 'subdomain'),
param('const Dune::Pymor::Parameter', 'mu')
],
is_const=True,
throw=exceptions)
Example.add_method(
'solve_for_local_correction',
retval(VectorType), [
param('std::vector< ' + VectorType + ' >&', 'local_vectors'),
param('const ' + CONFIG_H['DUNE_STUFF_SSIZE_T'], 'subdomain')
],
is_const=True,
throw=exceptions)
Example.add_method(
'solve_oversampled',
retval(VectorType), [
param('const ' + ssize_t, 'subdomain'),
param('const std::string&', 'boundary_value_type'),
param('const ' + VectorType + '&', 'boundary_value_type'),
param('const Dune::Pymor::Parameter', 'mu')
],
is_const=True,
throw=exceptions)
Example.add_method(
'solve_oversampled',
retval(VectorType), [
param('const ' + ssize_t, 'subdomain'),
param('const std::string&', 'boundary_value_type'),
param('const ' + VectorType + '&', 'boundary_value_type')
],
is_const=True,
throw=exceptions)
Example.add_method('visualize_on_coarse_grid',
None, [
param('const std::vector< double >&', 'vector'),
param('const std::string&', 'filename'),
param('const std::string&', 'name')
],
is_const=True,
throw=exceptions)
Example.add_method('alpha',
retval(RangeFieldType), [
param('const Dune::Pymor::Parameter&', 'mu_1'),
param('const Dune::Pymor::Parameter&', 'mu_2')
],
is_const=True,
throw=exceptions)
Example.add_method('gamma',
retval(RangeFieldType), [
param('const Dune::Pymor::Parameter&', 'mu_1'),
param('const Dune::Pymor::Parameter&', 'mu_2')
],
is_const=True,
throw=exceptions)
add_example('OS2015MultiscaleExample')
add_example('OS2015AcademicExample')
def add_example(GridType, space_backend, la_backend, name):
# build all types needed for the discretization
dimRange = '1'
polOrder = '1'
grid_layer = 'Dune::Stuff::Grid::ChooseLayer::leaf'
MatrixType = 'Dune::Stuff::LA::'
VectorType = 'Dune::Stuff::LA::'
if 'eigen_sparse' in la_backend:
MatrixType += 'EigenRowMajorSparseMatrix'
VectorType += 'EigenDenseVector'
elif 'istl_sparse' in la_backend:
MatrixType += 'IstlRowMajorSparseMatrix'
VectorType += 'IstlDenseVector'
MatrixType += '< ' + RangeFieldType + ' >'
VectorType += '< ' + RangeFieldType + ' >'
OperatorType = 'Dune::Pymor::Operators::LinearAffinelyDecomposedContainerBased< ' + MatrixType + ', ' + VectorType + ' >'
ProductType = OperatorType
FunctionalType = 'Dune::Pymor::Functionals::LinearAffinelyDecomposedVectorBased< ' + VectorType + ' >'
LocalDiscretizationName = 'Dune::HDD::LinearElliptic::Discretizations::SWIPDG'
local_layer = 'Dune::Stuff::Grid::ChooseLayer::local'
LocalDiscretizationFullName = (LocalDiscretizationName + '< ' +
GridType + ', ' + local_layer + ', ' +
RangeFieldType + ', ' + dimRange +
', ' + polOrder + ', ' + space_backend +
', ' + la_backend + ' >')
inject_StationaryDiscretizationImplementation(
module,
exceptions,
interfaces,
CONFIG_H,
LocalDiscretizationName,
Traits={
'VectorType': VectorType,
'OperatorType': OperatorType,
'FunctionalType': FunctionalType,
'ProductType': ProductType
},
template_parameters=[
GridType, local_layer, RangeFieldType, dimRange, polOrder,
space_backend, la_backend
])
oversampled_layer = 'Dune::Stuff::Grid::ChooseLayer::local_oversampled'
OversampledDiscretizationFullName = (
LocalDiscretizationName + '< ' + GridType + ', ' +
oversampled_layer + ', ' + RangeFieldType + ', ' + dimRange +
', ' + polOrder + ', ' + space_backend + ', ' + la_backend + ' >')
inject_StationaryDiscretizationImplementation(
module,
exceptions,
interfaces,
CONFIG_H,
LocalDiscretizationName,
Traits={
'VectorType': VectorType,
'OperatorType': OperatorType,
'FunctionalType': FunctionalType,
'ProductType': ProductType
},
template_parameters=[
GridType, oversampled_layer, RangeFieldType, dimRange,
polOrder, space_backend, la_backend
])
DiscretizationName = 'Dune::HDD::LinearElliptic::Discretizations::BlockSWIPDG'
DiscretizationType = (
DiscretizationName + '< ' + GridType + ', '
# + GridType + ', ' + grid_layer + ', '
+ RangeFieldType + ', ' + dimRange + ', ' + polOrder + ', ' +
la_backend + '>')
inject_StationaryMultiscaleDiscretizationImplementation(
module,
exceptions,
interfaces,
CONFIG_H,
DiscretizationName,
Traits={
'VectorType': VectorType,
'OperatorType': OperatorType,
'FunctionalType': FunctionalType,
'ProductType': ProductType,
'LocalDiscretizationType': LocalDiscretizationFullName,
'OversampledDiscretizationType':
OversampledDiscretizationFullName
},
template_parameters=[
GridType, RangeFieldType, dimRange, polOrder, la_backend
])
# then create the example
Example = module.add_class(
'PbGenericLinearellipticMultiscaleExample',
template_parameters=[GridType, space_backend, la_backend],
custom_name=name)
Example.add_method('logger_options',
retval('Dune::Stuff::Common::Configuration'), [],
is_static=True,
throw=exceptions)
Example.add_method('grid_options',
retval('std::vector< std::string >'), [],
is_static=True,
throw=exceptions)
Example.add_method('grid_options',
retval('Dune::Stuff::Common::Configuration'),
[param('const std::string&', 'type')],
is_static=True,
throw=exceptions)
Example.add_method('boundary_options',
retval('std::vector< std::string >'), [],
is_static=True,
throw=exceptions)
Example.add_method('boundary_options',
retval('Dune::Stuff::Common::Configuration'),
[param('const std::string&', 'type')],
is_static=True,
throw=exceptions)
Example.add_method('problem_options',
retval('std::vector< std::string >'), [],
is_static=True,
throw=exceptions)
Example.add_method('problem_options',
retval('Dune::Stuff::Common::Configuration'),
[param('const std::string&', 'type')],
is_static=True,
throw=exceptions)
Example.add_method('solver_options',
retval('std::vector< std::string >'), [],
is_static=True,
throw=exceptions)
Example.add_method('solver_options',
retval('Dune::Stuff::Common::Configuration'),
[param('const std::string&', 'type')],
is_static=True,
throw=exceptions)
Example.add_constructor([
param('const Dune::Stuff::Common::Configuration&', 'logger_cfg'),
param('const Dune::Stuff::Common::Configuration&', 'grid_cfg'),
param('const Dune::Stuff::Common::Configuration&', 'boundary_cfg'),
param('const Dune::Stuff::Common::Configuration&', 'problem_cfg'),
param('const std::vector<std::string>&', 'only_these_products')
],
throw=exceptions)
Example.add_method('pb_discretization_and_return_ptr',
retval(DiscretizationType + ' *',
caller_owns_return=True), [],
is_const=True,
throw=exceptions,
custom_name='discretization')
Example.add_method('visualize',
None,
[param('const std::string&', 'filename_prefix')],
is_const=True,
throw=exceptions)
Example.add_method('project',
retval(VectorType),
[param('const std::string', 'expression')],
is_const=True,
throw=exceptions)
Example.add_method(
'prolong',
retval(VectorType), [
param('const ' + DiscretizationType + '&', 'source_disc'),
param('const ' + VectorType + '&', 'source_vec')
],
is_const=True,
throw=exceptions)
Example.add_method('oswald_interpolate',
retval(VectorType),
[param('const ' + VectorType + '&', 'vector')],
is_const=True,
throw=exceptions)
Example.add_method('alpha',
retval('double'), [
param('const Dune::Pymor::Parameter&', 'mu_1'),
param('const Dune::Pymor::Parameter&', 'mu_2')
],
is_const=True,
throw=exceptions)
Example.add_method('gamma',
retval('double'), [
param('const Dune::Pymor::Parameter&', 'mu_1'),
param('const Dune::Pymor::Parameter&', 'mu_2')
],
is_const=True,
throw=exceptions)
Example.add_method('min_diffusion_ev',
retval('double'),
[param('const Dune::Pymor::Parameter&', 'mu')],
is_const=True,
throw=exceptions)
Example.add_method('max_diffusion_ev',
retval('double'),
[param('const Dune::Pymor::Parameter&', 'mu')],
is_const=True,
throw=exceptions)
Example.add_method(
'elliptic_reconstruction_estimate',
retval('double'), [
param('const ' + VectorType + '&', 'p_h'),
param('const Dune::Pymor::Parameter&', 'mu_min'),
param('const Dune::Pymor::Parameter&', 'mu_max'),
param('const Dune::Pymor::Parameter&', 'mu_hat'),
param('const Dune::Pymor::Parameter&', 'mu_bar'),
param('const Dune::Pymor::Parameter&', 'mu'),
param('const std::string', 'visualize')
],
is_const=True,
throw=exceptions)
Example.add_method(
'elliptic_reconstruction_estimate',
retval('double'), [
param('const ' + VectorType + '&', 'p_h'),
param('const ' + VectorType + '&', 'w_h'),
param('const ' + VectorType + '&', 'f_h'),
param('const Dune::Pymor::Parameter&', 'mu_min'),
param('const Dune::Pymor::Parameter&', 'mu_max'),
param('const Dune::Pymor::Parameter&', 'mu_hat'),
param('const Dune::Pymor::Parameter&', 'mu_bar'),
param('const Dune::Pymor::Parameter&', 'mu'),
param('const std::string', 'visualize')
],
is_const=True,
throw=exceptions)
Example.add_method('max_grid_width',
retval('double'), [],
is_const=True,
throw=exceptions)
Example.add_method('domain_diameter',
retval('double'), [],
is_const=True,
throw=exceptions)