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 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
