Thanks for replying,the shape function is different on each cell.
But how can I make a Triangulation/DoFhandler object from the reference cell?
Do I have to use the
create_triangulation() function?
> From: [email protected]
> To: [email protected]
> Subject: Re: [deal.II] Solving a local problem
> Date: Fri, 16 Jul 2010 07:17:04 -0500
> CC: [email protected]
>
>
> Travis,
>
> > Sorry for the inconvenience but I still cannot figure it out.The
> > algorithm my professor is working on is an implementation of a multiscale
> > finite element method.The key step of the multiscale finite element method
> > is to construct the basis function in each element. To construct the basis
> > function, we need to solve a local problem in each element. That is to say,
> > when we carry on the finite element method, the basis function in each
> > element is obtained by solving the local problem. I have a triangulation
> > object and it has been refined.Now I only need to solve an equation on one
> > of its cells and the solution will be regarded as the basis function of an
> > element in that cell. Is there an easier way other than associating
> > FE_Nothing elements to other cells?I wonder if active cell_iterator could
> > be somhow transformed into a triangulation object such that I can solve the
> > equation using the
> > techniques step-6 provide with only changes in the Coeffeicient
> > function.
>
> The question is whether you need to solve for different shape functions on
> each cell, or whether the local problem you need to solve is the same for all
> cells?
>
> In either case, the programs I know of where people have tried this before
> use
> two Triangulation/DoFHandler objects. One describes the reference cell and
> fine scale problem, and the other one describes the real problem and uses the
> shape functions computed on the reference cell.
>
> Best
> W.
>
> -------------------------------------------------------------------------
> Wolfgang Bangerth email: [email protected]
> www: http://www.math.tamu.edu/~bangerth/
>
_________________________________________________________________
约会说不清地方？来试试微软地图最新msn互动功能！
http://ditu.live.com/?form=TL&swm=1_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
