### Problem with visualization of Results in PARAVIEW

24
views
0
10 days ago by
Hi
I have problem with visualizing the results in Paraview. I want to solve an equation in a part of my domain (NOT the whole domain).
F = dot(grad(PHI), grad(v_4)) * dx(2)+...

n = File("PHI.pvd")

while t <= T:

J = derivative(F, z, dz)
problem = NonlinearVariationalProblem(F, z, bcs, J)
solver = NonlinearVariationalSolver(problem)
solver.solve()
t += dt

n << PHI​
The problem is when I want to see the results in Paraview I get this error:

ERROR: In /build/paraview-arsa8T/paraview-5.0.1+dfsg1/VTK/IO/XML/vtkXMLDataReader.cxx, line 449
vtkXMLUnstructuredGridReader (0x424c6b0): Cannot read point data array "f_7-3" from PointData in piece 0.  The data array in the element may be too short.​

The thing is when I solve the same equation in the whole domain (Using dx instead of dx(2)) the file is opened perfectly! Does anybody know of any solution for this weird problem?

Community: FEniCS Project