 |
|
|
|
|
|
| |
| |
|
|
|
|
| |
| |
|
|
Le 16/05/2017 à 13:41, Bald Eagle a écrit :
> I didn't have as much time as I would have liked to explore this,
> but after fiddling with the Dupin cyclide in both isosurface (implicit) and
> parametric form, I found the parameters to be highly unintuitive, the desired
> shape very difficult to achieve and control, and the constraints on the
> parameters too complex to be easily implemented.
>
> I hadn't gotten around to unraveling the implicit equation to fit the syntax for
> a polynomial.
I asked the other Internets about that, they had a round tuit left so I
got a nice answer.
> It is, however, SLLLLllllllllllooooooooooowwwwwwww.
It is not slow with a polynomial.
True Duplin Cyclide !
#version 3.7;
global_settings{ assumed_gamma 1.0}
#default { finish { ambient 0.2 diffuse 0.8 } }
#declare HEIGHT=9;
camera { orthographic
location <00,00,30>
up HEIGHT*y
right HEIGHT*image_width/image_height*x
look_at 0
rotate clock*x*90
}
light_source { <20,30,40>,1 }
#declare TC=texture { pigment { color srgb <1,1,1> } };
sphere { 0, 1 texture { TC }}
#declare A=3;
#declare B=A*0.98;
#declare C=sqrt(A*A-B*B);
#declare D=1.3*C;
#declare T=texture { pigment { color srgb <1,0.75,0.5> } };
polynomial{
4,
/* expand (x^2+y^2+z^2+b^2-d^2)^2-4*(a*x-c*d)^2-4*b^2*y^2=0 */
/* -4 a^2 x^2 + 2 b^2 x^2 - 2 d^2 x^2
+ 8 a c d x
+ b^4 - 2 b^2 d^2 - 4 c^2 d^2 + d^4
- 2 b^2 y^2 - 2 d^2 y^2
+ 2 b^2 z^2 - 2 d^2 z^2
+ 2 x^2 y^2
+ 2 x^2 z^2
+ 2 y^2 z^2
+ x^4
+ y^4
+ z^4
= 0
*/
xyz(2,0,0):(-4*A*A)+(2*B*B)-(2*D*D),
xyz(1,0,0):8*A*C*D,
xyz(0,0,0):(B*B*B*B)-(2*B*B*D*D)-(4*C*C*D*D)+(D*D*D*D),
xyz(0,2,0):(-2*B*B)-(2*D*D),
xyz(0,0,2):(2*B*B)-(2*D*D),
xyz(2,2,0):2,
xyz(2,0,2):2,
xyz(0,2,2):2,
xyz(4,0,0):1,
xyz(0,4,0):1,
xyz(0,0,4):1
texture { T }}
Post a reply to this message
Attachments:
Download 'cyclide4.png' (15 KB)
Download 'cyclide3.png' (39 KB)
Download 'cyclide2.png' (59 KB)
Download 'cyclide1.png' (73 KB)
Download 'cyclide0.png' (80 KB)
Preview of image 'cyclide4.png'
Preview of image 'cyclide3.png'
Preview of image 'cyclide2.png'
Preview of image 'cyclide1.png'
Preview of image 'cyclide0.png'
|
|
| |
| |
|
|
|
|
| |
| |
|
|
Le_Forgeron <jgr### [at] free fr> wrote:
> Le 16/05/2017 à 13:41, Bald Eagle a écrit :
> > I didn't have as much time as I would have liked to explore this,
> > but after fiddling with the Dupin cyclide in both isosurface (implicit) and
> > parametric form, I found the parameters to be highly unintuitive, the desired
> > shape very difficult to achieve and control, and the constraints on the
> > parameters too complex to be easily implemented.
> >
> > I hadn't gotten around to unraveling the implicit equation to fit the syntax for
> > a polynomial.
> I asked the other Internets about that, they had a round tuit left so I
> got a nice answer.
>
> > It is, however, SLLLLllllllllllooooooooooowwwwwwww.
> It is not slow with a polynomial.
I had wondered / suspected / hoped that were the case.
> True Duplin Cyclide !
That, sir, is a very beautiful thing :)
You did indeed get a very nice answer.
I'm wondering how a, b, c, and d relate to the radii of the shape, and how to
know what values to plug in. The large end always seems to be VERY large
compared to the smaller side, though IIRC, with certain values one can get a
normal torus, so I suppose that might not be inherent in the object.
I found a very interesting paper by Langevin: "Geometry with two screens and
computational graphics" (2014) where he points out some very interesting and
useful properties of the true Dupin cyclide.
Thank you as always, Jerome - and please thank your friends for me as well.
I'm sure I will have some fun playing with this new algebraic surface and
finding ways to harness it in some future scenes.
Post a reply to this message
|
|
| |
| |
|
|
|
|
| |
| |
|
|
In my quest to parameterize this, I found:
http://or.nsfc.gov.cn/bitstream/00001903-5/173475/1/1000009340012.pdf
Which has some useful information; however, it's stated that:
a,b > 0 and c,f >= 0 are constants.
it goes on to state that it's a ring cyclide if f < c < a.
How can f be less than c if for a torus c=0 and f=r?
(unless of course f < 0)...
it's a normal torus if a=b=R, c=0, and f=r
Currently in the process of expanding the polynomial and grouping to see how a&b
behave when equal, eliminating the c terms, and seeing how f affects the minor
radii of the Dupin cyclide.
Post a reply to this message
|
|
| |
| |
|
|
|
|
| |
| |
|
|
Le 17/05/2017 à 19:05, Bald Eagle a écrit :
> In my quest to parameterize this, I found:
>
> http://or.nsfc.gov.cn/bitstream/00001903-5/173475/1/1000009340012.pdf
>
> Which has some useful information; however, it's stated that:
>
> a,b > 0 and c,f >= 0 are constants.
> it goes on to state that it's a ring cyclide if f < c < a.
>
> How can f be less than c if for a torus c=0 and f=r?
> (unless of course f < 0)...
when c=0, the sign of f( aka d) is irrelevant
(x^2+y^2+z^2+b^2-d^2)^2-4*(a*x-c*d)^2-4*b^2*y^2=0
>
> it's a normal torus if a=b=R, c=0, and f=r
>
>
> Currently in the process of expanding the polynomial and grouping to see how a&b
> behave when equal, eliminating the c terms, and seeing how f affects the minor
> radii of the Dupin cyclide.
>
you already answered that: d(aka f) is then the minor radius, and a=b is
the major radius.
Comparing to the traditional torus equation:
(x^2+y^2+z^2+R^2-r^2)^2-4R^2(x^2+y^2) = 0
the change of radius is due to ax-cd instead of Rx.
Notice that b <= a, and they should be positive, yet c is amongst the
two roots of c^2 = a^2-b^2.
the sign of b is irrelevant, but the sign of a can have an impact with
the sign of c.
Post a reply to this message
|
|
| |
| |
|
|
|
|
| |
| |
|
|
I did a bit of work on this last night, and thanks to
figures 12.21 and 12.22 here:
http://what-when-how.com/computer-graphics-and-geometric-modeling/surfaces-in-computer-graphics-geometric-modeling-part
-7/
I was able to work out some of the circles that lie on the face of the algebraic
surface, the inner ellipse, and the hyperboloid that the surface is tangent to.
It's not perfect - I think there's a little bit of difference due to the scaling
in the z direction due to B=A*0.98, but most of the hard stuff is worked out.
[The x-z plane tori were the hardest to work out. I had a few renders that
looked more like a Hopf Fibration than a Dupin Cyclide. I'll post when I get
some free time ;) ]
I think once I get the Yvon-Vilarceau circles worked out, and everything looking
nice, it will make a nice little macro package.
*
There's a lot of information out there on the cyclides - very interesting work
related to the Appolonian Gasket Soddy Circles, and Soddy's Hexlet - as well as
research into blending surfaces in computer modeling packages - which ought to
be of great interest to those working on creating fillets and smooth transitions
with mathematical primitives.
*
Post a reply to this message
Attachments:
Download 'dupincyclideparameterization.png' (354 KB)
Preview of image 'dupincyclideparameterization.png'
|
|
| |
| |
|
|
|
|
| |
| |
|
|
Le 18/05/2017 à 13:57, Bald Eagle a écrit :
> I did a bit of work on this last night, and thanks to
>
> figures 12.21 and 12.22 here:
>
http://what-when-how.com/computer-graphics-and-geometric-modeling/surfaces-in-computer-graphics-geometric-modeling-part-7/
>
> I was able to work out some of the circles that lie on the face of the algebraic
> surface, the inner ellipse, and the hyperboloid that the surface is tangent to.
>
> It's not perfect - I think there's a little bit of difference due to the scaling
> in the z direction due to B=A*0.98, but most of the hard stuff is worked out.
Thanks to your link, I would say:
A = R
C = abs(r-rho)/2
B = sqrt(A^2-C^2)
D = (r+rho)/2
with R the traditional radius (major) of the torus
r & rho the two different minor radius
and voila!
Here with 3,2 & 1 (and still the white sphere of radius 1 at origin)
Post a reply to this message
Attachments:
Download 'cyclide0.png' (107 KB)
Download 'cyclide2.png' (76 KB)
Preview of image 'cyclide0.png'
Preview of image 'cyclide2.png'
|
|
| |
| |
|
|
|
|
| |
| |
|
|
Le_Forgeron <jgr### [at] freefr> wrote:
> Thanks to your link, I would say:
>
> A = R
> C = abs(r-rho)/2
> B = sqrt(A^2-C^2)
> D = (r+rho)/2
>
> with R the traditional radius (major) of the torus
> r & rho the two different minor radius
>
> and voila!
>
> Here with 3,2 & 1 (and still the white sphere of radius 1 at origin)
Yes, excellent - that's the same answer that I came up with, which is also in
agreement with the formulas I came up with for the parametric version.
Hopefully it works without any difficulty. I'll see if I can apply it to the
planetary and satellite orbitals tonight or this weekend. :)
Thanks again!
Post a reply to this message
|
|
| |
| |
|
|
|
|
| |
| |
|
|
"Bald Eagle" <cre### [at] netscapenet> wrote:
> [The x-z plane tori were the hardest to work out. I had a few renders that
> looked more like a Hopf Fibration than a Dupin Cyclide. I'll post when I get
> some free time ;) ]
Post a reply to this message
Attachments:
Download 'dupincyclideparameterization - coolmistake.png' (349 KB)
Preview of image 'dupincyclideparameterization - coolmistake.png'
|
|
| |
| |
|
|
|
|
| |
| |
|
|
2nd:
Post a reply to this message
Attachments:
Download 'dupincyclideparameterization - coolmistake2.png' (336 KB)
Preview of image 'dupincyclideparameterization - coolmistake2.png'
|
|
| |
| |
|
|
|
|
| |
| |
|
|
3rd:
Post a reply to this message
Attachments:
Download 'dupincyclideparameterization - coolmistake3.png' (294 KB)
Preview of image 'dupincyclideparameterization - coolmistake3.png'
|
|
| |
| |
|
|
|
|
| |
|
|
