Subversion Repositories gelsvn

Rev

Rev 635 | Rev 647 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 635 Rev 643
Line 4... Line 4...
4
 * For license and list of authors, see ../../doc/intro.pdf
4
 * For license and list of authors, see ../../doc/intro.pdf
5
 * ----------------------------------------------------------------------- */
5
 * ----------------------------------------------------------------------- */
6
 
6
 
7
#include "ManifoldRenderer.h"
7
#include "ManifoldRenderer.h"
8
 
8
 
-
 
9
#include <GLUT/GLUT.h>
-
 
10
 
9
#include <algorithm>
11
#include <algorithm>
10
#include <string>
12
#include <string>
11
#include <cstdlib>
13
#include <cstdlib>
12
#include "../CGLA/Mat3x3d.h"
14
#include "../CGLA/Mat3x3d.h"
13
#include "../GLGraphics/glsl_shader.h"
15
#include "../GLGraphics/glsl_shader.h"
Line 112... Line 114...
112
        
114
        
113
    }
115
    }
114
    
116
    
115
    void SimpleShaderRenderer::compile_display_list(const Manifold& m, bool smooth)
117
    void SimpleShaderRenderer::compile_display_list(const Manifold& m, bool smooth)
116
    {
118
    {
-
 
119
        GLint old_prog;
-
 
120
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
-
 
121
        glUseProgram(prog);
117
        glNewList(display_list,GL_COMPILE);
122
        glNewList(display_list,GL_COMPILE);
118
        GLGraphics::draw(m, smooth);
123
        GLGraphics::draw(m, smooth);
119
        glEndList();	
124
        glEndList();	
-
 
125
        glUseProgram(old_prog);
120
    }
126
    }
121
    
127
    
122
    void SimpleShaderRenderer::draw()
128
    void SimpleShaderRenderer::draw()
123
    {
129
    {
124
        GLint old_prog;
130
        GLint old_prog;
Line 126... Line 132...
126
        glUseProgram(prog);
132
        glUseProgram(prog);
127
        glCallList(display_list);
133
        glCallList(display_list);
128
        glUseProgram(old_prog);
134
        glUseProgram(old_prog);
129
    }
135
    }
130
    
136
    
131
    const string NormalRenderer::vss = 
137
    const string NormalRenderer::vss =
132
    "varying vec3 _n;\n"
138
    "varying vec3 _n;\n"
133
    "varying vec3 v;\n"
139
    "varying vec3 v;\n"
134
    "\n"
140
    "\n"
135
    "void main(void)\n"
141
    "void main(void)\n"
136
    "{\n"
142
    "{\n"
137
    "	gl_Position = ftransform();\n"
143
    "	gl_Position = ftransform();\n"
138
    "	v = vec3(gl_ModelViewMatrix * gl_Vertex);\n"
144
    "	v = vec3(gl_ModelViewMatrix * gl_Vertex);\n"
139
    "	_n = normalize(gl_NormalMatrix * gl_Normal);\n"
145
    "	_n = normalize(gl_NormalMatrix * gl_Normal);\n"
140
    "}\n";
146
    "}\n";
141
    
147
    
142
    const string NormalRenderer::fss = 
148
    const string NormalRenderer::fss =
143
    "varying vec3 _n;\n"
149
    "varying vec3 _n;\n"
144
    "varying vec3 v;\n"
150
    "varying vec3 v;\n"
145
    "\n"
151
    "\n"
146
    "void main(void)\n"
152
    "void main(void)\n"
147
    "{\n"
153
    "{\n"
Line 156... Line 162...
156
    "	vec4 s = vec4(0.3)*smoothstep(0.98,0.9999,dot(r, e));\n"
162
    "	vec4 s = vec4(0.3)*smoothstep(0.98,0.9999,dot(r, e));\n"
157
    "	\n"
163
    "	\n"
158
    "	gl_FragColor =  d+s;\n"
164
    "	gl_FragColor =  d+s;\n"
159
    "}\n";
165
    "}\n";
160
    
166
    
-
 
167
    const string DebugRenderer::vss =
-
 
168
    "varying vec3 _n;\n"
-
 
169
    "varying vec3 v;\n"
-
 
170
    "varying vec3 c;\n"
-
 
171
    "\n"
-
 
172
    "void main(void)\n"
-
 
173
    "{\n"
-
 
174
    "	gl_Position = ftransform();\n"
-
 
175
    "	v = vec3(gl_ModelViewMatrix * gl_Vertex);\n"
-
 
176
    "	_n = normalize(gl_NormalMatrix * gl_Normal);\n"
-
 
177
    "   c = gl_Color.rgb;\n"
-
 
