Find the global nodes index of a boundary mesh

6 months ago by
I am trying to find the index value of a boundary mesh with respect to the original mesh ordering.
from dolfin import *
import time
import numpy as np

#define domain
x_min = y_min = z_min = 0.0
x_max = y_max = z_max = 1.0

#define number of elements
nx = ny = 3
nz = 3

#create mesh
mesh = BoxMesh(Point(x_min, y_min, z_min), Point(x_max, y_max, z_max), nx, ny, nz)

# define function space
V = FunctionSpace(mesh, 'P', 1)

bmesh = BoundaryMesh(mesh, 'exterior')
print('connectivity of mesh',mesh.cells())
print('connectivity of bmesh', bmesh.cells())

is there anyway to find the nodal values of the bmesh nodes as they were defined in mesh?

Thanks in advance,

Community: FEniCS Project
Does this link help maybe?
written 6 months ago by klunkean  
So i figured I would be able to calculate the xyz coordinates of the dofs in my bmesh and then find the node in the mesh that has the same coordinates. But I am having trouble getting the xyz coordinates in the bmesh
written 6 months ago by Terrence  

1 Answer

6 months ago by
I had to write a function correlating the xyz coordinates of the boundary mesh with the coordinates of the global mesh nodes.
Please login to add an answer/comment or follow this question.

Similar posts:
Search »