crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 05:35:49 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Add support for curvilinear elements when exporting cylindrical and spherical meshes to VTK (#2533)

Browse source 
Co-authored-by: Paul Romano <paul.k.romano@gmail.com>
This commit is contained in:
Patrick Shriwise 2023-06-16 15:17:33 -05:00 • committed by GitHub
parent 96d150a64e
commit 0fcb174b2a
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23
12 changed files with 4012 additions and 227 deletions
Download patch file Download diff file Expand all files Collapse all files

486
openmc/mesh.py
View file

@ -177,7 +177,65 @@ class StructuredMesh(MeshBase):
unpacked along the first dimension with xx, yy, zz = mesh.vertices.
"""
return np.stack(np.meshgrid(*self._grids, indexing='ij'), axis=0)
return self._generate_vertices(*self._grids)
@staticmethod
def _generate_vertices(i_grid, j_grid, k_grid):
"""Returns an array with shape (3, i_grid.size+1, j_grid.size+1, k_grid.size+1)
containing the corner vertices of mesh elements.
"""
return np.stack(np.meshgrid(i_grid, j_grid, k_grid, indexing='ij'), axis=0)
@staticmethod
def _generate_edge_midpoints(grids):
"""Generates the midpoints of mesh element edges for each dimension of the mesh.
Parameters
----------
grids : numpy.ndarray
The vertex grids along each dimension of the mesh.
Returns
-------
midpoint_grids : list of numpy.ndarray
The edge midpoints for the i, j, and k midpoints of each element in
i, j, k ordering. The shapes of the resulting grids are
[(3, ni-1, nj, nk), (3, ni, nj-1, nk), (3, ni, nj, nk-1)]
"""
# generate a set of edge midpoints for each dimension
midpoint_grids = []
# generate the element edge midpoints in order s.t.
# the epxected element ordering is preserved with respect to the corner vertices
# each grid is comprised of the mid points for one dimension and the
# corner vertices of the other two
for dims in ((0, 1, 2), (1, 0, 2), (2, 0, 1)):
# compute the midpoints along the first dimension
midpoints = grids[dims[0]][:-1] + 0.5 * np.diff(grids[dims[0]])
coords = (midpoints, grids[dims[1]], grids[dims[2]])
i_grid, j_grid, k_grid = [coords[dims.index(i)] for i in range(3)]
# re-use the generate vertices method to create the full mesh grid
# transpose to get (i, j, k) ordering of the gridpoints
midpoint_grid = StructuredMesh._generate_vertices(i_grid, j_grid, k_grid)
midpoint_grids.append(midpoint_grid)
return midpoint_grids
@property
def midpoint_vertices(self):
"""Create vertices that lie on the midpoint of element edges
"""
# generate edge midpoints needed for curvilinear element definition
midpoint_vertices = self._generate_edge_midpoints(self._grids)
# convert each of the midpoint grids to cartesian coordinates
for vertices in midpoint_vertices:
self._convert_to_cartesian(vertices, self.origin)
return midpoint_vertices
@property
def centroids(self):
@ -202,13 +260,15 @@ class StructuredMesh(MeshBase):
def num_mesh_cells(self):
return np.prod(self.dimension)
def write_data_to_vtk(self, points, filename, datasets, volume_normalization=True):
def write_data_to_vtk(self,
filename,
datasets=None,
volume_normalization=True,
curvilinear=False):
"""Creates a VTK object of the mesh
Parameters
----------
points : list or np.array
List of (X,Y,Z) tuples.
filename : str
Name of the VTK file to write.
datasets : dict
@ -217,6 +277,9 @@ class StructuredMesh(MeshBase):
volume_normalization : bool, optional
Whether or not to normalize the data by
the volume of the mesh elements.
curvilinear : bool
Whether or not to write curvilinear elements. Only applies to
``SphericalMesh`` and ``CylindricalMesh``.
Raises
------
@ -225,14 +288,180 @@ class StructuredMesh(MeshBase):
Returns
-------
vtk.vtkStructuredGrid
the VTK object
vtk.StructuredGrid or vtk.UnstructuredGrid
a VTK grid object representing the mesh
"""
import vtk
from vtk.util import numpy_support as nps
# check that the data sets are appropriately sized
if datasets is not None:
self._check_vtk_datasets(datasets)
# write linear elements using a structured grid
if not curvilinear or isinstance(self, (RegularMesh, RectilinearMesh)):
vtk_grid = self._create_vtk_structured_grid()
writer = vtk.vtkStructuredGridWriter()
# write curvilinear elements using an unstructured grid
else:
vtk_grid = self._create_vtk_unstructured_grid()
writer = vtk.vtkUnstructuredGridWriter()
if datasets is not None:
# maintain a list of the datasets as added
# to the VTK arrays to ensure they persist
# in memory until the file is written
datasets_out = []
for label, dataset in datasets.items():
dataset = np.asarray(dataset).flatten()
datasets_out.append(dataset)
if volume_normalization:
dataset /= self.volumes.T.flatten()
dataset_array = vtk.vtkDoubleArray()
dataset_array.SetName(label)
dataset_array.SetArray(nps.numpy_to_vtk(dataset),
dataset.size,
True)
vtk_grid.GetCellData().AddArray(dataset_array)
writer.SetFileName(str(filename))
writer.SetInputData(vtk_grid)
writer.Write()
return vtk_grid
def _create_vtk_structured_grid(self):
"""Create a structured grid
Returns
-------
vtk.vtkStructuredGrid
a VTK structured grid object representing the mesh
"""
import vtk
from vtk.util import numpy_support as nps
vertices = self.cartesian_vertices.T.reshape(-1, 3)
vtkPts = vtk.vtkPoints()
vtkPts.SetData(nps.numpy_to_vtk(vertices, deep=True))
vtk_grid = vtk.vtkStructuredGrid()
vtk_grid.SetPoints(vtkPts)
vtk_grid.SetDimensions(*[dim + 1 for dim in self.dimension])
return vtk_grid
def _create_vtk_unstructured_grid(self):
"""Create an unstructured grid of curvilinear elements
representing the mesh
Returns
-------
vtk.vtkUnstructuredGrid
a VTK unstructured grid object representing the mesh
"""
import vtk
from vtk.util import numpy_support as nps
corner_vertices = self.cartesian_vertices.T.reshape(-1, 3)
vtkPts = vtk.vtkPoints()
vtk_grid = vtk.vtkUnstructuredGrid()
vtk_grid.SetPoints(vtkPts)
# add corner vertices to the point set for the unstructured grid
# only insert unique points, we'll get their IDs in the point set to
# define element connectivity later
vtkPts.SetData(nps.numpy_to_vtk(np.unique(corner_vertices, axis=0), deep=True))
# create a locator to assist with duplicate points
locator = vtk.vtkPointLocator()
locator.SetDataSet(vtk_grid)
locator.AutomaticOn() # autmoatically adds points to locator
locator.InitPointInsertion(vtkPts, vtkPts.GetBounds())
locator.BuildLocator()
# this function is used to add new points to the unstructured
# grid. It will return an existing point ID if the point is alread present
def _insert_point(pnt):
result = locator.IsInsertedPoint(pnt)
if result == -1:
point_id = vtkPts.InsertNextPoint(pnt)
locator.InsertPoint(point_id, pnt)
return point_id
else:
return result
### Add all points to the unstructured grid, maintaining a flat list of IDs as we go ###
# flat array storind point IDs for a given vertex
# in the grid
point_ids = []
# add element corner vertices to array
for pnt in corner_vertices:
point_ids.append(_insert_point(pnt))
# get edge midpoints and add them to the
# list of point IDs
midpoint_vertices = self.midpoint_vertices
for edge_grid in midpoint_vertices:
for pnt in edge_grid.T.reshape(-1, 3):
point_ids.append(_insert_point(pnt))
# determine how many elements in each dimension
# and how many points in each grid
n_elem = np.asarray(self.dimension)
n_pnts = n_elem + 1
# create hexes and set points for corner
# vertices
for i, j, k in self.indices:
# handle indices indexed from one
i -= 1
j -= 1
k -= 1
# create a new vtk hex
hex = vtk.vtkQuadraticHexahedron()
# set connectivity the hex corners
for n, (di, dj, dk) in enumerate(_HEX_VERTEX_CONN):
# compute flat index into the point ID list based on i, j, k
# of the vertex
flat_idx = np.ravel_multi_index((i+di, j+dj, k+dk), n_pnts, order='F')
# set corner vertices
hex.GetPointIds().SetId(n, point_ids[flat_idx])
# set connectivity of the hex midpoints
n_midpoint_vertices = [v.size // 3 for v in midpoint_vertices]
for n, (dim, (di, dj, dk)) in enumerate(_HEX_MIDPOINT_CONN):
# initial offset for corner vertices and midpoint dimension
flat_idx = corner_vertices.shape[0] + sum(n_midpoint_vertices[:dim])
# generate a flat index into the table of point IDs
midpoint_shape = midpoint_vertices[dim].shape[1:]
flat_idx += np.ravel_multi_index((i+di, j+dj, k+dk),
midpoint_shape,
order='F')
# set hex midpoint connectivity
hex.GetPointIds().SetId(_N_HEX_VERTICES + n, point_ids[flat_idx])
# add the hex to the grid
vtk_grid.InsertNextCell(hex.GetCellType(), hex.GetPointIds())
return vtk_grid
def _check_vtk_datasets(self, datasets):
"""Perform some basic checks that the datasets are valid for this mesh
Parameters
----------
datasets : dict
Dictionary whose keys are the data labels
and values are the data sets.
"""
for label, dataset in datasets.items():
errmsg = (
f"The size of the dataset '{label}' ({dataset.size}) should be"
@ -244,44 +473,7 @@ class StructuredMesh(MeshBase):
else:
if len(dataset) == self.num_mesh_cells:
raise ValueError(errmsg)
cv.check_type('label', label, str)
vtk_grid = vtk.vtkStructuredGrid()
vtk_grid.SetDimensions(*[dim + 1 for dim in self.dimension])
vtkPts = vtk.vtkPoints()
vtkPts.SetData(nps.numpy_to_vtk(points, deep=True))
vtk_grid.SetPoints(vtkPts)
# create VTK arrays for each of
# the data sets
# maintain a list of the datasets as added
# to the VTK arrays to ensure they persist
# in memory until the file is written
datasets_out = []
for label, dataset in datasets.items():
dataset = np.asarray(dataset).flatten()
datasets_out.append(dataset)
if volume_normalization:
dataset /= self.volumes.T.flatten()
dataset_array = vtk.vtkDoubleArray()
dataset_array.SetName(label)
dataset_array.SetArray(nps.numpy_to_vtk(dataset),
dataset.size,
True)
vtk_grid.GetCellData().AddArray(dataset_array)
# write the .vtk file
writer = vtk.vtkStructuredGridWriter()
writer.SetFileName(str(filename))
writer.SetInputData(vtk_grid)
writer.Write()
return vtk_grid
cv.check_type('data label', label, str)
class RegularMesh(StructuredMesh):
@ -363,6 +555,12 @@ class RegularMesh(StructuredMesh):
dims = self._dimension
return [(u - l) / d for u, l, d in zip(us, ls, dims)]
@property
def cartesian_vertices(self):
"""Returns vertices in cartesian coordiantes. Identical to ``vertices`` for RegularMesh and RectilinearMesh
"""
return self.vertices
@property
def volumes(self):
"""Return Volumes for every mesh cell
@ -763,42 +961,6 @@ class RegularMesh(StructuredMesh):
return root_cell, cells
def write_data_to_vtk(self, filename, datasets, volume_normalization=True):
"""Creates a VTK object of the mesh
Parameters
----------
filename : str or pathlib.Path
Name of the VTK file to write.
datasets : dict
Dictionary whose keys are the data labels
and values are the data sets.
volume_normalization : bool, optional
Whether or not to normalize the data by
the volume of the mesh elements.
Defaults to True.
Returns
-------
vtk.vtkStructuredGrid
the VTK object
"""
ll, ur = self.lower_left, self.upper_right
x_vals = np.linspace(ll[0], ur[0], num=self.dimension[0] + 1)
y_vals = np.linspace(ll[1], ur[1], num=self.dimension[1] + 1)
z_vals = np.linspace(ll[2], ur[2], num=self.dimension[2] + 1)
# create points
pts_cartesian = np.array([[x, y, z] for z in z_vals for y in y_vals for x in x_vals])
return super().write_data_to_vtk(
points=pts_cartesian,
filename=filename,
datasets=datasets,
volume_normalization=volume_normalization
)
def Mesh(*args, **kwargs):
warnings.warn("Mesh has been renamed RegularMesh. Future versions of "
"OpenMC will not accept the name Mesh.")
@ -870,6 +1032,12 @@ class RectilinearMesh(StructuredMesh):
def _grids(self):
return (self.x_grid, self.y_grid, self.z_grid)
@property
def cartesian_vertices(self):
"""Returns vertices in cartesian coordiantes. Identical to ``vertices`` for RegularMesh and RectilinearMesh
"""
return self.vertices
@property
def volumes(self):
"""Return Volumes for every mesh cell
@ -997,36 +1165,6 @@ class RectilinearMesh(StructuredMesh):
return element
def write_data_to_vtk(self, filename, datasets, volume_normalization=True):
"""Creates a VTK object of the mesh
Parameters
----------
filename : str or pathlib.Path
Name of the VTK file to write.
datasets : dict
Dictionary whose keys are the data labels
and values are the data sets.
volume_normalization : bool, optional
Whether or not to normalize the data by
the volume of the mesh elements.
Defaults to True.
Returns
-------
vtk.vtkStructuredGrid
the VTK object
"""
# create points
pts_cartesian = self.vertices.T.reshape(-1, 3)
return super().write_data_to_vtk(
points=pts_cartesian,
filename=filename,
datasets=datasets,
volume_normalization=volume_normalization
)
class CylindricalMesh(StructuredMesh):
"""A 3D cylindrical mesh
@ -1310,41 +1448,22 @@ class CylindricalMesh(StructuredMesh):
return np.multiply.outer(np.outer(V_r, V_p), V_z)
def write_data_to_vtk(self, filename, datasets, volume_normalization=True):
"""Creates a VTK object of the mesh
@property
def cartesian_vertices(self):
return self._convert_to_cartesian(self.vertices, self.origin)
Parameters
----------
filename : str or pathlib.Path
Name of the VTK file to write.
datasets : dict
Dictionary whose keys are the data labels
and values are the data sets.
volume_normalization : bool, optional
Whether or not to normalize the data by
the volume of the mesh elements.
Defaults to True.
Returns
-------
vtk.vtkStructuredGrid
the VTK object
@staticmethod
def _convert_to_cartesian(arr, origin):
"""Converts an array with xyz values in the first dimension (shape (3, ...))
to Cartesian coordinates.
"""
# create points
pts_cylindrical = self.vertices.T.reshape(-1, 3)
pts_cartesian = np.copy(pts_cylindrical)
x = arr[0, ...] * np.cos(arr[1, ...]) + origin[0]
y = arr[0, ...] * np.sin(arr[1, ...]) + origin[1]
arr[0, ...] = x
arr[1, ...] = y
arr[2, ...] += origin[2]
return arr
r, phi = pts_cylindrical[:, 0], pts_cylindrical[:, 1]
pts_cartesian[:, 0] = r * np.cos(phi) + self.origin[0]
pts_cartesian[:, 1] = r * np.sin(phi) + self.origin[1]
pts_cartesian[:, 2] += self.origin[2]
return super().write_data_to_vtk(
points=pts_cartesian,
filename=filename,
datasets=datasets,
volume_normalization=volume_normalization
)
class SphericalMesh(StructuredMesh):
"""A 3D spherical mesh
@ -1532,7 +1651,6 @@ class SphericalMesh(StructuredMesh):
Spherical mesh object
"""
mesh_id = int(get_text(elem, 'id'))
mesh = cls(mesh_id)
mesh.r_grid = [float(x) for x in get_text(elem, "r_grid").split()]
@ -1559,42 +1677,23 @@ class SphericalMesh(StructuredMesh):
return np.multiply.outer(np.outer(V_r, V_t), V_p)
def write_data_to_vtk(self, filename, datasets, volume_normalization=True):
"""Creates a VTK object of the mesh
@property
def cartesian_vertices(self):
return self._convert_to_cartesian(self.vertices, self.origin)
Parameters
----------
filename : str or pathlib.Path
Name of the VTK file to write.
datasets : dict
Dictionary whose keys are the data labels
and values are the data sets.
volume_normalization : bool, optional
Whether or not to normalize the data by
the volume of the mesh elements.
Defaults to True.
Returns
-------
vtk.vtkStructuredGrid
the VTK object
@staticmethod
def _convert_to_cartesian(arr, origin):
"""Converts an array with xyz values in the first dimension (shape (3, ...))
to Cartesian coordinates.
"""
# create points
pts_spherical = self.vertices.T.reshape(-1, 3)
pts_cartesian = np.copy(pts_spherical)
r, theta, phi = pts_spherical[:, 0], pts_spherical[:, 1], pts_spherical[:, 2]
pts_cartesian[:, 0] = r * np.sin(theta) * np.cos(phi) + self.origin[0]
pts_cartesian[:, 1] = r * np.sin(theta) * np.sin(phi) + self.origin[1]
pts_cartesian[:, 2] = r * np.cos(theta) + self.origin[2]
return super().write_data_to_vtk(
points=pts_cartesian,
filename=filename,
datasets=datasets,
volume_normalization=volume_normalization
)
r_xy = arr[0, ...] * np.sin(arr[1, ...])
x = r_xy * np.cos(arr[2, ...])
y = r_xy * np.sin(arr[2, ...])
z = arr[0, ...] * np.cos(arr[1, ...])
arr[0, ...] = x + origin[0]
arr[1, ...] = y + origin[1]
arr[2, ...] = z + origin[2]
return arr
class UnstructuredMesh(MeshBase):
@ -2005,3 +2104,34 @@ def _read_meshes(elem):
out[mesh.id] = mesh
return out
# hexahedron element connectivity
# lower-k connectivity offsets
_HEX_VERTEX_CONN = ((0, 0, 0),
(1, 0, 0),
(1, 1, 0),
(0, 1, 0))
# upper-k connectivity offsets
_HEX_VERTEX_CONN += ((0, 0, 1),
(1, 0, 1),
(1, 1, 1),
(0, 1, 1))
_N_HEX_VERTICES = 8
# lower-k connectivity offsets
_HEX_MIDPOINT_CONN = ((0, (0, 0, 0)),
(1, (1, 0, 0)),
(0, (0, 1, 0)),
(1, (0, 0, 0)))
# upper-k connectivity offsets
_HEX_MIDPOINT_CONN += ((0, (0, 0, 1)),
(1, (1, 0, 1)),
(0, (0, 1, 1)),
(1, (0, 0, 1)))
# mid-plane k connectivity
_HEX_MIDPOINT_CONN += ((2, (0, 0, 0)),
(2, (1, 0, 0)),
(2, (1, 1, 0)),
(2, (0, 1, 0)))

41
tests/unit_tests/mesh_to_vtk/cyl-data-actual.vtk Normal file
View file

@ -0,0 +1,41 @@
# vtk DataFile Version 5.1
vtk output
ASCII
DATASET STRUCTURED_GRID
DIMENSIONS 5 4 4
POINTS 80 double
10 10 -10 11.25 10 -10 12.5 10 -10
13.75 10 -10 15 10 -10 10 10 -10
9.375 11.082531755 -10 8.75 12.165063509 -10 8.125 13.247595264 -10
7.5 14.330127019 -10 10 10 -10 9.375 8.9174682453 -10
8.75 7.8349364905 -10 8.125 6.7524047358 -10 7.5 5.6698729811 -10
10 10 -10 11.25 10 -10 12.5 10 -10
13.75 10 -10 15 10 -10 10 10 -9.3333333333
11.25 10 -9.3333333333 12.5 10 -9.3333333333 13.75 10 -9.3333333333
15 10 -9.3333333333 10 10 -9.3333333333 9.375 11.082531755 -9.3333333333
8.75 12.165063509 -9.3333333333 8.125 13.247595264 -9.3333333333 7.5 14.330127019 -9.3333333333
10 10 -9.3333333333 9.375 8.9174682453 -9.3333333333 8.75 7.8349364905 -9.3333333333
8.125 6.7524047358 -9.3333333333 7.5 5.6698729811 -9.3333333333 10 10 -9.3333333333
11.25 10 -9.3333333333 12.5 10 -9.3333333333 13.75 10 -9.3333333333
15 10 -9.3333333333 10 10 -8.6666666667 11.25 10 -8.6666666667
12.5 10 -8.6666666667 13.75 10 -8.6666666667 15 10 -8.6666666667
10 10 -8.6666666667 9.375 11.082531755 -8.6666666667 8.75 12.165063509 -8.6666666667
8.125 13.247595264 -8.6666666667 7.5 14.330127019 -8.6666666667 10 10 -8.6666666667
9.375 8.9174682453 -8.6666666667 8.75 7.8349364905 -8.6666666667 8.125 6.7524047358 -8.6666666667
7.5 5.6698729811 -8.6666666667 10 10 -8.6666666667 11.25 10 -8.6666666667
12.5 10 -8.6666666667 13.75 10 -8.6666666667 15 10 -8.6666666667
10 10 -8 11.25 10 -8 12.5 10 -8
13.75 10 -8 15 10 -8 10 10 -8
9.375 11.082531755 -8 8.75 12.165063509 -8 8.125 13.247595264 -8
7.5 14.330127019 -8 10 10 -8 9.375 8.9174682453 -8
8.75 7.8349364905 -8 8.125 6.7524047358 -8 7.5 5.6698729811 -8
10 10 -8 11.25 10 -8 12.5 10 -8
13.75 10 -8 15 10 -8
CELL_DATA 36
FIELD FieldData 1
ascending_data 1 36 double
0 100 200 300 400 500 600 700 800
900 1000 1100 1200 1300 1400 1500 1600 1700
1800 1900 2000 2100 2200 2300 2400 2500 2600
2700 2800 2900 3000 3100 3200 3300 3400 3500

329
tests/unit_tests/mesh_to_vtk/cyl-data-curvilinear.vtk Normal file
View file

@ -0,0 +1,329 @@
# vtk DataFile Version 5.1
vtk output
ASCII
DATASET UNSTRUCTURED_GRID
POINTS 323 double
-4.0450849719 -2.9389262615 0 -4.0450849719 -2.9389262615 0.66666666667 -4.0450849719 -2.9389262615 1.3333333333
-4.0450849719 -2.9389262615 2 -4.0450849719 2.9389262615 0 -4.0450849719 2.9389262615 0.66666666667
-4.0450849719 2.9389262615 1.3333333333 -4.0450849719 2.9389262615 2 -3.0338137289 -2.2041946961 0
-3.0338137289 -2.2041946961 0.66666666667 -3.0338137289 -2.2041946961 1.3333333333 -3.0338137289 -2.2041946961 2
-3.0338137289 2.2041946961 0 -3.0338137289 2.2041946961 0.66666666667 -3.0338137289 2.2041946961 1.3333333333
-3.0338137289 2.2041946961 2 -2.0225424859 -1.4694631307 0 -2.0225424859 -1.4694631307 0.66666666667
-2.0225424859 -1.4694631307 1.3333333333 -2.0225424859 -1.4694631307 2 -2.0225424859 1.4694631307 0
-2.0225424859 1.4694631307 0.66666666667 -2.0225424859 1.4694631307 1.3333333333 -2.0225424859 1.4694631307 2
-1.011271243 -0.73473156537 0 -1.011271243 -0.73473156537 0.66666666667 -1.011271243 -0.73473156537 1.3333333333
-1.011271243 -0.73473156537 2 -1.011271243 0.73473156537 0 -1.011271243 0.73473156537 0.66666666667
-1.011271243 0.73473156537 1.3333333333 -1.011271243 0.73473156537 2 0 0 0
0 0 0.66666666667 0 0 1.3333333333 0 0 2
0.38627124297 -1.1888206454 0 0.38627124297 -1.1888206454 0.66666666667 0.38627124297 -1.1888206454 1.3333333333
0.38627124297 -1.1888206454 2 0.38627124297 1.1888206454 0 0.38627124297 1.1888206454 0.66666666667
0.38627124297 1.1888206454 1.3333333333 0.38627124297 1.1888206454 2 0.77254248594 -2.3776412907 0
0.77254248594 -2.3776412907 0.66666666667 0.77254248594 -2.3776412907 1.3333333333 0.77254248594 -2.3776412907 2
0.77254248594 2.3776412907 0 0.77254248594 2.3776412907 0.66666666667 0.77254248594 2.3776412907 1.3333333333
0.77254248594 2.3776412907 2 1.1588137289 -3.5664619361 0 1.1588137289 -3.5664619361 0.66666666667
1.1588137289 -3.5664619361 1.3333333333 1.1588137289 -3.5664619361 2 1.1588137289 3.5664619361 0
1.1588137289 3.5664619361 0.66666666667 1.1588137289 3.5664619361 1.3333333333 1.1588137289 3.5664619361 2
1.25 -3.0616169979e-16 0 1.25 -3.0616169979e-16 0.66666666667 1.25 -3.0616169979e-16 1.3333333333
1.25 -3.0616169979e-16 2 1.25 0 0 1.25 0 0.66666666667
1.25 0 1.3333333333 1.25 0 2 1.5450849719 -4.7552825815 0
1.5450849719 -4.7552825815 0.66666666667 1.5450849719 -4.7552825815 1.3333333333 1.5450849719 -4.7552825815 2
1.5450849719 4.7552825815 0 1.5450849719 4.7552825815 0.66666666667 1.5450849719 4.7552825815 1.3333333333
1.5450849719 4.7552825815 2 2.5 -6.1232339957e-16 0 2.5 -6.1232339957e-16 0.66666666667
2.5 -6.1232339957e-16 1.3333333333 2.5 -6.1232339957e-16 2 2.5 0 0
2.5 0 0.66666666667 2.5 0 1.3333333333 2.5 0 2
3.75 -9.1848509936e-16 0 3.75 -9.1848509936e-16 0.66666666667 3.75 -9.1848509936e-16 1.3333333333
3.75 -9.1848509936e-16 2 3.75 0 0 3.75 0 0.66666666667
3.75 0 1.3333333333 3.75 0 2 5 -1.2246467991e-15 0
5 -1.2246467991e-15 0.66666666667 5 -1.2246467991e-15 1.3333333333 5 -1.2246467991e-15 2
5 0 0 5 0 0.66666666667 5 0 1.3333333333
5 0 2 0.625 0 0 1.875 0 0
3.125 0 0 4.375 0 0 0.19313562148 0.59441032268 0
0.57940686445 1.7832309681 0 0.96567810742 2.9720516134 0 1.3519493504 4.1608722588 0
-0.50563562148 0.36736578268 0 -1.5169068645 1.102097348 0 -2.5281781074 1.8368289134 0
-3.5394493504 2.5715604788 0 -0.50563562148 -0.36736578268 0 -1.5169068645 -1.102097348 0
-2.5281781074 -1.8368289134 0 -3.5394493504 -2.5715604788 0 0.19313562148 -0.59441032268 0
0.57940686445 -1.7832309681 0 0.96567810742 -2.9720516134 0 1.3519493504 -4.1608722588 0
0.625 0 0.66666666667 1.875 0 0.66666666667 3.125 0 0.66666666667
4.375 0 0.66666666667 0.19313562148 0.59441032268 0.66666666667 0.57940686445 1.7832309681 0.66666666667
0.96567810742 2.9720516134 0.66666666667 1.3519493504 4.1608722588 0.66666666667 -0.50563562148 0.36736578268 0.66666666667
-1.5169068645 1.102097348 0.66666666667 -2.5281781074 1.8368289134 0.66666666667 -3.5394493504 2.5715604788 0.66666666667
-0.50563562148 -0.36736578268 0.66666666667 -1.5169068645 -1.102097348 0.66666666667 -2.5281781074 -1.8368289134 0.66666666667
-3.5394493504 -2.5715604788 0.66666666667 0.19313562148 -0.59441032268 0.66666666667 0.57940686445 -1.7832309681 0.66666666667
0.96567810742 -2.9720516134 0.66666666667 1.3519493504 -4.1608722588 0.66666666667 0.625 0 1.3333333333
1.875 0 1.3333333333 3.125 0 1.3333333333 4.375 0 1.3333333333
0.19313562148 0.59441032268 1.3333333333 0.57940686445 1.7832309681 1.3333333333 0.96567810742 2.9720516134 1.3333333333
1.3519493504 4.1608722588 1.3333333333 -0.50563562148 0.36736578268 1.3333333333 -1.5169068645 1.102097348 1.3333333333
-2.5281781074 1.8368289134 1.3333333333 -3.5394493504 2.5715604788 1.3333333333 -0.50563562148 -0.36736578268 1.3333333333
-1.5169068645 -1.102097348 1.3333333333 -2.5281781074 -1.8368289134 1.3333333333 -3.5394493504 -2.5715604788 1.3333333333
0.19313562148 -0.59441032268 1.3333333333 0.57940686445 -1.7832309681 1.3333333333 0.96567810742 -2.9720516134 1.3333333333
1.3519493504 -4.1608722588 1.3333333333 0.625 0 2 1.875 0 2
3.125 0 2 4.375 0 2 0.19313562148 0.59441032268 2
0.57940686445 1.7832309681 2 0.96567810742 2.9720516134 2 1.3519493504 4.1608722588 2
-0.50563562148 0.36736578268 2 -1.5169068645 1.102097348 2 -2.5281781074 1.8368289134 2
-3.5394493504 2.5715604788 2 -0.50563562148 -0.36736578268 2 -1.5169068645 -1.102097348 2
-2.5281781074 -1.8368289134 2 -3.5394493504 -2.5715604788 2 0.19313562148 -0.59441032268 2
0.57940686445 -1.7832309681 2 0.96567810742 -2.9720516134 2 1.3519493504 -4.1608722588 2
1.011271243 0.73473156537 0 2.0225424859 1.4694631307 0 3.0338137289 2.2041946961 0
4.0450849719 2.9389262615 0 -0.38627124297 1.1888206454 0 -0.77254248594 2.3776412907 0
-1.1588137289 3.5664619361 0 -1.5450849719 4.7552825815 0 -1.25 1.5308084989e-16 0
-2.5 3.0616169979e-16 0 -3.75 4.5924254968e-16 0 -5 6.1232339957e-16 0
-0.38627124297 -1.1888206454 0 -0.77254248594 -2.3776412907 0 -1.1588137289 -3.5664619361 0
-1.5450849719 -4.7552825815 0 1.011271243 -0.73473156537 0 2.0225424859 -1.4694631307 0
3.0338137289 -2.2041946961 0 4.0450849719 -2.9389262615 0 1.011271243 0.73473156537 0.66666666667
2.0225424859 1.4694631307 0.66666666667 3.0338137289 2.2041946961 0.66666666667 4.0450849719 2.9389262615 0.66666666667
-0.38627124297 1.1888206454 0.66666666667 -0.77254248594 2.3776412907 0.66666666667 -1.1588137289 3.5664619361 0.66666666667
-1.5450849719 4.7552825815 0.66666666667 -1.25 1.5308084989e-16 0.66666666667 -2.5 3.0616169979e-16 0.66666666667
-3.75 4.5924254968e-16 0.66666666667 -5 6.1232339957e-16 0.66666666667 -0.38627124297 -1.1888206454 0.66666666667
-0.77254248594 -2.3776412907 0.66666666667 -1.1588137289 -3.5664619361 0.66666666667 -1.5450849719 -4.7552825815 0.66666666667
1.011271243 -0.73473156537 0.66666666667 2.0225424859 -1.4694631307 0.66666666667 3.0338137289 -2.2041946961 0.66666666667
4.0450849719 -2.9389262615 0.66666666667 1.011271243 0.73473156537 1.3333333333 2.0225424859 1.4694631307 1.3333333333
3.0338137289 2.2041946961 1.3333333333 4.0450849719 2.9389262615 1.3333333333 -0.38627124297 1.1888206454 1.3333333333
-0.77254248594 2.3776412907 1.3333333333 -1.1588137289 3.5664619361 1.3333333333 -1.5450849719 4.7552825815 1.3333333333
-1.25 1.5308084989e-16 1.3333333333 -2.5 3.0616169979e-16 1.3333333333 -3.75 4.5924254968e-16 1.3333333333
-5 6.1232339957e-16 1.3333333333 -0.38627124297 -1.1888206454 1.3333333333 -0.77254248594 -2.3776412907 1.3333333333
-1.1588137289 -3.5664619361 1.3333333333 -1.5450849719 -4.7552825815 1.3333333333 1.011271243 -0.73473156537 1.3333333333
2.0225424859 -1.4694631307 1.3333333333 3.0338137289 -2.2041946961 1.3333333333 4.0450849719 -2.9389262615 1.3333333333
1.011271243 0.73473156537 2 2.0225424859 1.4694631307 2 3.0338137289 2.2041946961 2
4.0450849719 2.9389262615 2 -0.38627124297 1.1888206454 2 -0.77254248594 2.3776412907 2
-1.1588137289 3.5664619361 2 -1.5450849719 4.7552825815 2 -1.25 1.5308084989e-16 2
-2.5 3.0616169979e-16 2 -3.75 4.5924254968e-16 2 -5 6.1232339957e-16 2
-0.38627124297 -1.1888206454 2 -0.77254248594 -2.3776412907 2 -1.1588137289 -3.5664619361 2
-1.5450849719 -4.7552825815 2 1.011271243 -0.73473156537 2 2.0225424859 -1.4694631307 2
3.0338137289 -2.2041946961 2 4.0450849719 -2.9389262615 2 0 0 0.33333333333
1.25 0 0.33333333333 2.5 0 0.33333333333 3.75 0 0.33333333333
5 0 0.33333333333 0.38627124297 1.1888206454 0.33333333333 0.77254248594 2.3776412907 0.33333333333
1.1588137289 3.5664619361 0.33333333333 1.5450849719 4.7552825815 0.33333333333 -1.011271243 0.73473156537 0.33333333333
-2.0225424859 1.4694631307 0.33333333333 -3.0338137289 2.2041946961 0.33333333333 -4.0450849719 2.9389262615 0.33333333333
-1.011271243 -0.73473156537 0.33333333333 -2.0225424859 -1.4694631307 0.33333333333 -3.0338137289 -2.2041946961 0.33333333333
-4.0450849719 -2.9389262615 0.33333333333 0.38627124297 -1.1888206454 0.33333333333 0.77254248594 -2.3776412907 0.33333333333
1.1588137289 -3.5664619361 0.33333333333 1.5450849719 -4.7552825815 0.33333333333 0 0 1
1.25 0 1 2.5 0 1 3.75 0 1
5 0 1 0.38627124297 1.1888206454 1 0.77254248594 2.3776412907 1
1.1588137289 3.5664619361 1 1.5450849719 4.7552825815 1 -1.011271243 0.73473156537 1
-2.0225424859 1.4694631307 1 -3.0338137289 2.2041946961 1 -4.0450849719 2.9389262615 1
-1.011271243 -0.73473156537 1 -2.0225424859 -1.4694631307 1 -3.0338137289 -2.2041946961 1
-4.0450849719 -2.9389262615 1 0.38627124297 -1.1888206454 1 0.77254248594 -2.3776412907 1
1.1588137289 -3.5664619361 1 1.5450849719 -4.7552825815 1 0 0 1.6666666667
1.25 0 1.6666666667 2.5 0 1.6666666667 3.75 0 1.6666666667
5 0 1.6666666667 0.38627124297 1.1888206454 1.6666666667 0.77254248594 2.3776412907 1.6666666667
1.1588137289 3.5664619361 1.6666666667 1.5450849719 4.7552825815 1.6666666667 -1.011271243 0.73473156537 1.6666666667
-2.0225424859 1.4694631307 1.6666666667 -3.0338137289 2.2041946961 1.6666666667 -4.0450849719 2.9389262615 1.6666666667
-1.011271243 -0.73473156537 1.6666666667 -2.0225424859 -1.4694631307 1.6666666667 -3.0338137289 -2.2041946961 1.6666666667
-4.0450849719 -2.9389262615 1.6666666667 0.38627124297 -1.1888206454 1.6666666667 0.77254248594 -2.3776412907 1.6666666667
1.1588137289 -3.5664619361 1.6666666667 1.5450849719 -4.7552825815 1.6666666667
CELLS 61 1200
OFFSETS vtktypeint64
0 20 40 60 80 100 120 140 160
180 200 220 240 260 280 300 320 340
360 380 400 420 440 460 480 500 520
540 560 580 600 620 640 660 680 700
720 740 760 780 800 820 840 860 880
900 920 940 960 980 1000 1020 1040 1060
1080 1100 1120 1140 1160 1180 1200
CONNECTIVITY vtktypeint64
32 60 40 32 33 61 41 33 100
180 104 32 120 200 124 33 260 261
265 260 60 76 48 40 61 77 49
41 101 181 105 180 121 201 125 200
261 262 266 265 76 84 56 48 77
85 57 49 102 182 106 181 122 202
126 201 262 263 267 266 84 92 72
56 85 93 73 57 103 183 107 182
123 203 127 202 263 264 268 267 32
40 28 32 33 41 29 33 104 184
108 32 124 204 128 33 260 265 269
260 40 48 20 28 41 49 21 29
105 185 109 184 125 205 129 204 265
266 270 269 48 56 12 20 49 57
13 21 106 186 110 185 126 206 130
205 266 267 271 270 56 72 4 12
57 73 5 13 107 187 111 186 127
207 131 206 267 268 272 271 32 28
24 32 33 29 25 33 108 188 112
32 128 208 132 33 260 269 273 260
28 20 16 24 29 21 17 25 109
189 113 188 129 209 133 208 269 270
274 273 20 12 8 16 21 13 9
17 110 190 114 189 130 210 134 209
270 271 275 274 12 4 0 8 13
5 1 9 111 191 115 190 131 211
135 210 271 272 276 275 32 24 36
32 33 25 37 33 112 192 116 32
132 212 136 33 260 273 277 260 24
16 44 36 25 17 45 37 113 193
117 192 133 213 137 212 273 274 278
277 16 8 52 44 17 9 53 45
114 194 118 193 134 214 138 213 274
275 279 278 8 0 68 52 9 1
69 53 115 195 119 194 135 215 139
214 275 276 280 279 32 36 60 32
33 37 61 33 116 196 100 32 136
216 120 33 260 277 261 260 36 44
76 60 37 45 77 61 117 197 101
196 137 217 121 216 277 278 262 261
44 52 84 76 45 53 85 77 118
198 102 197 138 218 122 217 278 279
263 262 52 68 92 84 53 69 93
85 119 199 103 198 139 219 123 218
279 280 264 263 33 61 41 33 34
62 42 34 120 200 124 33 140 220
144 34 281 282 286 281 61 77 49
41 62 78 50 42 121 201 125 200
141 221 145 220 282 283 287 286 77
85 57 49 78 86 58 50 122 202
126 201 142 222 146 221 283 284 288
287 85 93 73 57 86 94 74 58
123 203 127 202 143 223 147 222 284
285 289 288 33 41 29 33 34 42
30 34 124 204 128 33 144 224 148
34 281 286 290 281 41 49 21 29
42 50 22 30 125 205 129 204 145
225 149 224 286 287 291 290 49 57
13 21 50 58 14 22 126 206 130
205 146 226 150 225 287 288 292 291
57 73 5 13 58 74 6 14 127
207 131 206 147 227 151 226 288 289
293 292 33 29 25 33 34 30 26
34 128 208 132 33 148 228 152 34
281 290 294 281 29 21 17 25 30
22 18 26 129 209 133 208 149 229
153 228 290 291 295 294 21 13 9
17 22 14 10 18 130 210 134 209
150 230 154 229 291 292 296 295 13
5 1 9 14 6 2 10 131 211
135 210 151 231 155 230 292 293 297
296 33 25 37 33 34 26 38 34
132 212 136 33 152 232 156 34 281
294 298 281 25 17 45 37 26 18
46 38 133 213 137 212 153 233 157
232 294 295 299 298 17 9 53 45
18 10 54 46 134 214 138 213 154
234 158 233 295 296 300 299 9 1
69 53 10 2 70 54 135 215 139
214 155 235 159 234 296 297 301 300
33 37 61 33 34 38 62 34 136
216 120 33 156 236 140 34 281 298
282 281 37 45 77 61 38 46 78
62 137 217 121 216 157 237 141 236
298 299 283 282 45 53 85 77 46
54 86 78 138 218 122 217 158 238
142 237 299 300 284 283 53 69 93
85 54 70 94 86 139 219 123 218
159 239 143 238 300 301 285 284 34
62 42 34 35 63 43 35 140 220
144 34 160 240 164 35 302 303 307
302 62 78 50 42 63 79 51 43
141 221 145 220 161 241 165 240 303
304 308 307 78 86 58 50 79 87
59 51 142 222 146 221 162 242 166
241 304 305 309 308 86 94 74 58
87 95 75 59 143 223 147 222 163
243 167 242 305 306 310 309 34 42
30 34 35 43 31 35 144 224 148
34 164 244 168 35 302 307 311 302
42 50 22 30 43 51 23 31 145
225 149 224 165 245 169 244 307 308
312 311 50 58 14 22 51 59 15
23 146 226 150 225 166 246 170 245
308 309 313 312 58 74 6 14 59
75 7 15 147 227 151 226 167 247
171 246 309 310 314 313 34 30 26
34 35 31 27 35 148 228 152 34
168 248 172 35 302 311 315 302 30
22 18 26 31 23 19 27 149 229
153 228 169 249 173 248 311 312 316
315 22 14 10 18 23 15 11 19
150 230 154 229 170 250 174 249 312
313 317 316 14 6 2 10 15 7
3 11 151 231 155 230 171 251 175
250 313 314 318 317 34 26 38 34
35 27 39 35 152 232 156 34 172
252 176 35 302 315 319 302 26 18
46 38 27 19 47 39 153 233 157
232 173 253 177 252 315 316 320 319
18 10 54 46 19 11 55 47 154
234 158 233 174 254 178 253 316 317
321 320 10 2 70 54 11 3 71
55 155 235 159 234 175 255 179 254
317 318 322 321 34 38 62 34 35
39 63 35 156 236 140 34 176 256
160 35 302 319 303 302 38 46 78
62 39 47 79 63 157 237 141 236
177 257 161 256 319 320 304 303 46
54 86 78 47 55 87 79 158 238
142 237 178 258 162 257 320 321 305
304 54 70 94 86 55 71 95 87
159 239 143 238 179 259 163 258 321
322 306 305
CELL_TYPES 60
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
CELL_DATA 60
FIELD FieldData 1
ascending_data 1 60 double
0 100 200 300 400 500 600 700 800
900 1000 1100 1200 1300 1400 1500 1600 1700
1800 1900 2000 2100 2200 2300 2400 2500 2600
2700 2800 2900 3000 3100 3200 3300 3400 3500
3600 3700 3800 3900 4000 4100 4200 4300 4400
4500 4600 4700 4800 4900 5000 5100 5200 5300
5400 5500 5600 5700 5800 5900

41
tests/unit_tests/mesh_to_vtk/cyl-data.vtk Normal file
View file

@ -0,0 +1,41 @@
# vtk DataFile Version 5.1
vtk output
ASCII
DATASET STRUCTURED_GRID
DIMENSIONS 5 4 4
POINTS 80 double
10 10 -10 11.25 10 -10 12.5 10 -10
13.75 10 -10 15 10 -10 10 10 -10
9.375 11.082531755 -10 8.75 12.165063509 -10 8.125 13.247595264 -10
7.5 14.330127019 -10 10 10 -10 9.375 8.9174682453 -10
8.75 7.8349364905 -10 8.125 6.7524047358 -10 7.5 5.6698729811 -10
10 10 -10 11.25 10 -10 12.5 10 -10
13.75 10 -10 15 10 -10 10 10 -9.3333333333
11.25 10 -9.3333333333 12.5 10 -9.3333333333 13.75 10 -9.3333333333
15 10 -9.3333333333 10 10 -9.3333333333 9.375 11.082531755 -9.3333333333
8.75 12.165063509 -9.3333333333 8.125 13.247595264 -9.3333333333 7.5 14.330127019 -9.3333333333
10 10 -9.3333333333 9.375 8.9174682453 -9.3333333333 8.75 7.8349364905 -9.3333333333
8.125 6.7524047358 -9.3333333333 7.5 5.6698729811 -9.3333333333 10 10 -9.3333333333
11.25 10 -9.3333333333 12.5 10 -9.3333333333 13.75 10 -9.3333333333
15 10 -9.3333333333 10 10 -8.6666666667 11.25 10 -8.6666666667
12.5 10 -8.6666666667 13.75 10 -8.6666666667 15 10 -8.6666666667
10 10 -8.6666666667 9.375 11.082531755 -8.6666666667 8.75 12.165063509 -8.6666666667
8.125 13.247595264 -8.6666666667 7.5 14.330127019 -8.6666666667 10 10 -8.6666666667
9.375 8.9174682453 -8.6666666667 8.75 7.8349364905 -8.6666666667 8.125 6.7524047358 -8.6666666667
7.5 5.6698729811 -8.6666666667 10 10 -8.6666666667 11.25 10 -8.6666666667
12.5 10 -8.6666666667 13.75 10 -8.6666666667 15 10 -8.6666666667
10 10 -8 11.25 10 -8 12.5 10 -8
13.75 10 -8 15 10 -8 10 10 -8
9.375 11.082531755 -8 8.75 12.165063509 -8 8.125 13.247595264 -8
7.5 14.330127019 -8 10 10 -8 9.375 8.9174682453 -8
8.75 7.8349364905 -8 8.125 6.7524047358 -8 7.5 5.6698729811 -8
10 10 -8 11.25 10 -8 12.5 10 -8
13.75 10 -8 15 10 -8
CELL_DATA 36
FIELD FieldData 1
ascending_data 1 36 double
0 100 200 300 400 500 600 700 800
900 1000 1100 1200 1300 1400 1500 1600 1700
1800 1900 2000 2100 2200 2300 2400 2500 2600
2700 2800 2900 3000 3100 3200 3300 3400 3500

198
tests/unit_tests/mesh_to_vtk/cylindrical-curvilinear.vtk Normal file
View file

@ -0,0 +1,198 @@
# vtk DataFile Version 5.1
vtk output
ASCII
DATASET UNSTRUCTURED_GRID
POINTS 195 double
7.5 5.6698729811 -10 7.5 5.6698729811 -9.3333333333 7.5 5.6698729811 -8.6666666667
7.5 5.6698729811 -8 7.5 14.330127019 -10 7.5 14.330127019 -9.3333333333
7.5 14.330127019 -8.6666666667 7.5 14.330127019 -8 8.125 6.7524047358 -10
8.125 6.7524047358 -9.3333333333 8.125 6.7524047358 -8.6666666667 8.125 6.7524047358 -8
8.125 13.247595264 -10 8.125 13.247595264 -9.3333333333 8.125 13.247595264 -8.6666666667
8.125 13.247595264 -8 8.75 7.8349364905 -10 8.75 7.8349364905 -9.3333333333
8.75 7.8349364905 -8.6666666667 8.75 7.8349364905 -8 8.75 12.165063509 -10
8.75 12.165063509 -9.3333333333 8.75 12.165063509 -8.6666666667 8.75 12.165063509 -8
9.375 8.9174682453 -10 9.375 8.9174682453 -9.3333333333 9.375 8.9174682453 -8.6666666667
9.375 8.9174682453 -8 9.375 11.082531755 -10 9.375 11.082531755 -9.3333333333
9.375 11.082531755 -8.6666666667 9.375 11.082531755 -8 10 10 -10
10 10 -9.3333333333 10 10 -8.6666666667 10 10 -8
11.25 10 -10 11.25 10 -9.3333333333 11.25 10 -8.6666666667
11.25 10 -8 12.5 10 -10 12.5 10 -9.3333333333
12.5 10 -8.6666666667 12.5 10 -8 13.75 10 -10
13.75 10 -9.3333333333 13.75 10 -8.6666666667 13.75 10 -8
13.75 10 -10 13.75 10 -9.3333333333 13.75 10 -8.6666666667
13.75 10 -8 15 10 -10 15 10 -9.3333333333
15 10 -8.6666666667 15 10 -8 15 10 -10
15 10 -9.3333333333 15 10 -8.6666666667 15 10 -8
10.625 10 -10 11.875 10 -10 13.125 10 -10
14.375 10 -10 9.6875 10.541265877 -10 9.0625 11.623797632 -10
8.4375 12.706329387 -10 7.8125 13.788861142 -10 9.6875 9.4587341226 -10
9.0625 8.3762023679 -10 8.4375 7.2936706132 -10 7.8125 6.2111388584 -10
10.625 10 -9.3333333333 11.875 10 -9.3333333333 13.125 10 -9.3333333333
14.375 10 -9.3333333333 9.6875 10.541265877 -9.3333333333 9.0625 11.623797632 -9.3333333333
8.4375 12.706329387 -9.3333333333 7.8125 13.788861142 -9.3333333333 9.6875 9.4587341226 -9.3333333333
9.0625 8.3762023679 -9.3333333333 8.4375 7.2936706132 -9.3333333333 7.8125 6.2111388584 -9.3333333333
10.625 10 -8.6666666667 11.875 10 -8.6666666667 13.125 10 -8.6666666667
14.375 10 -8.6666666667 9.6875 10.541265877 -8.6666666667 9.0625 11.623797632 -8.6666666667
8.4375 12.706329387 -8.6666666667 7.8125 13.788861142 -8.6666666667 9.6875 9.4587341226 -8.6666666667
9.0625 8.3762023679 -8.6666666667 8.4375 7.2936706132 -8.6666666667 7.8125 6.2111388584 -8.6666666667
10.625 10 -8 11.875 10 -8 13.125 10 -8
14.375 10 -8 9.6875 10.541265877 -8 9.0625 11.623797632 -8
8.4375 12.706329387 -8 7.8125 13.788861142 -8 9.6875 9.4587341226 -8
9.0625 8.3762023679 -8 8.4375 7.2936706132 -8 7.8125 6.2111388584 -8
10.625 11.082531755 -10 11.25 12.165063509 -10 11.875 13.247595264 -10
12.5 14.330127019 -10 8.75 10 -10 7.5 10 -10
6.25 10 -10 5 10 -10 10.625 8.9174682453 -10
11.25 7.8349364905 -10 11.875 6.7524047358 -10 12.5 5.6698729811 -10
10.625 11.082531755 -9.3333333333 11.25 12.165063509 -9.3333333333 11.875 13.247595264 -9.3333333333
12.5 14.330127019 -9.3333333333 8.75 10 -9.3333333333 7.5 10 -9.3333333333
6.25 10 -9.3333333333 5 10 -9.3333333333 10.625 8.9174682453 -9.3333333333
11.25 7.8349364905 -9.3333333333 11.875 6.7524047358 -9.3333333333 12.5 5.6698729811 -9.3333333333
10.625 11.082531755 -8.6666666667 11.25 12.165063509 -8.6666666667 11.875 13.247595264 -8.6666666667
12.5 14.330127019 -8.6666666667 8.75 10 -8.6666666667 7.5 10 -8.6666666667
6.25 10 -8.6666666667 5 10 -8.6666666667 10.625 8.9174682453 -8.6666666667
11.25 7.8349364905 -8.6666666667 11.875 6.7524047358 -8.6666666667 12.5 5.6698729811 -8.6666666667
10.625 11.082531755 -8 11.25 12.165063509 -8 11.875 13.247595264 -8
12.5 14.330127019 -8 8.75 10 -8 7.5 10 -8
6.25 10 -8 5 10 -8 10.625 8.9174682453 -8
11.25 7.8349364905 -8 11.875 6.7524047358 -8 12.5 5.6698729811 -8
10 10 -9.6666666667 11.25 10 -9.6666666667 12.5 10 -9.6666666667
13.75 10 -9.6666666667 15 10 -9.6666666667 9.375 11.082531755 -9.6666666667
8.75 12.165063509 -9.6666666667 8.125 13.247595264 -9.6666666667 7.5 14.330127019 -9.6666666667
9.375 8.9174682453 -9.6666666667 8.75 7.8349364905 -9.6666666667 8.125 6.7524047358 -9.6666666667
7.5 5.6698729811 -9.6666666667 10 10 -9 11.25 10 -9
12.5 10 -9 13.75 10 -9 15 10 -9
9.375 11.082531755 -9 8.75 12.165063509 -9 8.125 13.247595264 -9
7.5 14.330127019 -9 9.375 8.9174682453 -9 8.75 7.8349364905 -9
8.125 6.7524047358 -9 7.5 5.6698729811 -9 10 10 -8.3333333333
11.25 10 -8.3333333333 12.5 10 -8.3333333333 13.75 10 -8.3333333333
15 10 -8.3333333333 9.375 11.082531755 -8.3333333333 8.75 12.165063509 -8.3333333333
8.125 13.247595264 -8.3333333333 7.5 14.330127019 -8.3333333333 9.375 8.9174682453 -8.3333333333
8.75 7.8349364905 -8.3333333333 8.125 6.7524047358 -8.3333333333 7.5 5.6698729811 -8.3333333333
CELLS 37 720
OFFSETS vtktypeint64
0 20 40 60 80 100 120 140 160
180 200 220 240 260 280 300 320 340
360 380 400 420 440 460 480 500 520
540 560 580 600 620 640 660 680 700
720
CONNECTIVITY vtktypeint64
32 36 28 32 33 37 29 33 60
108 64 32 72 120 76 33 156 157
161 156 36 40 20 28 37 41 21
29 61 109 65 108 73 121 77 120
157 158 162 161 40 44 12 20 41
45 13 21 62 110 66 109 74 122
78 121 158 159 163 162 44 52 4
12 45 53 5 13 63 111 67 110
75 123 79 122 159 160 164 163 32
28 24 32 33 29 25 33 64 112
68 32 76 124 80 33 156 161 165
156 28 20 16 24 29 21 17 25
65 113 69 112 77 125 81 124 161
162 166 165 20 12 8 16 21 13
9 17 66 114 70 113 78 126 82
125 162 163 167 166 12 4 0 8
13 5 1 9 67 115 71 114 79
127 83 126 163 164 168 167 32 24
36 32 33 25 37 33 68 116 60
32 80 128 72 33 156 165 157 156
24 16 40 36 25 17 41 37 69
117 61 116 81 129 73 128 165 166
158 157 16 8 44 40 17 9 45
41 70 118 62 117 82 130 74 129
166 167 159 158 8 0 52 44 9
1 53 45 71 119 63 118 83 131
75 130 167 168 160 159 33 37 29
33 34 38 30 34 72 120 76 33
84 132 88 34 169 170 174 169 37
41 21 29 38 42 22 30 73 121
77 120 85 133 89 132 170 171 175
174 41 45 13 21 42 46 14 22
74 122 78 121 86 134 90 133 171
172 176 175 45 53 5 13 46 54
6 14 75 123 79 122 87 135 91
134 172 173 177 176 33 29 25 33
34 30 26 34 76 124 80 33 88
136 92 34 169 174 178 169 29 21
17 25 30 22 18 26 77 125 81
124 89 137 93 136 174 175 179 178
21 13 9 17 22 14 10 18 78
126 82 125 90 138 94 137 175 176
180 179 13 5 1 9 14 6 2
10 79 127 83 126 91 139 95 138
176 177 181 180 33 25 37 33 34
26 38 34 80 128 72 33 92 140
84 34 169 178 170 169 25 17 41
37 26 18 42 38 81 129 73 128
93 141 85 140 178 179 171 170 17
9 45 41 18 10 46 42 82 130
74 129 94 142 86 141 179 180 172
171 9 1 53 45 10 2 54 46
83 131 75 130 95 143 87 142 180
181 173 172 34 38 30 34 35 39
31 35 84 132 88 34 96 144 100
35 182 183 187 182 38 42 22 30
39 43 23 31 85 133 89 132 97
145 101 144 183 184 188 187 42 46
14 22 43 47 15 23 86 134 90
133 98 146 102 145 184 185 189 188
46 54 6 14 47 55 7 15 87
135 91 134 99 147 103 146 185 186
190 189 34 30 26 34 35 31 27
35 88 136 92 34 100 148 104 35
182 187 191 182 30 22 18 26 31
23 19 27 89 137 93 136 101 149
105 148 187 188 192 191 22 14 10
18 23 15 11 19 90 138 94 137
102 150 106 149 188 189 193 192 14
6 2 10 15 7 3 11 91 139
95 138 103 151 107 150 189 190 194
193 34 26 38 34 35 27 39 35
92 140 84 34 104 152 96 35 182
191 183 182 26 18 42 38 27 19
43 39 93 141 85 140 105 153 97
152 191 192 184 183 18 10 46 42
19 11 47 43 94 142 86 141 106
154 98 153 192 193 185 184 10 2
54 46 11 3 55 47 95 143 87
142 107 155 99 154 193 194 186 185
CELL_TYPES 36
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25

33
tests/unit_tests/mesh_to_vtk/cylindrical-linear.vtk Normal file
View file

@ -0,0 +1,33 @@
# vtk DataFile Version 5.1
vtk output
ASCII
DATASET STRUCTURED_GRID
DIMENSIONS 5 4 4
POINTS 80 double
10 10 -10 11.25 10 -10 12.5 10 -10
13.75 10 -10 15 10 -10 10 10 -10
9.375 11.082531755 -10 8.75 12.165063509 -10 8.125 13.247595264 -10
7.5 14.330127019 -10 10 10 -10 9.375 8.9174682453 -10
8.75 7.8349364905 -10 8.125 6.7524047358 -10 7.5 5.6698729811 -10
10 10 -10 11.25 10 -10 12.5 10 -10
13.75 10 -10 15 10 -10 10 10 -9.3333333333
11.25 10 -9.3333333333 12.5 10 -9.3333333333 13.75 10 -9.3333333333
15 10 -9.3333333333 10 10 -9.3333333333 9.375 11.082531755 -9.3333333333
8.75 12.165063509 -9.3333333333 8.125 13.247595264 -9.3333333333 7.5 14.330127019 -9.3333333333
10 10 -9.3333333333 9.375 8.9174682453 -9.3333333333 8.75 7.8349364905 -9.3333333333
8.125 6.7524047358 -9.3333333333 7.5 5.6698729811 -9.3333333333 10 10 -9.3333333333
11.25 10 -9.3333333333 12.5 10 -9.3333333333 13.75 10 -9.3333333333
15 10 -9.3333333333 10 10 -8.6666666667 11.25 10 -8.6666666667
12.5 10 -8.6666666667 13.75 10 -8.6666666667 15 10 -8.6666666667
10 10 -8.6666666667 9.375 11.082531755 -8.6666666667 8.75 12.165063509 -8.6666666667
8.125 13.247595264 -8.6666666667 7.5 14.330127019 -8.6666666667 10 10 -8.6666666667
9.375 8.9174682453 -8.6666666667 8.75 7.8349364905 -8.6666666667 8.125 6.7524047358 -8.6666666667
7.5 5.6698729811 -8.6666666667 10 10 -8.6666666667 11.25 10 -8.6666666667
12.5 10 -8.6666666667 13.75 10 -8.6666666667 15 10 -8.6666666667
10 10 -8 11.25 10 -8 12.5 10 -8
13.75 10 -8 15 10 -8 10 10 -8
9.375 11.082531755 -8 8.75 12.165063509 -8 8.125 13.247595264 -8
7.5 14.330127019 -8 10 10 -8 9.375 8.9174682453 -8
8.75 7.8349364905 -8 8.125 6.7524047358 -8 7.5 5.6698729811 -8
10 10 -8 11.25 10 -8 12.5 10 -8
13.75 10 -8 15 10 -8

340
tests/unit_tests/mesh_to_vtk/rectilinear.vtk Normal file
View file

@ -0,0 +1,340 @@
# vtk DataFile Version 5.1
vtk output
ASCII
DATASET STRUCTURED_GRID
DIMENSIONS 5 10 20
POINTS 1000 double
0 5 1 2.5 5 1 5 5 1
7.5 5 1 10 5 1 0 5.8326451979 1
2.5 5.8326451979 1 5 5.8326451979 1 7.5 5.8326451979 1
10 5.8326451979 1 0 6.8039500009 1 2.5 6.8039500009 1
5 6.8039500009 1 7.5 6.8039500009 1 10 6.8039500009 1
0 7.9370052598 1 2.5 7.9370052598 1 5 7.9370052598 1
7.5 7.9370052598 1 10 7.9370052598 1 0 9.2587471229 1
2.5 9.2587471229 1 5 9.2587471229 1 7.5 9.2587471229 1
10 9.2587471229 1 0 10.800597389 1 2.5 10.800597389 1
5 10.800597389 1 7.5 10.800597389 1 10 10.800597389 1
0 12.599210499 1 2.5 12.599210499 1 5 12.599210499 1
7.5 12.599210499 1 10 12.599210499 1 0 14.697344923 1
2.5 14.697344923 1 5 14.697344923 1 7.5 14.697344923 1
10 14.697344923 1 0 17.144879657 1 2.5 17.144879657 1
5 17.144879657 1 7.5 17.144879657 1 10 17.144879657 1
0 20 1 2.5 20 1 5 20 1
7.5 20 1 10 20 1 0 5 6.2105263158
2.5 5 6.2105263158 5 5 6.2105263158 7.5 5 6.2105263158
10 5 6.2105263158 0 5.8326451979 6.2105263158 2.5 5.8326451979 6.2105263158
5 5.8326451979 6.2105263158 7.5 5.8326451979 6.2105263158 10 5.8326451979 6.2105263158
0 6.8039500009 6.2105263158 2.5 6.8039500009 6.2105263158 5 6.8039500009 6.2105263158
7.5 6.8039500009 6.2105263158 10 6.8039500009 6.2105263158 0 7.9370052598 6.2105263158
2.5 7.9370052598 6.2105263158 5 7.9370052598 6.2105263158 7.5 7.9370052598 6.2105263158
10 7.9370052598 6.2105263158 0 9.2587471229 6.2105263158 2.5 9.2587471229 6.2105263158
5 9.2587471229 6.2105263158 7.5 9.2587471229 6.2105263158 10 9.2587471229 6.2105263158
0 10.800597389 6.2105263158 2.5 10.800597389 6.2105263158 5 10.800597389 6.2105263158
7.5 10.800597389 6.2105263158 10 10.800597389 6.2105263158 0 12.599210499 6.2105263158
2.5 12.599210499 6.2105263158 5 12.599210499 6.2105263158 7.5 12.599210499 6.2105263158
10 12.599210499 6.2105263158 0 14.697344923 6.2105263158 2.5 14.697344923 6.2105263158
5 14.697344923 6.2105263158 7.5 14.697344923 6.2105263158 10 14.697344923 6.2105263158
0 17.144879657 6.2105263158 2.5 17.144879657 6.2105263158 5 17.144879657 6.2105263158
7.5 17.144879657 6.2105263158 10 17.144879657 6.2105263158 0 20 6.2105263158
2.5 20 6.2105263158 5 20 6.2105263158 7.5 20 6.2105263158
10 20 6.2105263158 0 5 11.421052632 2.5 5 11.421052632
5 5 11.421052632 7.5 5 11.421052632 10 5 11.421052632
0 5.8326451979 11.421052632 2.5 5.8326451979 11.421052632 5 5.8326451979 11.421052632
7.5 5.8326451979 11.421052632 10 5.8326451979 11.421052632 0 6.8039500009 11.421052632
2.5 6.8039500009 11.421052632 5 6.8039500009 11.421052632 7.5 6.8039500009 11.421052632
10 6.8039500009 11.421052632 0 7.9370052598 11.421052632 2.5 7.9370052598 11.421052632
5 7.9370052598 11.421052632 7.5 7.9370052598 11.421052632 10 7.9370052598 11.421052632
0 9.2587471229 11.421052632 2.5 9.2587471229 11.421052632 5 9.2587471229 11.421052632
7.5 9.2587471229 11.421052632 10 9.2587471229 11.421052632 0 10.800597389 11.421052632
2.5 10.800597389 11.421052632 5 10.800597389 11.421052632 7.5 10.800597389 11.421052632
10 10.800597389 11.421052632 0 12.599210499 11.421052632 2.5 12.599210499 11.421052632
5 12.599210499 11.421052632 7.5 12.599210499 11.421052632 10 12.599210499 11.421052632
0 14.697344923 11.421052632 2.5 14.697344923 11.421052632 5 14.697344923 11.421052632
7.5 14.697344923 11.421052632 10 14.697344923 11.421052632 0 17.144879657 11.421052632
2.5 17.144879657 11.421052632 5 17.144879657 11.421052632 7.5 17.144879657 11.421052632
10 17.144879657 11.421052632 0 20 11.421052632 2.5 20 11.421052632
5 20 11.421052632 7.5 20 11.421052632 10 20 11.421052632
0 5 16.631578947 2.5 5 16.631578947 5 5 16.631578947
7.5 5 16.631578947 10 5 16.631578947 0 5.8326451979 16.631578947
2.5 5.8326451979 16.631578947 5 5.8326451979 16.631578947 7.5 5.8326451979 16.631578947
10 5.8326451979 16.631578947 0 6.8039500009 16.631578947 2.5 6.8039500009 16.631578947
5 6.8039500009 16.631578947 7.5 6.8039500009 16.631578947 10 6.8039500009 16.631578947
0 7.9370052598 16.631578947 2.5 7.9370052598 16.631578947 5 7.9370052598 16.631578947
7.5 7.9370052598 16.631578947 10 7.9370052598 16.631578947 0 9.2587471229 16.631578947
2.5 9.2587471229 16.631578947 5 9.2587471229 16.631578947 7.5 9.2587471229 16.631578947
10 9.2587471229 16.631578947 0 10.800597389 16.631578947 2.5 10.800597389 16.631578947
5 10.800597389 16.631578947 7.5 10.800597389 16.631578947 10 10.800597389 16.631578947
0 12.599210499 16.631578947 2.5 12.599210499 16.631578947 5 12.599210499 16.631578947
7.5 12.599210499 16.631578947 10 12.599210499 16.631578947 0 14.697344923 16.631578947
2.5 14.697344923 16.631578947 5 14.697344923 16.631578947 7.5 14.697344923 16.631578947
10 14.697344923 16.631578947 0 17.144879657 16.631578947 2.5 17.144879657 16.631578947
5 17.144879657 16.631578947 7.5 17.144879657 16.631578947 10 17.144879657 16.631578947
0 20 16.631578947 2.5 20 16.631578947 5 20 16.631578947
7.5 20 16.631578947 10 20 16.631578947 0 5 21.842105263
2.5 5 21.842105263 5 5 21.842105263 7.5 5 21.842105263
10 5 21.842105263 0 5.8326451979 21.842105263 2.5 5.8326451979 21.842105263
5 5.8326451979 21.842105263 7.5 5.8326451979 21.842105263 10 5.8326451979 21.842105263
0 6.8039500009 21.842105263 2.5 6.8039500009 21.842105263 5 6.8039500009 21.842105263
7.5 6.8039500009 21.842105263 10 6.8039500009 21.842105263 0 7.9370052598 21.842105263
2.5 7.9370052598 21.842105263 5 7.9370052598 21.842105263 7.5 7.9370052598 21.842105263
10 7.9370052598 21.842105263 0 9.2587471229 21.842105263 2.5 9.2587471229 21.842105263
5 9.2587471229 21.842105263 7.5 9.2587471229 21.842105263 10 9.2587471229 21.842105263
0 10.800597389 21.842105263 2.5 10.800597389 21.842105263 5 10.800597389 21.842105263
7.5 10.800597389 21.842105263 10 10.800597389 21.842105263 0 12.599210499 21.842105263
2.5 12.599210499 21.842105263 5 12.599210499 21.842105263 7.5 12.599210499 21.842105263
10 12.599210499 21.842105263 0 14.697344923 21.842105263 2.5 14.697344923 21.842105263
5 14.697344923 21.842105263 7.5 14.697344923 21.842105263 10 14.697344923 21.842105263
0 17.144879657 21.842105263 2.5 17.144879657 21.842105263 5 17.144879657 21.842105263
7.5 17.144879657 21.842105263 10 17.144879657 21.842105263 0 20 21.842105263
2.5 20 21.842105263 5 20 21.842105263 7.5 20 21.842105263
10 20 21.842105263 0 5 27.052631579 2.5 5 27.052631579
5 5 27.052631579 7.5 5 27.052631579 10 5 27.052631579
0 5.8326451979 27.052631579 2.5 5.8326451979 27.052631579 5 5.8326451979 27.052631579
7.5 5.8326451979 27.052631579 10 5.8326451979 27.052631579 0 6.8039500009 27.052631579
2.5 6.8039500009 27.052631579 5 6.8039500009 27.052631579 7.5 6.8039500009 27.052631579
10 6.8039500009 27.052631579 0 7.9370052598 27.052631579 2.5 7.9370052598 27.052631579
5 7.9370052598 27.052631579 7.5 7.9370052598 27.052631579 10 7.9370052598 27.052631579
0 9.2587471229 27.052631579 2.5 9.2587471229 27.052631579 5 9.2587471229 27.052631579
7.5 9.2587471229 27.052631579 10 9.2587471229 27.052631579 0 10.800597389 27.052631579
2.5 10.800597389 27.052631579 5 10.800597389 27.052631579 7.5 10.800597389 27.052631579
10 10.800597389 27.052631579 0 12.599210499 27.052631579 2.5 12.599210499 27.052631579
5 12.599210499 27.052631579 7.5 12.599210499 27.052631579 10 12.599210499 27.052631579
0 14.697344923 27.052631579 2.5 14.697344923 27.052631579 5 14.697344923 27.052631579
7.5 14.697344923 27.052631579 10 14.697344923 27.052631579 0 17.144879657 27.052631579
2.5 17.144879657 27.052631579 5 17.144879657 27.052631579 7.5 17.144879657 27.052631579
10 17.144879657 27.052631579 0 20 27.052631579 2.5 20 27.052631579
5 20 27.052631579 7.5 20 27.052631579 10 20 27.052631579
0 5 32.263157895 2.5 5 32.263157895 5 5 32.263157895
7.5 5 32.263157895 10 5 32.263157895 0 5.8326451979 32.263157895
2.5 5.8326451979 32.263157895 5 5.8326451979 32.263157895 7.5 5.8326451979 32.263157895
10 5.8326451979 32.263157895 0 6.8039500009 32.263157895 2.5 6.8039500009 32.263157895
5 6.8039500009 32.263157895 7.5 6.8039500009 32.263157895 10 6.8039500009 32.263157895
0 7.9370052598 32.263157895 2.5 7.9370052598 32.263157895 5 7.9370052598 32.263157895
7.5 7.9370052598 32.263157895 10 7.9370052598 32.263157895 0 9.2587471229 32.263157895
2.5 9.2587471229 32.263157895 5 9.2587471229 32.263157895 7.5 9.2587471229 32.263157895
10 9.2587471229 32.263157895 0 10.800597389 32.263157895 2.5 10.800597389 32.263157895
5 10.800597389 32.263157895 7.5 10.800597389 32.263157895 10 10.800597389 32.263157895
0 12.599210499 32.263157895 2.5 12.599210499 32.263157895 5 12.599210499 32.263157895
7.5 12.599210499 32.263157895 10 12.599210499 32.263157895 0 14.697344923 32.263157895
2.5 14.697344923 32.263157895 5 14.697344923 32.263157895 7.5 14.697344923 32.263157895
10 14.697344923 32.263157895 0 17.144879657 32.263157895 2.5 17.144879657 32.263157895
5 17.144879657 32.263157895 7.5 17.144879657 32.263157895 10 17.144879657 32.263157895
0 20 32.263157895 2.5 20 32.263157895 5 20 32.263157895
7.5 20 32.263157895 10 20 32.263157895 0 5 37.473684211
2.5 5 37.473684211 5 5 37.473684211 7.5 5 37.473684211
10 5 37.473684211 0 5.8326451979 37.473684211 2.5 5.8326451979 37.473684211
5 5.8326451979 37.473684211 7.5 5.8326451979 37.473684211 10 5.8326451979 37.473684211
0 6.8039500009 37.473684211 2.5 6.8039500009 37.473684211 5 6.8039500009 37.473684211
7.5 6.8039500009 37.473684211 10 6.8039500009 37.473684211 0 7.9370052598 37.473684211
2.5 7.9370052598 37.473684211 5 7.9370052598 37.473684211 7.5 7.9370052598 37.473684211
10 7.9370052598 37.473684211 0 9.2587471229 37.473684211 2.5 9.2587471229 37.473684211
5 9.2587471229 37.473684211 7.5 9.2587471229 37.473684211 10 9.2587471229 37.473684211
0 10.800597389 37.473684211 2.5 10.800597389 37.473684211 5 10.800597389 37.473684211
7.5 10.800597389 37.473684211 10 10.800597389 37.473684211 0 12.599210499 37.473684211
2.5 12.599210499 37.473684211 5 12.599210499 37.473684211 7.5 12.599210499 37.473684211
10 12.599210499 37.473684211 0 14.697344923 37.473684211 2.5 14.697344923 37.473684211
5 14.697344923 37.473684211 7.5 14.697344923 37.473684211 10 14.697344923 37.473684211
0 17.144879657 37.473684211 2.5 17.144879657 37.473684211 5 17.144879657 37.473684211
7.5 17.144879657 37.473684211 10 17.144879657 37.473684211 0 20 37.473684211
2.5 20 37.473684211 5 20 37.473684211 7.5 20 37.473684211
10 20 37.473684211 0 5 42.684210526 2.5 5 42.684210526
5 5 42.684210526 7.5 5 42.684210526 10 5 42.684210526
0 5.8326451979 42.684210526 2.5 5.8326451979 42.684210526 5 5.8326451979 42.684210526
7.5 5.8326451979 42.684210526 10 5.8326451979 42.684210526 0 6.8039500009 42.684210526
2.5 6.8039500009 42.684210526 5 6.8039500009 42.684210526 7.5 6.8039500009 42.684210526
10 6.8039500009 42.684210526 0 7.9370052598 42.684210526 2.5 7.9370052598 42.684210526
5 7.9370052598 42.684210526 7.5 7.9370052598 42.684210526 10 7.9370052598 42.684210526
0 9.2587471229 42.684210526 2.5 9.2587471229 42.684210526 5 9.2587471229 42.684210526
7.5 9.2587471229 42.684210526 10 9.2587471229 42.684210526 0 10.800597389 42.684210526
2.5 10.800597389 42.684210526 5 10.800597389 42.684210526 7.5 10.800597389 42.684210526
10 10.800597389 42.684210526 0 12.599210499 42.684210526 2.5 12.599210499 42.684210526
5 12.599210499 42.684210526 7.5 12.599210499 42.684210526 10 12.599210499 42.684210526
0 14.697344923 42.684210526 2.5 14.697344923 42.684210526 5 14.697344923 42.684210526
7.5 14.697344923 42.684210526 10 14.697344923 42.684210526 0 17.144879657 42.684210526
2.5 17.144879657 42.684210526 5 17.144879657 42.684210526 7.5 17.144879657 42.684210526
10 17.144879657 42.684210526 0 20 42.684210526 2.5 20 42.684210526
5 20 42.684210526 7.5 20 42.684210526 10 20 42.684210526
0 5 47.894736842 2.5 5 47.894736842 5 5 47.894736842
7.5 5 47.894736842 10 5 47.894736842 0 5.8326451979 47.894736842
2.5 5.8326451979 47.894736842 5 5.8326451979 47.894736842 7.5 5.8326451979 47.894736842
10 5.8326451979 47.894736842 0 6.8039500009 47.894736842 2.5 6.8039500009 47.894736842
5 6.8039500009 47.894736842 7.5 6.8039500009 47.894736842 10 6.8039500009 47.894736842
0 7.9370052598 47.894736842 2.5 7.9370052598 47.894736842 5 7.9370052598 47.894736842
7.5 7.9370052598 47.894736842 10 7.9370052598 47.894736842 0 9.2587471229 47.894736842
2.5 9.2587471229 47.894736842 5 9.2587471229 47.894736842 7.5 9.2587471229 47.894736842
10 9.2587471229 47.894736842 0 10.800597389 47.894736842 2.5 10.800597389 47.894736842
5 10.800597389 47.894736842 7.5 10.800597389 47.894736842 10 10.800597389 47.894736842
0 12.599210499 47.894736842 2.5 12.599210499 47.894736842 5 12.599210499 47.894736842
7.5 12.599210499 47.894736842 10 12.599210499 47.894736842 0 14.697344923 47.894736842
2.5 14.697344923 47.894736842 5 14.697344923 47.894736842 7.5 14.697344923 47.894736842
10 14.697344923 47.894736842 0 17.144879657 47.894736842 2.5 17.144879657 47.894736842
5 17.144879657 47.894736842 7.5 17.144879657 47.894736842 10 17.144879657 47.894736842
0 20 47.894736842 2.5 20 47.894736842 5 20 47.894736842
7.5 20 47.894736842 10 20 47.894736842 0 5 53.105263158
2.5 5 53.105263158 5 5 53.105263158 7.5 5 53.105263158
10 5 53.105263158 0 5.8326451979 53.105263158 2.5 5.8326451979 53.105263158
5 5.8326451979 53.105263158 7.5 5.8326451979 53.105263158 10 5.8326451979 53.105263158
0 6.8039500009 53.105263158 2.5 6.8039500009 53.105263158 5 6.8039500009 53.105263158
7.5 6.8039500009 53.105263158 10 6.8039500009 53.105263158 0 7.9370052598 53.105263158
2.5 7.9370052598 53.105263158 5 7.9370052598 53.105263158 7.5 7.9370052598 53.105263158
10 7.9370052598 53.105263158 0 9.2587471229 53.105263158 2.5 9.2587471229 53.105263158
5 9.2587471229 53.105263158 7.5 9.2587471229 53.105263158 10 9.2587471229 53.105263158
0 10.800597389 53.105263158 2.5 10.800597389 53.105263158 5 10.800597389 53.105263158
7.5 10.800597389 53.105263158 10 10.800597389 53.105263158 0 12.599210499 53.105263158
2.5 12.599210499 53.105263158 5 12.599210499 53.105263158 7.5 12.599210499 53.105263158
10 12.599210499 53.105263158 0 14.697344923 53.105263158 2.5 14.697344923 53.105263158
5 14.697344923 53.105263158 7.5 14.697344923 53.105263158 10 14.697344923 53.105263158
0 17.144879657 53.105263158 2.5 17.144879657 53.105263158 5 17.144879657 53.105263158
7.5 17.144879657 53.105263158 10 17.144879657 53.105263158 0 20 53.105263158
2.5 20 53.105263158 5 20 53.105263158 7.5 20 53.105263158
10 20 53.105263158 0 5 58.315789474 2.5 5 58.315789474
5 5 58.315789474 7.5 5 58.315789474 10 5 58.315789474
0 5.8326451979 58.315789474 2.5 5.8326451979 58.315789474 5 5.8326451979 58.315789474
7.5 5.8326451979 58.315789474 10 5.8326451979 58.315789474 0 6.8039500009 58.315789474
2.5 6.8039500009 58.315789474 5 6.8039500009 58.315789474 7.5 6.8039500009 58.315789474
10 6.8039500009 58.315789474 0 7.9370052598 58.315789474 2.5 7.9370052598 58.315789474
5 7.9370052598 58.315789474 7.5 7.9370052598 58.315789474 10 7.9370052598 58.315789474
0 9.2587471229 58.315789474 2.5 9.2587471229 58.315789474 5 9.2587471229 58.315789474
7.5 9.2587471229 58.315789474 10 9.2587471229 58.315789474 0 10.800597389 58.315789474
2.5 10.800597389 58.315789474 5 10.800597389 58.315789474 7.5 10.800597389 58.315789474
10 10.800597389 58.315789474 0 12.599210499 58.315789474 2.5 12.599210499 58.315789474
5 12.599210499 58.315789474 7.5 12.599210499 58.315789474 10 12.599210499 58.315789474
0 14.697344923 58.315789474 2.5 14.697344923 58.315789474 5 14.697344923 58.315789474
7.5 14.697344923 58.315789474 10 14.697344923 58.315789474 0 17.144879657 58.315789474
2.5 17.144879657 58.315789474 5 17.144879657 58.315789474 7.5 17.144879657 58.315789474
10 17.144879657 58.315789474 0 20 58.315789474 2.5 20 58.315789474
5 20 58.315789474 7.5 20 58.315789474 10 20 58.315789474
0 5 63.526315789 2.5 5 63.526315789 5 5 63.526315789
7.5 5 63.526315789 10 5 63.526315789 0 5.8326451979 63.526315789
2.5 5.8326451979 63.526315789 5 5.8326451979 63.526315789 7.5 5.8326451979 63.526315789
10 5.8326451979 63.526315789 0 6.8039500009 63.526315789 2.5 6.8039500009 63.526315789
5 6.8039500009 63.526315789 7.5 6.8039500009 63.526315789 10 6.8039500009 63.526315789
0 7.9370052598 63.526315789 2.5 7.9370052598 63.526315789 5 7.9370052598 63.526315789
7.5 7.9370052598 63.526315789 10 7.9370052598 63.526315789 0 9.2587471229 63.526315789
2.5 9.2587471229 63.526315789 5 9.2587471229 63.526315789 7.5 9.2587471229 63.526315789
10 9.2587471229 63.526315789 0 10.800597389 63.526315789 2.5 10.800597389 63.526315789
5 10.800597389 63.526315789 7.5 10.800597389 63.526315789 10 10.800597389 63.526315789
0 12.599210499 63.526315789 2.5 12.599210499 63.526315789 5 12.599210499 63.526315789
7.5 12.599210499 63.526315789 10 12.599210499 63.526315789 0 14.697344923 63.526315789
2.5 14.697344923 63.526315789 5 14.697344923 63.526315789 7.5 14.697344923 63.526315789
10 14.697344923 63.526315789 0 17.144879657 63.526315789 2.5 17.144879657 63.526315789
5 17.144879657 63.526315789 7.5 17.144879657 63.526315789 10 17.144879657 63.526315789
0 20 63.526315789 2.5 20 63.526315789 5 20 63.526315789
7.5 20 63.526315789 10 20 63.526315789 0 5 68.736842105
2.5 5 68.736842105 5 5 68.736842105 7.5 5 68.736842105
10 5 68.736842105 0 5.8326451979 68.736842105 2.5 5.8326451979 68.736842105
5 5.8326451979 68.736842105 7.5 5.8326451979 68.736842105 10 5.8326451979 68.736842105
0 6.8039500009 68.736842105 2.5 6.8039500009 68.736842105 5 6.8039500009 68.736842105
7.5 6.8039500009 68.736842105 10 6.8039500009 68.736842105 0 7.9370052598 68.736842105
2.5 7.9370052598 68.736842105 5 7.9370052598 68.736842105 7.5 7.9370052598 68.736842105
10 7.9370052598 68.736842105 0 9.2587471229 68.736842105 2.5 9.2587471229 68.736842105
5 9.2587471229 68.736842105 7.5 9.2587471229 68.736842105 10 9.2587471229 68.736842105
0 10.800597389 68.736842105 2.5 10.800597389 68.736842105 5 10.800597389 68.736842105
7.5 10.800597389 68.736842105 10 10.800597389 68.736842105 0 12.599210499 68.736842105
2.5 12.599210499 68.736842105 5 12.599210499 68.736842105 7.5 12.599210499 68.736842105
10 12.599210499 68.736842105 0 14.697344923 68.736842105 2.5 14.697344923 68.736842105
5 14.697344923 68.736842105 7.5 14.697344923 68.736842105 10 14.697344923 68.736842105
0 17.144879657 68.736842105 2.5 17.144879657 68.736842105 5 17.144879657 68.736842105
7.5 17.144879657 68.736842105 10 17.144879657 68.736842105 0 20 68.736842105
2.5 20 68.736842105 5 20 68.736842105 7.5 20 68.736842105
10 20 68.736842105 0 5 73.947368421 2.5 5 73.947368421
5 5 73.947368421 7.5 5 73.947368421 10 5 73.947368421
0 5.8326451979 73.947368421 2.5 5.8326451979 73.947368421 5 5.8326451979 73.947368421
7.5 5.8326451979 73.947368421 10 5.8326451979 73.947368421 0 6.8039500009 73.947368421
2.5 6.8039500009 73.947368421 5 6.8039500009 73.947368421 7.5 6.8039500009 73.947368421
10 6.8039500009 73.947368421 0 7.9370052598 73.947368421 2.5 7.9370052598 73.947368421
5 7.9370052598 73.947368421 7.5 7.9370052598 73.947368421 10 7.9370052598 73.947368421
0 9.2587471229 73.947368421 2.5 9.2587471229 73.947368421 5 9.2587471229 73.947368421
7.5 9.2587471229 73.947368421 10 9.2587471229 73.947368421 0 10.800597389 73.947368421
2.5 10.800597389 73.947368421 5 10.800597389 73.947368421 7.5 10.800597389 73.947368421
10 10.800597389 73.947368421 0 12.599210499 73.947368421 2.5 12.599210499 73.947368421
5 12.599210499 73.947368421 7.5 12.599210499 73.947368421 10 12.599210499 73.947368421
0 14.697344923 73.947368421 2.5 14.697344923 73.947368421 5 14.697344923 73.947368421
7.5 14.697344923 73.947368421 10 14.697344923 73.947368421 0 17.144879657 73.947368421
2.5 17.144879657 73.947368421 5 17.144879657 73.947368421 7.5 17.144879657 73.947368421
10 17.144879657 73.947368421 0 20 73.947368421 2.5 20 73.947368421
5 20 73.947368421 7.5 20 73.947368421 10 20 73.947368421
0 5 79.157894737 2.5 5 79.157894737 5 5 79.157894737
7.5 5 79.157894737 10 5 79.157894737 0 5.8326451979 79.157894737
2.5 5.8326451979 79.157894737 5 5.8326451979 79.157894737 7.5 5.8326451979 79.157894737
10 5.8326451979 79.157894737 0 6.8039500009 79.157894737 2.5 6.8039500009 79.157894737
5 6.8039500009 79.157894737 7.5 6.8039500009 79.157894737 10 6.8039500009 79.157894737
0 7.9370052598 79.157894737 2.5 7.9370052598 79.157894737 5 7.9370052598 79.157894737
7.5 7.9370052598 79.157894737 10 7.9370052598 79.157894737 0 9.2587471229 79.157894737
2.5 9.2587471229 79.157894737 5 9.2587471229 79.157894737 7.5 9.2587471229 79.157894737
10 9.2587471229 79.157894737 0 10.800597389 79.157894737 2.5 10.800597389 79.157894737
5 10.800597389 79.157894737 7.5 10.800597389 79.157894737 10 10.800597389 79.157894737
0 12.599210499 79.157894737 2.5 12.599210499 79.157894737 5 12.599210499 79.157894737
7.5 12.599210499 79.157894737 10 12.599210499 79.157894737 0 14.697344923 79.157894737
2.5 14.697344923 79.157894737 5 14.697344923 79.157894737 7.5 14.697344923 79.157894737
10 14.697344923 79.157894737 0 17.144879657 79.157894737 2.5 17.144879657 79.157894737
5 17.144879657 79.157894737 7.5 17.144879657 79.157894737 10 17.144879657 79.157894737
0 20 79.157894737 2.5 20 79.157894737 5 20 79.157894737
7.5 20 79.157894737 10 20 79.157894737 0 5 84.368421053
2.5 5 84.368421053 5 5 84.368421053 7.5 5 84.368421053
10 5 84.368421053 0 5.8326451979 84.368421053 2.5 5.8326451979 84.368421053
5 5.8326451979 84.368421053 7.5 5.8326451979 84.368421053 10 5.8326451979 84.368421053
0 6.8039500009 84.368421053 2.5 6.8039500009 84.368421053 5 6.8039500009 84.368421053
7.5 6.8039500009 84.368421053 10 6.8039500009 84.368421053 0 7.9370052598 84.368421053
2.5 7.9370052598 84.368421053 5 7.9370052598 84.368421053 7.5 7.9370052598 84.368421053
10 7.9370052598 84.368421053 0 9.2587471229 84.368421053 2.5 9.2587471229 84.368421053
5 9.2587471229 84.368421053 7.5 9.2587471229 84.368421053 10 9.2587471229 84.368421053
0 10.800597389 84.368421053 2.5 10.800597389 84.368421053 5 10.800597389 84.368421053
7.5 10.800597389 84.368421053 10 10.800597389 84.368421053 0 12.599210499 84.368421053
2.5 12.599210499 84.368421053 5 12.599210499 84.368421053 7.5 12.599210499 84.368421053
10 12.599210499 84.368421053 0 14.697344923 84.368421053 2.5 14.697344923 84.368421053
5 14.697344923 84.368421053 7.5 14.697344923 84.368421053 10 14.697344923 84.368421053
0 17.144879657 84.368421053 2.5 17.144879657 84.368421053 5 17.144879657 84.368421053
7.5 17.144879657 84.368421053 10 17.144879657 84.368421053 0 20 84.368421053
2.5 20 84.368421053 5 20 84.368421053 7.5 20 84.368421053
10 20 84.368421053 0 5 89.578947368 2.5 5 89.578947368
5 5 89.578947368 7.5 5 89.578947368 10 5 89.578947368
0 5.8326451979 89.578947368 2.5 5.8326451979 89.578947368 5 5.8326451979 89.578947368
7.5 5.8326451979 89.578947368 10 5.8326451979 89.578947368 0 6.8039500009 89.578947368
2.5 6.8039500009 89.578947368 5 6.8039500009 89.578947368 7.5 6.8039500009 89.578947368
10 6.8039500009 89.578947368 0 7.9370052598 89.578947368 2.5 7.9370052598 89.578947368
5 7.9370052598 89.578947368 7.5 7.9370052598 89.578947368 10 7.9370052598 89.578947368
0 9.2587471229 89.578947368 2.5 9.2587471229 89.578947368 5 9.2587471229 89.578947368
7.5 9.2587471229 89.578947368 10 9.2587471229 89.578947368 0 10.800597389 89.578947368
2.5 10.800597389 89.578947368 5 10.800597389 89.578947368 7.5 10.800597389 89.578947368
10 10.800597389 89.578947368 0 12.599210499 89.578947368 2.5 12.599210499 89.578947368
5 12.599210499 89.578947368 7.5 12.599210499 89.578947368 10 12.599210499 89.578947368
0 14.697344923 89.578947368 2.5 14.697344923 89.578947368 5 14.697344923 89.578947368
7.5 14.697344923 89.578947368 10 14.697344923 89.578947368 0 17.144879657 89.578947368
2.5 17.144879657 89.578947368 5 17.144879657 89.578947368 7.5 17.144879657 89.578947368
10 17.144879657 89.578947368 0 20 89.578947368 2.5 20 89.578947368
5 20 89.578947368 7.5 20 89.578947368 10 20 89.578947368
0 5 94.789473684 2.5 5 94.789473684 5 5 94.789473684
7.5 5 94.789473684 10 5 94.789473684 0 5.8326451979 94.789473684
2.5 5.8326451979 94.789473684 5 5.8326451979 94.789473684 7.5 5.8326451979 94.789473684
10 5.8326451979 94.789473684 0 6.8039500009 94.789473684 2.5 6.8039500009 94.789473684
5 6.8039500009 94.789473684 7.5 6.8039500009 94.789473684 10 6.8039500009 94.789473684
0 7.9370052598 94.789473684 2.5 7.9370052598 94.789473684 5 7.9370052598 94.789473684
7.5 7.9370052598 94.789473684 10 7.9370052598 94.789473684 0 9.2587471229 94.789473684
2.5 9.2587471229 94.789473684 5 9.2587471229 94.789473684 7.5 9.2587471229 94.789473684
10 9.2587471229 94.789473684 0 10.800597389 94.789473684 2.5 10.800597389 94.789473684
5 10.800597389 94.789473684 7.5 10.800597389 94.789473684 10 10.800597389 94.789473684
0 12.599210499 94.789473684 2.5 12.599210499 94.789473684 5 12.599210499 94.789473684
7.5 12.599210499 94.789473684 10 12.599210499 94.789473684 0 14.697344923 94.789473684
2.5 14.697344923 94.789473684 5 14.697344923 94.789473684 7.5 14.697344923 94.789473684
10 14.697344923 94.789473684 0 17.144879657 94.789473684 2.5 17.144879657 94.789473684
5 17.144879657 94.789473684 7.5 17.144879657 94.789473684 10 17.144879657 94.789473684
0 20 94.789473684 2.5 20 94.789473684 5 20 94.789473684
7.5 20 94.789473684 10 20 94.789473684 0 5 100
2.5 5 100 5 5 100 7.5 5 100
10 5 100 0 5.8326451979 100 2.5 5.8326451979 100
5 5.8326451979 100 7.5 5.8326451979 100 10 5.8326451979 100
0 6.8039500009 100 2.5 6.8039500009 100 5 6.8039500009 100
7.5 6.8039500009 100 10 6.8039500009 100 0 7.9370052598 100
2.5 7.9370052598 100 5 7.9370052598 100 7.5 7.9370052598 100
10 7.9370052598 100 0 9.2587471229 100 2.5 9.2587471229 100
5 9.2587471229 100 7.5 9.2587471229 100 10 9.2587471229 100
0 10.800597389 100 2.5 10.800597389 100 5 10.800597389 100
7.5 10.800597389 100 10 10.800597389 100 0 12.599210499 100
2.5 12.599210499 100 5 12.599210499 100 7.5 12.599210499 100
10 12.599210499 100 0 14.697344923 100 2.5 14.697344923 100
5 14.697344923 100 7.5 14.697344923 100 10 14.697344923 100
0 17.144879657 100 2.5 17.144879657 100 5 17.144879657 100
7.5 17.144879657 100 10 17.144879657 100 0 20 100
2.5 20 100 5 20 100 7.5 20 100
10 20 100

2394
tests/unit_tests/mesh_to_vtk/regular.vtk Normal file
View file

File diff suppressed because it is too large Load diff

210
tests/unit_tests/mesh_to_vtk/spherical-curvilinear.vtk Normal file
View file

@ -0,0 +1,210 @@
# vtk DataFile Version 5.1
vtk output
ASCII
DATASET UNSTRUCTURED_GRID
POINTS 175 double
5.4951556605 7.8305813044 -10 5.4951556605 12.169418696 -10 6.0986903619 8.1212282982 -7.5
6.0986903619 11.878771702 -7.5 7.7475778302 8.9152906522 -10 7.7475778302 8.9152906522 -5.6698729811
7.7475778302 11.084709348 -10 7.7475778302 11.084709348 -5.6698729811 8.0493451809 9.0606141491 -8.75
8.0493451809 10.939385851 -8.75 8.8737889151 9.4576453261 -7.8349364905 8.8737889151 10.542354674 -7.8349364905
8.8873953302 5.1253604391 -10 8.8873953302 14.874639561 -10 9.0364560916 5.778438306 -7.5
9.0364560916 14.221561694 -7.5 9.4436976651 7.5626802195 -10 9.4436976651 7.5626802195 -5.6698729811
9.4436976651 12.43731978 -10 9.4436976651 12.43731978 -5.6698729811 9.5182280458 7.889219153 -8.75
9.5182280458 12.110780847 -8.75 9.7218488326 8.7813401098 -7.8349364905 9.7218488326 11.21865989 -7.8349364905
10 10 -10 10 10 -7.5 10 10 -5
10.779362252 9.0227106469 -7.8349364905 10.779362252 10.977289353 -7.8349364905 11.25 10 -7.8349364905
11.349895019 8.3072851868 -8.75 11.349895019 11.692714813 -8.75 11.558724505 8.0454212938 -10
11.558724505 8.0454212938 -5.6698729811 11.558724505 11.954578706 -10 11.558724505 11.954578706 -5.6698729811
12.165063509 10 -8.75 12.5 10 -10 12.5 10 -5.6698729811
12.699790037 6.6145703735 -7.5 12.699790037 13.385429626 -7.5 13.117449009 6.0908425877 -10
13.117449009 13.909157412 -10 14.330127019 10 -7.5 14.330127019 10 -7.5
15 10 -10 15 10 -10 10 10 -8.75
10 10 -6.25 10.625 10 -8.9174682453 11.875 10 -6.7524047358
11.082531755 10 -9.375 13.247595264 10 -8.125 11.25 10 -10
13.75 10 -10 10.389681126 10.488644677 -8.9174682453 11.169043378 11.46593403 -6.7524047358
10.674947509 10.846357407 -9.375 12.024842528 12.53907222 -8.125 10.779362252 10.977289353 -10
12.338086757 12.931868059 -10 9.8609244163 10.609329945 -8.9174682453 9.5827732488 11.827989835 -6.7524047358
9.7591140229 11.055390424 -9.375 9.2773420687 13.166171271 -8.125 9.7218488326 11.21865989 -10
9.1655464977 13.655979671 -10 9.4368944576 10.271177337 -8.9174682453 8.3106833727 10.813532011 -6.7524047358
9.0246725905 10.469692925 -9.375 7.0740177714 11.409078776 -8.125 8.8737889151 10.542354674 -10
6.6213667454 11.627064022 -10 9.4368944576 9.7288226631 -8.9174682453 8.3106833727 9.1864679892 -6.7524047358
9.0246725905 9.5303070745 -9.375 7.0740177714 8.5909212236 -8.125 8.8737889151 9.4576453261 -10
6.6213667454 8.3729359783 -10 9.8609244163 9.3906700549 -8.9174682453 9.5827732488 8.1720101647 -6.7524047358
9.7591140229 8.9446095765 -9.375 9.2773420687 6.8338287295 -8.125 9.7218488326 8.7813401098 -10
9.1655464977 6.3440203293 -10 10.389681126 9.5113553235 -8.9174682453 11.169043378 8.5340659704 -6.7524047358
10.674947509 9.1536425934 -9.375 12.024842528 7.4609277801 -8.125 10.779362252 9.0227106469 -10
12.338086757 7.0681319407 -10 10.647047613 10 -7.5851854343 11.294095226 10 -5.1703708686
11.767766953 10 -8.232233047 13.535533906 10 -6.4644660941 12.414814566 10 -9.3529523872
14.829629131 10 -8.7059047745 10.403427588 10.505882194 -7.5851854343 10.806855176 11.011764389 -5.1703708686
11.102184667 11.382095857 -8.232233047 12.204369334 12.764191715 -6.4644660941 11.505612255 11.887978052 -9.3529523872
13.01122451 13.775956104 -8.7059047745 9.8560183609 10.630824778 -7.5851854343 9.7120367218 11.261649556 -5.1703708686
9.6066348466 11.723445345 -8.232233047 9.2132696932 13.446890689 -6.4644660941 9.4626532075 12.354270123 -9.3529523872
8.925306415 14.708540246 -8.7059047745 9.4170302449 10.280743438 -7.5851854343 8.8340604897 10.561486875 -5.1703708686
8.4072970097 10.767005335 -8.232233047 6.8145940193 11.534010671 -6.4644660941 7.8243272545 11.047748773 -9.3529523872
5.6486545091 12.095497546 -8.7059047745 9.4170302449 9.7192565624 -7.5851854343 8.8340604897 9.4385131248 -5.1703708686
8.4072970097 9.2329946646 -8.232233047 6.8145940193 8.4659893291 -6.4644660941 7.8243272545 8.9522512269 -9.3529523872
5.6486545091 7.9045024539 -8.7059047745 9.8560183609 9.3691752218 -7.5851854343 9.7120367218 8.7383504436 -5.1703708686
9.6066348466 8.2765546553 -8.232233047 9.2132696932 6.5531093106 -6.4644660941 9.4626532075 7.6457298771 -9.3529523872
8.925306415 5.2914597543 -8.7059047745 10.403427588 9.4941178057 -7.5851854343 10.806855176 8.9882356114 -5.1703708686
11.102184667 8.6179041425 -8.232233047 12.204369334 7.235808285 -6.4644660941 11.505612255 8.1120219482 -9.3529523872
13.01122451 6.2240438964 -8.7059047745 11.126211085 10.542354674 -7.8349364905 12.25242217 11.084709348 -5.6698729811
11.950654819 10.939385851 -8.75 13.901309638 11.878771702 -7.5 12.25242217 11.084709348 -10
14.50484434 12.169418696 -10 10.278151167 11.21865989 -7.8349364905 10.556302335 12.43731978 -5.6698729811
10.481771954 12.110780847 -8.75 10.963543908 14.221561694 -7.5 10.556302335 12.43731978 -10
11.11260467 14.874639561 -10 9.2206377477 10.977289353 -7.8349364905 8.4412754954 11.954578706 -5.6698729811
8.6501049815 11.692714813 -8.75 7.3002099629 13.385429626 -7.5 8.4412754954 11.954578706 -10
6.8825509907 13.909157412 -10 8.75 10 -7.8349364905 7.5 10 -5.6698729811
7.8349364905 10 -8.75 5.6698729811 10 -7.5 7.5 10 -10
5 10 -10 9.2206377477 9.0227106469 -7.8349364905 8.4412754954 8.0454212938 -5.6698729811
8.6501049815 8.3072851868 -8.75 7.3002099629 6.6145703735 -7.5 8.4412754954 8.0454212938 -10
6.8825509907 6.0908425877 -10 10.278151167 8.7813401098 -7.8349364905 10.556302335 7.5626802195 -5.6698729811
10.481771954 7.889219153 -8.75 10.963543908 5.778438306 -7.5 10.556302335 7.5626802195 -10
11.11260467 5.1253604391 -10 11.126211085 9.4576453261 -7.8349364905 12.25242217 8.9152906522 -5.6698729811
11.950654819 9.0606141491 -8.75 13.901309638 8.1212282982 -7.5 12.25242217 8.9152906522 -10
14.50484434 7.8305813044 -10
CELLS 43 840
OFFSETS vtktypeint64
0 20 40 60 80 100 120 140 160
180 200 220 240 260 280 300 320 340
360 380 400 420 440 460 480 500 520
540 560 580 600 620 640 660 680 700
720 740 760 780 800 820 840
CONNECTIVITY vtktypeint64
24 25 29 24 24 25 28 24 47
91 49 24 47 97 55 24 24 25
133 24 25 26 38 29 25 26 35
28 48 92 50 91 48 98 56 97
25 26 134 133 24 29 36 24 24
28 31 24 49 93 51 24 55 99
57 24 24 133 135 24 29 38 43
36 28 35 40 31 50 94 52 93
56 100 58 99 133 134 136 135 24
36 37 24 24 31 34 24 51 95
53 24 57 101 59 24 24 135 137
24 36 43 45 37 31 40 42 34
52 96 54 95 58 102 60 101 135
136 138 137 24 25 28 24 24 25
23 24 47 97 55 24 47 103 61
24 24 25 139 24 25 26 35 28
25 26 19 23 48 98 56 97 48
104 62 103 25 26 140 139 24 28
31 24 24 23 21 24 55 99 57
24 61 105 63 24 24 139 141 24
28 35 40 31 23 19 15 21 56
100 58 99 62 106 64 105 139 140
142 141 24 31 34 24 24 21 18
24 57 101 59 24 63 107 65 24
24 141 143 24 31 40 42 34 21
15 13 18 58 102 60 101 64 108
66 107 141 142 144 143 24 25 23
24 24 25 11 24 47 103 61 24
47 109 67 24 24 25 145 24 25
26 19 23 25 26 7 11 48 104
62 103 48 110 68 109 25 26 146
145 24 23 21 24 24 11 9 24
61 105 63 24 67 111 69 24 24
145 147 24 23 19 15 21 11 7
3 9 62 106 64 105 68 112 70
111 145 146 148 147 24 21 18 24
24 9 6 24 63 107 65 24 69
113 71 24 24 147 149 24 21 15
13 18 9 3 1 6 64 108 66
107 70 114 72 113 147 148 150 149
24 25 11 24 24 25 10 24 47
109 67 24 47 115 73 24 24 25
151 24 25 26 7 11 25 26 5
10 48 110 68 109 48 116 74 115
25 26 152 151 24 11 9 24 24
10 8 24 67 111 69 24 73 117
75 24 24 151 153 24 11 7 3
9 10 5 2 8 68 112 70 111
74 118 76 117 151 152 154 153 24
9 6 24 24 8 4 24 69 113
71 24 75 119 77 24 24 153 155
24 9 3 1 6 8 2 0 4
70 114 72 113 76 120 78 119 153
154 156 155 24 25 10 24 24 25
22 24 47 115 73 24 47 121 79
24 24 25 157 24 25 26 5 10
25 26 17 22 48 116 74 115 48
122 80 121 25 26 158 157 24 10
8 24 24 22 20 24 73 117 75
24 79 123 81 24 24 157 159 24
10 5 2 8 22 17 14 20 74
118 76 117 80 124 82 123 157 158
160 159 24 8 4 24 24 20 16
24 75 119 77 24 81 125 83 24
24 159 161 24 8 2 0 4 20
14 12 16 76 120 78 119 82 126
84 125 159 160 162 161 24 25 22
24 24 25 27 24 47 121 79 24
47 127 85 24 24 25 163 24 25
26 17 22 25 26 33 27 48 122
80 121 48 128 86 127 25 26 164
163 24 22 20 24 24 27 30 24
79 123 81 24 85 129 87 24 24
163 165 24 22 17 14 20 27 33
39 30 80 124 82 123 86 130 88
129 163 164 166 165 24 20 16 24
24 30 32 24 81 125 83 24 87
131 89 24 24 165 167 24 20 14
12 16 30 39 41 32 82 126 84
125 88 132 90 131 165 166 168 167
24 25 27 24 24 25 29 24 47
127 85 24 47 91 49 24 24 25
169 24 25 26 33 27 25 26 38
29 48 128 86 127 48 92 50 91
25 26 170 169 24 27 30 24 24
29 36 24 85 129 87 24 49 93
51 24 24 169 171 24 27 33 39
30 29 38 43 36 86 130 88 129
50 94 52 93 169 170 172 171 24
30 32 24 24 36 37 24 87 131
89 24 51 95 53 24 24 171 173
24 30 39 41 32 36 43 45 37
88 132 90 131 52 96 54 95 171
172 174 173
CELL_TYPES 42
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25

39
tests/unit_tests/mesh_to_vtk/spherical-linear.vtk Normal file
View file

@ -0,0 +1,39 @@
# vtk DataFile Version 5.1
vtk output
ASCII
DATASET STRUCTURED_GRID
DIMENSIONS 3 4 8
POINTS 96 double
10 10 -10 10 10 -7.5 10 10 -5
10 10 -10 11.25 10 -7.8349364905 12.5 10 -5.6698729811
10 10 -10 12.165063509 10 -8.75 14.330127019 10 -7.5
10 10 -10 12.5 10 -10 15 10 -10
10 10 -10 10 10 -7.5 10 10 -5
10 10 -10 10.779362252 10.977289353 -7.8349364905 11.558724505 11.954578706 -5.6698729811
10 10 -10 11.349895019 11.692714813 -8.75 12.699790037 13.385429626 -7.5
10 10 -10 11.558724505 11.954578706 -10 13.117449009 13.909157412 -10
10 10 -10 10 10 -7.5 10 10 -5
10 10 -10 9.7218488326 11.21865989 -7.8349364905 9.4436976651 12.43731978 -5.6698729811
10 10 -10 9.5182280458 12.110780847 -8.75 9.0364560916 14.221561694 -7.5
10 10 -10 9.4436976651 12.43731978 -10 8.8873953302 14.874639561 -10
10 10 -10 10 10 -7.5 10 10 -5
10 10 -10 8.8737889151 10.542354674 -7.8349364905 7.7475778302 11.084709348 -5.6698729811
10 10 -10 8.0493451809 10.939385851 -8.75 6.0986903619 11.878771702 -7.5
10 10 -10 7.7475778302 11.084709348 -10 5.4951556605 12.169418696 -10
10 10 -10 10 10 -7.5 10 10 -5
10 10 -10 8.8737889151 9.4576453261 -7.8349364905 7.7475778302 8.9152906522 -5.6698729811
10 10 -10 8.0493451809 9.0606141491 -8.75 6.0986903619 8.1212282982 -7.5
10 10 -10 7.7475778302 8.9152906522 -10 5.4951556605 7.8305813044 -10
10 10 -10 10 10 -7.5 10 10 -5
10 10 -10 9.7218488326 8.7813401098 -7.8349364905 9.4436976651 7.5626802195 -5.6698729811
10 10 -10 9.5182280458 7.889219153 -8.75 9.0364560916 5.778438306 -7.5
10 10 -10 9.4436976651 7.5626802195 -10 8.8873953302 5.1253604391 -10
10 10 -10 10 10 -7.5 10 10 -5
10 10 -10 10.779362252 9.0227106469 -7.8349364905 11.558724505 8.0454212938 -5.6698729811
10 10 -10 11.349895019 8.3072851868 -8.75 12.699790037 6.6145703735 -7.5
10 10 -10 11.558724505 8.0454212938 -10 13.117449009 6.0908425877 -10
10 10 -10 10 10 -7.5 10 10 -5
10 10 -10 11.25 10 -7.8349364905 12.5 10 -5.6698729811
10 10 -10 12.165063509 10 -8.75 14.330127019 10 -7.5
10 10 -10 12.5 10 -10 15 10 -10

126
tests/unit_tests/mesh_to_vtk/test.py
View file

@ -1,69 +1,99 @@
import difflib
import filecmp
from itertools import product
import numpy as np
from pathlib import Path
import numpy as np
import openmc
import openmc.lib
import pytest
from tests.regression_tests import config
pytest.importorskip('vtk')
def full_path(f):
return Path(__file__).parent.absolute() / f
def ids_func(param):
return f"{param['library']}_{param['elem_type']}"
def diff_file(file1, file2):
with open(file1) as fh:
f1_text = fh.readlines()
with open(file2) as fh:
f2_text = fh.readlines()
diff_lines = difflib.unified_diff(f1_text, f2_text)
return ''.join(diff_lines)
test_params = (['libmesh', 'moab'],
['tets', 'hexes'])
# test meshes
reg_mesh = openmc.RegularMesh()
reg_mesh.lower_left = (0, 0, 0)
reg_mesh.upper_right = (20, 50, 50)
reg_mesh.dimension = (10, 20, 30)
test_cases = [
{'library' : library, 'elem_type' : elem_type}
for library, elem_type in product(*test_params)
]
rect_mesh = openmc.RectilinearMesh()
rect_mesh.x_grid = np.linspace(0, 10, 5)
rect_mesh.y_grid = np.geomspace(5., 20., 10)
rect_mesh.z_grid = np.linspace(1, 100, 20)
@pytest.mark.parametrize("test_opts", test_cases, ids=ids_func)
def test_unstructured_mesh_to_vtk(run_in_tmpdir, request, test_opts):
cyl_mesh = openmc.CylindricalMesh()
cyl_mesh.origin = (10, 10, -10)
cyl_mesh.r_grid = np.linspace(0, 5, 5)
cyl_mesh.phi_grid = np.linspace(0, 2 * np.pi, 4)
cyl_mesh.z_grid = np.linspace(0, 2, 4)
if test_opts['library'] == 'moab' and test_opts['elem_type'] == 'hexes':
pytest.skip('Hexes are not supported with MOAB')
sphere_mesh = openmc.SphericalMesh()
sphere_mesh.origin = (10, 10, -10)
sphere_mesh.r_grid = np.linspace(0, 5, 3)
sphere_mesh.theta_grid = np.linspace(0, 0.5 * np.pi, 4)
sphere_mesh.phi_grid = np.linspace(0, 2*np.pi, 8)
if test_opts['library'] == 'libmesh' and not openmc.lib._libmesh_enabled():
pytest.skip('LibMesh is not enabled in this build.')
if test_opts['library'] == 'moab' and not openmc.lib._dagmc_enabled():
pytest.skip('DAGMC (and MOAB) mesh not enabled in this build.')
def mesh_data(mesh_dims):
data = 100 * np.arange(np.prod(mesh_dims), dtype=float)
return data.reshape(*mesh_dims)
# pull in a simple model -- just need to create the statepoint file
openmc.reset_auto_ids()
model = openmc.examples.pwr_pin_cell()
test_data = ((reg_mesh, False, 'regular'),
(rect_mesh, False, 'rectilinear'),
(cyl_mesh, False, 'cylindrical-linear'),
(cyl_mesh, True, 'cylindrical-curvilinear'),
(sphere_mesh, False, 'spherical-linear'),
(sphere_mesh, True, 'spherical-curvilinear'))
if test_opts['elem_type'] == 'tets':
filename = Path('tets.exo')
else:
filename = Path('hexes.exo')
@pytest.mark.parametrize('mesh_params',
test_data,
ids=lambda params: params[2])
def test_mesh_write_vtk(mesh_params, run_in_tmpdir):
mesh, curvilinear, filename = mesh_params
# create a basic tally using the unstructured mesh
umesh = openmc.UnstructuredMesh(request.node.path.parent / filename,
test_opts['library'])
umesh.output = False
mesh_filter = openmc.MeshFilter(umesh)
tally = openmc.Tally()
tally.filters = [mesh_filter]
tally.scores = ['flux']
tally.estimator = 'collision'
model.tallies = openmc.Tallies([tally])
sp_file = model.run()
test_data = full_path(filename + ".vtk")
kwargs = {}
if curvilinear:
kwargs['curvilinear'] = curvilinear
# set output filename based on test configuration
filename = test_data if config['update'] else filename + "-actual.vtk"
# write the mesh file and compare to the expected version
mesh.write_data_to_vtk(filename, **kwargs)
try:
assert filecmp.cmp(test_data, filename)
except AssertionError as e:
diff = diff_file(test_data, filename)
raise AssertionError(diff) from e
# check data writing
def test_mesh_write_vtk_data(run_in_tmpdir):
data = {'ascending_data': mesh_data(cyl_mesh.dimension)}
filename_expected = full_path('cyl-data.vtk')
filename_actual = full_path('cyl-data-actual.vtk')
# update the test file if requested
filename = filename_expected if config['update'] else filename_actual
cyl_mesh.write_data_to_vtk(filename, datasets=data, volume_normalization=False)
try:
assert filecmp.cmp(filename, filename_expected)
except AssertionError as e:
diff = diff_file(filename_expected, filename)
raise AssertionError(diff) from e
# check VTK output after reading mesh from statepoint file
with openmc.StatePoint(sp_file) as sp:
umesh = sp.meshes[umesh.id]
test_data = {'ids': np.arange(umesh.n_elements)}
umesh.write_data_to_vtk('umesh.vtk',
datasets=test_data,
volume_normalization=False)
# compare file content with reference file
ref_file = Path(f"{test_opts['library']}_{test_opts['elem_type']}_ref.vtk")
assert filecmp.cmp('umesh.vtk', request.node.path.parent / ref_file)

2
tests/unit_tests/test_mesh_to_vtk.py → tests/unit_tests/mesh_to_vtk/test_vtk_dims.py
View file

@ -212,7 +212,7 @@ def mesh_surf_id(param):
@pytest.mark.parametrize("mesh,surface", product(MESHES, SURFS), ids=mesh_surf_id)
def test_vtk_write_ordering(model, mesh, surface):
def test_vtk_write_ordering(run_in_tmpdir, model, mesh, surface):
tally = openmc.Tally()
tally.scores = ['flux']