Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-11-24 Thread Ronghong Fan
Hi Konrad, It's OK. I will follow your research. Thanks again anyway. Cheers, Ronghong On Wed, Nov 21, 2018 at 1:53 AM Konrad wrote: > Hi Ronghong, > > I have a version of the MsFEM for systems but this project is not > published yet so I hesitate to share code at this stage (also the code is

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-11-20 Thread Konrad
Hi Ronghong, I have a version of the MsFEM for systems but this project is not published yet so I hesitate to share code at this stage (also the code is not very clean yet). I will be happy though to guide you with the code design and even with some helpful code snippets. If you are interested

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-11-19 Thread Ronghong Fan
Dear Konrad, Thanks for your reply. Now I just want to solve two separate scales, which is just a start. I read the book "Multiscale Finite Element Methods. Theory and Applications by Thomas Hou" written by Yalchin Efendiev and Thomas Y. Hou. Maybe you have read it. I know this multiscale finite

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-11-16 Thread Konrad
Dear Ronghong Fan, Yes, I implemented several versions. The question I asked above got resolved as explained above. If you have a specific question feel free to ask, I will help if I can. Keep in mind that Multiscale FEMs are not generic, only work in certain (scaling) scenarios and do not

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-11-13 Thread Ronghong Fan
Dear Konrad, Do you have some progress about Multiscale Finite Element Method with dealii? I am very interested in mfem too and want to use dealii to implement. On Thursday, August 30, 2018 at 7:27:21 PM UTC+8, Konrad wrote: > > On Wednesday, August 29, 2018 at 5:25:29 PM UTC+2, Wolfgang

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-08-30 Thread Konrad
On Wednesday, August 29, 2018 at 5:25:29 PM UTC+2, Wolfgang Bangerth wrote: > > > I think that's actually an important point you are making here, maybe > inadvertently: You have as many *meshes* as there are coarse cells. > > Now, just refining the global mesh a number of more times gives you >

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-08-29 Thread Wolfgang Bangerth
Many thanks, Wolfgang. Yes, I need to solve for the basis in each coarse cell separately, i.e., I have as many local fine meshes as coarse cells. I think that's actually an important point you are making here, maybe inadvertently: You have as many *meshes* as there are coarse cells. Now,

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-08-29 Thread Konrad
On Wednesday, August 29, 2018 at 3:51:35 PM UTC+2, Jean-Paul Pelteret wrote: > > Hi Kendra, > > > Thank you, Jean-Paul, for the hint. Unfortunately this is not really > what I am looking for. I found a solution (that I don't like) by iterating > through two meshes on different levels

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-08-29 Thread Jean-Paul Pelteret
Hi Kendra, > Thank you, Jean-Paul, for the hint. Unfortunately this is not really what I > am looking for. I found a solution (that I don't like) by iterating through > two meshes on different levels simultaneously. Once I find a good solution I > will post it here and link to the code. I

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-08-29 Thread Konrad
On Wednesday, August 29, 2018 at 11:50:57 AM UTC+2, Jean-Paul Pelteret wrote: > > Hi Konrad, > > Do you think it would be possible to (ab)use the multigrid classes for > that? As far as I understand there are dofhandlers on each level of the > mesh. Can I get dofhandlers for all active subcells

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-08-29 Thread Jean-Paul Pelteret
Hi Konrad, > Do you think it would be possible to (ab)use the multigrid classes for that? > As far as I understand there are dofhandlers on each level of the mesh. Can I > get dofhandlers for all active subcells of say coarse cell number n in level > 5? The DoFHandler::begin_active()

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-08-29 Thread Konrad
On Wednesday, August 29, 2018 at 4:23:32 AM UTC+2, Wolfgang Bangerth wrote: > > Since you really only want to interpolate between two meshes, I would just > have two DoFHandlers and two Triangulations and use the class InterGridMap > to > translate between the two of them. > > Or do you want

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-08-28 Thread Wolfgang Bangerth
On 08/28/2018 04:47 AM, Konrad wrote: Is it possible to have system matrices, dofs and solution vectors etc on two different refinement levels of a mesh? Say my coarse FEM holds data on refinement level 5 and my basis function(s) on each coarse cell live on the leaves of each coarse cell on