Hello- I'm attempting to do a refinement based on an exact solution. My system has 3 variables in it.
I believe that the exact_error_estimator.estimate_error function has a bug in it, or at least I don't understand the current idea of the function. It has an outer loop for the variables and an inner loop for the elements. For each variable it seems to do as expected. However, the error_per_cell is essentially getting overwritten for each new variable. Thus, the returned error vector corresponds only to the last variable. For my case, the third variable is exactly interpolated on the mesh, and so this results in uniform refinement. I see that I can put the weight to zero for this variable, but it seems quite restrictive to assume only the error for only one variable should be considered. Or is this really the core idea? Thanks! Andrea ------------------------------------------------------------------------------ Download Intel® Parallel Studio Eval Try the new software tools for yourself. Speed compiling, find bugs proactively, and fine-tune applications for parallel performance. See why Intel Parallel Studio got high marks during beta. http://p.sf.net/sfu/intel-sw-dev _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
