More wierd behaviour from MatIsSymmetric?

2007-10-17 Thread Toby D. Young
I am having a little trouble, or confusion, with the MatIsSymmetric function. I am getting the reverse behaviour from this routine that I expect. I am convinced I have a symmetric matrix, but the answer I get from this function is always zero - ie, PetscTruth = PetscFalse. Why, if PETSC_YES is a

Non repeatability issue

2007-10-17 Thread Aldo Bonfiglioli
crubbed... Name: README.nonrepeatibility URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20071017/26ed5232/attachment.diff>

More wierd behaviour from MatIsSymmetric?

2007-10-17 Thread Matthew Knepley
Is there a reason for not using 2.3.3? I suspect this to be a problem with finite precision, which is corrected in the newest release.' Thanks, Matt On 10/17/07, Toby D. Young wrote: > > > I am having a little trouble, or confusion, with the MatIsSymmetric > function. I am getting the rev

Non repeatability issue

2007-10-17 Thread Barry Smith
The 2d results, as noted, are reasonable and expected behavior. The 3d results are not acceptable, when sometimes it converges and sometimes it does not. I would use a 1e-12 KSP rtol and a similar 1e-12 SNES rtol; does it now always converge? If not then I suspect there may be some error in t