25962e3b633e81722c9f96412911016597f8bffc
[occt.git] / src / OpenGl / OpenGl_telem_util.cxx
1 /***********************************************************************
2
3 FONCTION :
4 ----------
5 File OpenGl_telem_util :
6
7
8 REMARQUES:
9 ---------- 
10
11
12 HISTORIQUE DES MODIFICATIONS   :
13 --------------------------------
14 xx-xx-xx : xxx ; Creation.
15 07-02-96 : FMN ; - Ajout trace
16 - Suppression code inutile
17 08-03-96 : FMN ; - Ajout include manquant 
18 01-04-96 : CAL ; Integration MINSK portage WNT
19 15-04-96 : CAL ; Integration travail PIXMAP de Jim ROTH
20 22-04-96 : FMN ; Ajout TelReadImage TelDrawImage
21 10-05-96 : CAL ; Ajout d'un nouveau delta dans les copies
22 de pixels (voir CALL_DEF_DELTA)
23 25-06-96 : FMN ; Suppression utilisation de glScissor.
24 02-07-96 : FMN ; Suppression WSWSHeight et WSWSWidth
25 Suppression glViewport inutile.
26 18-07-96 : FMN ; Suppression TelFlush inutile.
27 08-07-96 : FMN ; Suppression de OPENGL_DEBUG inutile avec la nouvelle
28 version de ogldebug.
29 24-10-96 : CAL ; Portage WNT
30 23-01-97 : CAL ; Suppression de TelClearViews dans TelCopyBuffers
31 30-01-97 : FMN ; Ajout commentaires + WNT.
32 12-02-97 : FMN ; Suppression TelEnquireFacilities()
33 22-04-97 : FMN ; Ajout affichage du cadre pour la copie de buffer
34 30-06-97 : FMN ; Suppression OpenGl_telem_light.h
35 18-07-97 : FMN ; Utilisation de la toolkit sur les lights
36 07-10-97 : FMN ; Simplification WNT + correction Transient
37 05-12-97 : FMN ; PRO11168: Suppression TglActiveWs pour project/unproject
38 23-12-97 : FMN ; Suppression TelSetFrontFaceAttri et TelSetBackFaceAttri 
39 30-12-97 : FMN ; CTS18312: Correction back material
40 04-05-98 : CAL ; Contournement bug SGI octane bavure de pixels (PRO12899)
41 30-09-98 : CAL ; Optimisation pour eviter de charger inutilement
42 les matrices de la vue.
43 19-10-98 : FMN ; Suppression de glPixelTransferi dans TelEnable() car cela
44 rentre en conflit avec l'utilisation d'une image de fond.
45 02.14.100 : JR : Warnings on WNT truncations from double to float
46 08-03-01 : GG  ; BUC60823 Avoid crash in the normal computation method
47 when a face has confused or aligned points.
48
49 ************************************************************************/
50
51 #define IMP190100 /*GG To avoid too many REDRAW in immediat mode,
52 //      Add TelMakeFrontAndBackBufCurrent() function
53 */
54 #define QTOCC_PATCH
55
56 /*----------------------------------------------------------------------*/
57 /*
58 * Includes
59 */
60
61 #include <string.h>
62 #include <stdio.h>
63 #include <stdlib.h>
64
65 #ifndef WNT
66 # include <X11/Xlib.h>
67 #else
68 # define STRICT
69 # include <windows.h>
70 #endif  /* WNT */
71
72
73 #include <GL/gl.h>
74 #include <GL/glu.h>
75 #ifndef WNT
76 #include <GL/glx.h>
77 #endif /* WNT */
78
79 #include <OpenGl_tgl_all.hxx>
80 #include <OpenGl_telem_util.hxx>
81 #include <InterfaceGraphic_Graphic3d.hxx>
82 #include <InterfaceGraphic_Visual3d.hxx>
83
84 #define GPRECIS 0.000001
85 Tint TelGetPolygonNormal(tel_point pnts, Tint* indexs, Tint npnt, Tfloat *norm ) {
86   Tint status=0;
87
88   norm[0] = norm[1] = norm[2] = 0.;
89   if( npnt > 2 ) { 
90     Tfloat a[3], b[3], c[3];
91     Tint i,j,i0,ii=0,jj;
92
93     i0 = 0; if( indexs ) i0 = indexs[0];
94     for( i=1 ; i<npnt ; i++ ) {
95       ii = i; if( indexs ) ii = indexs[i];
96       vecsub( a, pnts[ii].xyz, pnts[i0].xyz );
97       if( vecmg2(a) > GPRECIS ) break;
98     }
99     if( i < npnt-1 ) {
100       for( j=i+1 ; j<npnt ; j++ ) {
101         jj = j; if( indexs ) jj = indexs[j];
102         vecsub( b, pnts[jj].xyz, pnts[i0].xyz );
103         vecsub( c, pnts[jj].xyz, pnts[ii].xyz );
104         if( (vecmg2(b) > GPRECIS) && (vecmg2(c) > GPRECIS) ) break;
105       }
106       if( j < npnt ) {
107         Tfloat d;
108         veccrs( norm, a, b );
109         d = vecnrmd( norm, d );
110         status = (d > 0.) ? 1 : 0;
111       }
112     }
113   }
114 #ifdef DEB
115   if( !status )
116     printf(" *** OpenGl_TelGetPolygonNormal.has found confused or aligned points\n");
117 #endif
118
119   return status;
120 }
121
122 Tint TelGetNormal(Tfloat *data1, Tfloat *data2, Tfloat *data3, Tfloat *norm ) {
123   Tfloat a[3], b[3];
124   Tint status=0;
125
126   norm[0] = norm[1] = norm[2] = 0.;
127   vecsub( a, data2, data1 );
128   vecsub( b, data3, data2 );
129   if( (vecmg2(a) > GPRECIS) && (vecmg2(b) > GPRECIS) ) {
130     Tfloat d;
131     veccrs( norm, a, b );
132     d = vecnrmd( norm, d );
133     status = (d > 0.) ? 1 : 0;
134   }
135 #ifdef DEB
136   if( !status )
137     printf(" *** OpenGl_TelGetNormal.has found confused or aligned points\n");
138 #endif
139
140   return status;
141 }
142
143 void TelMultiplymat3 (Tmatrix3 c, Tmatrix3 a, Tmatrix3 b)
144 {
145   Tint row, col, i;
146   Tmatrix3 res;
147   Tint dim = 4;
148
149   /* on multiplie d'abord les 2 matrices dim x dim */
150   for (row = 0; row < dim; row++) {
151     for (col = 0; col < dim; col++) {
152       Tfloat sum = ( float )0.0;
153       for (i = 0; i < dim; i++)
154         sum += a[row][i] * b[i][col];
155       res[row][col] = sum;
156     }
157   }
158
159   /* on copie ensuite le resultat */
160   matcpy (c, res);
161
162   return;
163 }