[dune-pdelab] Periodic Yasp-grid

Christian Engwer christi at uni-hd.de
Sat Sep 4 22:20:16 CEST 2010


Hi,

the peridocity in yaspgrid is build upon the parallel interface. You
have a grid with overlap and the overlap is wrpped around. You then
need to run a parallel program to use the periodicity.

Judging from your code, I assume that you didn't do the above. As
PDELab is able to do parallel computations, it should also be able to
do periodic computations. I hope these changes will solve your
problem.

Cheers
Christian

On Sat, Sep 04, 2010 at 04:00:27PM +0200, Christoph GrĂ¼ninger wrote:
> Ladies and Gentlemen,
> I would like to solve an equation with a periodic one-dimensional grid. I have not found any example using a periodic grid. Using Yasp-grid I tried
> 
>       const int dim = 1;
>       Dune::FieldVector<double,dim> L(1.0);
>       Dune::FieldVector<int,dim> gridsize(1);
>       Dune::FieldVector<bool,dim> periodic(true);
>       Dune::YaspGrid<dim> grid(L,gridsize,periodic,0);
>       grid.globalRefine(5);
> 
> with PDELab. But alpha_boundary from my local operator is called anyway. PDELab seems to ignore the periodicity of my grid. Have I overseen something or is it Yasp-grid / PDELab?
> 
> Thanks in advance,
> Christoph 
> 
> _______________________________________________
> dune-pdelab mailing list
> dune-pdelab at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune-pdelab
> 




More information about the dune-pdelab mailing list