def test_rdpostop2_jsc_cla_model():
post = op2.rdpostop2("tests/nas2cam_extseout/cla_test_model.op2")
assert len(post["geom1"]) == len(post["geom1_list"]) == 2
grid = np.array([
[1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0],
[2.0, 0.0, 1.0, 0.0, 0.0, 3.0, 0.0, 0.0],
[3.0, 0.0, 2.0, 0.0, 0.0, 4.0, 0.0, 0.0],
[4.0, 0.0, 3.0, 0.0, 0.0, 5.0, 0.0, 0.0],
[5.0, 0.0, 3.0, 0.0, 0.0, 0.0, 0.0, 0.0],
])
assert (post["geom1"][0]["grid"] == grid).all()
extrn = np.array([
[5, 123456],
[90560001, 1],
[90560002, 1],
[90560003, 1],
[90560004, 1],
[90560005, 1],
[90560006, 1],
[-1, -1],
])
assert (post["geom1"][56]["extrn"] == extrn).all()
def test_rdpostop2_assemble():
post = op2.rdpostop2("tests/nas2cam_extseout/assemble.op2", 1, 1)
assert np.all(post["selist"] == [[101, 0], [102, 0], [0, 0]])
sebulk = np.array([[101, 5, -1, 2, 0.0, 1, 3, 101],
[102, 5, -1, 2, 0.0, 1, 3, 102]])
assert np.all(post["sebulk"] == sebulk)
with op2.OP2("tests/nas2cam_extseout/assemble.op2") as o2:
o2._rowsCutoff = 0
# fpos = o2.dbnames['GEOM1S'][2][0][0]
# o2._fileh.seek(fpos)
o2.set_position("GEOM1S", 2)
name, trailer, dbtype = o2.rdop2nt()
# (cords, sebulk, selist, seload, seconct) = o2._rdop2geom1cord2()
dct = o2._rdop2geom1cord2()
assert np.all(post["selist"] == dct["selist"])
assert np.all(post["sebulk"] == dct["sebulk"])
def test_rdpostop2():
post = op2.rdpostop2("tests/nas2cam_extseout/inboard.op2", 1, 1, 1, 1)
dof = post["mats"]["ougv1"][0]["dof"]
lam = post["mats"]["ougv1"][0]["lambda"]
ougv1 = post["mats"]["ougv1"][0]["ougv1"]
o4 = op4.load("tests/nas2cam_extseout/inboard.op4")
mug1 = o4["mug1"][0]
tug1 = nastran.rddtipch("tests/nas2cam_extseout/inboard.pch")
tef1 = nastran.rddtipch("tests/nas2cam_extseout/inboard.pch", "TEF1")
tes1 = nastran.rddtipch("tests/nas2cam_extseout/inboard.pch", "TES1")
# ougv1, oef1, oes1 ... they don't have the constraint modes
# or the resflex modes! How can they be useful? Anyway, this
# checks the values present:
# mug1 has 24 b-set ... get first 3 modes (rest are resflex):
modes = mug1[:, 24:27]
pv = locate.mat_intersect(tug1, dof)[0]
assert np.allclose(modes[pv], ougv1)
assert np.allclose(o4["mef1"][0][:, 24:27], post["mats"]["oef1"][0][0])
assert np.all(post["mats"]["oef1"][0][1][:, 0] == tef1[:, 0])
assert np.all(post["mats"]["oef1"][0][1][:, 1] == 34)
pv = np.ones(15, bool)
pv[5:7] = False
pv[11:15] = False
pv = np.hstack((pv, pv, pv))
assert np.allclose(o4["mes1"][0][:, 24:27], post["mats"]["oes1"][0][0][pv])
assert np.all(post["mats"]["oes1"][0][1][pv, 0] == tes1[:, 0])
assert np.all(post["mats"]["oes1"][0][1][:, 1] == 34)
with op2.OP2("tests/nas2cam_extseout/inboard.op2") as o2:
o2._rowsCutoff = 0
fpos = o2.dbnames["OUGV1"][0][0][0]
o2._fileh.seek(fpos)
name, trailer, dbtype = o2.rdop2nt()
oug = o2._rdop2ougv1("OUGV1")
assert np.all(oug["ougv1"] == ougv1)
assert np.all(oug["dof"] == dof)
assert np.all(oug["lambda"] == lam)
def test_rdpostop2_many_nodes():
post = op2.rdpostop2("tests/nas2cam_extseout/cant_beam.op2")
assert (post["geom1"][0]["grid"][:, 0].astype(int) == np.arange(8222) +
1).all()
