Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
ElmerCSC
GitHub Repository: ElmerCSC/elmerfem
 Path: blob/devel/post/src/glaux/teapot.c
3203 views
1
/* teapot.c */

2


3
#include <GL/gl.h>

4
#include "glaux.h"

5


6


7
/* Rim, body, lid, and bottom data must be reflected in x and y; 

8
 * handle and spout data across the y axis only.

9
 */

10


11
static long patchdata[][16] = {

12
    {102,103,104,105,4,5,6,7,8,9,10,11,12,13,14,15},		   /* rim */

13
    {12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27},		   /* body */

14
    {24,25,26,27,29,30,31,32,33,34,35,36,37,38,39,40},		   /* body */

15
    {96,96,96,96,97,98,99,100,101,101,101,101,0,1,2,3,},	   /* lid */

16
    {0,1,2,3,106,107,108,109,110,111,112,113,114,115,116,117},	   /* lid */

17
    {118,118,118,118,124,122,119,121,123,126,125,120,40,39,38,37}, /* bottom */

18
    {41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56},		   /* handle */

19
    {53,54,55,56,57,58,59,60,61,62,63,64,28,65,66,67},		   /* handle */

20
    {68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83},		   /* spout */

21
    {80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95}		   /* spout */

22
};

23


24
static float cpdata[][3] = {

25
{0.2,0,2.7},{0.2,-0.112,2.7},{0.112,-0.2,2.7},{0,-0.2,2.7},{1.3375,0,2.53125},

26
{1.3375,-0.749,2.53125},{0.749,-1.3375,2.53125},{0,-1.3375,2.53125},

27
{1.4375,0,2.53125},{1.4375,-0.805,2.53125},{0.805,-1.4375,2.53125},

28
{0,-1.4375,2.53125},{1.5,0,2.4},{1.5,-0.84,2.4},{0.84,-1.5,2.4},{0,-1.5,2.4},

29
{1.75,0,1.875},{1.75,-0.98,1.875},{0.98,-1.75,1.875},{0,-1.75,1.875},{2,0,1.35},

30
{2,-1.12,1.35},{1.12,-2,1.35},{0,-2,1.35},{2,0,0.9},{2,-1.12,0.9},{1.12,-2,0.9},

31
{0,-2,0.9},{-2,0,0.9},{2,0,0.45},{2,-1.12,0.45},{1.12,-2,0.45},{0,-2,0.45},

32
{1.5,0,0.225},{1.5,-0.84,0.225},{0.84,-1.5,0.225},{0,-1.5,0.225},{1.5,0,0.15},

33
{1.5,-0.84,0.15},{0.84,-1.5,0.15},{0,-1.5,0.15},{-1.6,0,2.025},{-1.6,-0.3,2.025},

34
{-1.5,-0.3,2.25},{-1.5,0,2.25},{-2.3,0,2.025},{-2.3,-0.3,2.025},{-2.5,-0.3,2.25},

35
{-2.5,0,2.25},{-2.7,0,2.025},{-2.7,-0.3,2.025},{-3,-0.3,2.25},{-3,0,2.25},

36
{-2.7,0,1.8},{-2.7,-0.3,1.8},{-3,-0.3,1.8},{-3,0,1.8},{-2.7,0,1.575},

37
{-2.7,-0.3,1.575},{-3,-0.3,1.35},{-3,0,1.35},{-2.5,0,1.125},{-2.5,-0.3,1.125},

38
{-2.65,-0.3,0.9375},{-2.65,0,0.9375},{-2,-0.3,0.9},{-1.9,-0.3,0.6},{-1.9,0,0.6},

39
{1.7,0,1.425},{1.7,-0.66,1.425},{1.7,-0.66,0.6},{1.7,0,0.6},{2.6,0,1.425},

40
{2.6,-0.66,1.425},{3.1,-0.66,0.825},{3.1,0,0.825},{2.3,0,2.1},{2.3,-0.25,2.1},

41
{2.4,-0.25,2.025},{2.4,0,2.025},{2.7,0,2.4},{2.7,-0.25,2.4},{3.3,-0.25,2.4},

42
{3.3,0,2.4},{2.8,0,2.475},{2.8,-0.25,2.475},{3.525,-0.25,2.49375},

43
{3.525,0,2.49375},{2.9,0,2.475},{2.9,-0.15,2.475},{3.45,-0.15,2.5125},

44
{3.45,0,2.5125},{2.8,0,2.4},{2.8,-0.15,2.4},{3.2,-0.15,2.4},{3.2,0,2.4},

45
{0,0,3.15},{0.8,0,3.15},{0.8,-0.45,3.15},{0.45,-0.8,3.15},{0,-0.8,3.15},

46
{0,0,2.85},{1.4,0,2.4},{1.4,-0.784,2.4},{0.784,-1.4,2.4},{0,-1.4,2.4},

47
{0.4,0,2.55},{0.4,-0.224,2.55},{0.224,-0.4,2.55},{0,-0.4,2.55},{1.3,0,2.55},

48
{1.3,-0.728,2.55},{0.728,-1.3,2.55},{0,-1.3,2.55},{1.3,0,2.4},{1.3,-0.728,2.4},

49
{0.728,-1.3,2.4},{0,-1.3,2.4},{0,0,0},{1.425,-0.798,0},{1.5,0,0.075},{1.425,0,0},

50
{0.798,-1.425,0},{0,-1.5,0.075},{0,-1.425,0},{1.5,-0.84,0.075},{0.84,-1.5,0.075}

51
};