178
    "}\n";
-
 
179
    
-
 
180
    const string DebugRenderer::fss =
-
 
181
    "varying vec3 _n;\n"
-
 
182
    "varying vec3 v;\n"
-
 
183
    "varying vec3 c;\n"
-
 
184
    "\n"
-
 
185
    "void main(void)\n"
-
 
186
    "{\n"
-
 
187
    "   vec3 n = normalize(_n);\n"
-
 
188
    "	vec3 l = normalize(-v);\n"
-
 
189
    "	vec3 e = l;\n"
-
 
190
    "	vec3 r = normalize(2.0*dot(l, n)*n - l);\n"
-
 
191
    "	\n"
-
 
192
    "	vec4 a = vec4(0.0,0.1,.3,1.0);\n"
-
 
193
    "   float dot_ln = abs(dot(l, n));\n"
-
 
194
    "	vec4 d = vec4(c,1) * 0.7 * dot_ln;\n"
-
 
195
    "	vec4 s = vec4(c,1) * 0.3 * smoothstep(0.98,0.9999,dot(r, e));\n"
-
 
196
    "	\n"
-
 
197
    "	gl_FragColor =  d+s;\n"
-
 
198
    "}\n";
-
 
199
    
-
 
200
    HMesh::VertexAttributeVector<CGLA::Vec3f> DebugRenderer::vertex_colors;
-
 
201
    HMesh::FaceAttributeVector<CGLA::Vec3f> DebugRenderer::face_colors;
-
 
202
 
-
 
203
    
-
 
204
    void DebugRenderer::compile_display_list(const HMesh::Manifold& m, bool smooth)
-
 
205
    {
-
 
206
        GLint old_prog;
-
 
207
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
-
 
208
        glUseProgram(prog);
-
 
209
        glNewList(display_list,GL_COMPILE);
-
 
210
        for(FaceID f: m.faces()){
-
 
211
            Vec3f c = face_colors[f];
-
 
212
            glColor3f(c[0], c[1], c[2]);
-
 
213
            if(!smooth)
-
 
214
                glNormal3dv(normal(m, f).get());
-
 
215
            if(no_edges(m, f)== 3)
-
 
216
                glBegin(GL_TRIANGLES);
-
 
217
            else
-
 
218
                glBegin(GL_POLYGON);
-
 
219
            
-
 
220
            for(Walker w = m.walker(f); !w.full_circle(); w = w.circulate_face_ccw()){
-
 
221
                Vec3d n = normal(m, w.vertex());
-
 
222
                if(smooth)
-
 
223
                    glNormal3dv(n.get());
-
 
224
                glVertex3dv(m.pos(w.vertex()).get());
-
 
225
            }
-
 
226
            glEnd();
-
 
227
        }
-
 
228
 
-
 
229
        Vec3d c;
-
 
230
        float r;
-
 
231
        bsphere(m, c, r);
-
 
232
        for(auto vid : m.vertices())
-
 
233
        {
-
 
234
            Vec3d p = m.pos(vid);
-
 
235
            Vec3f c = vertex_colors[vid];
-
 
236
            glColor3f(c[0], c[1], c[2]);
-
 
237
            glPushMatrix();
-
 
238
            glTranslated(p[0], p[1], p[2]);
-
 
239
            glutSolidSphere(0.01*r, 7, 7);
-
 
240
            glPopMatrix();
-
 
241
        }
-
 
242
        glEndList();
-
 
243
        glUseProgram(old_prog);
-
 
244
    }
161
    
245
    
162
    const string ReflectionLineRenderer::vss = 
246
    const string ReflectionLineRenderer::vss = 
163
    "varying vec3 _n;\n"
247
    "varying vec3 _n;\n"
164
    "varying vec3 v;\n"
248
    "varying vec3 v;\n"
165
    "\n"
249
    "\n"
Line 252... Line 336...
252
    "	\n"
336
    "	\n"
253
    "	gl_FragColor =  d+s;\n"
337
    "	gl_FragColor =  d+s;\n"
254
    "}\n";
338
    "}\n";
255
    
339
    
256
    
340
    
257
    void GlazedRenderer::draw()
341
    void GlazedRenderer::compile_display_list(const HMesh::Manifold& m, bool smooth, float bsphere_rad)
