pttype=CGK.PointList_s)
g = CGL.newGridLocation(n, value=CGK.FaceCenter_s)
return (T, b, z, zbc, n, g)
vertexsize = 20
cellsize = 7
ntris = 12
(T, b, z, zbc, n, g) = makeCorrectTree(vertexsize, cellsize, ntris)
TESTS.append((tag, T, diag))
# -------------------------------------------------------------------------
tag = 'bc bad location'
diag = False
(T, b, z, zbc, n, g) = makeCorrectTree(vertexsize, cellsize, ntris)
CGL.newPointRange(n, value=NPY.array([[1, cellsize]], 'i'))
zbc[2] = []
n = CGL.newBoundary(zbc,
'BC', [range(1, cellsize + 1)],
btype=CGK.Null_s,
family=None,
pttype=CGK.PointList_s)
g = CGL.newGridLocation(n, value=CGK.CellCenter_s)
TESTS.append((tag, T, diag))
# -------------------------------------------------------------------------
tag = 'bc both PointList and PointRange'
diag = False
(T, b, z, zbc, n, g) = makeCorrectTree(vertexsize, cellsize, ntris)
CGL.newPointRange(n, value=NPY.array([[cellsize + 1, cellsize + ntris]], 'i'))
TESTS.append((tag, T, diag))
d = CGL.newDataClass(a, CGK.NondimensionalParameter_s)
i = CGL.newBaseIterativeData(b, '{BaseIterativeData}')
d = CGL.newDataClass(i, CGK.NondimensionalParameter_s)
i = CGL.newZoneIterativeData(z, '{ZoneIterativeData}')
d = CGL.newDataClass(i, CGK.NondimensionalParameter_s)
m = CGL.newRigidGridMotion(z,
'{RigidGridMotion}',
vector=NPY.array([[0.0, 0.0], [0.0, 0.0],
[0.0, 0.0]]))
d = CGL.newDataClass(m, CGK.NondimensionalParameter_s)
m = CGL.newArbitraryGridMotion(z, '{ArbitraryGridMotion}')
d = CGL.newDataClass(m, CGK.NondimensionalParameter_s)
x = CGL.newZoneGridConnectivity(z)
x = CGL.newGridConnectivity(x, '{GridConnectivity}', z[0])
CGL.newPointRange(x,
value=NPY.array([[1, 1], [1, 1], [1, 1]],
dtype=NPY.int32,
order='Fortran'))
p = CGL.newGridConnectivityProperty(x)
m = CGL.newPeriodic(p)
d = CGL.newDataClass(m, CGK.NondimensionalParameter_s)
w = CGL.newConvergenceHistory(b)
d = CGL.newDataClass(w, CGK.NondimensionalParameter_s)
i = CGL.newIntegralData(b, '{IntegralData}')
d = CGL.newDataClass(i, CGK.NondimensionalParameter_s)
i = CGL.newUserDefinedData(b, '{UserDefinedData}')
d = CGL.newDataClass(i, CGK.NondimensionalParameter_s)
i = CGL.newGravity(b)
d = CGL.newDataClass(i, CGK.NondimensionalParameter_s)
f = CGL.newFlowEquationSet(b)
d = CGL.newDataClass(f, CGK.NondimensionalParameter_s)
m = CGL.newGasModel(f)
NPY.array([[1, cellsize]], 'i', order='F'))
tris = CGL.newElements(
z2, 'TRIS', CGK.TRI_3_s, NPY.ones((ntris * 3), dtype='int32'),
NPY.array([[cellsize + 1, cellsize + ntris]], 'i', order='F'))
z3 = CGU.copyNode(z1, 'Zone3')
b[2].append(z3)
z4 = CGU.copyNode(z2, 'Zone4')
b[2].append(z4)
z = [z1, z2, z3, z4]
return (T, b, z)
(T, b, z) = makeCorrectTree()
zgc = CGL.newZoneGridConnectivity(z[0])
gc = CGL.newGridConnectivity(zgc, 'join1_3', 'Zone3', ctype=CGK.Abutting1to1_s)
CGL.newPointRange(gc, value=NPY.array([[1, 1], [1, 7], [1, 9]], order='F'))
CGL.newPointRange(gc,
name=CGK.PointRangeDonor_s,
value=NPY.array([[1, 1], [1, 7], [1, 9]], order='F'))
gc = CGL.newGridConnectivity(zgc, 'join1_2', 'Zone2', ctype=CGK.Abutting1to1_s)
CGL.newPointRange(gc, value=NPY.array([[1, 5], [1, 1], [1, 9]], order='F'))
CGL.newIndexArray(gc,
CGK.PointListDonor_s,
value=NPY.array([range(cellsize + 1, cellsize + ntris + 1)],
order='F'))
CGL.newGridLocation(gc, value=CGK.FaceCenter_s)
TESTS.append((tag, T, diag))
# -------------------------------------------------------------------------
tag = 'gridconnectivity bad datatype'
diag = False
# pyCGNS - Python package for CFD General Notation System -
# See license.txt file in the root directory of this Python module source
# ---------------------------------------------------------------------------
#
import CGNS.PAT.cgnslib as C
import CGNS.PAT.cgnserrors as E
import CGNS.PAT.cgnskeywords as K
import numpy as N
import copy
#
import BCData_t
#
data = C.newBCDataSet(None, '{BCDataSet}')
C.newGridLocation(data)
C.newPointRange(data)
C.newPointList(data)
C.newDescriptor(data, '{Descriptor}')
C.newDataClass(data)
C.newDimensionalUnits(data)
C.newReferenceState(data)
C.newUserDefinedData(data, '{UserDefinedData}')
#
d1 = copy.deepcopy(BCData_t.pattern[0])
d1[0] = K.NeumannData_s
data[2].append(d1)
#
d2 = copy.deepcopy(BCData_t.pattern[0])
d2[0] = K.DirichletData_s
data[2].append(d2)
#
