/***** * drawsurface.cc * * Stores a surface that has been added to a picture. *****/ #include "drawsurface.h" #include "arrayop.h" #include #include #include using namespace prc; namespace camp { const triple drawElement::zero; //double Tx[3]; // x-component of current transform //double Ty[3]; // y-component of current transform double Tz[3]; // z-component of current transform using vm::array; #ifdef HAVE_GL BezierCurve drawSurface::C; BezierPatch drawBezierPatch::S; BezierTriangle drawBezierTriangle::S; void storecolor(GLfloat *colors, int i, const vm::array &pens, int j) { pen p=vm::read(pens,j); p.torgb(); colors[i]=p.red(); colors[i+1]=p.green(); colors[i+2]=p.blue(); colors[i+3]=p.opacity(); } void storecolor(GLfloat *colors, int i, const RGBAColour& p) { colors[i]=p.R; colors[i+1]=p.G; colors[i+2]=p.B; colors[i+3]=p.A; } void setcolors(bool colors, bool lighton, const RGBAColour& diffuse, const RGBAColour& ambient, const RGBAColour& emissive, const RGBAColour& specular, double shininess) { static prc::RGBAColour lastdiffuse; static prc::RGBAColour lastambient; static prc::RGBAColour lastemissive; static prc::RGBAColour lastspecular; static double lastshininess=0.0; static bool lastcolors=false; if(colors != lastcolors) { drawBezierPatch::S.draw(); lastcolors=colors; } if(!colors && (diffuse != lastdiffuse || ambient != lastambient || emissive != lastemissive || specular != lastspecular || shininess != lastshininess)) { drawBezierPatch::S.draw(); lastdiffuse=diffuse; lastambient=ambient; lastemissive=emissive; lastspecular=specular; lastshininess=shininess; } if(colors) { if(!lighton) glColorMaterial(GL_FRONT_AND_BACK,GL_EMISSION); GLfloat Black[]={0,0,0,(GLfloat) diffuse.A}; glMaterialfv(GL_FRONT_AND_BACK,GL_DIFFUSE,Black); glMaterialfv(GL_FRONT_AND_BACK,GL_AMBIENT,Black); glMaterialfv(GL_FRONT_AND_BACK,GL_EMISSION,Black); } else { GLfloat Diffuse[]={(GLfloat) diffuse.R,(GLfloat) diffuse.G, (GLfloat) diffuse.B,(GLfloat) diffuse.A}; glMaterialfv(GL_FRONT_AND_BACK,GL_DIFFUSE,Diffuse); GLfloat Ambient[]={(GLfloat) ambient.R,(GLfloat) ambient.G, (GLfloat) ambient.B,(GLfloat) ambient.A}; glMaterialfv(GL_FRONT_AND_BACK,GL_AMBIENT,Ambient); GLfloat Emissive[]={(GLfloat) emissive.R,(GLfloat) emissive.G, (GLfloat) emissive.B,(GLfloat) emissive.A}; glMaterialfv(GL_FRONT_AND_BACK,GL_EMISSION,Emissive); } if(lighton) { GLfloat Specular[]={(GLfloat) specular.R,(GLfloat) specular.G, (GLfloat) specular.B,(GLfloat) specular.A}; glMaterialfv(GL_FRONT_AND_BACK,GL_SPECULAR,Specular); glMaterialf(GL_FRONT_AND_BACK,GL_SHININESS,128.0*shininess); } } #endif void drawBezierPatch::bounds(const double* t, bbox3& b) { double x,y,z; double X,Y,Z; if(straight) { triple Vertices[3]; if(t == NULL) { Vertices[0]=controls[0]; Vertices[1]=controls[3]; Vertices[2]=controls[12]; Vertices[3]=controls[15]; } else { Vertices[0]=t*controls[0]; Vertices[1]=t*controls[3]; Vertices[2]=t*controls[12]; Vertices[3]=t*controls[15]; } boundstriples(x,y,z,X,Y,Z,4,Vertices); } else { double cx[16]; double cy[16]; double cz[16]; if(t == NULL) { for(int i=0; i < 16; ++i) { triple v=controls[i]; cx[i]=v.getx(); cy[i]=v.gety(); cz[i]=v.getz(); } } else { for(int i=0; i < 16; ++i) { triple v=t*controls[i]; cx[i]=v.getx(); cy[i]=v.gety(); cz[i]=v.getz(); } } double c0=cx[0]; double fuzz=sqrtFuzz*run::norm(cx,16); x=bound(cx,min,b.empty ? c0 : min(c0,b.left),fuzz,maxdepth); X=bound(cx,max,b.empty ? c0 : max(c0,b.right),fuzz,maxdepth); c0=cy[0]; fuzz=sqrtFuzz*run::norm(cy,16); y=bound(cy,min,b.empty ? c0 : min(c0,b.bottom),fuzz,maxdepth); Y=boundtri(cy,max,b.empty ? c0 : max(c0,b.top),fuzz,maxdepth); c0=cz[0]; fuzz=sqrtFuzz*run::norm(cz,16); z=bound(cz,min,b.empty ? c0 : min(c0,b.lower),fuzz,maxdepth); Z=bound(cz,max,b.empty ? c0 : max(c0,b.upper),fuzz,maxdepth); } b.add(x,y,z); b.add(X,Y,Z); if(t == NULL) { Min=triple(x,y,z); Max=triple(X,Y,Z); } } void drawBezierPatch::ratio(const double* t, pair &b, double (*m)(double, double), double fuzz, bool &first) { triple buf[16]; triple* Controls; if(straight) { if(t == NULL) Controls=controls; else { Controls=buf; Controls[0]=t*controls[0]; Controls[3]=t*controls[3]; Controls[12]=t*controls[12]; Controls[15]=t*controls[15]; } triple v=Controls[0]; double x=xratio(v); double y=yratio(v); if(first) { first=false; b=pair(x,y); } else { x=m(b.getx(),x); y=m(b.gety(),y); } v=Controls[3]; x=m(x,xratio(v)); y=m(y,yratio(v)); v=Controls[12]; x=m(x,xratio(v)); y=m(y,yratio(v)); v=Controls[15]; x=m(x,xratio(v)); y=m(y,yratio(v)); b=pair(x,y); } else { if(t == NULL) Controls=controls; else { Controls=buf; for(unsigned int i=0; i < 16; ++i) Controls[i]=t*controls[i]; } if(first) { triple v=Controls[0]; b=pair(xratio(v),yratio(v)); first=false; } b=pair(bound(Controls,m,xratio,b.getx(),fuzz,maxdepth), bound(Controls,m,yratio,b.gety(),fuzz,maxdepth)); } } bool drawBezierPatch::write(prcfile *out, unsigned int *, double, groupsmap&) { if(invisible || !prc) return true; PRCmaterial m(ambient,diffuse,emissive,specular,opacity,PRCshininess); if(straight) { triple vertices[]={controls[0],controls[12],controls[3],controls[15]}; if(colors) out->addQuad(vertices,colors); else out->addRectangle(vertices,m); } else out->addPatch(controls,m); return true; } void drawBezierPatch::render(GLUnurbs *nurb, double size2, const triple& b, const triple& B, double perspective, bool lighton, bool transparent) { #ifdef HAVE_GL if(invisible || ((colors ? colors[0].A+colors[1].A+colors[2].A+colors[3].A < 4.0 : diffuse.A < 1.0) ^ transparent)) return; const bool billboard=interaction == BILLBOARD && !settings::getSetting("offscreen"); triple m,M; double f,F,s; if(perspective) { f=Min.getz()*perspective; F=Max.getz()*perspective; m=triple(min(f*b.getx(),F*b.getx()),min(f*b.gety(),F*b.gety()),b.getz()); M=triple(max(f*B.getx(),F*B.getx()),max(f*B.gety(),F*B.gety()),B.getz()); s=max(f,F); } else { m=b; M=B; s=1.0; } const pair size3(s*(B.getx()-b.getx()),s*(B.gety()-b.gety())); double t[16]; // current transform glGetDoublev(GL_MODELVIEW_MATRIX,t); // Like Fortran, OpenGL uses transposed (column-major) format! run::transpose(t,4); /* Tx[0]=t[0]; Tx[1]=t[1]; Tx[2]=t[2]; Ty[0]=t[4]; Ty[1]=t[5]; Ty[2]=t[6]; */ Tz[0]=t[8]; Tz[1]=t[9]; Tz[2]=t[10]; run::inverse(t,4); bbox3 box(m,M); box.transform(t); m=box.Min(); M=box.Max(); if(!billboard && (Max.getx() < m.getx() || Min.getx() > M.getx() || Max.gety() < m.gety() || Min.gety() > M.gety() || Max.getz() < m.getz() || Min.getz() > M.getz())) return; setcolors(colors,lighton,diffuse,ambient,emissive,specular,shininess); if(billboard) BB.init(center); GLfloat v[16]; if(colors) for(size_t i=0; i < 4; ++i) storecolor(v,4*i,colors[i]); triple *Controls; triple Controls0[16]; if(billboard) { Controls=Controls0; for(size_t i=0; i < 16; i++) Controls[i]=BB.transform(controls[i]); } else Controls=controls; if(gl::outlinemode) { triple edge0[]={Controls[0],Controls[4],Controls[8],Controls[12]}; C.queue(edge0,straight,size3.length()/size2,m,M); triple edge1[]={Controls[12],Controls[13],Controls[14],Controls[15]}; C.queue(edge1,straight,size3.length()/size2,m,M); triple edge2[]={Controls[15],Controls[11],Controls[7],Controls[3]}; C.queue(edge2,straight,size3.length()/size2,m,M); triple edge3[]={Controls[3],Controls[2],Controls[1],Controls[0]}; C.queue(edge3,straight,size3.length()/size2,m,M); C.draw(); } else { S.queue(Controls,straight,size3.length()/size2,m,M,transparent, colors ? v : NULL); if(!transparent) S.draw(); } #endif } drawElement *drawBezierPatch::transformed(const double* t) { return new drawBezierPatch(t,this); } void drawBezierTriangle::bounds(const double* t, bbox3& b) { double x,y,z; double X,Y,Z; if(straight) { triple Vertices[3]; if(t == NULL) { Vertices[0]=controls[0]; Vertices[1]=controls[6]; Vertices[2]=controls[9]; } else { Vertices[0]=t*controls[0]; Vertices[1]=t*controls[6]; Vertices[2]=t*controls[9]; } boundstriples(x,y,z,X,Y,Z,3,Vertices); } else { double cx[10]; double cy[10]; double cz[10]; if(t == NULL) { for(unsigned int i=0; i < 10; ++i) { triple v=controls[i]; cx[i]=v.getx(); cy[i]=v.gety(); cz[i]=v.getz(); } } else { for(unsigned int i=0; i < 10; ++i) { triple v=t*controls[i]; cx[i]=v.getx(); cy[i]=v.gety(); cz[i]=v.getz(); } } double c0=cx[0]; double fuzz=sqrtFuzz*run::norm(cx,10); x=boundtri(cx,min,b.empty ? c0 : min(c0,b.left),fuzz,maxdepth); X=boundtri(cx,max,b.empty ? c0 : max(c0,b.right),fuzz,maxdepth); c0=cy[0]; fuzz=sqrtFuzz*run::norm(cy,10); y=boundtri(cy,min,b.empty ? c0 : min(c0,b.bottom),fuzz,maxdepth); Y=boundtri(cy,max,b.empty ? c0 : max(c0,b.top),fuzz,maxdepth); c0=cz[0]; fuzz=sqrtFuzz*run::norm(cz,10); z=boundtri(cz,min,b.empty ? c0 : min(c0,b.lower),fuzz,maxdepth); Z=boundtri(cz,max,b.empty ? c0 : max(c0,b.upper),fuzz,maxdepth); } b.add(x,y,z); b.add(X,Y,Z); if(t == NULL) { Min=triple(x,y,z); Max=triple(X,Y,Z); } } void drawBezierTriangle::ratio(const double* t, pair &b, double (*m)(double, double), double fuzz, bool &first) { triple buf[10]; triple* Controls; if(straight) { if(t == NULL) Controls=controls; else { Controls=buf; Controls[0]=t*controls[0]; Controls[6]=t*controls[6]; Controls[9]=t*controls[9]; } triple v=Controls[0]; double x=xratio(v); double y=yratio(v); if(first) { first=false; b=pair(x,y); } else { x=m(b.getx(),x); y=m(b.gety(),y); } v=Controls[6]; x=m(x,xratio(v)); y=m(y,yratio(v)); v=Controls[9]; x=m(x,xratio(v)); y=m(y,yratio(v)); b=pair(x,y); } else { if(t == NULL) Controls=controls; else { Controls=buf; for(unsigned int i=0; i < 10; ++i) Controls[i]=t*controls[i]; } if(first) { triple v=Controls[0]; b=pair(xratio(v),yratio(v)); first=false; } b=pair(boundtri(Controls,m,xratio,b.getx(),fuzz,maxdepth), boundtri(Controls,m,yratio,b.gety(),fuzz,maxdepth)); } } bool drawBezierTriangle::write(prcfile *out, unsigned int *, double, groupsmap&) { if(invisible) return true; PRCmaterial m(ambient,diffuse,emissive,specular,opacity,PRCshininess); static const double third=1.0/3.0; static const double third2=2.0/3.0; triple Controls[]={controls[0],controls[0],controls[0],controls[0], controls[1],third2*controls[1]+third*controls[2], third*controls[1]+third2*controls[2], controls[2],controls[3], third*controls[3]+third2*controls[4], third2*controls[4]+third*controls[5], controls[5],controls[6],controls[7], controls[8],controls[9]}; out->addPatch(Controls,m); return true; } void drawBezierTriangle::render(GLUnurbs *nurb, double size2, const triple& b, const triple& B, double perspective, bool lighton, bool transparent) { #ifdef HAVE_GL if(invisible || ((colors ? colors[0].A+colors[1].A+colors[2].A < 3.0 : diffuse.A < 1.0) ^ transparent)) return; const bool billboard=interaction == BILLBOARD && !settings::getSetting("offscreen"); triple m,M; double f,F,s; if(perspective) { f=Min.getz()*perspective; F=Max.getz()*perspective; m=triple(min(f*b.getx(),F*b.getx()),min(f*b.gety(),F*b.gety()),b.getz()); M=triple(max(f*B.getx(),F*B.getx()),max(f*B.gety(),F*B.gety()),B.getz()); s=max(f,F); } else { m=b; M=B; s=1.0; } const pair size3(s*(B.getx()-b.getx()),s*(B.gety()-b.gety())); double t[16]; // current transform glGetDoublev(GL_MODELVIEW_MATRIX,t); // Like Fortran, OpenGL uses transposed (column-major) format! run::transpose(t,4); /* Tx[0]=t[0]; Tx[1]=t[1]; Tx[2]=t[2]; Ty[0]=t[4]; Ty[1]=t[5]; Ty[2]=t[6]; */ Tz[0]=t[8]; Tz[1]=t[9]; Tz[2]=t[10]; run::inverse(t,4); bbox3 box(m,M); box.transform(t); m=box.Min(); M=box.Max(); if(!billboard && (Max.getx() < m.getx() || Min.getx() > M.getx() || Max.gety() < m.gety() || Min.gety() > M.gety() || Max.getz() < m.getz() || Min.getz() > M.getz())) return; setcolors(colors,lighton,diffuse,ambient,emissive,specular,shininess); if(billboard) BB.init(center); GLfloat v[12]; if(colors) for(size_t i=0; i < 3; ++i) storecolor(v,4*i,colors[i]); triple *Controls; triple Controls0[10]; if(billboard) { Controls=Controls0; for(size_t i=0; i < 10; i++) Controls[i]=BB.transform(controls[i]); } else Controls=controls; if(gl::outlinemode) { triple edge0[]={Controls[0],Controls[1],Controls[3],Controls[6]}; C.queue(edge0,straight,size3.length()/size2,m,M); triple edge1[]={Controls[6],Controls[7],Controls[8],Controls[9]}; C.queue(edge1,straight,size3.length()/size2,m,M); triple edge2[]={Controls[9],Controls[5],Controls[2],Controls[0]}; C.queue(edge2,straight,size3.length()/size2,m,M); C.draw(); } else { S.queue(Controls,straight,size3.length()/size2,m,M,transparent, colors ? v : NULL); if(!transparent) S.draw(); } #endif } drawElement *drawBezierTriangle::transformed(const double* t) { return new drawBezierTriangle(t,this); } bool drawNurbs::write(prcfile *out, unsigned int *, double, groupsmap&) { if(invisible) return true; PRCmaterial m(ambient,diffuse,emissive,specular,opacity,PRCshininess); out->addSurface(udegree,vdegree,nu,nv,controls,uknots,vknots,m,weights); return true; } // Approximate bounds by bounding box of control polyhedron. void drawNurbs::bounds(const double* t, bbox3& b) { double x,y,z; double X,Y,Z; const size_t n=nu*nv; triple* Controls; if(t == NULL) Controls=controls; else { Controls=new triple[n]; for(size_t i=0; i < n; i++) Controls[i]=t*controls[i]; } boundstriples(x,y,z,X,Y,Z,n,Controls); b.add(x,y,z); b.add(X,Y,Z); if(t == NULL) { Min=triple(x,y,z); Max=triple(X,Y,Z); } else delete[] Controls; } drawElement *drawNurbs::transformed(const double* t) { return new drawNurbs(t,this); } void drawNurbs::ratio(const double *t, pair &b, double (*m)(double, double), double, bool &first) { const size_t n=nu*nv; triple* Controls; if(t == NULL) Controls=controls; else { Controls=new triple[n]; for(unsigned int i=0; i < n; ++i) Controls[i]=t*controls[i]; } if(first) { first=false; triple v=Controls[0]; b=pair(xratio(v),yratio(v)); } double x=b.getx(); double y=b.gety(); for(size_t i=0; i < n; ++i) { triple v=Controls[i]; x=m(x,xratio(v)); y=m(y,yratio(v)); } b=pair(x,y); if(t != NULL) delete[] Controls; } void drawNurbs::displacement() { #ifdef HAVE_GL size_t n=nu*nv; size_t nuknots=udegree+nu+1; size_t nvknots=vdegree+nv+1; if(Controls == NULL) { Controls=new(UseGC) GLfloat[(weights ? 4 : 3)*n]; uKnots=new(UseGC) GLfloat[nuknots]; vKnots=new(UseGC) GLfloat[nvknots]; } if(weights) for(size_t i=0; i < n; ++i) store(Controls+4*i,controls[i],weights[i]); else for(size_t i=0; i < n; ++i) store(Controls+3*i,controls[i]); for(size_t i=0; i < nuknots; ++i) uKnots[i]=uknots[i]; for(size_t i=0; i < nvknots; ++i) vKnots[i]=vknots[i]; #endif } void drawNurbs::render(GLUnurbs *nurb, double size2, const triple& Min, const triple& Max, double perspective, bool lighton, bool transparent) { #ifdef HAVE_GL if(invisible || ((colors ? colors[3]+colors[7]+colors[11]+colors[15] < 4.0 : diffuse.A < 1.0) ^ transparent)) return; double t[16]; // current transform glGetDoublev(GL_MODELVIEW_MATRIX,t); run::transpose(t,4); bbox3 B(this->Min,this->Max); B.transform(t); triple m=B.Min(); triple M=B.Max(); if(perspective) { double f=m.getz()*perspective; double F=M.getz()*perspective; if(M.getx() < min(f*Min.getx(),F*Min.getx()) || m.getx() > max(f*Max.getx(),F*Max.getx()) || M.gety() < min(f*Min.gety(),F*Min.gety()) || m.gety() > max(f*Max.gety(),F*Max.gety()) || M.getz() < Min.getz() || m.getz() > Max.getz()) return; } else { if(M.getx() < Min.getx() || m.getx() > Max.getx() || M.gety() < Min.gety() || m.gety() > Max.gety() || M.getz() < Min.getz() || m.getz() > Max.getz()) return; } setcolors(colors,lighton,diffuse,ambient,emissive,specular,shininess); gluBeginSurface(nurb); int uorder=udegree+1; int vorder=vdegree+1; size_t stride=weights ? 4 : 3; gluNurbsSurface(nurb,uorder+nu,uKnots,vorder+nv,vKnots,stride*nv,stride, Controls,uorder,vorder, weights ? GL_MAP2_VERTEX_4 : GL_MAP2_VERTEX_3); if(lighton) gluNurbsCallback(nurb,GLU_NURBS_NORMAL,(_GLUfuncptr) glNormal3fv); if(colors) { static GLfloat linear[]={0.0,0.0,1.0,1.0}; gluNurbsSurface(nurb,4,linear,4,linear,8,4,colors,2,2,GL_MAP2_COLOR_4); } gluEndSurface(nurb); if(colors) glDisable(GL_COLOR_MATERIAL); #endif } void drawSphere::P(triple& t, double x, double y, double z) { if(half) { double temp=z; z=x; x=-temp; } if(T == NULL) { t=triple(x,y,z); return; } double f=T[12]*x+T[13]*y+T[14]*z+T[15]; if(f == 0.0) run::dividebyzero(); f=1.0/f; t=triple((T[0]*x+T[1]*y+T[2]*z+T[3])*f,(T[4]*x+T[5]*y+T[6]*z+T[7])*f, (T[8]*x+T[9]*y+T[10]*z+T[11])*f); } bool drawSphere::write(prcfile *out, unsigned int *, double, groupsmap&) { if(invisible) return true; PRCmaterial m(ambient,diffuse,emissive,specular,opacity,shininess); switch(type) { case 0: // PRCsphere { if(half) out->addHemisphere(1.0,m,NULL,NULL,NULL,1.0,T); else out->addSphere(1.0,m,NULL,NULL,NULL,1.0,T); break; } case 1: // NURBSsphere { static double uknot[]={0.0,0.0,1.0/3.0,0.5,1.0,1.0}; static double vknot[]={0.0,0.0,0.0,0.0,1.0,1.0,1.0,1.0}; static double Weights[12]={2.0/3.0,2.0/9.0,2.0/9.0,2.0/3.0, 1.0/3.0,1.0/9.0,1.0/9.0,1.0/3.0, 1.0,1.0/3.0,1.0/3.0,1.0}; // NURBS representation of a sphere using 10 distinct control points // K. Qin, J. Comp. Sci. and Tech. 12, 210-216 (1997). triple N,S,P1,P2,P3,P4,P5,P6,P7,P8; P(N,0.0,0.0,1.0); P(P1,-2.0,-2.0,1.0); P(P2,-2.0,-2.0,-1.0); P(S,0.0,0.0,-1.0); P(P3,2.0,-2.0,1.0); P(P4,2.0,-2.0,-1.0); P(P5,2.0,2.0,1.0); P(P6,2.0,2.0,-1.0); P(P7,-2.0,2.0,1.0); P(P8,-2.0,2.0,-1.0); triple p0[]={N,P1,P2,S, N,P3,P4,S, N,P5,P6,S, N,P7,P8,S, N,P1,P2,S, N,P3,P4,S}; out->addSurface(2,3,3,4,p0,uknot,vknot,m,Weights); out->addSurface(2,3,3,4,p0+4,uknot,vknot,m,Weights); if(!half) { out->addSurface(2,3,3,4,p0+8,uknot,vknot,m,Weights); out->addSurface(2,3,3,4,p0+12,uknot,vknot,m,Weights); } break; } default: reportError("Invalid sphere type"); } return true; } bool drawCylinder::write(prcfile *out, unsigned int *, double, groupsmap&) { if(invisible) return true; PRCmaterial m(ambient,diffuse,emissive,specular,opacity,shininess); out->addCylinder(1.0,1.0,m,NULL,NULL,NULL,1.0,T); return true; } bool drawDisk::write(prcfile *out, unsigned int *, double, groupsmap&) { if(invisible) return true; PRCmaterial m(ambient,diffuse,emissive,specular,opacity,shininess); out->addDisk(1.0,m,NULL,NULL,NULL,1.0,T); return true; } bool drawTube::write(prcfile *out, unsigned int *, double, groupsmap&) { if(invisible) return true; PRCmaterial m(ambient,diffuse,emissive,specular,opacity,shininess); Int n=center.length(); if(center.piecewisestraight()) { triple *centerControls=new(UseGC) triple[n+1]; for(Int i=0; i <= n; ++i) centerControls[i]=center.point(i); size_t N=n+1; triple *controls=new(UseGC) triple[N]; for(Int i=0; i <= n; ++i) controls[i]=g.point(i); out->addTube(N,centerControls,controls,true,m); } else { size_t N=3*n+1; triple *centerControls=new(UseGC) triple[N]; centerControls[0]=center.point((Int) 0); centerControls[1]=center.postcontrol((Int) 0); size_t k=1; for(Int i=1; i < n; ++i) { centerControls[++k]=center.precontrol(i); centerControls[++k]=center.point(i); centerControls[++k]=center.postcontrol(i); } centerControls[++k]=center.precontrol(n); centerControls[++k]=center.point(n); triple *controls=new(UseGC) triple[N]; controls[0]=g.point((Int) 0); controls[1]=g.postcontrol((Int) 0); k=1; for(Int i=1; i < n; ++i) { controls[++k]=g.precontrol(i); controls[++k]=g.point(i); controls[++k]=g.postcontrol(i); } controls[++k]=g.precontrol(n); controls[++k]=g.point(n); out->addTube(N,centerControls,controls,false,m); } return true; } bool drawPixel::write(prcfile *out, unsigned int *, double, groupsmap&) { if(invisible) return true; out->addPoint(v,c,width); return true; } void drawPixel::render(GLUnurbs *nurb, double size2, const triple& Min, const triple& Max, double perspective, bool lighton, bool transparent) { #ifdef HAVE_GL if(invisible) return; GLfloat V[4]; glEnable(GL_COLOR_MATERIAL); glColorMaterial(GL_FRONT_AND_BACK,GL_EMISSION); static GLfloat Black[]={0,0,0,1}; glMaterialfv(GL_FRONT_AND_BACK,GL_DIFFUSE,Black); glMaterialfv(GL_FRONT_AND_BACK,GL_AMBIENT,Black); glMaterialfv(GL_FRONT_AND_BACK,GL_SPECULAR,Black); glMaterialf(GL_FRONT_AND_BACK,GL_SHININESS,0.0); glPointSize(1.0+width); glBegin(GL_POINT); storecolor(V,0,c); glColor4fv(V); store(V,v); glVertex3fv(V); glEnd(); glPointSize(1.0); glDisable(GL_COLOR_MATERIAL); #endif } const string drawBaseTriangles::wrongsize= "triangle indices require 3 components"; const string drawBaseTriangles::outofrange="index out of range"; void drawBaseTriangles::bounds(const double* t, bbox3& b) { double x,y,z; double X,Y,Z; triple* tP; if(t == NULL) tP=P; else { tP=new triple[nP]; for(size_t i=0; i < nP; i++) tP[i]=t*P[i]; } boundstriples(x,y,z,X,Y,Z,nP,tP); b.add(x,y,z); b.add(X,Y,Z); if(t == NULL) { Min=triple(x,y,z); Max=triple(X,Y,Z); } else delete[] tP; } void drawBaseTriangles::ratio(const double* t, pair &b, double (*m)(double, double), double fuzz, bool &first) { triple* tP; if(t == NULL) tP=P; else { tP=new triple[nP]; for(size_t i=0; i < nP; i++) tP[i]=t*P[i]; } ratiotriples(b,m,first,nP,tP); if(t != NULL) delete[] tP; } bool drawTriangles::write(prcfile *out, unsigned int *, double, groupsmap&) { if(invisible) return true; if (nC) { const RGBAColour white(1,1,1,opacity); const RGBAColour black(0,0,0,opacity); const PRCmaterial m(black,white,black,specular,opacity,PRCshininess); out->addTriangles(nP,P,nI,PI,m,nN,N,NI,0,NULL,NULL,nC,C,CI,0,NULL,NULL,30); } else { const PRCmaterial m(ambient,diffuse,emissive,specular,opacity,PRCshininess); out->addTriangles(nP,P,nI,PI,m,nN,N,NI,0,NULL,NULL,0,NULL,NULL,0,NULL,NULL,30); } return true; } void drawTriangles::render(GLUnurbs *nurb, double size2, const triple& Min, const triple& Max, double perspective, bool lighton, bool transparent) { #ifdef HAVE_GL if(invisible) return; if(invisible || ((diffuse.A < 1.0) ^ transparent)) return; triple m,M; double t[16]; // current transform glGetDoublev(GL_MODELVIEW_MATRIX,t); run::transpose(t,4); bbox3 B(this->Min,this->Max); B.transform(t); m=B.Min(); M=B.Max(); if(perspective) { const double f=m.getz()*perspective; const double F=M.getz()*perspective; if((M.getx() < min(f*Min.getx(),F*Min.getx()) || m.getx() > max(f*Max.getx(),F*Max.getx()) || M.gety() < min(f*Min.gety(),F*Min.gety()) || m.gety() > max(f*Max.gety(),F*Max.gety()) || M.getz() < Min.getz() || m.getz() > Max.getz())) return; } else { if((M.getx() < Min.getx() || m.getx() > Max.getx() || M.gety() < Min.gety() || m.gety() > Max.gety() || M.getz() < Min.getz() || m.getz() > Max.getz())) return; } setcolors(nC,!nC,diffuse,ambient,emissive,specular,shininess); if(!nN) lighton=false; glBegin(GL_TRIANGLES); for(size_t i=0; i < nI; i++) { const uint32_t *pi=PI[i]; const uint32_t *ni=NI[i]; const uint32_t *ci=nC ? CI[i] : 0; if(lighton) glNormal3f(N[ni[0]].getx(),N[ni[0]].gety(),N[ni[0]].getz()); if(nC) glColor4f(C[ci[0]].R,C[ci[0]].G,C[ci[0]].B,C[ci[0]].A); glVertex3f(P[pi[0]].getx(),P[pi[0]].gety(),P[pi[0]].getz()); if(lighton) glNormal3f(N[ni[1]].getx(),N[ni[1]].gety(),N[ni[1]].getz()); if(nC) glColor4f(C[ci[1]].R,C[ci[1]].G,C[ci[1]].B,C[ci[1]].A); glVertex3f(P[pi[1]].getx(),P[pi[1]].gety(),P[pi[1]].getz()); if(lighton) glNormal3f(N[ni[2]].getx(),N[ni[2]].gety(),N[ni[2]].getz()); if(nC) glColor4f(C[ci[2]].R,C[ci[2]].G,C[ci[2]].B,C[ci[2]].A); glVertex3f(P[pi[2]].getx(),P[pi[2]].gety(),P[pi[2]].getz()); } glEnd(); if(nC) glDisable(GL_COLOR_MATERIAL); #endif } } //namespace camp