void circle( GLVApp& app ){
static Cir c = CXY(1);
DRAWANDTOUCH3(c,1,0,0);
static Cir src = CXZ(1).trs(1,0,0);
DRAWANDTOUCH3(src,0,1,0);
Cir nc = src.re(c);
DRAW3(nc,0,0,1);
}
virtual void onDraw(){
getMouse();
static Cir cir = CXY(1);
Touch(interface,cir);
Draw(cir);
Dlp dlp(1,0,0,adj);
//Touch(interface,dlp);
//Euclidean vector of circle at theta
Vec v = Ro::vec(cir,phi);
Draw( v );
for(int i=0;i<100;++i){
Pnt p = Ro::pnt_cir(cir,(float)i/100 * TWOPI);
auto np = Conic::Transform(p, v, amt );
Draw(np,0,1,0);
}
auto tnp = Conic::Transform( Ro::pnt_cir(cir,adj*TWOPI), v, amt ); //tangent
Draw( Ro::dls(tnp,.2),1,0,0);
// auto itnp = Conic::I
Dlp idlp = Conic::ITransform(dlp, v, amt);
Par par = (idlp ^ cir.dual() ).dual();
// Pnt ipnt = Pnt(idlp) + Ori(1);
auto tpnt = Ta::at(cir, Pnt(idlp));
Draw( tpnt, 1,1,0);
cout << Ro::size(par,false) << endl;
Draw(dlp,0,0,1);
Draw(idlp);
/* Dlp idlp2 = idlp.rot(Biv::xy * adj); */
/* Draw( idlp2,1,0,0 ); */
/* Draw( Conic::Transform(idlp2,v,amt) ); */
for(auto i : Ro::split(par)){
Draw( Conic::Transform(i,v,amt),1,0,0 );
}
}
virtual void onDraw(Graphics& gl){
//Model Transform
Rot t = Gen::aa( scene().model.rot() );
GL::rotate( t.w() );
//A Circle in the XY plane with radius 1.5
static Cir c = CXY(1.5);
DRAW(c); //<-- DRAW is a macro that draws an element
interface.touch(c);
//A Plane (with normal along y axis, at distance -.5 units from the origin)
static Dlp d(0,1,0,-.5);
DRAW3(d,1,0,0); //<-- DRAW3 is a macro that draws an element a certain color (red here)
interface.touch(d);
//Intersection of Circle and Plane (two points)
DRAW3( (d ^ c.dual() ).dual(), 0,1,0 );
}
void hopf(al::VsrApp& app){
HopfFiber hf;
static int time = 0; time++; double tphase = PI * time/180.0;
//GUI Controlled Parameters
static int num = 2;
static int res = 3;
static bool handed, bDrawCir, bDrawPnt,bReset;
static double phase,amt,minDistance;
//SET UP Gui
static bool bSet = 0;
if (!bSet) {
bSet = 1;
app.glv.gui(num,"num",0,1000)(res,"res",0,1000)(phase,"phase",-100,100)(amt,"amt",0,10)(minDistance,"minDistance",0,1000);
app.glv.gui(handed)(bDrawCir, "drawcir")(bDrawPnt,"drawpnt")(bReset,"reset");
}
vector<Pnt> vp;
for (int i = 0; i < num; ++i){
double u = 1.0 * i/num;
for (int j = 0; j < res; ++j){
double v = 1.0 * j/res;
Cir tc = hf.fiber(-1.0 + 2*u, -.5 + v);
Pnt tp = Ro::pnt_cir( tc, tphase * Ro::cur(tc) * phase);
if(bDrawPnt)DRAW3(tp,u,v,1-u);
if(bDrawCir)DRAW3(tc,v,u,1-u);
vp.push_back( Ro::pnt_cir( tc, v * u * PI ) );
}
}
static Dll dll = DLN(0,1,0); app.interface.touch(dll); DRAW(dll);
static Point point = PT(2,0,0); app.interface.touch(point); DRAW3(point,1,0,0);
for (int i = 0; i <= 1; ++i){
double t = i/1.0;
Dll tdll = dll.sp( Gen::rot(Biv::xy * t*PIOVERTWO) );
Pnt pa = Ro::split( ( tdll ^ Ro::sur(hf.cir) ).dual(), true);
Pnt pb = Ro::split( ( tdll ^ Ro::sur(hf.cir) ).dual(), false);
DRAW3(pa,1,1,0); DRAW3(pb,1,1,0);
Coord::Sph spha = Coord::vec2sph( Vector(pa).unit() );
Coord::Sph sphb = Coord::vec2sph( Vector(pb).unit() );
Cir cca = hf.fiber( spha.theta / PI, spha.phi / PI );
Cir ccb = hf.fiber( sphb.theta / PI, sphb.phi / PI );
Bst bst = Gen::bst( (cca.dual() * PI/3.0 + ccb.dual() * PI/2.0 )* amt );
DRAW3(cca,0,1,0); DRAW3(ccb,0,0,1);
vector<Pnt> pp;
Pnt tp = point;//.trs(-t,0,0);
for (int i = 0; i<1000; ++i){
tp = Ro::loc( tp.sp( bst ) );
pp.push_back(tp);
}
glColor3f(1,0,0);
glBegin(GL_LINE_STRIP);
for (int i = 0; i<pp.size(); ++i ){
GL::vertex( pp[i].w() );
}
glEnd();
}
int tnum = 10;
// for (int j = 0; j < tnum; ++j){
//
// double t = 1.0 * j/tnum;
//
// Bst bst2 = hf.mtt( );
//
// for (int i = 0; i<pp.size(); ++i ){
// for(int k = 0; k < 10; ++k){
//
// }
// }
//
// }
}