Hello,
I am trying to render a non-uniform rectilinear grid with catalyst.
The code I use to generate the grid is:
mesh["coordsets/coords/type"] = "rectilinear";
mesh["coordsets/coords/values/x"] = grid.g_x;
mesh["coordsets/coords/values/y"] = grid.g_y;
mesh["coordsets/coords/values/z"] = grid.g_z;
mesh["topologies/mesh/coordset"] = "coords";
mesh["topologies/mesh/type"] = "rectilinear";
mesh["topologies/mesh/elements/origin/i"] = grid.g_x[0];
mesh["topologies/mesh/elements/origin/j"] = grid.g_y[0];
mesh["topologies/mesh/elements/origin/k"] = grid.g_z[0];
calling mesh[“coordsets/coords/values”].print(); gives the following:
{
"x":
{
"dtype":"float64",
"number_of_elements": 31,
"offset": 0,
"stride": 8,
"element_bytes": 8,
"endianness": "little",
"value": [0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0, 21.0, 22.0, 23.0, 24.0, 25.0, 26.0, 27.0, 28.0, 29.0, 30.0]
},
"y":
{
"dtype":"float64",
"number_of_elements": 8,
"offset": 0,
"stride": 8,
"element_bytes": 8,
"endianness": "little",
"value": [-1.0, -0.6, 0.0, 1.0, 2.0, 3.0, 3.6, 4.0]
},
"z":
{
"dtype":"float64",
"number_of_elements": 8,
"offset": 0,
"stride": 8,
"element_bytes": 8,
"endianness": "little",
"value": [-1.0, -0.6, 0.0, 1.0, 2.0, 3.0, 3.6, 4.0]
}
}
This for some reason results in a uniform grid in Paraview, with an origin at {0,0,0}.
How can I properly render the non-uniform rectilinear grid?
The file I use for catalyst is the one from Lulesh:
savedata.py (578 Bytes)
Thanks in advance!