> >> One could set up the link to determine the jump between the opposite faces >> using a link-topology matrix in a similar way that you determine the >> matching dofs for periodic boundary constraints, but this seems like a lot >> of work. > > Right. You can't ask the cell for its neighbor -- it will just say it's at > the > boundary. > >> Wolfgang's suggestion of using hp_collector and fe_nothing combinations >> makes sense as you also have double dof at the interface. How would one >> determine the jumps across the interface in this approach? > > Well, on the other side you've got a FENothing which is constant zero. So the > jump equals the value on this side. In essence, using a FE_Q in one subdomain > and a FE_Nothing in the other yields a discontinuous trial space. You would > evaluate the jump in the same way as you do for DG spaces.
I think my question is slightly different to Thomas's. I have the same fields on either side of the interface but they are discontinuous over the interface but continuous on either side. So I would like the interface to be DG and away from that I would use CG. (again, apologies if I'm hijacking the original post) > > W. > > ------------------------------------------------------------------------- > Wolfgang Bangerth email: [email protected] > www: http://www.math.tamu.edu/~bangerth/ > _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
