Hi everybody I started experimenting with mesh adaptation in order to apply it in my model. Some time ago I asked about the way that Neumann boundaries should be indicated in the estimate function of the error estimation class.I was getting some run time exceptions with the way I was implementing it. I just wanted to ask the following:
Assume that we have a class that describes these Neumann b.c say class Nbc, which inherits from Function<dim>, I did the following in the refine_mesh routine of the driver code Nbc traction; Nbc* pt=&traction; typename FunctionMap<dim>::type Neumann_boundary say the Neumann bc are imposed at face with boundary indicator 0 Neumann_boundary['0']=pt; And then pass this FunctionMap to the estimate function. As a result I did not get any run time exceptions and the code run normally although as far as I can understand there seems to be some problems with the convergence of the solver.Anyway, is this the right way that I am describing above that the implementation of Neumann bc should be implemented with mesh adaptation or am I doing something really stupid? All the best Alex This message has been checked for viruses but the contents of an attachment may still contain software viruses, which could damage your computer system: you are advised to perform your own checks. Email communications with the University of Nottingham may be monitored as permitted by UK legislation.
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
