FEMSystem users, I've attached a patch that fixes a problem when numerical Jacobians are forced on the user side and verify_analytic_jacobians = 0.0. In particular, elem_jacobian got zeroed without the stored residual getting added back in during the side loop. A private review by Roy also picked up the case for when element Jacobians are computed, but side Jacobians are not. This patch should fix both cases. I also updated the adjoint examples to test this. In particular, I set verify_analytic_jacobians = 0.0 and analytic_jacobians = false in adjoints_ex1 and then enabled verify_analtyic_jacobians in adjoints_ex2 to still exercise that capability. However, none of the examples currently exercise the case where element Jacobians are analytically computed while the side Jacobians are not.
Feedback before committing? Best, Paul
femsystem.patch
Description: Binary data
------------------------------------------------------------------------------ Live Security Virtual Conference Exclusive live event will cover all the ways today's security and threat landscape has changed and how IT managers can respond. Discussions will include endpoint security, mobile security and the latest in malware threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
_______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
