Hi all,
I have this small patch to avoid inverse_map() in nodal_solution(...) and 
change the output of ex21 from gmv to ExodusII.
I applied the change only to MONOMIALS but it should work with every reference 
frame basis function.
Tell me if you are interested and I'll do it.

I also noticed that the new L2_HIERARCHIC basis doesn't work with ex21. Is 
there any reason why?
Thanks
Lorenzo

Attachment: patch
Description: Binary data



On Nov 9, 2011, at 4:46 AM, Derek Gaston wrote:

> Holy crap.... that is scary as hell!  We'll see if it modifies our
> regression tests at all...
> 
> Derek
> 
> Sent from my iPhone
> 
> On Nov 8, 2011, at 7:16 PM, Roy Stogner <[email protected]> wrote:
> 
>> 
>> On Tue, 8 Nov 2011, Kirk, Benjamin (JSC-EG311) wrote:
>> 
>>> None here. I was pretty much sequestered today - what's the gist?
>> 
>> The short version is that it's a miracle I manage to dress myself and
>> feed myself despite gross mental handicaps.
>> 
>> The long version is that my "has_affine_map" optimization, which
>> obtains the correct init_shape_functions result much faster on
>> triangles and tets, obtains incorrect dphidxi etc. maps on rectangluar
>> quads, rectangular prism hexes, and rectangular prisms.  As far as I
>> can tell, the only reason we managed to avoid noticing the problem for
>> so long is that for most codes the errors can act as a self-consistent
>> weighting function in the weak formulation integrals, so e.g.
>> manufactured solution tests still converge at the right order and
>> exact solutions in the FEM space are still recovered.  After the
>> change all our FIN-S regression test results are still coming in under
>> the 1e-15 tolerance, even!
>> 
>> I think Boyce only managed to catch things by using different order
>> quadrature rules to integrate the matrix vs forcing function for an L2
>> projection.
>> ---
>> Roy
>> 
>> ------------------------------------------------------------------------------
>> RSA(R) Conference 2012
>> Save $700 by Nov 18
>> Register now
>> http://p.sf.net/sfu/rsa-sfdev2dev1
>> _______________________________________________
>> Libmesh-devel mailing list
>> [email protected]
>> https://lists.sourceforge.net/lists/listinfo/libmesh-devel
> 
> ------------------------------------------------------------------------------
> RSA(R) Conference 2012
> Save $700 by Nov 18
> Register now
> http://p.sf.net/sfu/rsa-sfdev2dev1
> _______________________________________________
> Libmesh-devel mailing list
> [email protected]
> https://lists.sourceforge.net/lists/listinfo/libmesh-devel

------------------------------------------------------------------------------
RSA(R) Conference 2012
Save $700 by Nov 18
Register now
http://p.sf.net/sfu/rsa-sfdev2dev1
_______________________________________________
Libmesh-devel mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-devel

Reply via email to