def test_05SEScalculations():
print("___________________")
print("test_05SEScalculations")
print("___________________")
from mslib import MSMS, msms
m = MSMS(coords=coords, radii=rad)
m.compute_rs(allComponents=1)
m.compute_ses()
print(m.info())
msms.MS_reset_RSR(m.this)
m.compute_rs(allComponents=1)
m.compute_ses(0)
m.compute_ses(1)
print(m.info())
def test_05SEScalculations():
print("___________________")
print("test_05SEScalculations")
print("___________________")
from mslib import MSMS, msms
m = MSMS(coords=coords, radii = rad)
m.compute_rs( allComponents=1 )
m.compute_ses()
print(m.info())
msms.MS_reset_RSR(m.this)
m.compute_rs( allComponents=1 )
m.compute_ses( 0 )
m.compute_ses( 1 )
print(m.info())
def test_06Triangulation():
print("___________________")
print("test_06Triangulation")
print("___________________")
from mslib import MSMS, msms
# triangulate with default parameters, component = 0
m = MSMS(coords=coords, radii=rad)
m.compute_rs(allComponents=1)
m.compute_ses(0)
m.compute_ses(1)
print("trinagulate 1")
m.triangulate(component=0)
print(m.info())
# triangulate with default parameters, component = 1
print("trinagulate 2")
m.triangulate(component=1)
print(m.info())
# re-triangulate at density = 3.0
print("trinagulate 3")
m.triangulate(component=0, density=3.0)
print(m.info())
# triangulate all components
print("trinagulate 4")
m.triangulate(density=1.1)
# write triangulation
#m.write_triangulation('testAll')
#m.write_triangulation('testAllnoh', no_header=1)
#m.write_triangulation('testASES', component=1, format=msms.MS_ASES_ASCII)
#m.write_triangulation('testTAVS', component=0, format=msms.MS_TSES_ASCII_AVS)
#m.write_triangulation('testAAVS', component=0, format=msms.MS_ASES_ASCII_AVS)
vfloat, vint, tri = m.getTriangles()
print("getTriangles(): len(vfloat): %d, len(vint): %d, len(tri): %d" %
(len(vfloat), len(vint), len(tri)))
atomindices = [
131, 120, 143, 147, 128, 146, 142, 137, 139, 130, 135, 136, 134, 129,
145, 119, 133, 132, 141, 144, 138, 140
]
vfloat, vint, tri = m.getTriangles(atomIndices=atomindices)
print(
"getTriangles(atomIndices=atomindices): len(vfloat): %d, len(vint): %d, len(tri): %d"
% (len(vfloat), len(vint), len(tri)))
# reset all after triangulation
msms.MS_reset_RSR(m.this)
print(m.info())
def test_06Triangulation():
print("___________________")
print("test_06Triangulation")
print("___________________")
from mslib import MSMS, msms
# triangulate with default parameters, component = 0
m = MSMS(coords=coords, radii = rad)
m.compute_rs(allComponents=1)
m.compute_ses(0 )
m.compute_ses(1)
print("trinagulate 1")
m.triangulate( component=0 )
print(m.info())
# triangulate with default parameters, component = 1
print("trinagulate 2")
m.triangulate( component=1 )
print(m.info())
# re-triangulate at density = 3.0
print("trinagulate 3")
m.triangulate( component=0, density=3.0 )
print(m.info())
# triangulate all components
print("trinagulate 4")
m.triangulate( density=1.1 )
# write triangulation
#m.write_triangulation('testAll')
#m.write_triangulation('testAllnoh', no_header=1)
#m.write_triangulation('testASES', component=1, format=msms.MS_ASES_ASCII)
#m.write_triangulation('testTAVS', component=0, format=msms.MS_TSES_ASCII_AVS)
#m.write_triangulation('testAAVS', component=0, format=msms.MS_ASES_ASCII_AVS)
vfloat, vint, tri = m.getTriangles()
print("getTriangles(): len(vfloat): %d, len(vint): %d, len(tri): %d" % (len(vfloat), len(vint), len(tri)))
atomindices = [131, 120, 143, 147, 128, 146, 142, 137, 139, 130, 135,
136, 134, 129, 145, 119, 133, 132, 141, 144, 138, 140]
vfloat, vint, tri = m.getTriangles(atomIndices=atomindices)
print("getTriangles(atomIndices=atomindices): len(vfloat): %d, len(vint): %d, len(tri): %d" % (len(vfloat), len(vint), len(tri)))
# reset all after triangulation
msms.MS_reset_RSR(m.this)
print(m.info())
def test_04RScalculations():
print("___________________")
print("test_04RScalculations")
print("___________________")
from mslib import MSMS, readxyzr, msms
# compute all components with default probe radius 1.5
global coords, rad
m1 = MSMS(coords=coords, radii=rad)
m1.compute_rs()
ar = m1.rsr.fst.far0
print("printing -n, -ct...")
print(ar._n())
print(ar._ct())
print(ar._cp())
print(ar._cc())
print(ar._rcc())
print(ar._ps())
print(ar._s())
# recompute reduced surface with probe sphere radius = 2.0
m1.compute_rs(probe_radius=2.0)
# Multiple components
output = readxyzr('Data/1crn.xyzr')
new_coords = output[0]
names = output[1]
#new_coords, names = readxyzr('Data/1crn.xyzr')
coords = new_coords[:, :3]
rad = new_coords[:, 3]
m2 = MSMS(coords=coords, radii=rad)
# "***** COMPUTE EXTERNAL COMPONENT"
m2.compute_rs()
print("Info1:")
print(m2.info())
# " ADD OTHER COMPONENTS"
m2.compute_rs(allComponents=1)
print("Info2:")
print(m2.info())
# " CLEAR RSR"
# " COMPUTE ONE CAVITY"
msms.MS_reset_RSR(m2.this)
m2.compute_rs(atoms=(11, -1, -1))
print("Info3:")
print(m2.info())
msms.MS_reset_RSR(m2.this)
m2.compute_rs(atoms=(13, -1, -1))
print("Info4:")
print(m2.info())
msms.MS_reset_RSR(m2.this)
m2.compute_rs(atoms=(93, -1, -1))
print("Info5:")
print(m2.info())
msms.MS_reset_RSR(m2.this)
m2.compute_rs(atoms=(11, 13, -1))
print("Info6:")
print(m2.info())
msms.MS_reset_RSR(m2.this)
m2.compute_rs(atoms=(13, 11, -1))
print("Info7:")
print(m2.info())
msms.MS_reset_RSR(m2.this)
m2.compute_rs(atoms=(13, 93, -1))
print("Info8:")
print(m2.info())
msms.MS_reset_RSR(m2.this)
m2.compute_rs(atoms=(13, 93, 11))
print("Info9:")
print(m2.info())
def test_04RScalculations():
print("___________________")
print("test_04RScalculations")
print("___________________")
from mslib import MSMS, readxyzr, msms
# compute all components with default probe radius 1.5
global coords, rad
m1 = MSMS(coords=coords, radii=rad)
m1.compute_rs()
ar = m1.rsr.fst.far0
print("printing -n, -ct...")
print(ar._n())
print(ar._ct())
print(ar._cp())
print(ar._cc())
print(ar._rcc())
print(ar._ps())
print(ar._s())
# recompute reduced surface with probe sphere radius = 2.0
m1.compute_rs(probe_radius=2.0)
# Multiple components
output = readxyzr('Data/1crn.xyzr')
new_coords = output[0]
names = output[1]
#new_coords, names = readxyzr('Data/1crn.xyzr')
coords = new_coords[:,:3]
rad = new_coords[:,3]
m2 = MSMS(coords=coords, radii = rad)
# "***** COMPUTE EXTERNAL COMPONENT"
m2.compute_rs()
print("Info1:")
print(m2.info())
# " ADD OTHER COMPONENTS"
m2.compute_rs( allComponents=1 )
print("Info2:")
print(m2.info())
# " CLEAR RSR"
# " COMPUTE ONE CAVITY"
msms.MS_reset_RSR(m2.this)
m2.compute_rs( atoms = (11, -1, -1) )
print("Info3:")
print(m2.info())
msms.MS_reset_RSR(m2.this)
m2.compute_rs( atoms = (13, -1, -1) )
print("Info4:")
print(m2.info())
msms.MS_reset_RSR(m2.this)
m2.compute_rs( atoms = (93, -1, -1) )
print("Info5:")
print(m2.info())
msms.MS_reset_RSR(m2.this)
m2.compute_rs( atoms = (11, 13, -1) )
print("Info6:")
print(m2.info())
msms.MS_reset_RSR(m2.this)
m2.compute_rs( atoms = (13, 11, -1) )
print("Info7:")
print(m2.info())
msms.MS_reset_RSR(m2.this)
m2.compute_rs( atoms = (13, 93, -1) )
print("Info8:")
print(m2.info())
msms.MS_reset_RSR(m2.this)
m2.compute_rs( atoms = (13, 93, 11) )
print("Info9:")
print(m2.info())