52


53
long GRD;

54


55
#define TEAPOTSOLID 0

56
#define TEAPOTWIRE 1

57


58
static GLuint teapots[2] = {0, 0};

59


60
static float tex[2][2][2] = {{{0, 0},{1, 0}},{{0, 1},{1, 1}}};

61


62
static void solidTeapot(long grid, GLdouble scale)

63
{

64
    float p[4][4][3], q[4][4][3], r[4][4][3], s[4][4][3];

65
    long i, j, k, l;

66


67
    if (grid < 2) grid = 7;

68
    GRD = grid;

69
    teapots[TEAPOTSOLID] = glGenLists (1);

70
    glNewList(teapots[TEAPOTSOLID], GL_COMPILE);

71
    glPushMatrix ();

72
    glRotatef (270.0, 1.0, 0.0, 0.0);

73
    glScalef (0.5*scale, 0.5*scale, 0.5*scale);

74
    glTranslatef (0.0, 0.0, -1.5);

75
    for (i = 0; i < 10; i++) {

76
	for (j = 0; j < 4; j++)

77
	    for (k = 0; k < 4; k++) 

78
		for (l = 0; l < 3; l++) {

79
		    p[j][k][l] = cpdata[patchdata[i][j*4+k]][l];

80
		    q[j][k][l] = cpdata[patchdata[i][j*4+(3-k)]][l];

81
		    if (l == 1) q[j][k][l] *= -1.0;

82
		    if (i < 6) {

83
			r[j][k][l] = cpdata[patchdata[i][j*4+(3-k)]][l];

84
			if (l == 0) r[j][k][l] *= -1.0;

85
			s[j][k][l] = cpdata[patchdata[i][j*4+k]][l];

86
			if (l == 0) s[j][k][l] *= -1.0;

87
			if (l == 1) s[j][k][l] *= -1.0;

88
		    }

89
		}

90
	glMap2f(GL_MAP2_TEXTURE_COORD_2, 0, 1, 2, 2, 0, 1, 4, 2, &tex[0][0][0]);

91
	glMap2f(GL_MAP2_VERTEX_3, 0, 1, 3, 4, 0, 1, 12, 4, &p[0][0][0]);

92
	glEnable(GL_MAP2_VERTEX_3); glEnable(GL_MAP2_TEXTURE_COORD_2);

93
	glMapGrid2f(GRD, 0.0, 1.0, GRD, 0.0, 1.0);

94
	glEvalMesh2(GL_FILL, 0, GRD, 0, GRD);

95
	glMap2f(GL_MAP2_VERTEX_3, 0, 1, 3, 4, 0, 1, 12, 4, &q[0][0][0]);

96
	glEvalMesh2(GL_FILL, 0, GRD, 0, GRD);

97
	if (i < 6) {

98
	    glMap2f(GL_MAP2_VERTEX_3, 0, 1, 3, 4, 0, 1, 12, 4, &r[0][0][0]);

99
	    glEvalMesh2(GL_FILL, 0, GRD, 0, GRD);

100
	    glMap2f(GL_MAP2_VERTEX_3, 0, 1, 3, 4, 0, 1, 12, 4, &s[0][0][0]);

101
	    glEvalMesh2(GL_FILL, 0, GRD, 0, GRD);

102
	}

103
    }

104
    glDisable(GL_MAP2_VERTEX_3); glDisable(GL_MAP2_TEXTURE_COORD_2);

105
    glPopMatrix ();

106
    glEndList();

107
}

