> 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. W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