258
    {
342
    {
259
        GLint old_prog;
343
        GLint old_prog;
260
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
344
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
261
        glUseProgram(prog);
345
        glUseProgram(prog);
-
 
346
        glNewList(display_list,GL_COMPILE);
262
        glBindTexture(GL_TEXTURE_3D, get_noise_texture_id());
347
        glBindTexture(GL_TEXTURE_3D, get_noise_texture_id());
263
        glUniform1iARB(glGetUniformLocationARB(prog, "noise_tex"),0);
348
        glUniform1iARB(glGetUniformLocationARB(prog, "noise_tex"),0);
264
        glUniform1fARB(glGetUniformLocationARB(prog, "noise_scale"),12.0/bsphere_rad);
349
        glUniform1fARB(glGetUniformLocationARB(prog, "noise_scale"),12.0/bsphere_rad);
-
 
350
        GLGraphics::draw(m, smooth);
265
        glCallList(display_list);
351
        glEndList();
266
        glUseProgram(old_prog);
352
        glUseProgram(old_prog);
-
 
353
 
267
    }
354
    }
268
    
-
 
269
    
355
 
270
    
356
    
271
    const string GlazedRenderer::vss = 
357
    const string GlazedRenderer::vss = 
272
    "varying vec3 _n;\n"
358
    "varying vec3 _n;\n"
273
    "varying vec3 v;\n"
359
    "varying vec3 v;\n"
274
    "varying vec3 v_obj;\n"
360
    "varying vec3 v_obj;\n"
Line 300... Line 386...
300
    "	return 0.15*refl + diff;\n"
386
    "	return 0.15*refl + diff;\n"
301
    "}\n"
387
    "}\n"
302
    "\n"
388
    "\n"
303
    "void main(void)\n"
389
    "void main(void)\n"
304
    "{\n"
390
    "{\n"
305
    "	vec4 mat_col = vec4(0.9,1.0,0.4,1.0) + vec4(-0.1,-0.1,0.12,0.0) * texture3D(noise_tex, noise_scale*v_obj).x\n"
391
    "	vec4 mat_col = vec4(0.9,1.0,0.4,1.0) +  vec4(-0.1,-0.1,0.12,0.0) * texture3D(noise_tex, noise_scale*v_obj).x\n"
306
    " + vec4(0.05) * texture3D(noise_tex, 500.0*v_obj).x;\n"
392
    " + vec4(0.05) * texture3D(noise_tex, 500.0*v_obj).x;\n"
307
    "	\n"
393
    "	\n"
308
    "	vec3 light0_dir = vec3(0.0,1.0,0.0);\n"
394
    "	vec3 light0_dir = vec3(0.0,1.0,0.0);\n"
309
    "	vec4 light0_col = vec4(0.7,0.9,1.0,1.0);\n"
395
    "	vec4 light0_col = vec4(0.7,0.9,1.0,1.0);\n"
310
    "	\n"
396
    "	\n"
Line 353... Line 439...
353
    "       gl_FragColor.g = pow(gl_FragColor.g, 1.0/gamma);\n"
439
    "       gl_FragColor.g = pow(gl_FragColor.g, 1.0/gamma);\n"
354
    "       gl_FragColor.b = pow(gl_FragColor.b, 1.0/gamma);\n"
440
    "       gl_FragColor.b = pow(gl_FragColor.b, 1.0/gamma);\n"
355
    "		//gl_FragColor += stripe_col * smoothstep(0.8,1.0,cos(stripe_signal));\n"
441
    "		//gl_FragColor += stripe_col * smoothstep(0.8,1.0,cos(stripe_signal));\n"
356
    "	}\n";
442
    "	}\n";
357
    
443
    
358
    ScalarFieldRenderer::ScalarFieldRenderer(const Manifold& m, 
444
    void ScalarFieldRenderer::compile_display_list(const HMesh::Manifold& m, bool smooth,
359
                                             bool smooth, 
-
 
360
                                             VertexAttributeVector<double>& field, 
445
                                                   HMesh::VertexAttributeVector<double>& field, double max_val, float gamma)
361
                                             double max_val, float gamma): SimpleShaderRenderer(vss, fss)
-
 
362
    {
446
    {
363
        
447
        
364
        GLint old_prog;
448
        GLint old_prog;
365
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
449
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
366
        glUseProgram(prog);
450
        glUseProgram(prog);
Line 393... Line 477...
393
        glEndList();	
477
        glEndList();	
394
        glUseProgram(old_prog);
478
        glUseProgram(old_prog);
395
        
479
        
396
    }
480
    }
397
 
481
 
398
    const string PeriodicScalarFieldRenderer::vss =
-
 
399
    "	attribute float scalar;\n"
-
 
400
    "	varying vec3 _normal;\n"
-
 
401
    "	varying float s;\n"
-
 
402
    "	\n"
-
 
403
    "	void main(void)\n"
-
 
404
    "	{\n"
-
 
405
    "		gl_Position =  ftransform();\n"
-
 
406
    "		_normal = normalize(gl_NormalMatrix * gl_Normal);\n"
