PageRenderTime 57ms CodeModel.GetById 0ms RepoModel.GetById 0ms app.codeStats 0ms

/indra/llmath/tests/v3dmath_test.cpp

https://bitbucket.org/lindenlab/viewer-beta/
C++ | 526 lines | 472 code | 27 blank | 27 comment | 108 complexity | 72071fa6258f0e102d82154e94d45e9f MD5 | raw file
Possible License(s): LGPL-2.1 
    

  1. /**
    2. 

  2.  * @file v3dmath_test.cpp
    3. 

  3.  * @author Adroit
    4. 

  4.  * @date 2007-03
    5. 

  5.  * @brief v3dmath test cases.
    6. 

  6.  *
    7. 

  7.  * $LicenseInfo:firstyear=2007&license=viewerlgpl$
    8. 

  8.  * Second Life Viewer Source Code
    9. 

  9.  * Copyright (C) 2010, Linden Research, Inc.
    10. 

  10.  * 
    11. 

  11.  * This library is free software; you can redistribute it and/or
    12. 

  12.  * modify it under the terms of the GNU Lesser General Public
    13. 

  13.  * License as published by the Free Software Foundation;
    14. 

  14.  * version 2.1 of the License only.
    15. 

  15.  * 
    16. 

  16.  * This library is distributed in the hope that it will be useful,
    17. 

  17.  * but WITHOUT ANY WARRANTY; without even the implied warranty of
    18. 

  18.  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    19. 

  19.  * Lesser General Public License for more details.
    20. 

  20.  * 
    21. 

  21.  * You should have received a copy of the GNU Lesser General Public
    22. 

  22.  * License along with this library; if not, write to the Free Software
    23. 

  23.  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
    24. 

  24.  * 
    25. 

  25.  * Linden Research, Inc., 945 Battery Street, San Francisco, CA  94111  USA
    26. 

  26.  * $/LicenseInfo$
    27. 

  27.  */
    28. 

  28.  
    29. 

  29. #include "linden_common.h"
    30. 

  30. #include "llsd.h"
    31. 

  31. #include "../test/lltut.h"
    32. 

  32. 

  33. #include "../m3math.h"
    34. 

  34. #include "../v4math.h"
    35. 

  35. #include "../v3dmath.h"
    36. 

  36. #include "../v3dmath.h"
    37. 

  37. #include "../llquaternion.h"
    38. 

  38. 

  39. namespace tut
    40. 

  40. {
    41. 

  41. 	struct v3dmath_data
    42. 

  42. 	{
    43. 

  43. 	};
    44. 

  44. 	typedef test_group<v3dmath_data> v3dmath_test;
    45. 

  45. 	typedef v3dmath_test::object v3dmath_object;
    46. 

  46. 	tut::v3dmath_test v3dmath_testcase("v3dmath_h");
    47. 

  47. 

  48. 	template<> template<>
    49. 

  49. 	void v3dmath_object::test<1>()
    50. 

  50. 	{
    51. 

  51. 		LLVector3d vec3D;
    52. 

  52. 		ensure("1:LLVector3d:Fail to initialize ", ((0 == vec3D.mdV[VX]) && (0 == vec3D.mdV[VY]) && (0 == vec3D.mdV[VZ])));
    53. 

  53. 		F64 x = 2.32f, y = 1.212f, z = -.12f;
    54. 

  54. 		LLVector3d vec3Da(x,y,z);
    55. 

  55. 		ensure("2:LLVector3d:Fail to initialize ", ((2.32f == vec3Da.mdV[VX]) && (1.212f == vec3Da.mdV[VY]) && (-.12f == vec3Da.mdV[VZ])));
    56. 

  56. 		const F64 vec[3] = {1.2f ,3.2f, -4.2f};
    57. 

  57. 		LLVector3d vec3Db(vec);
    58. 

  58. 		ensure("3:LLVector3d:Fail to initialize ", ((1.2f == vec3Db.mdV[VX]) && (3.2f == vec3Db.mdV[VY]) && (-4.2f == vec3Db.mdV[VZ])));
    59. 

  59. 		LLVector3 vec3((F32)x,(F32)y,(F32)z);
    60. 

  60. 		LLVector3d vec3Dc(vec3);
    61. 

  61. 		ensure_equals("4:LLVector3d Fail to initialize",vec3Da,vec3Dc);
    62. 

  62. 	}
    63. 

  63. 

  64. 	template<> template<>
    65. 

  65. 	void v3dmath_object::test<2>()
    66. 

  66. 	{
    67. 

  67. 		S32 a = -235;
    68. 

  68. 		LLSD llsd(a);
    69. 

  69. 		LLVector3d vec3d(llsd);
    70. 

  70. 		LLSD sd = vec3d.getValue();
    71. 

  71. 		LLVector3d vec3da(sd);
    72. 

  72. 		ensure("1:getValue:Fail ", (vec3d == vec3da));
    73. 

  73. 	}
    74. 

  74. 

  75. 	template<> template<>
    76. 

  76. 	void v3dmath_object::test<3>()
    77. 

  77. 	{
    78. 

  78. 		F64 a = 232345521.411132;
    79. 

  79. 		LLSD llsd(a);
    80. 

  80. 		LLVector3d vec3d;
    81. 

  81. 		vec3d.setValue(llsd);
    82. 

  82. 		LLSD sd = vec3d.getValue();
    83. 

  83. 		LLVector3d vec3da(sd);
    84. 

  84. 		ensure("1:setValue:Fail to initialize ", (vec3d == vec3da));
    85. 

  85. 	}
    86. 

  86. 

  87. 	template<> template<>
    88. 

  88. 	void v3dmath_object::test<4>()
    89. 

  89. 	{
    90. 

  90. 		F64 a[3] = {222231.43222, 12345.2343, -434343.33222};
    91. 

  91. 		LLSD llsd;
    92. 

  92. 		llsd[0] = a[0];
    93. 

  93. 		llsd[1] = a[1];
    94. 

  94. 		llsd[2] = a[2];
    95. 

  95. 		LLVector3d vec3D;
    96. 

  96. 		vec3D = (LLVector3d)llsd;
    97. 

  97. 		ensure("1:operator=:Fail to initialize ", ((llsd[0].asReal()== vec3D.mdV[VX]) && (llsd[1].asReal() == vec3D.mdV[VY]) && (llsd[2].asReal() == vec3D.mdV[VZ])));
    98. 

  98. 	}
    99. 

  99. 	
    100. 

  100. 	template<> template<>
    101. 

  101. 	void v3dmath_object::test<5>()
    102. 

  102. 	{
    103. 

  103. 		F64 x = 2.32f, y = 1.212f, z = -.12f;
    104. 

  104. 		LLVector3d vec3D(x,y,z);
    105. 

  105. 		vec3D.clearVec();
    106. 

  106. 		ensure("1:clearVec:Fail to initialize ", ((0 == vec3D.mdV[VX]) && (0 == vec3D.mdV[VY]) && (0 == vec3D.mdV[VZ])));
    107. 

  107. 		vec3D.setVec(x,y,z);
    108. 

  108. 		ensure("2:setVec:Fail to initialize ", ((x == vec3D.mdV[VX]) && (y == vec3D.mdV[VY]) && (z == vec3D.mdV[VZ])));
    109. 

  109. 		vec3D.zeroVec();
    110. 

  110. 		ensure("3:zeroVec:Fail to initialize ", ((0 == vec3D.mdV[VX]) && (0 == vec3D.mdV[VY]) && (0 == vec3D.mdV[VZ])));
    111. 

  111. 		vec3D.clearVec();
    112. 

  112. 		LLVector3 vec3((F32)x,(F32)y,(F32)z);
    113. 

  113. 		vec3D.setVec(vec3);
    114. 

  114. 		ensure("4:setVec:Fail to initialize ", ((x == vec3D.mdV[VX]) && (y == vec3D.mdV[VY]) && (z == vec3D.mdV[VZ])));
    115. 

  115. 		vec3D.clearVec();
    116. 

  116. 		const F64 vec[3] = {x,y,z};
    117. 

  117. 		vec3D.setVec(vec);
    118. 

  118. 		ensure("5:setVec:Fail to initialize ", ((x == vec3D.mdV[VX]) && (y == vec3D.mdV[VY]) && (z == vec3D.mdV[VZ])));
    119. 

  119. 		LLVector3d vec3Da;
    120. 

  120. 		vec3Da.setVec(vec3D);
    121. 

  121. 		ensure_equals("6:setVec: Fail to initialize", vec3D, vec3Da);
    122. 

  122. 	}
    123. 

  123. 	
    124. 

  124. 	template<> template<>
    125. 

  125. 	void v3dmath_object::test<6>()
    126. 

  126. 	{
    127. 

  127. 		F64 x = -2.32, y = 1.212, z = -.12;
    128. 

  128. 		LLVector3d vec3D(x,y,z);
    129. 

  129. 		vec3D.abs();
    130. 

  130. 		ensure("1:abs:Fail  ", ((-x == vec3D.mdV[VX]) && (y == vec3D.mdV[VY]) && (-z == vec3D.mdV[VZ])));
    131. 

  131. 		ensure("2:isNull():Fail ", (FALSE == vec3D.isNull()));	
    132. 

  132. 		vec3D.clearVec();
    133. 

  133. 		x =.00000001, y = .000001001, z = .000001001;
    134. 

  134. 		vec3D.setVec(x,y,z);
    135. 

  135. 		ensure("3:isNull():Fail ", (TRUE == vec3D.isNull()));	
    136. 

  136. 		ensure("4:isExactlyZero():Fail ", (FALSE == vec3D.isExactlyZero()));	
    137. 

  137. 		x =.0000000, y = .00000000, z = .00000000;
    138. 

  138. 		vec3D.setVec(x,y,z);
    139. 

  139. 		ensure("5:isExactlyZero():Fail ", (TRUE == vec3D.isExactlyZero()));	
    140. 

  140. 	}
    141. 

  141. 

  142. 	template<> template<>
    143. 

  143. 	void v3dmath_object::test<7>()
    144. 

  144. 	{
    145. 

  145. 		F64 x = -2.32, y = 1.212, z = -.12;
    146. 

  146. 		LLVector3d vec3D(x,y,z);
    147. 

  147. 		
    148. 

  148. 		ensure("1:operator [] failed",( x ==  vec3D[0]));	
    149. 

  149. 		ensure("2:operator [] failed",( y ==  vec3D[1]));
    150. 

  150. 		ensure("3:operator [] failed",( z ==  vec3D[2]));
    151. 

  151. 		vec3D.clearVec();
    152. 

  152. 		x = 23.23, y = -.2361, z = 3.25;
    153. 

  153. 		vec3D.setVec(x,y,z);
    154. 

  154. 		F64 &ref1 = vec3D[0];
    155. 

  155. 		ensure("4:operator [] failed",( ref1 ==  vec3D[0]));
    156. 

  156. 		F64 &ref2 = vec3D[1];
    157. 

  157. 		ensure("5:operator [] failed",( ref2 ==  vec3D[1]));
    158. 

  158. 		F64 &ref3 = vec3D[2];
    159. 

  159. 		ensure("6:operator [] failed",( ref3 ==  vec3D[2]));
    160. 

  160. 	}
    161. 

  161. 

  162. 	template<> template<>
    163. 

  163. 	void v3dmath_object::test<8>()
    164. 

  164. 	{
    165. 

  165. 		F32 x = 1.f, y = 2.f, z = -1.f;
    166. 

  166. 		LLVector4 vec4(x,y,z);		
    167. 

  167. 		LLVector3d vec3D;
    168. 

  168. 		vec3D = vec4;
    169. 

  169. 		ensure("1:operator=:Fail to initialize ", ((vec4.mV[VX] == vec3D.mdV[VX]) && (vec4.mV[VY] == vec3D.mdV[VY]) && (vec4.mV[VZ] == vec3D.mdV[VZ])));
    170. 

  170. 	}
    171. 

  171. 

  172. 	template<> template<>
    173. 

  173. 	void v3dmath_object::test<9>()
    174. 

  174. 	{
    175. 

  175. 		F64 x1 = 1.78787878, y1 = 232322.2121, z1 = -12121.121212;
    176. 

  176. 		F64 x2 = 1.2, y2 = 2.5, z2 = 1.;
    177. 

  177. 		LLVector3d vec3D(x1,y1,z1),vec3Da(x2,y2,z2),vec3Db;
    178. 

  178. 		vec3Db = vec3Da+ vec3D;
    179. 

  179. 		ensure("1:operator+:Fail to initialize ", ((x1+x2 == vec3Db.mdV[VX]) && (y1+y2 == vec3Db.mdV[VY]) && (z1+z2 == vec3Db.mdV[VZ])));
    180. 

  180. 		x1 = -2.45, y1 = 2.1, z1 = 3.0;
    181. 

  181. 		vec3D.clearVec();
    182. 

  182. 		vec3Da.clearVec();
    183. 

  183. 		vec3D.setVec(x1,y1,z1);
    184. 

  184. 		vec3Da += vec3D;
    185. 

  185. 		ensure_equals("2:operator+=: Fail to initialize", vec3Da,vec3D);
    186. 

  186. 		vec3Da += vec3D;
    187. 

  187. 		ensure("3:operator+=:Fail to initialize ", ((2*x1 == vec3Da.mdV[VX]) && (2*y1 == vec3Da.mdV[VY]) && (2*z1 == vec3Da.mdV[VZ])));
    188. 

  188. 	}
    189. 

  189. 

  190. 	template<> template<>
    191. 

  191. 	void v3dmath_object::test<10>()
    192. 

  192. 	{
    193. 

  193. 		F64 x1 = 1., y1 = 2., z1 = -1.1;
    194. 

  194. 		F64 x2 = 1.2, y2 = 2.5, z2 = 1.;
    195. 

  195. 		LLVector3d vec3D(x1,y1,z1),vec3Da(x2,y2,z2),vec3Db;
    196. 

  196. 		vec3Db = vec3Da - vec3D;
    197. 

  197. 		ensure("1:operator-:Fail to initialize ", ((x2-x1 == vec3Db.mdV[VX]) && (y2-y1 == vec3Db.mdV[VY]) && (z2-z1 == vec3Db.mdV[VZ])));
    198. 

  198. 		x1 = -2.45, y1 = 2.1, z1 = 3.0;
    199. 

  199. 		vec3D.clearVec();
    200. 

  200. 		vec3Da.clearVec();
    201. 

  201. 		vec3D.setVec(x1,y1,z1);
    202. 

  202. 		vec3Da -=vec3D;
    203. 

  203. 		ensure("2:operator-=:Fail to initialize ", ((2.45 == vec3Da.mdV[VX]) && (-2.1 == vec3Da.mdV[VY]) && (-3.0 == vec3Da.mdV[VZ])));
    204. 

  204. 		vec3Da -= vec3D;
    205. 

  205. 		ensure("3:operator-=:Fail to initialize ", ((-2*x1 == vec3Da.mdV[VX]) && (-2*y1 == vec3Da.mdV[VY]) && (-2*z1 == vec3Da.mdV[VZ])));
    206. 

  206. 	}
    207. 

  207. 	template<> template<>
    208. 

  208. 	void v3dmath_object::test<11>()
    209. 

  209. 	{
    210. 

  210. 		F64 x1 = 1., y1 = 2., z1 = -1.1;
    211. 

  211. 		F64 x2 = 1.2, y2 = 2.5, z2 = 1.;
    212. 

  212. 		LLVector3d vec3D(x1,y1,z1),vec3Da(x2,y2,z2);
    213. 

  213. 		F64 res = vec3D * vec3Da;
    214. 

  214. 		ensure_approximately_equals(
    215. 

  215. 			"1:operator* failed",
    216. 

  216. 			res,
    217. 

  217. 			(x1*x2 + y1*y2 + z1*z2),
    218. 

  218. 			8);
    219. 

  219. 		vec3Da.clearVec();
    220. 

  220. 		F64 mulVal = 4.2;
    221. 

  221. 		vec3Da = vec3D * mulVal;
    222. 

  222. 		ensure_approximately_equals(
    223. 

  223. 			"2a:operator* failed",
    224. 

  224. 			vec3Da.mdV[VX],
    225. 

  225. 			x1*mulVal,
    226. 

  226. 			8);
    227. 

  227. 		ensure_approximately_equals(
    228. 

  228. 			"2b:operator* failed",
    229. 

  229. 			vec3Da.mdV[VY],
    230. 

  230. 			y1*mulVal,
    231. 

  231. 			8);
    232. 

  232. 		ensure_approximately_equals(
    233. 

  233. 			"2c:operator* failed",
    234. 

  234. 			vec3Da.mdV[VZ],
    235. 

  235. 			z1*mulVal,
    236. 

  236. 			8);
    237. 

  237. 		vec3Da.clearVec();
    238. 

  238. 		vec3Da = mulVal * vec3D;
    239. 

  239. 		ensure_approximately_equals(
    240. 

  240. 			"3a:operator* failed",
    241. 

  241. 			vec3Da.mdV[VX],
    242. 

  242. 			x1*mulVal,
    243. 

  243. 			8);
    244. 

  244. 		ensure_approximately_equals(
    245. 

  245. 			"3b:operator* failed",
    246. 

  246. 			vec3Da.mdV[VY],
    247. 

  247. 			y1*mulVal,
    248. 

  248. 			8);
    249. 

  249. 		ensure_approximately_equals(
    250. 

  250. 			"3c:operator* failed",
    251. 

  251. 			vec3Da.mdV[VZ],
    252. 

  252. 			z1*mulVal,
    253. 

  253. 			8);
    254. 

  254. 		vec3D *= mulVal;
    255. 

  255. 		ensure_approximately_equals(
    256. 

  256. 			"4a:operator*= failed",
    257. 

  257. 			vec3D.mdV[VX],
    258. 

  258. 			x1*mulVal,
    259. 

  259. 			8);
    260. 

  260. 		ensure_approximately_equals(
    261. 

  261. 			"4b:operator*= failed",
    262. 

  262. 			vec3D.mdV[VY],
    263. 

  263. 			y1*mulVal,
    264. 

  264. 			8);
    265. 

  265. 		ensure_approximately_equals(
    266. 

  266. 			"4c:operator*= failed",
    267. 

  267. 			vec3D.mdV[VZ],
    268. 

  268. 			z1*mulVal,
    269. 

  269. 			8);
    270. 

  270. 	}
    271. 

  271. 

  272. 	template<> template<>
    273. 

  273. 	void v3dmath_object::test<12>()
    274. 

  274. 	{
    275. 

  275. 		F64 x1 = 1., y1 = 2., z1 = -1.1;
    276. 

  276. 		F64 x2 = 1.2, y2 = 2.5, z2 = 1.;
    277. 

  277. 		F64 val1, val2, val3;
    278. 

  278. 		LLVector3d vec3D(x1,y1,z1),vec3Da(x2,y2,z2), vec3Db;
    279. 

  279. 		vec3Db = vec3D % vec3Da;
    280. 

  280. 		val1 = y1*z2 - y2*z1;
    281. 

  281. 		val2 = z1*x2 -z2*x1;
    282. 

  282. 		val3 = x1*y2-x2*y1;
    283. 

  283. 		ensure("1:operator% failed",(val1 == vec3Db.mdV[VX]) && (val2 == vec3Db.mdV[VY]) && (val3 == vec3Db.mdV[VZ])); 
    284. 

  284. 		vec3D %= vec3Da;
    285. 

  285. 		ensure("2:operator%= failed",
    286. 

  286. 		       is_approx_equal(vec3D.mdV[VX],vec3Db.mdV[VX]) &&
    287. 

  287. 		       is_approx_equal(vec3D.mdV[VY],vec3Db.mdV[VY]) &&
    288. 

  288. 		       is_approx_equal(vec3D.mdV[VZ],vec3Db.mdV[VZ]) ); 
    289. 

  289. 	}
    290. 

  290. 

  291. 	template<> template<>
    292. 

  292. 	void v3dmath_object::test<13>()
    293. 

  293. 	{
    294. 

  294. 		F64 x1 = 1., y1 = 2., z1 = -1.1,div = 4.2;
    295. 

  295. 		F64 t = 1.f / div;
    296. 

  296. 		LLVector3d vec3D(x1,y1,z1), vec3Da;
    297. 

  297. 		vec3Da = vec3D/div;
    298. 

  298. 		ensure_approximately_equals(
    299. 

  299. 			"1a:operator/ failed",
    300. 

  300. 			vec3Da.mdV[VX],
    301. 

  301. 			x1*t,
    302. 

  302. 			8);
    303. 

  303. 		ensure_approximately_equals(
    304. 

  304. 			"1b:operator/ failed",
    305. 

  305. 			vec3Da.mdV[VY],
    306. 

  306. 			y1*t,
    307. 

  307. 			8);
    308. 

  308. 		ensure_approximately_equals(
    309. 

  309. 			"1c:operator/ failed",
    310. 

  310. 			vec3Da.mdV[VZ],
    311. 

  311. 			z1*t,
    312. 

  312. 			8);
    313. 

  313. 		x1 = 1.23, y1 = 4., z1 = -2.32;
    314. 

  314. 		vec3D.clearVec();
    315. 

  315. 		vec3Da.clearVec();
    316. 

  316. 		vec3D.setVec(x1,y1,z1);
    317. 

  317. 		vec3Da = vec3D/div;
    318. 

  318. 		ensure_approximately_equals(
    319. 

  319. 			"2a:operator/ failed",
    320. 

  320. 			vec3Da.mdV[VX],
    321. 

  321. 			x1*t,
    322. 

  322. 			8);
    323. 

  323. 		ensure_approximately_equals(
    324. 

  324. 			"2b:operator/ failed",
    325. 

  325. 			vec3Da.mdV[VY],
    326. 

  326. 			y1*t,
    327. 

  327. 			8);
    328. 

  328. 		ensure_approximately_equals(
    329. 

  329. 			"2c:operator/ failed",
    330. 

  330. 			vec3Da.mdV[VZ],
    331. 

  331. 			z1*t,
    332. 

  332. 			8);
    333. 

  333. 		vec3D /= div;
    334. 

  334. 		ensure_approximately_equals(
    335. 

  335. 			"3a:operator/= failed",
    336. 

  336. 			vec3D.mdV[VX],
    337. 

  337. 			x1*t,
    338. 

  338. 			8);
    339. 

  339. 		ensure_approximately_equals(
    340. 

  340. 			"3b:operator/= failed",
    341. 

  341. 			vec3D.mdV[VY],
    342. 

  342. 			y1*t,
    343. 

  343. 			8);
    344. 

  344. 		ensure_approximately_equals(
    345. 

  345. 			"3c:operator/= failed",
    346. 

  346. 			vec3D.mdV[VZ],
    347. 

  347. 			z1*t,
    348. 

  348. 			8);
    349. 

  349. 	}
    350. 

  350. 

  351. 	template<> template<>
    352. 

  352. 	void v3dmath_object::test<14>()
    353. 

  353. 	{
    354. 

  354. 		F64 x1 = 1., y1 = 2., z1 = -1.1;
    355. 

  355. 		LLVector3d vec3D(x1,y1,z1), vec3Da;
    356. 

  356. 		ensure("1:operator!= failed",(TRUE == (vec3D !=vec3Da)));
    357. 

  357. 		vec3Da = vec3D;
    358. 

  358. 		ensure("2:operator== failed",(vec3D ==vec3Da)); 
    359. 

  359. 		vec3D.clearVec();
    360. 

  360. 		vec3Da.clearVec();
    361. 

  361. 		x1 = .211, y1 = 21.111, z1 = 23.22;
    362. 

  362. 		vec3D.setVec(x1,y1,z1);
    363. 

  363. 		vec3Da.setVec(x1,y1,z1);
    364. 

  364. 		ensure("3:operator== failed",(vec3D ==vec3Da)); 
    365. 

  365. 		ensure("4:operator!= failed",(FALSE == (vec3D !=vec3Da)));
    366. 

  366. 	}
    367. 

  367. 		
    368. 

  368. 	template<> template<>
    369. 

  369. 	void v3dmath_object::test<15>()
    370. 

  370. 	{
    371. 

  371. 		F64 x1 = 1., y1 = 2., z1 = -1.1;
    372. 

  372. 		LLVector3d vec3D(x1,y1,z1), vec3Da;
    373. 

  373. 		std::ostringstream stream1, stream2;
    374. 

  374. 		stream1 << vec3D;
    375. 

  375. 		vec3Da.setVec(x1,y1,z1);
    376. 

  376. 		stream2 << vec3Da;
    377. 

  377. 		ensure("1:operator << failed",(stream1.str() == stream2.str()));	
    378. 

  378. 	}
    379. 

  379. 

  380. 	template<> template<>
    381. 

  381. 	void v3dmath_object::test<16>()
    382. 

  382. 	{
    383. 

  383. 		F64 x1 = 1.23, y1 = 2.0, z1 = 4.;
    384. 

  384. 		std::string buf("1.23 2. 4");
    385. 

  385. 		LLVector3d vec3D, vec3Da(x1,y1,z1);
    386. 

  386. 		LLVector3d::parseVector3d(buf, &vec3D);
    387. 

  387. 		ensure_equals("1:parseVector3d: failed " , vec3D, vec3Da);	
    388. 

  388. 	}
    389. 

  389. 

  390. 	template<> template<>
    391. 

  391. 	void v3dmath_object::test<17>()
    392. 

  392. 	{
    393. 

  393. 		F64 x1 = 1., y1 = 2., z1 = -1.1;
    394. 

  394. 		LLVector3d vec3D(x1,y1,z1), vec3Da;
    395. 

  395. 		vec3Da = -vec3D;
    396. 

  396. 		ensure("1:operator- failed", (vec3D == - vec3Da));	
    397. 

  397. 	}
    398. 

  398. 

  399. 	template<> template<>
    400. 

  400. 	void v3dmath_object::test<18>()
    401. 

  401. 	{
    402. 

  402. 		F64 x = 1., y = 2., z = -1.1;
    403. 

  403. 		LLVector3d vec3D(x,y,z);
    404. 

  404. 		F64 res = (x*x + y*y + z*z) - vec3D.magVecSquared();
    405. 

  405. 		ensure("1:magVecSquared:Fail ", ((-F_APPROXIMATELY_ZERO <= res)&& (res <=F_APPROXIMATELY_ZERO)));
    406. 

  406. 		res = (F32) sqrt(x*x + y*y + z*z) - vec3D.magVec();
    407. 

  407. 		ensure("2:magVec: Fail ", ((-F_APPROXIMATELY_ZERO <= res)&& (res <=F_APPROXIMATELY_ZERO)));	
    408. 

  408. 	}
    409. 

  409. 

  410. 	template<> template<>
    411. 

  411. 	void v3dmath_object::test<19>()
    412. 

  412. 	{
    413. 

  413. 		F64 x = 1., y = 2., z = -1.1;
    414. 

  414. 		LLVector3d vec3D(x,y,z);
    415. 

  415. 		F64 mag = vec3D.normVec();
    416. 

  416. 		mag = 1.f/ mag;
    417. 

  417. 		ensure_approximately_equals(
    418. 

  418. 			"1a:normVec: Fail ",
    419. 

  419. 			vec3D.mdV[VX],
    420. 

  420. 			x * mag,
    421. 

  421. 			8);
    422. 

  422. 		ensure_approximately_equals(
    423. 

  423. 			"1b:normVec: Fail ",
    424. 

  424. 			vec3D.mdV[VY],
    425. 

  425. 			y * mag,
    426. 

  426. 			8);
    427. 

  427. 		ensure_approximately_equals(
    428. 

  428. 			"1c:normVec: Fail ",
    429. 

  429. 			vec3D.mdV[VZ],
    430. 

  430. 			z * mag,
    431. 

  431. 			8);
    432. 

  432. 		x = 0.000000001, y = 0.000000001, z = 0.000000001;
    433. 

  433. 		vec3D.clearVec();
    434. 

  434. 		vec3D.setVec(x,y,z);
    435. 

  435. 		mag = vec3D.normVec();
    436. 

  436. 		ensure_approximately_equals(
    437. 

  437. 			"2a:normVec: Fail ",
    438. 

  438. 			vec3D.mdV[VX],
    439. 

  439. 			x * mag,
    440. 

  440. 			8);
    441. 

  441. 		ensure_approximately_equals(
    442. 

  442. 			"2b:normVec: Fail ",
    443. 

  443. 			vec3D.mdV[VY],
    444. 

  444. 			y * mag,
    445. 

  445. 			8);
    446. 

  446. 		ensure_approximately_equals(
    447. 

  447. 			"2c:normVec: Fail ",
    448. 

  448. 			vec3D.mdV[VZ],
    449. 

  449. 			z * mag,
    450. 

  450. 			8);
    451. 

  451. 	}
    452. 

  452. 

  453. 	template<> template<>
    454. 

  454. 	void v3dmath_object::test<20>()
    455. 

  455. 	{
    456. 

  456. 		F64 x1 = 1111.232222;
    457. 

  457. 		F64 y1 = 2222222222.22;
    458. 

  458. 		F64 z1 = 422222222222.0;
    459. 

  459. 		std::string buf("1111.232222 2222222222.22 422222222222");
    460. 

  460. 		LLVector3d vec3Da, vec3Db(x1,y1,z1);
    461. 

  461. 		LLVector3d::parseVector3d(buf, &vec3Da);
    462. 

  462. 		ensure_equals("1:parseVector3 failed", vec3Da, vec3Db);	
    463. 

  463. 	}
    464. 

  464. 

  465. 	template<> template<>
    466. 

  466. 	void v3dmath_object::test<21>()
    467. 

  467. 	{
    468. 

  468. 		F64 x1 = 1., y1 = 2., z1 = -1.1;
    469. 

  469. 		F64 x2 = 1.2, y2 = 2.5, z2 = 1.;
    470. 

  470. 		F64 val = 2.3f,val1,val2,val3;
    471. 

  471. 		val1 = x1 + (x2 - x1)* val;
    472. 

  472. 		val2 = y1 + (y2 - y1)* val;
    473. 

  473. 		val3 = z1 + (z2 - z1)* val;
    474. 

  474. 		LLVector3d vec3Da(x1,y1,z1),vec3Db(x2,y2,z2);
    475. 

  475. 		LLVector3d vec3d = lerp(vec3Da,vec3Db,val);
    476. 

  476. 		ensure("1:lerp failed", ((val1 ==vec3d.mdV[VX])&& (val2 ==vec3d.mdV[VY]) && (val3 ==vec3d.mdV[VZ])));		
    477. 

  477. 	}
    478. 

  478. 

  479. 	template<> template<>
    480. 

  480. 	void v3dmath_object::test<22>()
    481. 

  481. 	{
    482. 

  482. 		F64 x = 2.32, y = 1.212, z = -.12;
    483. 

  483. 		F64 min = 0.0001, max = 3.0;
    484. 

  484. 		LLVector3d vec3d(x,y,z);		
    485. 

  485. 		ensure("1:clamp:Fail ", (TRUE == (vec3d.clamp(min, max))));
    486. 

  486. 		x = 0.000001f, z = 5.3f;
    487. 

  487. 		vec3d.setVec(x,y,z);
    488. 

  488. 		ensure("2:clamp:Fail ", (TRUE == (vec3d.clamp(min, max))));
    489. 

  489. 	}
    490. 

  490. 

  491. 	template<> template<>
    492. 

  492. 	void v3dmath_object::test<23>()
    493. 

  493. 	{
    494. 

  494. 		F64 x = 10., y = 20., z = -15.;
    495. 

  495. 		F64 epsilon = .23425;
    496. 

  496. 		LLVector3d vec3Da(x,y,z), vec3Db(x,y,z);
    497. 

  497. 		ensure("1:are_parallel: Fail ", (TRUE == are_parallel(vec3Da,vec3Db,epsilon)));
    498. 

  498. 		F64 x1 = -12., y1 = -20., z1 = -100.;
    499. 

  499. 		vec3Db.clearVec();
    500. 

  500. 		vec3Db.setVec(x1,y1,z1);
    501. 

  501. 		ensure("2:are_parallel: Fail ", (FALSE == are_parallel(vec3Da,vec3Db,epsilon)));
    502. 

  502. 	}
    503. 

  503. 

  504. 	template<> template<>
    505. 

  505. 	void v3dmath_object::test<24>()
    506. 

  506. 	{
    507. 

  507. #if LL_WINDOWS && _MSC_VER < 1400
    508. 

  508. 		skip("This fails on VS2003!");
    509. 

  509. #else
    510. 

  510. 		F64 x = 10., y = 20., z = -15.;
    511. 

  511. 		F64 angle1, angle2;
    512. 

  512. 		LLVector3d vec3Da(x,y,z), vec3Db(x,y,z);
    513. 

  513. 		angle1 = angle_between(vec3Da, vec3Db);
    514. 

  514. 		ensure("1:angle_between: Fail ", (0 == angle1));
    515. 

  515. 		F64 x1 = -1., y1 = -20., z1 = -1.;
    516. 

  516. 		vec3Da.clearVec();
    517. 

  517. 		vec3Da.setVec(x1,y1,z1);
    518. 

  518. 		angle2 = angle_between(vec3Da, vec3Db);
    519. 

  519. 		vec3Db.normVec();
    520. 

  520. 		vec3Da.normVec();
    521. 

  521. 		F64 angle = vec3Db*vec3Da;
    522. 

  522. 		angle = acos(angle);
    523. 

  523. 		ensure("2:angle_between: Fail ", (angle == angle2));
    524. 

  524. #endif
    525. 

  525. 	}
    526. 

  526. }
    527. 

  527.