On Tue, 6 Jan 2015, Kameeko Kiwi wrote: > I've not yet been able to replicate the phenomenon with a single-variable > system, though I am not sure whether that is a coincidence or related to > the problem.
I can replicate your results: ~*~*~*~*~*~*~*~*~ adjoint solve start ~*~*~*~*~*~*~*~*~ number of iterations to solve adjoint: 10 final residual of adjoint solve: 1.58046e-13 ~*~*~*~*~*~*~*~*~ adjoint solve end ~*~*~*~*~*~*~*~*~ ------------ herp derp ------------ adjoint system residual (discrete L2): 0.000308402 adjoint system residual (L2, all): 8.17061e-05 adjoint system residual (L2, 0): 1.07364e-14 adjoint system residual (L2, 1): 8.17061e-05 But I've still no idea what's causing them. You're sure your system is invertible? I get a zero pivot if I try "-pc_type lu", but that may just be because PETSc only does partial pivoting. --- Roy ------------------------------------------------------------------------------ Dive into the World of Parallel Programming! The Go Parallel Website, sponsored by Intel and developed in partnership with Slashdot Media, is your hub for all things parallel software development, from weekly thought leadership blogs to news, videos, case studies, tutorials and more. Take a look and join the conversation now. http://goparallel.sourceforge.net _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
