# possible bug in vtkRectilinearGrid

**URL:** https://discourse.paraview.org/t/possible-bug-in-vtkrectilineargrid/3451
**Category:** In Situ Support
**Created:** [January 29, 2020, 5:22pm UTC](https://discourse.paraview.org/t/possible-bug-in-vtkrectilineargrid/3451 "2020-01-29T17:22:49Z")
**Posts on this page:** 3
**Page:** 1

<div class="post-metadata">

### Author: ![Rigel](https://discourse.paraview.org/letter_avatar_proxy/v4/letter/r/8e7dd6/32.png) [@Rigel](https://discourse.paraview.org/u/Rigel)
#### Post date: [January 29, 2020, 5:22pm UTC](https://discourse.paraview.org/t/possible-bug-in-vtkrectilineargrid/3451/1 "2020-01-29T17:22:49Z")

</div>

Hi,

I have a simulation which runs on a simple, rectilinear grid (with 8 points in each Cartesian direction). When running in serial, Catalyst outputs the correct information, which can then be visualized in ParaView out of the box:

 ![image](https://discourse.paraview.org/uploads/default/original/2X/b/bf03dd9ab99d8948175d9782b336d38745f9e3ba.png)

[serial.pvsm](https://discourse.paraview.org/uploads/short-url/dBCXVVNsAmBBxG2ckRU8ToIaFId.pvsm) (166.1 KB) is the correspondent ParaView state file, whose data is [serial.zip](https://discourse.paraview.org/uploads/short-url/cEoXJXV8yGsqQRR3sQEbw2YY9ZC.zip) (13.2 KB).

However, when running in parallel (with 2 ranks per Cartesian direction, 8 in total), Catalyst seems to output the correct grid data, but when attempting to open the `pvtr` file series in ParaView (at the moment I am only doing post-mortem analysis), only the part of the grid pertaining to rank 0 is shown:

 ![image](https://discourse.paraview.org/uploads/default/original/2X/d/da1b3008c0bbc7e2206658e43a6b863f1ad32738.png)

[parallel.pvsm](https://discourse.paraview.org/uploads/short-url/1VvLSiprVL61rY0xC1tAyUp3mzc.pvsm) (146.8 KB) is the correspondent ParaView state file, whose data is [parallel.zip](https://discourse.paraview.org/uploads/short-url/dDNmcQeMYTJSaH10xfz0tiDVMvl.zip) (61.0 KB). If I take an arbitrary time-step folder and open separately the `vtr` files, I can see the remaining parts of the grid are there:

 ![separated](https://discourse.paraview.org/uploads/default/original/2X/c/cf74e395b419984a188d5247a4d193b9e9c0f24f.png)

In my simulation, all fields are point-based, so having an ‘empty’ cell layer between each subdomain should not be a problem.

Finally, trying to send to Catalyst also the coordinates of the halo points did not help: only the subdomain pertaining to rank 0 is shown (when opening the `pvtr` series), but the remaining ones are actually defined in the respective `vtr` files:

 ![separated](https://discourse.paraview.org/uploads/default/original/2X/3/3234796753698788058ed87a1bc64999fa5263da.png)

[separated.pvsm](https://discourse.paraview.org/uploads/short-url/pWqEsPkMi8vU3di3r4kDIk2Ne6T.pvsm) (780.4 KB) is the state file for this separate visualization; the data is [halos.zip](https://discourse.paraview.org/uploads/short-url/78pdz3qDQKqmB39PaWVUMwzQgNg.zip) (64.0 KB).

Does anyone spot an error on my implementation? If not, then there might be a bug within the vtkRectilinearGrid class.

Thank you very much,

---

<div class="post-metadata">

### Author: ![Andy\_Bauer](https://discourse.paraview.org/user_avatar/discourse.paraview.org/andy_bauer/32/5442_2.png) [@Andy\_Bauer](https://discourse.paraview.org/u/Andy_Bauer)
#### Post date: [January 29, 2020, 5:59pm UTC](https://discourse.paraview.org/t/possible-bug-in-vtkrectilineargrid/3451/2 "2020-01-29T17:59:15Z")

</div>

Looking at your grid\_2.pvtr file from your halos.zip your extents are incorrect. Check out Figure 3.6 in the Catalyst User’s Guide Version 5.8 for how to set it for each process and then also make sure to also set your vtkCPInputDataDescription::SetWholeExtent(). Additionally, even if you only have point data you still won’t want to have an ‘empty’ cell layer between each subdomain because certain filters will require that cell to exist.

---

<div class="post-metadata">

### Author: ![Rigel](https://discourse.paraview.org/letter_avatar_proxy/v4/letter/r/8e7dd6/32.png) [@Rigel](https://discourse.paraview.org/u/Rigel)
#### Post date: [January 30, 2020, 2:24pm UTC](https://discourse.paraview.org/t/possible-bug-in-vtkrectilineargrid/3451/3 "2020-01-30T14:24:10Z")

</div>

Thank you very much, Andy. You are right: I did as you mentioned and the issue has been solved.
