def test4(): n = 3 ; p = 2 geo = periodic_square(n=[n,n], p=[p,p]) geo_ref, list_lmatrices = geo.bezier_extract() from scipy.io import mmwrite M = list_lmatrices[0][0] mmwrite("M_conversion_test_4.mtx", M)
def test4(): n = 3 p = 2 geo = periodic_square(n=[n, n], p=[p, p]) geo_ref, list_lmatrices = geo.bezier_extract() from scipy.io import mmwrite M = list_lmatrices[0][0] mmwrite("M_conversion_test_4.mtx", M)
def test5(): # ... number of elemens n = 4 # ... k is the b-spline order k = 3 # ... print ("==== connectivity on a periodic square ====") geo2d = periodic_square(n=[n-1,n-1], p=[k-1,k-1]) con2d = connectivity(geo2d) con2d.init_data_structure() con2d.printinfo() print ("==========================================")
def test5(): # ... number of elemens n = 4 # ... k is the b-spline order k = 3 # ... print("==== connectivity on a periodic square ====") geo2d = periodic_square(n=[n - 1, n - 1], p=[k - 1, k - 1]) con2d = connectivity(geo2d) con2d.init_data_structure() con2d.printinfo() print("==========================================")
def test5(): print("====== test 5 =====") # nx = 3 ; px = 2 # ny = 3 ; py = 2 # nx = 7 ; px = 2 # ny = 7 ; py = 2 # nx = 15 ; px = 2 # ny = 15 ; py = 2 # nx = 31 ; px = 2 # ny = 31 ; py = 2 # nx = 63 ; px = 2 # ny = 63 ; py = 2 # nx = 3 ; px = 3 # ny = 3 ; py = 3 # nx = 7 ; px = 3 # ny = 7 ; py = 3 # nx = 15 ; px = 3 # ny = 15 ; py = 3 # nx = 31 ; px = 3 # ny = 31 ; py = 3 # nx = 63 ; px = 3 # ny = 63 ; py = 3 # nx = 3 ; px = 4 # ny = 3 ; py = 4 # nx = 7 ; px = 4 # ny = 7 ; py = 4 # nx = 15 ; px = 4 # ny = 15 ; py = 4 # nx = 31 ; px = 4 # ny = 31 ; py = 4 # nx = 63 ; px = 4 # ny = 63 ; py = 4 # nx = 3 ; px = 5 # ny = 3 ; py = 5 # nx = 7 ; px = 5 # ny = 7 ; py = 5 nx = 15 ; px = 5 ny = 15 ; py = 5 # nx = 31 ; px = 5 # ny = 31 ; py = 5 # nx = 63 ; px = 5 # ny = 63 ; py = 5 # geo = square(n=[nx, ny], p=[px, py]) geo = periodic_square(n=[nx, ny], p=[px, py]) basis_only = False #basis_only = True basename = "splines" #basename = None dirname = "tmp" #dirname = None import os os.system("mkdir -p " + dirname) from caid.io import BZR rw = BZR() rw.write(geo, fmt="txt", basename=basename, dirname=dirname, basis_only=basis_only)