108


109
static void wireTeapot(long grid, GLdouble scale)

110
{

111
    float p[4][4][3], q[4][4][3], r[4][4][3], s[4][4][3];

112
    long i, j, k, l;

113
    

114
    if (grid < 2) grid = 7;

115
    GRD = grid;

116
    teapots[TEAPOTWIRE] = glGenLists (1);

117
    glNewList(teapots[TEAPOTWIRE], GL_COMPILE);

118
    glPushMatrix ();

119
    glRotatef (270.0, 1.0, 0.0, 0.0);

120
    glScalef (0.5*scale, 0.5*scale, 0.5*scale);

121
    glTranslatef (0.0, 0.0, -1.5);

122
    for (i = 0; i < 10; i++) {

123
	for (j = 0; j < 4; j++)

124
	    for (k = 0; k < 4; k++) 

125
		for (l = 0; l < 3; l++) {

126
		    p[j][k][l] = cpdata[patchdata[i][j*4+k]][l];

127
		    q[j][k][l] = cpdata[patchdata[i][j*4+(3-k)]][l];

128
		    if (l == 1) q[j][k][l] *= -1.0;

129
		    if (i < 6) {

130
			r[j][k][l] = cpdata[patchdata[i][j*4+(3-k)]][l];

131
			if (l == 0) r[j][k][l] *= -1.0;

132
			s[j][k][l] = cpdata[patchdata[i][j*4+k]][l];

133
			if (l == 0) s[j][k][l] *= -1.0;

134
			if (l == 1) s[j][k][l] *= -1.0;

135
		    }

136
		}

137
	glMap2f(GL_MAP2_TEXTURE_COORD_2, 0, 1, 2, 2, 0, 1, 4, 2, &tex[0][0][0]);

138
	glMap2f(GL_MAP2_VERTEX_3, 0, 1, 3, 4, 0, 1, 12, 4, &p[0][0][0]);

139
	glEnable(GL_MAP2_VERTEX_3); glEnable(GL_MAP2_TEXTURE_COORD_2);

140
	glMapGrid2f(GRD, 0.0, 1.0, GRD, 0.0, 1.0);

141
	glEvalMesh2(GL_LINE, 0, GRD, 0, GRD);

142
	glMap2f(GL_MAP2_VERTEX_3, 0, 1, 3, 4, 0, 1, 12, 4, &q[0][0][0]);

143
	glEvalMesh2(GL_LINE, 0, GRD, 0, GRD);

144
	if (i < 6) {

145
	    glMap2f(GL_MAP2_VERTEX_3, 0, 1, 3, 4, 0, 1, 12, 4, &r[0][0][0]);

146
	    glEvalMesh2(GL_LINE, 0, GRD, 0, GRD);

147
	    glMap2f(GL_MAP2_VERTEX_3, 0, 1, 3, 4, 0, 1, 12, 4, &s[0][0][0]);

148
	    glEvalMesh2(GL_LINE, 0, GRD, 0, GRD);

149
	}

150
    }

151
    glDisable(GL_MAP2_VERTEX_3); glDisable(GL_MAP2_TEXTURE_COORD_2);

152
    glPopMatrix ();

153
    glEndList();

154
}

155


156
void auxSolidTeapot(GLdouble scale)

157
{

158
    if (glIsList(teapots[TEAPOTSOLID]) == 0)

159
	solidTeapot (14, scale);

160
    glCallList(teapots[TEAPOTSOLID]);

161
}

162


163
void auxWireTeapot(GLdouble scale)

164
{

165
    if (glIsList(teapots[TEAPOTWIRE]) == 0)

166
	wireTeapot (10, scale);

167
    glCallList(teapots[TEAPOTWIRE]);

168
}

169


170


171