-
 
407
    "		s=scalar;\n"
-
 
408
    "	}\n";
-
 
409
    
-
 
410
    const string PeriodicScalarFieldRenderer::fss = 	
-
 
411
    "	varying vec3 _normal;\n"
-
 
412
    "	varying float s;\n"
-
 
413
    "   uniform float gamma;\n"
-
 
414
    "	const vec3 light_dir = vec3(0,0,1);\n"
-
 
415
    "	\n"
-
 
416
    "	void main()\n"
-
 
417
    "	{\n"
-
 
418
    "       vec3 normal = normalize(_normal);\n"
-
 
419
    "		float dot_ln = max(0.0,dot(light_dir, normal));\n"
-
 
420
    "		\n"
-
 
421
    "		gl_FragColor = sin(s) * vec4(-1,0,1,0);\n"
-
 
422
    "       gl_FragColor *= dot_ln;\n"
-
 
423
    "       gl_FragColor.r = pow(gl_FragColor.r, 1.0/gamma);\n"
-
 
424
    "       gl_FragColor.g = pow(gl_FragColor.g, 1.0/gamma);\n"
-
 
425
    "       gl_FragColor.b = pow(gl_FragColor.b, 1.0/gamma);\n"
-
 
426
    "	}\n";
-
 
427
    
-
 
428
    PeriodicScalarFieldRenderer::PeriodicScalarFieldRenderer(const Manifold& m, 
-
 
429
                                             bool smooth, 
-
 
430
                                             VertexAttributeVector<double>& field, 
-
 
431
                                             float gamma): SimpleShaderRenderer(vss, fss)
-
 
432
    {
-
 
433
        
-
 
434
        GLint old_prog;
-
 
435
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
-
 
436
        glUseProgram(prog);
-
 
437
        
-
 
438
        GLuint scalar_attrib = glGetAttribLocation(prog, "scalar");
-
 
439
        
-
 
440
        //    static float& gamma = CreateCVar("display.scalar_field_renderer.gamma",2.2f);
-
 
441
        glUniform1fARB(glGetUniformLocationARB(prog, "gamma"), gamma);
-
 
442
        glNewList(display_list,GL_COMPILE);
-
 
443
        
-
 
444
        for(FaceIDIterator f = m.faces_begin(); f != m.faces_end(); ++f){      
-
 
445
            if(!smooth) 
-
 
446
                glNormal3dv(normal(m, *f).get());
-
 
447
            if(no_edges(m, *f)== 3) 
-
 
448
                glBegin(GL_TRIANGLES);
-
 
449
            else 
-
 
450
                glBegin(GL_POLYGON);
-
 
451
            
-
 
452
            
-
 
453
            for(Walker w = m.walker(*f); !w.full_circle(); w = w.circulate_face_ccw()){
-
 
454
                Vec3d n(normal(m, w.vertex()));
-
 
455
                if(smooth) 
-
 
456
                    glNormal3dv(n.get());
-
 
457
                glVertexAttrib1d(scalar_attrib, field[w.vertex()]);
-
 
458
                glVertex3dv(m.pos(w.vertex()).get());
-
 
459
            }
-
 
460
            glEnd();
-
 
461
        }
-
 
462
        glEndList();	
-
 
463
        glUseProgram(old_prog);
-
 
464
        
-
 
465
    }
-
 
466
    
-
 
467
    const string CircleFieldRenderer::vss =
-
 
468
    "	attribute vec2 circlepos;\n"
-
 
469
    "	varying vec3 _normal;\n"
-
 
470
    "	varying vec2 s;\n"
-
 
471
    "	\n"
-
 
472
    "	void main(void)\n"
-
 
473
    "	{\n"
-
 
474
    "		gl_Position =  ftransform();\n"
-
 
475
    "		_normal = normalize(gl_NormalMatrix * gl_Normal);\n"
-
 
476
    "		s=circlepos;\n"
-
 
477
    "	}\n";
-
 
478
    
-
 
479
    const string CircleFieldRenderer::fss =
-
 
480
    "	varying vec3 _normal;\n"
-
 
481
    "	varying vec2 s;\n"
-
 
482
    "   uniform float gamma;\n"
-
 
483
    "	const vec3 light_dir = vec3(0,0,1);\n"
-
 
484
    "	\n"
-
 
485
    "	void main()\n"
-
 
486
    "	{\n"
-
 
487
    "       vec3 normal = normalize(_normal);\n"
-
 
488
    "		float dot_ln = max(0.0,dot(light_dir, normal));\n"
-
 
489
    "		\n"
-
 
490
    "       vec2 _s = normalize(s);\n"
