> > But the question is: do you use *2* DoFHandler objects for your two > > finite elements, or do you use one that you re-initialize with a > > different finite element between steps 1 and 2? > > One dof_handler that I reinitialize. Is that likely to complcate things? > Do you think I should be using two dof_handlers?
You will have to. The DoFHandler object stores which degree of freedom lives where, and so makes the association between the elements of your solution vector and which shape function/node they belong to. If you re-initialize your DoFHandler object with another finite element, you forget this connection, and in effect you end up with a solution vector from the first step for which you don't know any more what vector elements correspond to the shape functions on your current cell. step-18 gets around this because it stores stresses with the individual cell, and doesn't use a DoFHandler for this in the first place. But that doesn't seem to be the case for you. Does that explain things a bit? Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
