def computeSurfaces(self): geometry = self.geometry oml0 = geometry.oml0 quad = self.quad premeshFaces = self.premeshFaces self.surfaceNames = [] B0 = [] quads0 = [] nnode0 = [0] iList = 0 for k in range(len(geometry.comps)): comp = geometry.comps[geometry.keys[k]] for f in range(len(comp.Ks)): verts,edges,edge_group,edgeLengths = premeshFaces[iList] iList += 1 nvert = verts.shape[0] nedge = edges.shape[0] ni, nj = comp.Ks[f].shape idims, jdims = self.faceDims[k][f] print 'Computing skin elements:', geometry.keys[k], f for i in range(ni): for j in range(nj): surf = comp.Ks[f][i,j] if oml0.visible[surf]: nedge1 = PSMlib.countsurfaceedges(nvert, nedge, idims[i], idims[i+1], jdims[j], jdims[j+1], verts, edges) edges1 = PSMlib.computesurfaceedges(nvert, nedge, nedge1, idims[i], idims[i+1], jdims[j], jdims[j+1], verts, edges) print geometry.keys[k], f, i, j quad.importEdges(edges1) if geometry.keys[k]=='fu' and f==1 and i==0 and j==0: nodes, quads = quad.mesh(self.maxL, self.surfEdgeLengths[surf,:,:])#,True,True) else: nodes, quads = quad.mesh(self.maxL, self.surfEdgeLengths[surf,:,:]) mu, mv = oml0.edgeProperty(surf,1) for u in range(mu): for v in range(mv): oml0.C[oml0.getIndex(surf,u,v,1),:3] = [u/(mu-1), v/(mv-1), 0] oml0.computePointsC() P0, Q = PSMlib.computesurfaceprojections(nodes.shape[0], nodes) s,u,v = oml0.evaluateProjection(P0, [surf], Q) B = oml0.evaluateBases(s,u,v) name = geometry.keys[k] + ':' + str(f) + '-' + str(i) + '-' + str(j) self.surfaceNames.append(name) B0.append(B) nnode0.append(nnode0[-1] + P0.shape[0]) quads0.append(quads) B0 = scipy.sparse.vstack(B0) self.meshS = [B0, quads0, nnode0]
def computeSurfaces(self): geometry = self.geometry oml0 = geometry.oml0 quad = self.quad premeshFaces = self.premeshFaces self.surfaceNames = [] B0 = [] quads0 = [] nnode0 = [0] iList = 0 for comp in geometry.comps.values(): for face in comp.faces.values(): verts,edges,edge_group,edgeLengths = premeshFaces[iList] iList += 1 nvert = verts.shape[0] nedge = edges.shape[0] ni, nj = face.num_surf idims, jdims = self.faceDims[comp.name][face.name] print 'Computing skin elements:', comp.name, face.name for i in range(ni): for j in range(nj): surf = face.surf_indices[i,j] if surf >= 0:# and oml0.visible[surf]: nedge1 = PSMlib.countsurfaceedges(nvert, nedge, idims[i], idims[i+1], jdims[j], jdims[j+1], verts, edges) edges1 = PSMlib.computesurfaceedges(nvert, nedge, nedge1, idims[i], idims[i+1], jdims[j], jdims[j+1], verts, edges) print comp.name, face.name, i, j quad.importEdges(edges1) output = False if comp.name=='lw' and face.name=='upp' and i==0 and j==2: output = True nodes, quads = quad.mesh(self.maxL, self.surfEdgeLengths[surf,:,:], output, output) mu, mv = oml0.edgeProperty(surf,1) for u in range(mu): for v in range(mv): oml0.C[oml0.getIndex(surf,u,v,1),:3] = [u/(mu-1), v/(mv-1), 0] oml0.computePointsC() P0, Q = PSMlib.computesurfaceprojections(nodes.shape[0], nodes) s,u,v = oml0.evaluateProjection(P0, [surf], Q) B = oml0.evaluateBases(s,u,v) name = comp.name + ':' + str(face.name) + ':' + str(i) + ':' + str(j) self.surfaceNames.append(name) B0.append(B) nnode0.append(nnode0[-1] + P0.shape[0]) quads0.append(quads) B0 = scipy.sparse.vstack(B0) self.meshS = [B0, quads0, nnode0]
def computeSurfaces(self): geometry = self.geometry bse = geometry._bse quad = self.quad premeshFaces = self.premeshFaces self.surfaceNames = [] B0 = [] quads0 = [] nnode0 = [0] mem0 = [] ucoord0 = [] vcoord0 = [] iList = 0 imem = 0 for comp in geometry.comps.values(): for face in comp.faces.values(): verts, edges, edge_group, edgeLengths = premeshFaces[iList] iList += 1 nvert = verts.shape[0] nedge = edges.shape[0] ni, nj = face._num_surf['u'], face._num_surf['v'] idims, jdims = self.faceDims[comp._name][face._name] print 'Computing skin elements:', comp._name, face._name for i in range(ni): for j in range(nj): surf = face._surf_indices[i, j] if surf >= 0: nedge1 = PSMlib.countsurfaceedges( nvert, nedge, idims[i], idims[i + 1], jdims[j], jdims[j + 1], verts, edges) edges1 = PSMlib.computesurfaceedges( nvert, nedge, nedge1, idims[i], idims[i + 1], jdims[j], jdims[j + 1], verts, edges) print comp._name, face._name, i, j quad.importEdges(edges1) output = False if comp._name == 'lw' and face._name == 'upp' and i == 0 and j == 2: output = True nodes, quads = quad.mesh( self.maxL, self.surfEdgeLengths[surf, :, :], output, output) mu = bse.get_bspline_option('num_cp', surf, 'u') mv = bse.get_bspline_option('num_cp', surf, 'v') nu = bse.get_bspline_option('num_pt', surf, 'u') nv = bse.get_bspline_option('num_pt', surf, 'v') for u in range(mu): for v in range(mv): bse.vec['cp_str'](surf)[u, v, :] = [ u / (mu - 1), v / (mv - 1), 0 ] for u in range(nu): for v in range(nv): bse.vec['pt_str'](surf)[u, v, :] = [ u / (nu - 1), v / (nv - 1), 0 ] P0, Q = PSMlib.computesurfaceprojections( nodes.shape[0], nodes) bse.compute_projection('temp', P0, [surf], ndim=3) B = bse.jac['d(temp)/d(cp_str)'] name = '{}:{}:{:03d}:{:03d}'.format( comp._name, face._name, i, j) self.surfaceNames.append(name) B0.append(B) nnode0.append(nnode0[-1] + P0.shape[0]) quads0.append(quads) mem0.append(imem * numpy.ones(P0.shape[0])) imem += 1 ucoord0.append(P0[:, 0]) vcoord0.append(P0[:, 1]) bse.apply_jacobian('cp_str', 'd(cp_str)/d(cp)', 'cp') B0 = scipy.sparse.vstack(B0) self.meshS = [B0, quads0, nnode0, mem0, ucoord0, vcoord0]
def computeSurfaces(self): geometry = self.geometry oml0 = geometry.oml0 quad = self.quad premeshFaces = self.premeshFaces self.surfaceNames = [] B0 = [] quads0 = [] nnode0 = [0] iList = 0 for comp in geometry.comps.values(): for face in comp.faces.values(): verts, edges, edge_group, edgeLengths = premeshFaces[iList] iList += 1 nvert = verts.shape[0] nedge = edges.shape[0] ni, nj = face.num_surf idims, jdims = self.faceDims[comp.name][face.name] print 'Computing skin elements:', comp.name, face.name for i in range(ni): for j in range(nj): surf = face.surf_indices[i, j] if surf >= 0: # and oml0.visible[surf]: nedge1 = PSMlib.countsurfaceedges( nvert, nedge, idims[i], idims[i + 1], jdims[j], jdims[j + 1], verts, edges) edges1 = PSMlib.computesurfaceedges( nvert, nedge, nedge1, idims[i], idims[i + 1], jdims[j], jdims[j + 1], verts, edges) print comp.name, face.name, i, j quad.importEdges(edges1) output = False if comp.name == 'lw' and face.name == 'upp' and i == 0 and j == 2: output = True nodes, quads = quad.mesh( self.maxL, self.surfEdgeLengths[surf, :, :], output, output) mu, mv = oml0.edgeProperty(surf, 1) for u in range(mu): for v in range(mv): oml0.C[ oml0.getIndex(surf, u, v, 1), :3] = [ u / (mu - 1), v / (mv - 1), 0 ] oml0.computePointsC() P0, Q = PSMlib.computesurfaceprojections( nodes.shape[0], nodes) s, u, v = oml0.evaluateProjection(P0, [surf], Q) B = oml0.evaluateBases(s, u, v) name = comp.name + ':' + str( face.name) + ':' + str(i) + ':' + str(j) self.surfaceNames.append(name) B0.append(B) nnode0.append(nnode0[-1] + P0.shape[0]) quads0.append(quads) B0 = scipy.sparse.vstack(B0) self.meshS = [B0, quads0, nnode0]
def computeSurfaces(self): geometry = self.geometry bse = geometry._bse quad = self.quad premeshFaces = self.premeshFaces self.surfaceNames = [] B0 = [] quads0 = [] nnode0 = [0] mem0 = [] ucoord0 = [] vcoord0 = [] iList = 0 imem = 0 for comp in geometry.comps.values(): for face in comp.faces.values(): verts,edges,edge_group,edgeLengths = premeshFaces[iList] iList += 1 nvert = verts.shape[0] nedge = edges.shape[0] ni, nj = face._num_surf['u'], face._num_surf['v'] idims, jdims = self.faceDims[comp._name][face._name] print 'Computing skin elements:', comp._name, face._name for i in range(ni): for j in range(nj): surf = face._surf_indices[i,j] if surf >= 0: nedge1 = PSMlib.countsurfaceedges(nvert, nedge, idims[i], idims[i+1], jdims[j], jdims[j+1], verts, edges) edges1 = PSMlib.computesurfaceedges(nvert, nedge, nedge1, idims[i], idims[i+1], jdims[j], jdims[j+1], verts, edges) print comp._name, face._name, i, j quad.importEdges(edges1) output = False if comp._name=='lw' and face._name=='upp' and i==0 and j==2: output = True nodes, quads = quad.mesh(self.maxL, self.surfEdgeLengths[surf,:,:], output, output) mu = bse.get_bspline_option('num_cp', surf, 'u') mv = bse.get_bspline_option('num_cp', surf, 'v') nu = bse.get_bspline_option('num_pt', surf, 'u') nv = bse.get_bspline_option('num_pt', surf, 'v') for u in range(mu): for v in range(mv): bse.vec['cp_str'](surf)[u, v, :] = [u/(mu-1), v/(mv-1), 0] for u in range(nu): for v in range(nv): bse.vec['pt_str'](surf)[u, v, :] = [u/(nu-1), v/(nv-1), 0] P0, Q = PSMlib.computesurfaceprojections(nodes.shape[0], nodes) bse.compute_projection('temp', P0, [surf], ndim=3) B = bse.jac['d(temp)/d(cp_str)'] name = comp._name + ':' + str(face._name) + ':' + str(i) + ':' + str(j) self.surfaceNames.append(name) B0.append(B) nnode0.append(nnode0[-1] + P0.shape[0]) quads0.append(quads) mem0.append(imem*numpy.ones(P0.shape[0])) imem += 1 ucoord0.append(P0[:,0]) vcoord0.append(P0[:,1]) bse.apply_jacobian('cp_str', 'd(cp_str)/d(cp)', 'cp') B0 = scipy.sparse.vstack(B0) self.meshS = [B0, quads0, nnode0, mem0, ucoord0, vcoord0]