intravenous

annotate src/vein.cc @ 3:94d4c60af435

some progress
author John Tsiombikas <nuclear@member.fsf.org>
date Sun, 22 Apr 2012 03:35:18 +0300
parents 3ea290d35984
children c6a6a64df6de
rev   line source
nuclear@3 1 #include <stdio.h>
nuclear@1 2 #ifndef _MSC_VER
nuclear@1 3 #include <alloca.h>
nuclear@1 4 #else
nuclear@1 5 #include <malloc.h>
nuclear@1 6 #endif
nuclear@1 7 #include "vein.h"
nuclear@1 8 #include "geom.h"
nuclear@3 9 #include "sdr.h"
nuclear@1 10
nuclear@1 11 Vein::Vein()
nuclear@1 12 {
nuclear@3 13 gen_dist = 16.0;
nuclear@3 14 subdiv = 16;
nuclear@1 15 ring_subdiv = 16;
nuclear@1 16 rad = 1.0;
nuclear@1 17
nuclear@1 18 idxbuf = 0;
nuclear@1 19 }
nuclear@1 20
nuclear@1 21 Vein::~Vein()
nuclear@1 22 {
nuclear@1 23 delete [] idxbuf;
nuclear@1 24 }
nuclear@1 25
nuclear@1 26 Vector3 Vein::calc_center(const Vector3 &ppos) const
nuclear@1 27 {
nuclear@3 28 Vector3 pt{0, 0, ppos.z};
nuclear@3 29 pt.x = sin(ppos.z * 0.75);
nuclear@3 30 pt.y = cos(ppos.z * 0.2) * 0.6;
nuclear@3 31 return pt;
nuclear@1 32 }
nuclear@1 33
nuclear@1 34 Vector3 Vein::calc_dir(const Vector3 &ppos) const
nuclear@1 35 {
nuclear@3 36 Vector3 dir = calc_center(ppos + Vector3(0, 0, 0.01)) - calc_center(ppos - Vector3(0, 0, 0.01));
nuclear@3 37 return dir.normalized();
nuclear@1 38 }
nuclear@1 39
nuclear@1 40 void Vein::build_idxbuf()
nuclear@1 41 {
nuclear@1 42 delete [] idxbuf;
nuclear@1 43
nuclear@1 44 int nfaces = subdiv * ring_subdiv;
nuclear@1 45 int nidx = nfaces * 4;
nuclear@1 46 idxbuf = new unsigned int[nidx];
nuclear@1 47 unsigned int *idxptr = idxbuf;
nuclear@1 48
nuclear@1 49 for(int i=0; i<subdiv; i++) {
nuclear@1 50 for(int j=0; j<ring_subdiv; j++) {
nuclear@1 51 idxptr[0] = i * ring_subdiv + j;
nuclear@1 52 idxptr[1] = i * ring_subdiv + ((j + 1) % ring_subdiv);
nuclear@1 53 idxptr[2] = idxptr[1] + ring_subdiv;
nuclear@1 54 idxptr[3] = idxptr[0] + ring_subdiv;
nuclear@1 55 idxptr += 4;
nuclear@1 56 }
nuclear@1 57 }
nuclear@1 58 }
nuclear@1 59
nuclear@3 60 bool Vein::init()
nuclear@3 61 {
nuclear@3 62 if(!(sdr = create_program_load("sdr/vein.v.glsl", "sdr/vein.p.glsl"))) {
nuclear@3 63 fprintf(stderr, "failed to load vein shaders\n");
nuclear@3 64 return false;
nuclear@3 65 }
nuclear@3 66 if((attr_tang_loc = get_attrib_loc(sdr, "attr_tang")) == -1) {
nuclear@3 67 fprintf(stderr, "can't find tangent attribute!\n");
nuclear@3 68 }
nuclear@3 69 return true;
nuclear@3 70 }
nuclear@3 71
nuclear@1 72 void Vein::draw(const Vector3 &ppos) const
nuclear@1 73 {
nuclear@1 74 float start_z = ppos.z - gen_dist / 2.0;
nuclear@1 75 float dz = gen_dist / subdiv;
nuclear@1 76
nuclear@1 77 int nslices = subdiv + 1;
nuclear@1 78 int nverts = nslices * ring_subdiv;
nuclear@1 79 int nfaces = subdiv * ring_subdiv;
nuclear@1 80 Vertex *vbuf = (Vertex*)alloca(nverts * sizeof *vbuf);
nuclear@1 81 Vertex *vptr = vbuf;
nuclear@1 82
nuclear@1 83 Vector3 pt = ppos;
nuclear@1 84 pt.z = start_z;
nuclear@1 85
nuclear@1 86 for(int i=0; i<nslices; i++) {
nuclear@1 87 Vector3 cent = calc_center(pt);
nuclear@1 88 Vector3 dir = calc_dir(pt);
nuclear@1 89 Vector3 up(0, 1, 0);
nuclear@1 90 Vector3 right = cross_product(up, dir);
nuclear@1 91 up = cross_product(dir, right);
nuclear@1 92
nuclear@1 93 Matrix3x3 vrot{right, up, dir};
nuclear@3 94 vrot.transpose();
nuclear@1 95
nuclear@1 96 float theta = 0.0, dtheta = 2.0 * M_PI / ring_subdiv;
nuclear@1 97 for(int j=0; j<ring_subdiv; j++) {
nuclear@3 98 Vector3 vec = Vector3{-cos(theta) * rad, sin(theta) * rad, 0.0};
nuclear@1 99 vec.transform(vrot);
nuclear@3 100 vec += cent;
nuclear@1 101
nuclear@1 102 vptr->pos = vec;
nuclear@1 103 vptr->norm = cent - vec;
nuclear@3 104 vptr->tang = dir;
nuclear@1 105 vptr->tc = Vector2(); // TODO
nuclear@1 106 vptr++;
nuclear@1 107
nuclear@1 108 theta += dtheta;
nuclear@1 109 }
nuclear@1 110
nuclear@1 111 pt.z += dz;
nuclear@1 112 }
nuclear@1 113
nuclear@1 114 // also create the index buffer if it's not valid
nuclear@1 115 if(!idxbuf) {
nuclear@1 116 ((Vein*)this)->build_idxbuf();
nuclear@1 117 }
nuclear@1 118
nuclear@1 119 // awesome, now draw it
nuclear@3 120 bind_program(sdr);
nuclear@3 121 draw_mesh(GL_QUADS, nfaces * 4, vbuf, idxbuf, attr_tang_loc);
nuclear@3 122 bind_program(0);
nuclear@1 123 }