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