On 02/17/2017 09:14 PM, Fred wrote:
I currently have a few codes written in the meshworker framework that solve
systems of PDE's. (Elasticity, Laplace/Poisson, Stokes) I am interested in
experimenting with Multigrid preconditioner for these. I have been unable to
get this to work using the
MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks function. I am
interested in seeing an example of a code where this strategy is used. If
there is any interest in seeing my code, I will be happy to produce a minimal
non-working example. For now, the compilation error that I get from clang is
as follows (g++ is similar).
git/dealiiBase/dealiiInstall/include/deal.II/lac/matrix_block.h:932:10: error:
no viable conversion from returned value of type 'const
MGLevelObject<dealii::SparseMatrix<double>>' to function return type
'const
MGLevelObject<MatrixBlock<dealii::SparseMatrix<double> >>'
return *matrices.read<const MGLevelObject<MatrixType>* >(i);
The line that triggers it is the following
dealii::MeshWorker::integration_loop<dim,dim>
( dof_handler.begin_mg(),
dof_handler.end_mg(),
dof_Info,
info_box,
LDGintegrator,
assembleSystem);
LDGintegrator is my custom integrator of a type inheriated from
MeshWorker::LocalIntegrator<dim>;
and assembleSystem is of type
MeshWorker
::Assembler
::MGMatrixLocalBlocksToGlobalBlock< SparseMatrix<double>,
double>;
Any help would be greatly appreciated. I would be happy to contribute any
changes required to get this to work.
Can you open an issue on the github forum at
https://github.com/dealii/dealii/issues
It would be great to have a minimal example that shows the problem.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.