def test3():
n = 3 ; p = 2
geo = periodic_line(n=[n], p=[p])
geo_ref, list_lmatrices = geo.bezier_extract()
from scipy.io import mmwrite
M = list_lmatrices[0][0]
mmwrite("M_conversion_test_3.mtx", M)
def test3():
n = 3
p = 2
geo = periodic_line(n=[n], p=[p])
geo_ref, list_lmatrices = geo.bezier_extract()
from scipy.io import mmwrite
M = list_lmatrices[0][0]
mmwrite("M_conversion_test_3.mtx", M)
def test4():
# ... number of elemens
n = 4
# ... k is the b-spline order
k = 3
#
geo1d = periodic_line(n=[n-1], p=[k-1])
# ...
print ("===== connectivity on a periodic line ====")
con1d = connectivity(geo1d)
con1d.init_data_structure()
con1d.printinfo()
print ("==========================================")
def test4():
# ... number of elemens
n = 4
# ... k is the b-spline order
k = 3
#
geo1d = periodic_line(n=[n - 1], p=[k - 1])
# ...
print("===== connectivity on a periodic line ====")
con1d = connectivity(geo1d)
con1d.init_data_structure()
con1d.printinfo()
print("==========================================")
def test4():
print("====== test 4 =====")
nx = 3 ; px = 3
geo = periodic_line(n=[nx], p=[px])
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)