Hallo Markus, I have already used SolutionTransfer during the forward model simulation to get grid2 from grid1 (coarsening & refinement), and to transfer the forward solution from grid1 to grid2. Now, in the discrete adjoint simulation, I need to transfer the discrete adjoint solution from grid2 to grid1 (reverse order), where grid1 and grid2 are fixed (they have been saved to a file during the forward model computation). It's a relatively simple inverse problem setup.
SolutionTransfer needs me to explicitly coarsen and refine grid2 to get to grid1... that's why I thought I cannot use it here. Can this class help me? Perhaps if I somehow record the coarsening/refining cell-wise operations that got me grid2 from grid1, and somehow mark each cell in grid2 for coarsening/refinement s.t. i get back to grid1 exactly... can that be done? Thanks, and best regards, Mihai ________________________________ Von: Markus Bürg <[email protected]> An: [email protected] Gesendet: Donnerstag, den 18. Februar 2010, 11:26:45 Uhr Betreff: Re: [deal.II] transfering solution between two different grids Hello Mihai, I think the class SolutionTransfer is what you want. Best Regards, Markus Am 18.02.10 17:15, schrieb mihai alexe: > >Hello all, > > >How can I transfer a discrete solution from one grid to another? >My two grids are fixed - I am reading them from a file, and grid #2 >holds my discontinuous Galerkin solution u^h. Grid #2 has been obtained >through refinement and coarsening directly from grid #1. I want to >transfer u^h back from grid #2 to grid #1. Is there an API function for >doing this? I came upon static >void VectorTools::interpolate, >and I am wondering how I can set up my "transfer" matrix beforehand >such that >I >successfully do the solution transfer. > > > > >Thank >you >~~ >Mihai > >template<int >dim, class InVector , class OutVector , int spacedim> > >static void >VectorTools::interpolate ( const DoFHandler< >dim, spacedim > & dof_1, > > > const DoFHandler< >dim, spacedim > & dof_2, > > > const FullMatrix< double > >& transfer, > > > const >InVector & data_1, > > > OutVector >& data_2 > > ) > > > > > >>__________________________________________________ >>Do You Yahoo!? >>Sie sind Spam leid? Yahoo! Mail verfügt über einen herausragenden >Schutz gegen Massenmails. >http://mail.yahoo.com > >_______________________________________________ >dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii __________________________________________________ Do You Yahoo!? Sie sind Spam leid? Yahoo! Mail verfügt über einen herausragenden Schutz gegen Massenmails. http://mail.yahoo.com
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
