ljj>dZddlZddlmZddlZejdddlmZddl Z ddl Z ddl m Z ejdeeejejddddlZddlmZmZmZmZmZmZejdZgd gd d Zeejd kr ejd nd Z e j!gde j"#gdz Z$ddZ%dZ&ddZ'dZ(e)dkr e(dSdS)uThe clean bowl_gso mesh at 256 layers: does its dimple texture (about 50 periods around, ~0.1 % of the radius) come through? Variants of N (harmonics per ring) and M (Chebyshev degree) with the atom cap raised for the experiment (the package's MAX_ATOMS is 6144). Per variant: atoms stored, Chamfer (‰ of the diagonal, against the mesh), the amplitude of harmonics 45-55 of a decoded mid-height section against GT's; renders: shaded (smooth shading from vertex normals; GT subdivided twice so its large triangles do not show) and coloured by distance to GT. report/texture_bowl_.html. env -u PYTHONPATH PYTHONWARNINGS=ignore .venv/bin/python ~/code/data/meshes/fcp_2710_tests/scripts/texture_bowl.py N)PathAgg)Poly3DCollectionz/home/joao/code/fcp_devel/api)CSSIMGOBJREPORTUPstatsz fcp.encode)))r @rr)r rr`r r0r)r)rr)rrr@ )rrrr)256512r)g333333?g?-8c |dd|fd tjSt fd jDt ddddf}|d}tj|dddf|dz |dddf|dz }tj ||z d }tj |}tj tj tj d d } tj | ||||dtj z } tjtj| | z t | z dz} t#tj| ||dz|z dzS)zdAmplitude (% of the mean radius) of harmonics lo..hi-1 of the largest loop of section z (canonical).rrrr) plane_origin plane_normalNc3<K|]}j|jVdSN)verticespoints).0ess B/home/joao/code/data/meshes/fcp_2710_tests/scripts/texture_bowl.py zband..$s, 6 6aQZ ! 6 6 6 6 6 6)keyr)axisiF)endpoint)periodd)sectionnpnanmaxentitieslenmeanarctan2linalgnormargsortlinspacepiinterpabsfftrfftfloatsqrtsum) GzlohiPcaroagrgAr's @r(bandrQs 1ayy AAAyv  6 6 6 61: 6 6 6C@@@BQBGA q A 1QQQT7QqT>1QQQT7QqT>22A q1u1%%A 1 A beVRUD5 9 9 9B 2qtQqT!be) 4 4 4B rv{{2 >**++c"gg59A !BrE(a-,,..//!&&((:S@ A AAr*g?c |jdddfdk}|j|}|j|d}|t jtj|dddz}ddt j|tzd dzz}|tj |||d zf}nqtj d t j||dtz d ddddd fd d|zzdddfz}|t|j||d ddd d|dd|dd|dd|d|dd|||ddS)zqSmooth shading: each face lit by the mean of its vertices' normals; err (per vertex, % of the diagonal): colours.Nr皙?T)r-keepdimsg-q=g333333?g333333?rg{Gz?magma_rg?g?noneF) facecolors edgecolors linewidths antialiased)rrri )fontsize) face_normalsfacesvertex_normalsr7r2maximumr9r:clipLIGHTc_pltget_cmapEMAXadd_collection3drr#set_xlimset_ylimset_zlimset_box_aspect view_init set_axis_off set_title) axmtitleerrfrontFnshadefcs r(showr{2s N111a4 4 'E A   ##ABINN11tN<,>,>i   ",,q#"6"6"68I8I8I2<|j>d?d@g}t~d5tdAz @dBA|tG|dS)DNbowl_gsomesh)forceiPr)seedrcuda)atomslayersupdeviceNM)rF)processc0g|]}t|SrQ)r%rFrEs r( zmain..Os!777$q!**777r*)g333333ӿgg?rzGT harmonics 45-55 = z.3fz % of rGT)variantrchamferrQc0g|]}t|Srr)r%rFCs r(rzmain..Ts!:::AT!QZZ:::r*g@@z layers, N=z, M=z, z atomsrVr0z Chamfer u‰, harmonics 45-55 = z %r\T)flushr)Imagec tjdd}t|ddddd||||jtj|j ddddddf }tj ||S) N)g@nfigsizedpir3dortho) projection proj_typerV) rgfigurer{ add_subplotcanvasdraw fromarrayr2asarray buffer_rgbacopyclose)rstrufigimrs r(tilezmain..tile^sjs333 S__Q1_ I I1aQTUUU  __RZ (>(>(@(@AA!!!QQQ(KPPRR S S # r*c0g|]\}}}||Srr)r%rsr_rs r(rzmain..fs) , , ,'!Q441:: , , ,r*cLg|] \}}}||dtd|nd!S)Nzdistance to GT, 0-z % of the diagonal)ri)r%rsrrurs r(rzmain..gsU $ $ $q!SKN/44@@@@# F F F_c $ $ $r*)rg?rr)rSrg?g?rU)cmapr: horizontalz!distance to GT, % of the diagonal)cax orientationlabelRGBr,white texture_bowl_z.pngzDBowl texturez"

bowl_gso: the dimple texture (z layers set)

z

The texture is about 50 periods around (harmonics 45-55 of a section, about 0.1 % of the radius). N harmonics per ring stores harmonics 1..N only. The atom cap is raised beyond the package's 6144 for the experiment.

)indexzzz.html r")Btrimeshloadrr2rsampleVARIANTSSETENC MAX_ATOMSfcpencoder getdecodeTrimeshpose to_canonicalr#ra subdividerBr9r:ptpappendr7dictr r6rroundr? proximity closest_pointprintPILrrgradd_axescolorbarcmScalarMappable Normalizerirrrrrrsizenew enumeratepastersavepandasr DataFrameto_htmlnamer write_textjoin)"grefoutpanelsrrrcapdrsgtGsdiagbrJlabrutopbotrrrbarwhsheetkrimgpdhtmlrrErrs" @@@@r(mainrCs  S_F333A *QXXf1X-- . .CbC A%c]##1c JqFrvj)7T7T]cghlm n n n Jq ( ( ( OAF// ;;QWe T T T 9 3 3AJ ? ?RWXXXA7777&6777B((**Bqz1(=(=>>??D MM2P PPPPRVW X X X JJtDBRWR[[QQQ R R R G::::)9::: ; ; !S!  Q # %GGAGG1GGAGAJGGG 43qwqz??E!QKKVWXXXYYYfW&44Q CCAFGG$NQTT qSSSASSS!SSSSUXYZZZ c"gT"""""CM - , , ,V , , ,C $ $ $ $" $ $ $C *X3 / / /C ,,, - -BLL&&ICM!T>qqq!!!RaRxHMMOO P PCIcNNN q6;DAq IIea#c((lAE2G < r s  u77777733ttH~~,--...2333 9999999999999999 k,iiiooo q qS]]Q&&chqkkE!!!""RY^^4D4D4D%E%EE B B B B  kkkk">>>B zDFFFFFr*