Hi,

I am preparing an example, and now there is better than even chance I 
figured it out.

I will put up the example on github only if I can't figure it out.

I really appreciate your help!

On Friday, February 17, 2017 at 11:51:07 PM UTC-5, Wolfgang Bangerth wrote:
>
> 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:                 bang...@colostate.edu 
> <javascript:> 
>                             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.

Reply via email to