-
 
491
    "       float a = atan(_s.y,_s.x);\n"
-
 
492
    "		gl_FragColor = sin(32.0*a) * vec4(-.7,0,.7,0)+vec4(0.3);\n"
-
 
493
    "       gl_FragColor *= dot_ln*length(s);\n"
-
 
494
    "       gl_FragColor.r = pow(gl_FragColor.r, 1.0/gamma);\n"
-
 
495
    "       gl_FragColor.g = pow(gl_FragColor.g, 1.0/gamma);\n"
-
 
496
    "       gl_FragColor.b = pow(gl_FragColor.b, 1.0/gamma);\n"
-
 
497
    "	}\n";
-
 
498
    
-
 
499
    CircleFieldRenderer::CircleFieldRenderer(const Manifold& m,
-
 
500
                                                             bool smooth,
-
 
501
                                                             VertexAttributeVector<Vec2d>& field,
-
 
502
                                                             float gamma): SimpleShaderRenderer(vss, fss)
-
 
503
    {
-
 
504
        
-
 
505
        GLint old_prog;
-
 
506
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
-
 
507
        glUseProgram(prog);
-
 
508
        
-
 
509
        GLuint scalar_attrib = glGetAttribLocation(prog, "circlepos");
-
 
510
        
-
 
511
        //    static float& gamma = CreateCVar("display.scalar_field_renderer.gamma",2.2f);
-
 
512
        glUniform1fARB(glGetUniformLocationARB(prog, "gamma"), gamma);
-
 
513
        glNewList(display_list,GL_COMPILE);
-
 
514
        
-
 
515
        for(FaceIDIterator f = m.faces_begin(); f != m.faces_end(); ++f){
-
 
516
            if(!smooth)
-
 
517
                glNormal3dv(normal(m, *f).get());
-
 
518
            if(no_edges(m, *f)== 3)
-
 
519
                glBegin(GL_TRIANGLES);
-
 
520
            else
-
 
521
                glBegin(GL_POLYGON);
-
 
522
            
-
 
523
            
-
 
524
            for(Walker w = m.walker(*f); !w.full_circle(); w = w.circulate_face_ccw()){
-
 
525
                Vec3d n(normal(m, w.vertex()));
-
 
526
                if(smooth)
-
 
527
                    glNormal3dv(n.get());
-
 
528
                glVertexAttrib2dv(scalar_attrib, field[w.vertex()].get());
-
 
529
                glVertex3dv(m.pos(w.vertex()).get());
-
 
530
            }
-
 
531
            glEnd();
-
 
532
        }
-
 
533
        glEndList();
-
 
534
        glUseProgram(old_prog);
-
 
535
        
-
 
536
    }
-
 
537
 
482
 
538
 
483
 
539
    const string AmbientOcclusionRenderer::vss =
484
    const string AmbientOcclusionRenderer::vss =
540
    "	attribute float scalar;\n"
485
    "	attribute float scalar;\n"
541
    "	varying vec3 _normal;\n"
486
    "	varying vec3 _normal;\n"
Line 566... Line 511...
566
    "       gl_FragColor.r = pow(gl_FragColor.r, 1.0);\n"
511
    "       gl_FragColor.r = pow(gl_FragColor.r, 1.0);\n"
567
    "       gl_FragColor.g = pow(gl_FragColor.g, 1.0);\n"
512
    "       gl_FragColor.g = pow(gl_FragColor.g, 1.0);\n"
568
    "       gl_FragColor.b = pow(gl_FragColor.b, 1.0);\n"
513
    "       gl_FragColor.b = pow(gl_FragColor.b, 1.0);\n"
569
    "	}\n";
514
    "	}\n";
570
    
515
    
571
    AmbientOcclusionRenderer::AmbientOcclusionRenderer(const Manifold& m, bool smooth, VertexAttributeVector<double>& field, double max_val):
516
    void AmbientOcclusionRenderer::compile_display_list(const HMesh::Manifold& m, HMesh::VertexAttributeVector<double>& field, double max_val)
572
    SimpleShaderRenderer(vss,fss)
-
 
