
私は、リーブ葉状構造. これまでのところ、コードと画像は以下のとおりです。問題は、1 つのサーフェスが、その前のサーフェス (最前面のサーフェス) からわずかに突き出ていることです。湾曲した植木鉢のように、各サーフェスを前のサーフェス内に完全にネストする必要があります。
settings.outformat="pdf";
settings.prc=false;
settings.render=4;
import graph;
import graph3;
import three;
size(12cm,0);
real rcutoff=.8;
real rinner=1;
real router=2;
currentprojection=orthographic(1,-2,1);
currentlight=(1,-1,0.5);
real f(pair z) {return 1/(1-abs(z)^2);}
// We change variables so that the unit cylinder x^2+y^2<1 gets mapped to the torus.
triple changevariables(triple X)
{
return (
(router+rinner*X.x)*cos(X.z),
(router+rinner*X.x)*sin(X.z),
rinner*X.y
);
}
// g returns the location of a point on a ``paraboloid'' inside the unit cylinder with given angle and radius above the unit disk in the (x,y) plane.
triple g(pair p){
real x=p.x*cos(p.y), y=p.x*sin(p.y);
return (x,y,f((x,y)));
}
triple h(pair p){
return changevariables(g(p));
}
real anglestep=10;
for (int z = 0; z <= 20; ++z) {
draw(
rotate(-anglestep*z,Z)*surface(
h,
(0,0),
(rcutoff,2*pi),Spline
),lightgray+opacity(1)
);
// real zz=f((rcutoff,0));
// triple cent=changevariables((0,0,zz));
// draw(
// rotate(-anglestep,Z)*
// circle(c=cent, r=rcutoff,normal=(-router*sin(zz),router*cos(zz),0)),
// white
// );
}
最後の 2 つの面だけを異なる色で描画すると、問題はさらに明確になるかもしれません。
答え1
これは、ベジェパッチによる表面の数値近似によるものです。デフォルトでは、10個のメッシュ(または10x10)があります。PCに応じて、2つの変数に関してメッシュの数を個別に増やすことができます。ここでは、最後の2つの最終表面をテストしました。
draw(
rotate(-anglestep*z,Z)*surface(
h,
(0,0),
(rcutoff,2*pi),nu=16,nv=10,Spline
),lightgreen+opacity(1)
);
大丈夫そうです。