Mi/.SSKrSSKJrJr SSKJrJr SSKJ r SSKJrJr SSKJr SS KJr SS KJr SS KJr SS jrSS jrSSjrSSjrSS\S\S\\4SjjrSrSr g!\ ar SSK J r \ "\ 5r\\\srr rSr C NsSr C ff=f)N) ArrayLikeOptional) coo_matrixeye)spsolve)ExceptionWrapper)graph triangles)Trimesh) index_sparsemass_properties)unitizectUc [U5nU(a URnURR5R [ R 5nURR5R [ R 5nU(a,URSn [U 5U[U 5U- --n [U5H`n U(dURU5U- n XqU -- nO [W U5nU(dM;[R"XxSS9Sn UWU - S--nMb XplU$)a* Smooth a mesh in-place using laplacian smoothing. Articles 1 - "Improved Laplacian Smoothing of Noisy Surface Meshes" J. Vollmer, R. Mencl, and H. Muller 2 - "Implicit Fairing of Irregular Meshes using Diffusion and Curvature Flow". M. Desbrun, M. Meyer, P. Schroder, A.H.B. Caltech Parameters ------------ mesh : trimesh.Trimesh Mesh to be smoothed in place lamb : float Diffusion speed constant If 0.0, no diffusion If > 0.0, diffusion occurs implicit_time_integration: boolean if False: explicit time integration -lamb <= 1.0 - Stability Limit (Article 1) if True: implicit time integration -lamb no limit (Article 2) iterations : int Number of passes to run filter laplacian_operator : None or scipy.sparse.coo.coo_matrix Sparse matrix laplacian operator Will be autogenerated if None rT skip_inertiavolumegUUUUUU?)laplacian_calculationrverticescopyviewnpndarrayfacesshaperrangedotrr r)meshlamb iterationsimplicit_time_integrationvolume_constraintlaplacian_operatorvol_inirrdlapAA_indexrvol_news k/var/www/eduai.edurigo.com/storigo/production/storigo_env/lib/python3.13/site-packages/trimesh/smoothing.pyfilter_laplacianr+s!J!248++}}!!#((4H JJOO  " "2:: .E!!''* YT-?!?@ @ #($((2X=C s "Hr8,H  //dSG 7* : :H#$(M KcUc [U5nURR5R[R 5nUR5n[ U5HSnUR5nURU5nXQU-SU- U--- n XRU -SU- URU 5---nMU XPlU$)az Smooth a mesh in-place using laplacian smoothing and Humphrey filtering. Articles "Improved Laplacian Smoothing of Noisy Surface Meshes" J. Vollmer, R. Mencl, and H. Muller Parameters ------------ mesh : trimesh.Trimesh Mesh to be smoothed in place alpha : float Controls shrinkage, range is 0.0 - 1.0 If 0.0, not considered If 1.0, no smoothing beta : float Controls how aggressive smoothing is If 0.0, no smoothing If 1.0, full aggressiveness iterations : int Number of passes to run filter laplacian_operator : None or scipy.sparse.coo.coo_matrix Sparse matrix laplacian operator Will be autogenerated if None ?rrrrrrrr) ralphabetar!r$roriginalr(vert_qvert_bs r*filter_humphreyr5cs4!248}}!!#((4H}}H #%))(3X-u0FFG6MS4Z3E3I3I&3Q$QQQ $M Kr,cUc [U5nURR5R[R 5n[ U5H0nURU5U- nUS-S:Xa XQU-- nM)XRU--nM2 XPlU$)aL Smooth a mesh in-place using laplacian smoothing and taubin filtering. Articles "Improved Laplacian Smoothing of Noisy Surface Meshes" J. Vollmer, R. Mencl, and H. Muller Parameters ------------ mesh : trimesh.Trimesh Mesh to be smoothed in place. lamb : float Controls shrinkage, range is 0.0 - 1.0 nu : float Controls dilation, range is 0.0 - 1.0 Nu shall be between 0.0 < 1.0/lambda - 1.0/nu < 0.1 iterations : int Number of passes to run the filter laplacian_operator : None or scipy.sparse.coo.coo_matrix Sparse matrix laplacian operator Will be autogenerated if None rr/)rr nur!r$rindexrs r* filter_taubinr:s.!248}}!!#((4Hz" $$X.9 19> s "H S H#M Kr,c Uc [U5nU(a URnURR5R [ R 5nURR5R [ R 5nS[ R"UR5S--n[U5GHn [U5n URU5n Xk- n [ R"X-R[ R"S555n S[ R"SU 5- n [ R"SU-[ R "SX-[ R""U 5- 55nURU5nXnX- -- nU(dM[%XgSS9S nU S :Xa['XgU UU5nXjW-WU- -- nGM X`lU$) a Smooth a mesh in-place using laplacian smoothing using a mutable diffusion laplacian. Articles Barroqueiro, B., Andrade-Campos, A., Dias-de-Oliveira, J., and Valente, R. (January 21, 2021). "Bridging between topology optimization and additive manufacturing via Laplacian smoothing." ASME. J. Mech. Des. Parameters ------------ mesh : trimesh.Trimesh Mesh to be smoothed in place lamb : float Diffusion speed constant If 0.0, no diffusion If > 0.0, diffusion occurs iterations : int Number of passes to run filter laplacian_operator : None or scipy.sparse.coo.coo_matrix Sparse matrix laplacian operator Will be autogenerated if None g{Gz??)rr.g-q=g?Trrr)rrrrrrrrmax area_facesrget_vertices_normalsrabsonesmaximumminimummeanr dilate_slope)rr r!r#r$v_inirrepsr(normalsqipi_qiadillamberrvolslopes r*filter_mut_dif_laplacianrPs<!248 }}!!#((4H JJOO  " "2:: .E "&&)S0 1C #&t,  # #H - vvw++BGGFO<=RZZt,,C$J 3 bggdm8S(TU!$$X.cn--  !(/EhOC{$XgsCH %53;7 7H'$,M Kr,r equal_weightpinned_verticesc&U(Ga [R"UR5nUbSUR[R "UR U5'[R"UR U45Ul[R"URU45Ul[R"UR[R"[U5[S945UlURSURSS9- 5nU$URnUbUHnU/XE'M URR![R"5n[R"U5n[R"[%U5VVs/sHupXU/[U5-PM snn5n [R"S5n [%U5VVs/sHRupXS[R&"S[R("[R*"XeXh- S-U 555- PMT n nn[R"U Vs/sHoUUR5- PM sn5n [-XU44[U5/S-S 9nU$s snnfs snnfs snf) a( Calculate a sparse matrix for laplacian operations. Note that setting equal_weight to False significantly hampers performance. Parameters ------------- mesh : trimesh.Trimesh Input geometry equal_weight : bool If True, all neighbors will be considered equally If False, all neighbors will be weighted by inverse distance pinned_vertices : None or list of ints If None, no vertices are pinned If list, vertices will be pinned, such that they will not be moved Returns ---------- laplacian : scipy.sparse.coo.coo_matrix Laplacian operator F)dtyper.r)axisr=gư>r7)r)r edges_to_cooedgesdatarisinrow concatenatecolrBlenboolmultiplysumvertex_neighborsrrr enumeraterCsqrtrr) rrQrR laplacian neighborsirr\nrZrBnormsrXs r*rrs2&&tzz2  &FKINN2779==/B CNNIMM?+KLIMNNIMM?+KLIM^^_)=T!JKIN&&sY]]]-B'BC D =))   &$ !s %==%%bjj1nnY'nny7KL7KtqqcCFl7KLMwwqz "), - jjrwwrvvx{X[/HQ.NPT'UVW X,  ~~E:Eq1557{E:;Cj1#h-19LM #M  ;s J AJJcURnURnURn[[ U5U5R U5n[ U5$)z Compute Vertex normals using equal weighting of neighbors faces. Parameters ------------- mesh : trimesh.Trimesh Input geometry Returns ---------- vertices_normals: array Vertices normals )rr face_normalsr r]rr)rrrrj vert_normalss r*r@r@QsK}}H JJE$$L H u599,GL <  r,c<XU--n[XQSS9SnXFU- - $)az Get the derivate of dilation scalar by the volume variation by finite differences Thus, Vertices += vertex_normals*dilate_slope*(Initial_Volume - Srinked_Volume) Parameters ------------- mesh : trimesh.Trimesh Input geometry vertices: mesh.vertices faces: mesh.faces normals: array vertices normals Returns ---------- dilate_slope: float derivative Trrr)rrrIvrH vertices2v2s r*rFrFks1&S=(I ) =h GB F r,)r< FTN)g?r<rpN)r<r<rpN)r<rpTN)TN)!numpyrtypedrr scipy.sparserrscipy.sparse.linalgr ImportErrorE exceptionsr wrapperr r baser geometryr rutilrr+r5r:rPr^rr@rFr,r*r~s&9,+"& # K\+\)ZOSCP+/K KKi(K\!4G 9,q!G&C* 9sA//B5BB