573
    {	
517
    {	
574
        GLint old_prog;
518
        GLint old_prog;
575
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
519
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
576
        glUseProgram(prog);
520
        glUseProgram(prog);
577
        
521
        
Line 579... Line 523...
579
        glUniform1fARB(glGetUniformLocationARB(prog, "scalar_max"), max_val);
523
        glUniform1fARB(glGetUniformLocationARB(prog, "scalar_max"), max_val);
580
        
524
        
581
        glNewList(display_list,GL_COMPILE);
525
        glNewList(display_list,GL_COMPILE);
582
        
526
        
583
        for(FaceIDIterator f = m.faces_begin(); f != m.faces_end(); ++f){
527
        for(FaceIDIterator f = m.faces_begin(); f != m.faces_end(); ++f){
584
            
528
 
585
            if(!smooth) 
-
 
586
                glNormal3dv(normal(m, *f).get());
-
 
587
            if(no_edges(m, *f)== 3) 
529
            if(no_edges(m, *f)== 3)
588
                glBegin(GL_TRIANGLES);
530
                glBegin(GL_TRIANGLES);
589
            else 
531
            else 
590
                glBegin(GL_POLYGON);
532
                glBegin(GL_POLYGON);
591
            
533
            
592
            for(Walker w = m.walker(*f); !w.full_circle(); w = w.circulate_face_ccw())
534
            for(Walker w = m.walker(*f); !w.full_circle(); w = w.circulate_face_ccw())
593
            {
535
            {
594
                Vec3d n(normal(m, w.vertex()));
536
                Vec3d n(normal(m, w.vertex()));
595
                if(smooth) 
-
 
596
                    glNormal3dv(n.get());
537
                glNormal3dv(n.get());
597
                glVertexAttrib1d(scalar_attrib, field[w.vertex()]);
538
                glVertexAttrib1d(scalar_attrib, field[w.vertex()]);
598
                glVertex3dv(m.pos(w.vertex()).get());
539
                glVertex3dv(m.pos(w.vertex()).get());
599
            }
540
            }
600
            glEnd();
541
            glEnd();
601
        }
542
        }
Line 603... Line 544...
603
        glUseProgram(old_prog);
544
        glUseProgram(old_prog);
604
        
545
        
605
    }
546
    }
606
    
547
    
607
    
548
    
608
    LineFieldRenderer::LineFieldRenderer(const Manifold& m, bool smooth, VertexAttributeVector<Vec3d>& lines, float _r): 
549
    void LineFieldRenderer::compile_display_list(const HMesh::Manifold& m,HMesh::VertexAttributeVector<CGLA::Vec3d>& lines, float _r)
609
    SimpleShaderRenderer(vss,fss), r(_r)
-
 
