def __initializePoints(self, P_arrays): """ Rearrange the list of surfaces into a single vector of unique points Input P_arrays: list of ndarrays (nu,nv,3) Each element of the list is an nu x nv array of x-y-z coordinates """ self.P = numpy.zeros((self.nP, 3), order='F') for k in range(self.nsurf): nu = P_arrays[k].shape[0] nv = P_arrays[k].shape[1] PUBSlib.populatep(self.nP, nu, nv, k + 1, self.nsurf, self.nedge, self.ngroup, self.nvert, self.surf_vert, self.surf_edge, self.edge_group, self.group_n, self.surf_index_P, self.edge_index_P, P_arrays[k], self.P) PUBSlib.avgboundaries(self.nP, self.nedge, self.ngroup, self.nvert, self.edge_group, self.group_n, self.vert_count, self.edge_count, self.edge_index_P, self.P) self.vert_symm, self.edge_symm = PUBSlib.determinesymm( self.symmPlane + 1, 1e-10, 1e-8, self.nP, self.nedge, self.ngroup, self.nvert, self.edge_group, self.group_n, self.P) self.P = numpy.hstack([self.P, numpy.zeros((self.nP, 3), order='F')])
def __initializePoints(self, P_arrays): """ Rearrange the list of surfaces into a single vector of unique points * Input * P_arrays: list of ndarrays (nu,nv,3) Each element of the list is an nu x nv array of x-y-z coordinates """ self.P = numpy.zeros((self.nP,3),order='F') for k in range(self.nsurf): nu = P_arrays[k].shape[0] nv = P_arrays[k].shape[1] PUBSlib.populatep(self.nP, nu, nv, k+1, self.nsurf, self.nedge, self.ngroup, self.nvert, self.surf_vert, self.surf_edge, self.edge_group, self.group_n, self.surf_index_P, self.edge_index_P, P_arrays[k], self.P) PUBSlib.avgboundaries(self.nP, self.nedge, self.ngroup, self.nvert, self.edge_group, self.group_n, self.vert_count, self.edge_count, self.edge_index_P, self.P) self.vert_symm, self.edge_symm = PUBSlib.determinesymm(self.symmPlane+1, 1e-10, 1e-8, self.nP, self.nedge, self.ngroup, self.nvert, self.edge_group, self.group_n, self.P) self.P = numpy.hstack([self.P,numpy.zeros((self.nP,3),order='F')])