Skip to content
Snippets Groups Projects
Commit 92db2f30 authored by Nima's avatar Nima
Browse files

Removing some line for tests from free-form.. file

parent 99126b67
Branches add/line_load
No related tags found
No related merge requests found
......@@ -37,11 +37,12 @@ surfs = sgs.SurfaceGeo(data, 0, thk)
############# Visualization
# sgs.visualization(surfs)
############# Visualization
p_1 = surfs.physical_crd(0., 0.)
p_2 = surfs.physical_crd(1., 0.)
p_3 = surfs.physical_crd(1., 1.)
print("p_1:", p_1, " p_2:", p_2, " p_3:", p_3)
# p_1 = surfs.physical_crd(0., 0.)
# p_2 = surfs.physical_crd(1., 0.)
# p_3 = surfs.physical_crd(1., 1.)
# print("p_1:", p_1, " p_2:", p_2, " p_3:", p_3)
print("\nSelecting the type of refinement: \nPlease enter 'p' for p_refinement and 'h' for h_refinement")
refine_type = input()
......@@ -63,9 +64,9 @@ if refine_type =="p":
mesh = gsm.mesh_func(surfs, u_manual, v_manual, c_order_u, c_order_v)
element_boundaries_u = mesh[0]
element_boundaries_v = mesh[1]
print('elment bc in u', element_boundaries_u)
print('\nelement bc in v', element_boundaries_v)
print("element boundaries: ", element_boundaries_u, element_boundaries_v)
# print('elment bc in u', element_boundaries_u)
# print('\nelement bc in v', element_boundaries_v)
# print("element boundaries: ", element_boundaries_u, element_boundaries_v)
order_displm_array = np.zeros((max_order_elem - min_order_elem + 1, 2))
time_assembling = np.zeros((max_order_elem - min_order_elem + 1, 2))
time_solver = np.zeros((max_order_elem - min_order_elem + 1, 2))
......@@ -81,8 +82,6 @@ if refine_type =="p":
bc = gsm.global_boundary_condition(lobatto_pw, bc_h_bott, bc_h_top,\
bc_v_left, bc_v_right, element_boundaries_u,\
element_boundaries_v)
print('bc', bc)
print("\n\n\nOrder of elements: {}".format(str(i_main)))
print("\nAssembling global stiffness matrix ...")
t_1_assembling = time.perf_counter()
......@@ -97,7 +96,6 @@ if refine_type =="p":
element_boundaries_v)
global_load_bc = np.delete(global_load, bc, 0)
# print('bc:', bc)
print(global_load)
print('resultant load',np.sum(global_load))
print("\nLinear solver in action! ...")
t_1_solver = time.perf_counter()
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment