Skip to content
Snippets Groups Projects
Commit 606aaa81 authored by David's avatar David
Browse files

delete five patch mesh definition from example script

parent 16ea3a3e
No related branches found
No related tags found
No related merge requests found
...@@ -162,103 +162,6 @@ subdomain_def_points = substructuring.subdomain_def_points ...@@ -162,103 +162,6 @@ subdomain_def_points = substructuring.subdomain_def_points
outer_boundary_def_points = substructuring.outer_boundary_def_points outer_boundary_def_points = substructuring.outer_boundary_def_points
# MODEL CONFIGURATION ######################################################### # MODEL CONFIGURATION #########################################################
# # subdomain1.
# sub_domain1_vertices = [interface23_vertices[0],
# interface23_vertices[1],
# interface12_vertices[1],
# sub_domain0_vertices[2],
# df.Point(0.0, 1.0)]
#
# # vertex coordinates of the outer boundaries. If it can not be specified as a
# # polygon, use an entry per boundary polygon. This information is used for
# # defining the Dirichlet boundary conditions. If a domain is completely inter-
# # nal, the dictionary entry should be 0: None
# subdomain1_outer_boundary_verts = {
# 0: [interface12_vertices[1],
# sub_domain0_vertices[2],
# df.Point(0.0, 1.0)]
# }
# # subdomain2
# sub_domain2_vertices = [interface24_vertices[1],
# sub_domain0_vertices[1],
# interface12_vertices[1],
# interface23_vertices[1],
# interface23_vertices[0]]
#
# subdomain2_outer_boundary_verts = {
# 0: [interface24_vertices[1],
# sub_domain0_vertices[1],
# interface12_vertices[1]]
# }
#
# sub_domain3_vertices = [interface23_vertices[0],
# interface23_vertices[1],
# interface13_vertices[0]]
#
# subdomain3_outer_boundary_verts = None
#
#
# sub_domain4_vertices = [sub_domain0_vertices[0],
# interface24_vertices[1],
# interface34_vertices[1],
# interface34_vertices[0],
# interface45_vertices[0]]
#
# subdomain4_outer_boundary_verts = {
# 0: [interface45_vertices[0],
# sub_domain0_vertices[0],
# interface24_vertices[1]]
# }
#
# sub_domain5_vertices = [interface45_vertices[0],
# interface15_vertices[0],
# interface15_vertices[1],
# sub_domain0_vertices[3]]
#
# subdomain5_outer_boundary_verts = {
# 0: [interface15_vertices[1],
# sub_domain0_vertices[3],
# interface45_vertices[0]]
# }
#
# # list of subdomains given by the boundary polygon vertices.
# # Subdomains are given as a list of dolfin points forming
# # a closed polygon, such that mshr.Polygon(subdomain_def_points[i]) can be used
# # to create the subdomain. subdomain_def_points[0] contains the
# # vertices of the global simulation domain and subdomain_def_points[i] contains
# # the vertices of the subdomain i.
# subdomain_def_points = [sub_domain0_vertices,
# sub_domain1_vertices,
# sub_domain2_vertices,
# sub_domain3_vertices,
# sub_domain4_vertices,
# sub_domain5_vertices]
# # in the below list, index 0 corresponds to the 12 interface which has global
# # marker value 1
# interface_def_points = [interface13_vertices,
# interface12_vertices,
# interface23_vertices,
# interface24_vertices,
# interface34_vertices,
# interface45_vertices,
# interface15_vertices,]
#
# # adjacent_subdomains[i] contains the indices of the subdomains sharing the
# # interface i (i.e. given by interface_def_points[i]).
# adjacent_subdomains = [[1, 3], [1, 2], [2, 3], [2, 4], [3, 4], [4, 5], [1, 5]]
#
# # if a subdomain has no outer boundary write None instead, i.e.
# # i: None
# # if i is the index of the inner subdomain.
# outer_boundary_def_points = {
# # subdomain number
# 1: subdomain1_outer_boundary_verts,
# 2: subdomain2_outer_boundary_verts,
# 3: subdomain3_outer_boundary_verts,
# 4: subdomain4_outer_boundary_verts,
# 5: subdomain5_outer_boundary_verts
# }
isRichards = { isRichards = {
1: True, 1: True,
2: False, 2: False,
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment