☰

Рисование фолиации Риба

tag-icon • tag-icon asymptote
Рисование фолиации Риба

Я пытаюсь нарисоватьшелушение Риба. Пока что у меня есть код и изображение ниже. Проблема в том, что одна поверхность немного выпирает через ту, что перед ней (через самую переднюю поверхность). Они должны быть вложены каждая полностью внутрь предыдущей, как изогнутые цветочные горшки.

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
//  );
}

введите описание изображения здесь

Возможно, проблема станет еще яснее, если я нарисую только две последние поверхности разными цветами:

введите описание изображения здесь

решение1

Это только из-за численного приближения вашей поверхности с помощью патчей Безье. По умолчанию у вас есть 10 сеток (или 10x10). В зависимости от вашего ПК вы можете независимо увеличить количество сеток относительно двух переменных. Здесь я протестировал последние две конечные поверхности с

draw(
     rotate(-anglestep*z,Z)*surface(
                                    h,
                                    (0,0),
                                    (rcutoff,2*pi),nu=16,nv=10,Spline
                                    ),lightgreen+opacity(1)
     );

и вроде бы все в порядке.

Связанный контент