Subversion Repositories gelsvn

/* ----------------------------------------------------------------------- *

 * ----------------------------------------------------------------------- */

    const string PeriodicScalarFieldRenderer::vss =

    "	attribute float scalar;\n"

    "	varying vec3 _normal;\n"

    "	varying float s;\n"

    "	\n"

    "	void main(void)\n"

    "	{\n"

    "		gl_Position =  ftransform();\n"

    "		_normal = normalize(gl_NormalMatrix * gl_Normal);\n"

    "		s=scalar;\n"

    "	}\n";

    const string PeriodicScalarFieldRenderer::fss = 	

    "	varying vec3 _normal;\n"

    "	varying float s;\n"

    "   uniform float gamma;\n"

    "	const vec3 light_dir = vec3(0,0,1);\n"

    "	\n"

    "	void main()\n"

    "	{\n"

    "       vec3 normal = normalize(_normal);\n"

    "		float dot_ln = max(0.0,dot(light_dir, normal));\n"

    "		\n"

    "		gl_FragColor = sin(s) * vec4(-1,0,1,0);\n"

    "       gl_FragColor *= dot_ln;\n"

    "       gl_FragColor.r = pow(gl_FragColor.r, 1.0/gamma);\n"

    "       gl_FragColor.g = pow(gl_FragColor.g, 1.0/gamma);\n"

    "       gl_FragColor.b = pow(gl_FragColor.b, 1.0/gamma);\n"

    "	}\n";

    PeriodicScalarFieldRenderer::PeriodicScalarFieldRenderer(const Manifold& m, 

                                             bool smooth, 

                                             VertexAttributeVector<double>& field, 

                                             float gamma): SimpleShaderRenderer(vss, fss)

    {

        GLint old_prog;

        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);

        glUseProgram(prog);

        GLuint scalar_attrib = glGetAttribLocation(prog, "scalar");

        //    static float& gamma = CreateCVar("display.scalar_field_renderer.gamma",2.2f);

        glUniform1fARB(glGetUniformLocationARB(prog, "gamma"), gamma);

        glNewList(display_list,GL_COMPILE);

        for(FaceIDIterator f = m.faces_begin(); f != m.faces_end(); ++f){      

            if(!smooth) 

                glNormal3fv(normal(m, *f).get());

            if(no_edges(m, *f)== 3) 

                glBegin(GL_TRIANGLES);

            else 

                glBegin(GL_POLYGON);

            for(HalfEdgeWalker w = m.halfedgewalker(*f); !w.full_circle(); w = w.circulate_face_ccw()){

                Vec3d n(normal(m, w.vertex()));

                if(smooth) 

                    glNormal3dv(n.get());

                glVertexAttrib1d(scalar_attrib, field[w.vertex()]);

                glVertex3fv(m.pos(w.vertex()).get());

            }

            glEnd();

        }

        glEndList();	

        glUseProgram(old_prog);

    }