def computeCindices(self):
""" Compute where each vertex, edge, and surface is located in the global list of Cs """
self.surf_index_C = PUBSlib.computesurfindices(self.nsurf, self.nedge, self.ngroup, self.surf_edge, self.edge_group, self.group_m)
self.edge_index_C = PUBSlib.computeedgeindices(self.nedge, self.ngroup, self.edge_group, self.group_m)
self.nC = self.nvert
self.nC += self.edge_index_C[-1,1]
self.nC += self.surf_index_C[-1,1]
if self.printInfo:
print '# Control points =',self.nC
def computePindices(self):
""" Compute where each vertex, edge, and surface is located in the global list of Ps """
self.surf_index_P = PUBSlib.computesurfindices(self.nsurf, self.nedge, self.ngroup, self.surf_edge, self.edge_group, self.group_n)
self.edge_index_P = PUBSlib.computeedgeindices(self.nedge, self.ngroup, self.edge_group, self.group_n)
self.nT = 0
self.nT += self.edge_index_P[-1,1]
self.nT += 2*self.surf_index_P[-1,1]
self.nP = self.nvert
self.nP += self.edge_index_P[-1,1]
self.nP += self.surf_index_P[-1,1]
if self.printInfo:
print '# Points =',self.nP
def computeQindices(self):
""" Compute where each vertex, edge, and surface is located in the global list of Qs """
self.surf_index_Q = PUBSlib.computesurfindices(self.nsurf, self.nedge, self.ngroup, self.surf_edge, self.edge_group, self.group_m)
self.edge_index_Q = PUBSlib.computeedgeindicesq(self.nsurf, self.nedge, self.ngroup, self.surf_edge, self.edge_group, self.group_m, self.surf_c1)
self.vert_index_Q = PUBSlib.computevertindicesq(self.nsurf, self.nedge, self.nvert, self.surf_vert, self.surf_edge, self.surf_c1, self.edge_c1)
self.nQ = 0
self.nQ += max(self.vert_index_Q)
self.nQ += max(self.edge_index_Q[:,1])
self.nQ += self.surf_index_Q[-1,1]
self.Q = numpy.zeros((self.nQ,self.nvar),order='F')
if self.printInfo:
print '# Degrees of freedom =',self.nQ
def computeCindices(self):
""" Compute where each vertex, edge, and surface is located in the global list of Cs """
self.surf_index_C = PUBSlib.computesurfindices(self.nsurf, self.nedge,
self.ngroup,
self.surf_edge,
self.edge_group,
self.group_m)
self.edge_index_C = PUBSlib.computeedgeindices(self.nedge, self.ngroup,
self.edge_group,
self.group_m)
self.nC = self.nvert
self.nC += self.edge_index_C[-1, 1]
self.nC += self.surf_index_C[-1, 1]
if self.printInfo:
print '# Control points =', self.nC
def computePindices(self):
""" Compute where each vertex, edge, and surface is located in the global list of Ps """
self.surf_index_P = PUBSlib.computesurfindices(self.nsurf, self.nedge,
self.ngroup,
self.surf_edge,
self.edge_group,
self.group_n)
self.edge_index_P = PUBSlib.computeedgeindices(self.nedge, self.ngroup,
self.edge_group,
self.group_n)
self.nT = 0
self.nT += self.edge_index_P[-1, 1]
self.nT += 2 * self.surf_index_P[-1, 1]
self.nP = self.nvert
self.nP += self.edge_index_P[-1, 1]
self.nP += self.surf_index_P[-1, 1]
if self.printInfo:
print '# Points =', self.nP
def computeQindices(self):
""" Compute where each vertex, edge, and surface is located in the global list of Qs """
self.surf_index_Q = PUBSlib.computesurfindices(self.nsurf, self.nedge,
self.ngroup,
self.surf_edge,
self.edge_group,
self.group_m)
self.edge_index_Q = PUBSlib.computeedgeindicesq(
self.nsurf, self.nedge, self.ngroup, self.surf_edge,
self.edge_group, self.group_m, self.surf_c1)
self.vert_index_Q = PUBSlib.computevertindicesq(
self.nsurf, self.nedge, self.nvert, self.surf_vert, self.surf_edge,
self.surf_c1, self.edge_c1)
self.nQ = 0
self.nQ += max(self.vert_index_Q)
self.nQ += max(self.edge_index_Q[:, 1])
self.nQ += self.surf_index_Q[-1, 1]
self.Q = numpy.zeros((self.nQ, self.nvar), order='F')
if self.printInfo:
print '# Degrees of freedom =', self.nQ
