Dharmendar Reddy <dharmareddy84 at gmail.com> writes: > Hello, > I am seeing nan's at boundary locations for a plot over line data > on a 2D mesh. There are no nan's in the data obtained from solver. > > For example, When i click plot over line along y-axis in paraview, it > selects the y-coordiantes to be from > > -21.141444145 to 26.426805182 , The boundary Node at the positive end is > showing a nan. If i change the end node to 26.426804 the nan goes away > and i get the expected value (which is zero, as per imposed boundary > condition). > > I have not seen this before with the ascii vtk writer in my fortran code. > > I am wondering if there is something i should be aware of petsc binary vtk > viewer. I run the code on a 64 bit linux cluster. I am visualizing the vtk > file on a 32 bit paraview 3.10.1 on windows xp.
32/64-bit values shouldn't matter, but it's possible there is a bug with 64-bit-indices (I don't remember testing it). Is the writing code valgrind-clean? Is it reproducible?
