> Moreover --I probably wasn't clear--, I'm solving the problem using *two* > meshes. First is a small rectangle, and second a large rectangle with a > hole (at the size of the first rectangle) in its center. > What I thought I should do is create two classes similar to LaplaceProblem > in Tutorial 3, one for each of the problems. Then call each of them to > solve their > respective problem (by first feeding to them the new boundary conditions), > and when they're done take the gradient from each of the solutions and > update the boundary conditions.
I think this is quite a workable approach. You'll need to use the FEFaceValues::get_function_gradients() function to extract gradients on the boundary of one mesh and feed it to the other. Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
