691ad3e7b3aa98db661087864cc63c3fd68ef76b
[occt.git] / src / OpenGl / OpenGl_telem_util.hxx
1 // Copyright (c) 1995-1999 Matra Datavision
2 // Copyright (c) 1999-2014 OPEN CASCADE SAS
3 //
4 // This file is part of Open CASCADE Technology software library.
5 //
6 // This library is free software; you can redistribute it and/or modify it under
7 // the terms of the GNU Lesser General Public License version 2.1 as published
8 // by the Free Software Foundation, with special exception defined in the file
9 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
10 // distribution for complete text of the license and disclaimer of any warranty.
11 //
12 // Alternatively, this file may be used under the terms of Open CASCADE
13 // commercial license or contractual agreement.
14
15 /***********************************************************************
16
17 FONCTION :
18 ----------
19 File Include OpenGl_telem_util :
20
21
22 REMARQUES:
23 ---------- 
24
25
26 HISTORIQUE DES MODIFICATIONS   :
27 --------------------------------
28 xx-xx-xx : xxx ; Creation.
29 22-04-96 : FMN ; Ajout TelReadImage TelDrawImage
30 23-12-96 : CAL ; Probleme dans la macro vecmag
31 12-02-97 : FMN ; Suppression TelEnquireFacilities()
32 07-10-97 : FMN ; Simplification WNT
33 23-12-97 : FMN ; Suppression TelBackInteriorStyle, TelBackInteriorStyleIndex
34
35 ************************************************************************/
36
37 #ifndef  OPENGL_TELEM_UTIL_H
38 #define  OPENGL_TELEM_UTIL_H
39
40 #include <InterfaceGraphic_telem.hxx>
41 #include <cmath>
42
43 /*
44 * ShortRealLast () = 3.40282346638528860e+38
45 * ShortRealFirst () = 3.40282346638528860e+38
46 */
47 #define  shortreallast()     (3.e+38)
48 #define  shortrealfirst()    (-3.e+38)
49
50 /* vector manipulation macros */
51 #define  square(a)      ((a)*(a))
52
53 /* add */
54 #define  vecadd(a,b,c)  { (a)[0]=(b)[0]+(c)[0]; \
55   (a)[1]=(b)[1]+(c)[1]; \
56   (a)[2]=(b)[2]+(c)[2]; }
57
58 /* subtract */
59 #define  vecsub(a,b,c)  { (a)[0]=(b)[0]-(c)[0]; \
60   (a)[1]=(b)[1]-(c)[1]; \
61   (a)[2]=(b)[2]-(c)[2]; }
62
63 /* dot product */
64 #define  vecdot(a,b)    ((a)[0]*(b)[0]+(a)[1]*(b)[1]+(a)[2]*(b)[2])
65
66 /* cross product */
67 #define  veccrs(a,b,c)  { Tfloat x, y, z; \
68   x = (b)[1]*(c)[2] - (b)[2]*(c)[1]; \
69   y = (b)[2]*(c)[0] - (b)[0]*(c)[2]; \
70   z = (b)[0]*(c)[1] - (b)[1]*(c)[0]; \
71   (a)[0] = x; \
72   (a)[1] = y; \
73   (a)[2] = z; }
74 /* scale */
75 #define  vecscl(a,b)     { (b)[0] *= (a); (b)[1] *= (a); (b)[2] *= (a); }
76
77 /* magnitude square */
78 #define  vecmg2(a)       (square((a)[0])+square((a)[1])+square((a)[2]))
79
80 /* magnitude */
81 #define  vecmag(a)       (std::sqrt((double)vecmg2(a)))
82
83 /* normalize */
84 #define  vecnrmd(a,d)    ( d = (Tfloat)vecmag(a), \
85   ( d > 1e-10 ? (a[0] /= d, a[1] /= d, a[2] /= d, d) : (Tfloat)0. ) )
86 #define  vecnrm(a)       { Tfloat d; vecnrmd(a,d); }
87
88 /* angle between two vectors */
89 #define  vecang(a,b,d)   { d = (Tfloat)(vecmag(a)*vecmag(b)); \
90   d = vecdot(a,b)/d; \
91   d = d < -1.0F ? -1.0F : d > 1.0F ? 1.0F : d; \
92   d = ( Tfloat )acos(d); }
93
94 /* point along a vector at a distance */
95 #define  vecgnd(a,b,c,d) { Tfloat w; w = d/vecmag(c); \
96   (a)[0] = (b)[0]+(c)[0]*w;  \
97   (a)[1] = (b)[1]+(c)[1]*w;  \
98   (a)[2] = (b)[2]+(c)[2]*w;  }
99 /* copy */
100 #define  veccpy(a,b)     ((a)[0]=(b)[0],(a)[1]=(b)[1],(a)[2]=(b)[2])
101
102 /* end vector macros */
103
104 typedef  struct
105 {
106   Tmatrix3  mat;
107 }
108 Tmatrix3Struct;
109
110 #define  matcpy(d,s) { *((Tmatrix3Struct*)(d)) = *((Tmatrix3Struct*)(s)); }
111 /*
112 (d)[0][0]=(s)[0][0],(d)[0][1]=(s)[0][1],(d)[0][2]=(s)[0][2],(d)[0][3]=(s)[0][3],
113 (d)[1][0]=(s)[1][0],(d)[1][1]=(s)[1][1],(d)[1][2]=(s)[1][2],(d)[1][3]=(s)[1][3],
114 (d)[2][0]=(s)[2][0],(d)[2][1]=(s)[2][1],(d)[2][2]=(s)[2][2],(d)[2][3]=(s)[2][3],
115 (d)[3][0]=(s)[3][0],(d)[3][1]=(s)[3][1],(d)[3][2]=(s)[3][2],(d)[3][3]=(s)[3][3]
116 */
117
118 #define matdump(m) { \
119   int i, j; \
120   for (i=0; i<4; i++) {\
121   printf ("\t"); \
122   for (j=0; j<4; j++) \
123   printf ("%f ", m[i][j]); \
124   printf ("\n"); \
125   } \
126 }
127
128 extern int TelGetPolygonNormal( tel_point, Tint*, Tint, Tfloat* );
129 extern int TelGetNormal( Tfloat*, Tfloat*, Tfloat*, Tfloat* );
130 extern void TelMultiplymat3( Tmatrix3, Tmatrix3, Tmatrix3 );
131
132 #endif