610
    {
550
    {
-
 
551
        glNewList(display_list,GL_COMPILE);
611
        float noise_scale = 10.0f/r;
552
        float noise_scale = 10.0f/r;
612
        float line_scale = 0.003f;
553
        float line_scale = 0.003f;
613
        
554
        
614
        GLint old_prog;
555
        GLint old_prog;
615
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
556
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
616
        glUseProgram(prog);	
557
        glUseProgram(prog);	
617
        glUniform1fARB(glGetUniformLocationARB(prog, "scale_line"),line_scale*noise_scale);
558
        glUniform1fARB(glGetUniformLocationARB(prog, "scale_line"),line_scale*noise_scale);
618
        glUniform1fARB(glGetUniformLocationARB(prog, "noise_scale"),noise_scale);
559
        glUniform1fARB(glGetUniformLocationARB(prog, "noise_scale"),noise_scale);
619
        glUniform1iARB(glGetUniformLocationARB(prog, "noise_tex"),0);
560
        glUniform1iARB(glGetUniformLocationARB(prog, "noise_tex"),0);
620
        GLuint direction = glGetAttribLocation(prog, "direction");	
561
        GLuint direction = glGetAttribLocation(prog, "direction");	
621
        glNewList(display_list,GL_COMPILE);
562
        glBindTexture(GL_TEXTURE_3D, get_noise_texture_id());
-
 
563
 
622
        for(FaceIDIterator f = m.faces_begin(); f != m.faces_end(); ++f){
564
        for(FaceIDIterator f = m.faces_begin(); f != m.faces_end(); ++f){
623
            if(!smooth) 
-
 
624
                glNormal3dv(normal(m, *f).get());
-
 
625
            if(no_edges(m, *f) == 3) 
565
            if(no_edges(m, *f) == 3)
626
                glBegin(GL_TRIANGLES);
566
                glBegin(GL_TRIANGLES);
627
            else 
567
            else 
628
                glBegin(GL_POLYGON);
568
                glBegin(GL_POLYGON);
629
            
569
            
630
            for(Walker w = m.walker(*f); !w.full_circle(); w = w.circulate_face_ccw()){
570
            for(Walker w = m.walker(*f); !w.full_circle(); w = w.circulate_face_ccw()){
631
                Vec3d n(normal(m, w.vertex()));
571
                Vec3d n(normal(m, w.vertex()));
632
                if(smooth) 
-
 
633
                    glNormal3dv(n.get());
572
                glNormal3dv(n.get());
634
                
573
                
635
                Vec3d d = lines[w.vertex()];
574
                Vec3d d = lines[w.vertex()];
636
                d = normalize(d-n*dot(n,d));
575
                d = normalize(d-n*dot(n,d));
637
                glVertexAttrib3dv(direction, d.get());
576
                glVertexAttrib3dv(direction, d.get());
638
                glVertex3dv(m.pos(w.vertex()).get());
577
                glVertex3dv(m.pos(w.vertex()).get());
Line 644... Line 583...
644
        glEndList();	
583
        glEndList();	
645
        glUseProgram(old_prog);
584
        glUseProgram(old_prog);
646
        
585
        
647
    }
586
    }
648
    
587
    
649
    void LineFieldRenderer::draw()
-
 
650
    {
-
 
651
        GLint old_prog;
-
 
652
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
-
 
653
        glUseProgram(prog);
-
 
654
        glBindTexture(GL_TEXTURE_3D, get_noise_texture_id());
-
 
655
        glCallList(display_list);
-
 
656
        glBindTexture(GL_TEXTURE_3D, 0);
-
 
657
        glUseProgram(old_prog);	
-
 
658
    }
-
 
659
    
588
    
660
    const string LineFieldRenderer::vss = 
589
    const string LineFieldRenderer::vss = 
661
    "attribute vec3 direction;\n"
590
    "attribute vec3 direction;\n"
662
    "varying vec3 _n;\n"
591
    "varying vec3 _n;\n"
663
    "varying vec3 dir_obj;\n"
592
    "varying vec3 dir_obj;\n"
Line 701... Line 630...
701
    "   float diff = max(0.0,dot(n,vec3(0.0, 0.0, 1.0)));\n"
630
    "   float diff = max(0.0,dot(n,vec3(0.0, 0.0, 1.0)));\n"
702
    "	gl_FragColor.rgb = vec3(diff*I*(1.0/13.0));\n"
631
    "	gl_FragColor.rgb = vec3(diff*I*(1.0/13.0));\n"
703
    "	gl_FragColor.a = 1.0;\n"
632
    "	gl_FragColor.a = 1.0;\n"
704
    "}\n";
633
    "}\n";
705
    
634
    
706
    const string DualVertexRenderer::vss = 
-
 
707
    "#version 120\n"
-
 
708
    "#extension GL_EXT_gpu_shader4 : enable\n"
-
 
709
    "varying vec4 diffuseIn;\n"
-
 
710
    "varying vec3 _normalIn;\n"
-
 
711
    "void main(void)\n"
-
 
712
    "{\n"
-
 
713
    "   diffuseIn = gl_Color;\n"
-
 
714
    "   _normalIn = normalize(gl_NormalMatrix*gl_Normal);\n"
-
 
715
    "   gl_Position =  ftransform();\n"
-
 
716
    "}\n";
-
 
717
    
-
 
718
    const string DualVertexRenderer::gss = 
-
 
719
    "#version 120\n"
-
 
720
    "#extension GL_EXT_gpu_shader4 : enable\n"
-
 
721
    "#extension GL_EXT_geometry_shader4 : enable\n"
-
 
722
    "\n"
-
 
723
    "varying in vec4 diffuseIn[3];\n"
-
 
724
    "varying in vec3 _normalIn[3];\n"
-
 
725
    "varying vec4 diffuse[3];\n"
-
 
726
    "varying float f;\n"
-
 
727
    "varying vec3 _normal;\n"
-
 
728
    "void main(void)\n"
-
 
729
    "{\n"
-
 
730
    "  diffuse[0] = diffuseIn[0];\n"
-
 
731
    "  diffuse[1] = diffuseIn[1];\n"
-
 
732
    "  diffuse[2] = diffuseIn[2];\n"
-
 
733
    "\n"
-
 
734
    "  f = diffuseIn[0].a;\n"
-
 
735
    "  gl_Position = gl_PositionIn[0];\n"
-
 
736
    "  _normal = _normalIn[0];\n"
-
 
737
    "  EmitVertex();\n"
-
 
738
    "	\n"
-
 
739
    "  f = diffuseIn[1].a;\n"
-
 
740
    "  gl_Position = gl_PositionIn[1];\n"
-
 
741
    "  _normal = _normalIn[1];\n"
-
 
742
    "  EmitVertex();\n"
-
 
743
    "\n"
-
 
744
    "  f = diffuseIn[2].a;\n"
-
 
745
    "  gl_Position = gl_PositionIn[2];\n"
-
 
746
    "  _normal = _normalIn[2];\n"
-
 
747
    "  EmitVertex();\n"
-
 
748
    "\n"
-
 
749
    "  EndPrimitive();\n"
-
 
750
    "}\n";
-
 
751
    
-
 
752
    const string DualVertexRenderer::fss =
-
 
753
    "#version 120\n"
-
 
754
    "#extension GL_EXT_gpu_shader4 : enable\n"
-
 
755
    "\n"
-
 
756
    "varying float f;\n"
-
 
757
    "varying vec4 diffuse[3];\n"
-
 
758
    "varying vec3 _normal;\n"
-
 
759
    "\n"
-
 
760
    "void main(void)\n"
-
 
761
    "{\n"
-
 
762
    "   vec3 normal = normalize(_normal);\n"
-
 
763
    "   float col_idx=0;\n"
-
 
764
    "   if(f>diffuse[0].g && f<diffuse[0].b)\n"
-
 
765
    "      col_idx = diffuse[0].r;\n"
-
 
766
    "   else if(f>diffuse[1].g && f<diffuse[1].b)\n"
-
 
767
    "      col_idx = diffuse[1].r;\n"
-
 
768
    "   else if(f>diffuse[2].g && f<diffuse[2].b)\n"
-
 
769
    "      col_idx = diffuse[2].r;\n"
-
 
770
    "   vec4 col = col_idx < .5 ? vec4(1,0,0,0) : vec4(0,0,1,0);\n"
-
 
771
    "\n"
-
 
772
    " 	gl_FragColor =col*dot(normal, vec3(0,0,1));\n"
-
 
773
    "}\n";
-
 
774
    
-
 
775
    
-
 
776
    
-
 
777
    
-
 
778
    DualVertexRenderer::DualVertexRenderer(const HMesh::Manifold& m, VertexAttributeVector<Vec4d>& field)
-
 
779
    {		
-
 
780
        // Create the program
-
 
781
        static GLuint prog = glCreateProgram();
-
 
782
        
-
 
783
        static bool was_here = false;
-
 
784
        if(!was_here)
-
 
785
        {
-
 
786
            was_here = true;
-
 
787
            // Create s	haders directly from file
-
 
788
            static GLuint vs = create_glsl_shader(GL_VERTEX_SHADER, vss);
-
 
789
            static GLuint gs = create_glsl_shader(GL_GEOMETRY_SHADER_EXT, gss);
-
 
790
            static GLuint fs = create_glsl_shader(GL_FRAGMENT_SHADER, fss);
-
 
791
            
-
 
792
            // Attach all shaders
-
 
793
            if(vs) glAttachShader(prog, vs);
-
 
794
            if(gs) glAttachShader(prog, gs);
-
 
795
            if(fs) glAttachShader(prog, fs);
-
 
796
            
-
 
797
            // Specify input and output for the geometry shader. Note that this must be
-
 
798
            // done before linking the program.
-
 
799
            glProgramParameteriEXT(prog,GL_GEOMETRY_INPUT_TYPE_EXT,GL_TRIANGLES);
-
 
800
            glProgramParameteriEXT(prog,GL_GEOMETRY_VERTICES_OUT_EXT,3);
-
 
801
            glProgramParameteriEXT(prog,GL_GEOMETRY_OUTPUT_TYPE_EXT,GL_TRIANGLE_STRIP);
-
 
802
            
-
 
803
            // Link the program object and print out the info log
-
 
804
            glLinkProgram(prog);
-
 
805
        }
-
 
806
        
-
 
807
        
-
 
808
        
-
 
809
        GLint old_prog;
-
 
810
        glGetIntegerv(GL_CURRENT_PROGRAM, &old_prog);
-
 
811
        
-
 
812
        glNewList(display_list,GL_COMPILE);
-
 
813
        glUseProgram(prog);
-
 
814
        for(FaceIDIterator f = m.faces_begin(); f != m.faces_end(); ++f){
-
 
815
            if(no_edges(m, *f) != 3) 
-
 
816
                continue;
-
 
817
            else 
-
 
818
                glBegin(GL_TRIANGLES);
-
 
819
            
-
 
820
            for(Walker w = m.walker(*f); !w.full_circle(); w = w.circulate_face_ccw()){
-
 
821
                Vec3d n(normal(m, w.vertex()));
-
 
822
                glNormal3dv(n.get());
-
 
823
                glColor4dv(field[w.vertex()].get());
-
 
824
                glVertex3dv(m.pos(w.vertex()).get());
-
 
825
            }
-
 
826
            glEnd();
-
 
827
        }
-
 
828
        glUseProgram(old_prog);
-
 
829
        glEndList();	
-
 
830
        
-
 
831
    }
-
 
832
    
635
    
833
}
636
}
834
 
637
 
835
 
638