orthanc-stone: UnitTestsSources/GeometryToolboxTests.cpp annotate

annotate UnitTestsSources/GeometryToolboxTests.cpp @ 1762:604fc328dc10

typo

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Mon, 10 May 2021 11:56:47 +0200
parents	9ac2a65d4172
children	177f93c8371c

rev	line source
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1 /**
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	2 * Stone of Orthanc
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	4 * Department, University Hospital of Liege, Belgium
1739 9ac2a65d4172 upgrade to year 2021 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1679 diff changeset	5 * Copyright (C) 2017-2021 Osimis S.A., Belgium
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	6 *
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	7 * This program is free software: you can redistribute it and/or
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	8 * modify it under the terms of the GNU Affero General Public License
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	9 * as published by the Free Software Foundation, either version 3 of
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	10 * the License, or (at your option) any later version.
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	11 *
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	12 * This program is distributed in the hope that it will be useful, but
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	13 * WITHOUT ANY WARRANTY; without even the implied warranty of
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	15 * Affero General Public License for more details.
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	16 *
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	17 * You should have received a copy of the GNU Affero General Public License
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	18 * along with this program. If not, see <http://www.gnu.org/licenses/>.
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	19 **/
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	20
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	21
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	22 #include <gtest/gtest.h>
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	23
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	24 #include "../OrthancStone/Sources/Toolbox/FiniteProjectiveCamera.h"
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	25 #include "../OrthancStone/Sources/Toolbox/GenericToolbox.h"
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	26 #include "../OrthancStone/Sources/Toolbox/GeometryToolbox.h"
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	27
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	28 #include <Logging.h>
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	29 #include <OrthancException.h>
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	30
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	31 #include <boost/math/special_functions/round.hpp>
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	32
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	33
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	34 TEST(GeometryToolbox, Interpolation)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	35 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	36 using namespace OrthancStone::GeometryToolbox;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	37
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	38 // https://en.wikipedia.org/wiki/Bilinear_interpolation#Application_in_image_processing
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	39 ASSERT_DOUBLE_EQ(146.1f, ComputeBilinearInterpolationUnitSquare(0.5f, 0.2f, 91, 210, 162, 95));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	40
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	41 ASSERT_DOUBLE_EQ(91, ComputeBilinearInterpolationUnitSquare(0, 0, 91, 210, 162, 95));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	42 ASSERT_DOUBLE_EQ(210, ComputeBilinearInterpolationUnitSquare(1, 0, 91, 210, 162, 95));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	43 ASSERT_DOUBLE_EQ(162, ComputeBilinearInterpolationUnitSquare(0, 1, 91, 210, 162, 95));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	44 ASSERT_DOUBLE_EQ(95, ComputeBilinearInterpolationUnitSquare(1, 1, 91, 210, 162, 95));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	45
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	46 ASSERT_DOUBLE_EQ(123.35f, ComputeTrilinearInterpolationUnitSquare
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	47 (0.5f, 0.2f, 0.7f,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	48 91, 210, 162, 95,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	49 51, 190, 80, 92));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	50
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	51 ASSERT_DOUBLE_EQ(ComputeBilinearInterpolationUnitSquare(0.5f, 0.2f, 91, 210, 162, 95),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	52 ComputeTrilinearInterpolationUnitSquare(0.5f, 0.2f, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	53 91, 210, 162, 95,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	54 51, 190, 80, 92));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	55
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	56 ASSERT_DOUBLE_EQ(ComputeBilinearInterpolationUnitSquare(0.5f, 0.2f, 51, 190, 80, 92),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	57 ComputeTrilinearInterpolationUnitSquare(0.5f, 0.2f, 1,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	58 91, 210, 162, 95,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	59 51, 190, 80, 92));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	60 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	61
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	62
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	63 static bool CompareMatrix(const OrthancStone::Matrix& a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	64 const OrthancStone::Matrix& b,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	65 double threshold = 0.00000001)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	66 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	67 if (a.size1() != b.size1() \|\|
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	68 a.size2() != b.size2())
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	69 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	70 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	71 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	72
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	73 for (size_t i = 0; i < a.size1(); i++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	74 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	75 for (size_t j = 0; j < a.size2(); j++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	76 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	77 if (fabs(a(i, j) - b(i, j)) > threshold)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	78 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	79 LOG(ERROR) << "Too large difference in component ("
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	80 << i << "," << j << "): " << a(i,j) << " != " << b(i,j);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	81 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	82 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	83 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	84 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	85
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	86 return true;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	87 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	88
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	89
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	90 static bool CompareVector(const OrthancStone::Vector& a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	91 const OrthancStone::Vector& b,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	92 double threshold = 0.00000001)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	93 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	94 if (a.size() != b.size())
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	95 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	96 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	97 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	98
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	99 for (size_t i = 0; i < a.size(); i++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	100 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	101 if (fabs(a(i) - b(i)) > threshold)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	102 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	103 LOG(ERROR) << "Too large difference in component "
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	104 << i << ": " << a(i) << " != " << b(i);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	105 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	106 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	107 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	108
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	109 return true;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	110 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	111
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	112
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	113
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	114 TEST(FiniteProjectiveCamera, Decomposition1)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	115 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	116 // Example 6.2 of "Multiple View Geometry in Computer Vision - 2nd
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	117 // edition" (page 163)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	118 const double p[12] = {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	119 3.53553e+2, 3.39645e+2, 2.77744e+2, -1.44946e+6,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	120 -1.03528e+2, 2.33212e+1, 4.59607e+2, -6.32525e+5,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	121 7.07107e-1, -3.53553e-1, 6.12372e-1, -9.18559e+2
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	122 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	123
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	124 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	125 ASSERT_EQ(3u, camera.GetMatrix().size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	126 ASSERT_EQ(4u, camera.GetMatrix().size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	127 ASSERT_EQ(3u, camera.GetIntrinsicParameters().size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	128 ASSERT_EQ(3u, camera.GetIntrinsicParameters().size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	129 ASSERT_EQ(3u, camera.GetRotation().size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	130 ASSERT_EQ(3u, camera.GetRotation().size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	131 ASSERT_EQ(3u, camera.GetCenter().size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	132
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	133 ASSERT_NEAR(1000.0, camera.GetCenter()[0], 0.01);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	134 ASSERT_NEAR(2000.0, camera.GetCenter()[1], 0.01);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	135 ASSERT_NEAR(1500.0, camera.GetCenter()[2], 0.01);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	136
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	137 ASSERT_NEAR(468.2, camera.GetIntrinsicParameters() (0, 0), 0.1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	138 ASSERT_NEAR(91.2, camera.GetIntrinsicParameters() (0, 1), 0.1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	139 ASSERT_NEAR(300.0, camera.GetIntrinsicParameters() (0, 2), 0.1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	140 ASSERT_NEAR(427.2, camera.GetIntrinsicParameters() (1, 1), 0.1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	141 ASSERT_NEAR(200.0, camera.GetIntrinsicParameters() (1, 2), 0.1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	142 ASSERT_NEAR(1.0, camera.GetIntrinsicParameters() (2, 2), 0.1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	143
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	144 ASSERT_NEAR(0, camera.GetIntrinsicParameters() (1, 0), 0.0000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	145 ASSERT_NEAR(0, camera.GetIntrinsicParameters() (2, 0), 0.0000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	146 ASSERT_NEAR(0, camera.GetIntrinsicParameters() (2, 1), 0.0000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	147
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	148 ASSERT_NEAR(0.41380, camera.GetRotation() (0, 0), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	149 ASSERT_NEAR(0.90915, camera.GetRotation() (0, 1), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	150 ASSERT_NEAR(0.04708, camera.GetRotation() (0, 2), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	151 ASSERT_NEAR(-0.57338, camera.GetRotation() (1, 0), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	152 ASSERT_NEAR(0.22011, camera.GetRotation() (1, 1), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	153 ASSERT_NEAR(0.78917, camera.GetRotation() (1, 2), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	154 ASSERT_NEAR(0.70711, camera.GetRotation() (2, 0), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	155 ASSERT_NEAR(-0.35355, camera.GetRotation() (2, 1), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	156 ASSERT_NEAR(0.61237, camera.GetRotation() (2, 2), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	157
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	158 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	159
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	160 OrthancStone::FiniteProjectiveCamera camera2(camera.GetIntrinsicParameters(),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	161 camera.GetRotation(),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	162 camera.GetCenter());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	163
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	164 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	165 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	166 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	167 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	168 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	169
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	170
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	171 TEST(FiniteProjectiveCamera, Decomposition2)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	172 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	173 const double p[] = { 1188.111986, 580.205341, -808.445330, 128000.000000, -366.466264, 1446.510501, 418.499736, 128000.000000, -0.487118, 0.291726, -0.823172, 500.000000 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	174 const double k[] = { -1528.494743, 0.000000, 256.000000, 0.000000, 1528.494743, 256.000000, 0.000000, 0.000000, 1.000000 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	175 const double r[] = { -0.858893, -0.330733, 0.391047, -0.158171, 0.897503, 0.411668, -0.487118, 0.291726, -0.823172 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	176 const double c[] = { 243.558936, -145.863085, 411.585964 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	177
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	178 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	179 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	180
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	181 OrthancStone::FiniteProjectiveCamera camera2(k, r, c);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	182 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix(), 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	183 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters(), 0.001));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	184 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation(), 0.000001));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	185 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter(), 0.0001));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	186 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	187
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	188
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	189 TEST(FiniteProjectiveCamera, Decomposition3)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	190 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	191 const double p[] = { 10, 0, 0, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	192 0, 20, 0, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	193 0, 0, 30, 0 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	194
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	195 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	196 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	197 ASSERT_DOUBLE_EQ(10, camera.GetIntrinsicParameters() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	198 ASSERT_DOUBLE_EQ(20, camera.GetIntrinsicParameters() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	199 ASSERT_DOUBLE_EQ(30, camera.GetIntrinsicParameters() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	200 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	201 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	202 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	203 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	204 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	205 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	206 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	207
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	208
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	209 TEST(FiniteProjectiveCamera, Decomposition4)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	210 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	211 const double p[] = { 1, 0, 0, 10,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	212 0, 1, 0, 20,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	213 0, 0, 1, 30 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	214
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	215 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	216 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	217 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	218 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	219 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	220 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	221 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	222 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	223 ASSERT_DOUBLE_EQ(-10, camera.GetCenter() (0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	224 ASSERT_DOUBLE_EQ(-20, camera.GetCenter() (1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	225 ASSERT_DOUBLE_EQ(-30, camera.GetCenter() (2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	226 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	227
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	228
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	229 TEST(FiniteProjectiveCamera, Decomposition5)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	230 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	231 const double p[] = { 0, 0, 10, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	232 0, 20, 0, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	233 30, 0, 0, 0 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	234
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	235 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	236 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	237 ASSERT_DOUBLE_EQ(-10, camera.GetIntrinsicParameters() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	238 ASSERT_DOUBLE_EQ(20, camera.GetIntrinsicParameters() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	239 ASSERT_DOUBLE_EQ(30, camera.GetIntrinsicParameters() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	240 ASSERT_DOUBLE_EQ(-1, camera.GetRotation() (0, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	241 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	242 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	243 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	244 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	245 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	246
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	247 OrthancStone::FiniteProjectiveCamera camera2(camera.GetIntrinsicParameters(),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	248 camera.GetRotation(),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	249 camera.GetCenter());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	250 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	251 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	252 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	253 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	254 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	255
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	256
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	257 static double GetCosAngle(const OrthancStone::Vector& a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	258 const OrthancStone::Vector& b)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	259 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	260 // Returns the cosine of the angle between two vectors
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	261 // https://en.wikipedia.org/wiki/Dot_product#Geometric_definition
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	262 return boost::numeric::ublas::inner_prod(a, b) /
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	263 (boost::numeric::ublas::norm_2(a) * boost::numeric::ublas::norm_2(b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	264 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	265
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	266
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	267 TEST(FiniteProjectiveCamera, Ray)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	268 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	269 const double pp[] = { -1499.650894, 2954.618773, -259.737419, 637891.819097,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	270 -2951.517707, -1501.019129, -285.785281, 637891.819097,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	271 0.008528, 0.003067, -0.999959, 2491.764918 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	272
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	273 const OrthancStone::FiniteProjectiveCamera camera(pp);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	274
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	275 ASSERT_NEAR(-21.2492, camera.GetCenter() (0), 0.0001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	276 ASSERT_NEAR(-7.64234, camera.GetCenter() (1), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	277 ASSERT_NEAR(2491.66, camera.GetCenter() (2), 0.01);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	278
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	279 // Image plane that led to these parameters, with principal point at
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	280 // (256,256). The image has dimensions 512x512.
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	281 OrthancStone::Vector o =
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	282 OrthancStone::LinearAlgebra::CreateVector(7.009620, 2.521030, -821.942000);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	283 OrthancStone::Vector ax =
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	284 OrthancStone::LinearAlgebra::CreateVector(-0.453219, 0.891399, -0.001131);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	285 OrthancStone::Vector ay =
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	286 OrthancStone::LinearAlgebra::CreateVector(-0.891359, -0.453210, -0.008992);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	287
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	288 OrthancStone::CoordinateSystem3D imagePlane(o, ax, ay);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	289
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	290 // Back-projection of the principal point
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	291 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	292 OrthancStone::Vector ray = camera.GetRayDirection(256, 256);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	293
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	294 // The principal axis vector is orthogonal to the image plane
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	295 // (i.e. parallel to the plane normal), in the opposite direction
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	296 // ("-1" corresponds to "cos(pi)").
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	297 ASSERT_NEAR(-1, GetCosAngle(ray, imagePlane.GetNormal()), 0.0000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	298
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	299 // Forward projection of principal axis, resulting in the principal point
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	300 double x, y;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	301 camera.ApplyFinite(x, y, camera.GetCenter() - ray);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	302
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	303 ASSERT_NEAR(256, x, 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	304 ASSERT_NEAR(256, y, 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	305 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	306
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	307 // Back-projection of the 4 corners of the image
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	308 std::vector<double> cx, cy;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	309 cx.push_back(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	310 cy.push_back(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	311 cx.push_back(512);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	312 cy.push_back(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	313 cx.push_back(512);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	314 cy.push_back(512);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	315 cx.push_back(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	316 cy.push_back(512);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	317
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	318 bool first = true;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	319 double angle;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	320
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	321 for (size_t i = 0; i < cx.size(); i++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	322 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	323 OrthancStone::Vector ray = camera.GetRayDirection(cx[i], cy[i]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	324
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	325 // Check that the angle wrt. principal axis is the same for all
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	326 // the 4 corners
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	327 double a = GetCosAngle(ray, imagePlane.GetNormal());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	328 if (first)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	329 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	330 first = false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	331 angle = a;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	332 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	333 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	334 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	335 ASSERT_NEAR(angle, a, 0.000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	336 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	337
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	338 // Forward projection of the ray, going back to the original point
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	339 double x, y;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	340 camera.ApplyFinite(x, y, camera.GetCenter() - ray);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	341
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	342 ASSERT_NEAR(cx[i], x, 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	343 ASSERT_NEAR(cy[i], y, 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	344
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	345 // Alternative construction, by computing the intersection of the
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	346 // ray with the image plane
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	347 OrthancStone::Vector p;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	348 ASSERT_TRUE(imagePlane.IntersectLine(p, camera.GetCenter(), -ray));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	349 imagePlane.ProjectPoint(x, y, p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	350 ASSERT_NEAR(cx[i], x + 256, 0.01);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	351 ASSERT_NEAR(cy[i], y + 256, 0.01);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	352 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	353 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	354
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	355
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	356 TEST(Matrix, Inverse1)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	357 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	358 OrthancStone::Matrix a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	359
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	360 a.resize(0, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	361 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	362 ASSERT_EQ(0u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	363 ASSERT_EQ(0u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	364
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	365 a.resize(2, 3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	366 ASSERT_THROW(OrthancStone::LinearAlgebra::InvertMatrix(b, a), Orthanc::OrthancException);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	367
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	368 a.resize(1, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	369 a(0, 0) = 45.0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	370
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	371 ASSERT_DOUBLE_EQ(45, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	372 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	373 ASSERT_EQ(1u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	374 ASSERT_EQ(1u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	375 ASSERT_DOUBLE_EQ(1.0 / 45.0, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	376
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	377 a(0, 0) = 0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	378 ASSERT_DOUBLE_EQ(0, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	379 ASSERT_THROW(OrthancStone::LinearAlgebra::InvertMatrix(b, a), Orthanc::OrthancException);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	380 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	381
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	382
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	383 TEST(Matrix, Inverse2)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	384 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	385 OrthancStone::Matrix a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	386 a.resize(2, 2);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	387 a(0, 0) = 4;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	388 a(0, 1) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	389 a(1, 0) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	390 a(1, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	391
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	392 ASSERT_DOUBLE_EQ(-1, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	393 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	394 ASSERT_EQ(2u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	395 ASSERT_EQ(2u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	396
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	397 ASSERT_DOUBLE_EQ(-2, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	398 ASSERT_DOUBLE_EQ(3, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	399 ASSERT_DOUBLE_EQ(3, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	400 ASSERT_DOUBLE_EQ(-4, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	401
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	402 a(0, 0) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	403 a(0, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	404 a(1, 0) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	405 a(1, 1) = 4;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	406
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	407 ASSERT_DOUBLE_EQ(-2, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	408 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	409
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	410 ASSERT_DOUBLE_EQ(-2, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	411 ASSERT_DOUBLE_EQ(1, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	412 ASSERT_DOUBLE_EQ(1.5, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	413 ASSERT_DOUBLE_EQ(-0.5, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	414 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	415
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	416
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	417 TEST(Matrix, Inverse3)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	418 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	419 OrthancStone::Matrix a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	420 a.resize(3, 3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	421 a(0, 0) = 7;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	422 a(0, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	423 a(0, 2) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	424 a(1, 0) = 0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	425 a(1, 1) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	426 a(1, 2) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	427 a(2, 0) = -3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	428 a(2, 1) = 4;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	429 a(2, 2) = -2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	430
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	431 ASSERT_DOUBLE_EQ(1, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	432 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	433 ASSERT_EQ(3u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	434 ASSERT_EQ(3u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	435
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	436 ASSERT_DOUBLE_EQ(-2, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	437 ASSERT_DOUBLE_EQ(8, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	438 ASSERT_DOUBLE_EQ(-5, b(0, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	439 ASSERT_DOUBLE_EQ(3, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	440 ASSERT_DOUBLE_EQ(-11, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	441 ASSERT_DOUBLE_EQ(7, b(1, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	442 ASSERT_DOUBLE_EQ(9, b(2, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	443 ASSERT_DOUBLE_EQ(-34, b(2, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	444 ASSERT_DOUBLE_EQ(21, b(2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	445
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	446
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	447 a(0, 0) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	448 a(0, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	449 a(0, 2) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	450 a(1, 0) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	451 a(1, 1) = 0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	452 a(1, 2) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	453 a(2, 0) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	454 a(2, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	455 a(2, 2) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	456
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	457 ASSERT_DOUBLE_EQ(2, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	458 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	459 ASSERT_EQ(3u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	460 ASSERT_EQ(3u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	461
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	462 ASSERT_DOUBLE_EQ(-1, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	463 ASSERT_DOUBLE_EQ(1, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	464 ASSERT_DOUBLE_EQ(1, b(0, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	465 ASSERT_DOUBLE_EQ(0, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	466 ASSERT_DOUBLE_EQ(-0.5, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	467 ASSERT_DOUBLE_EQ(0.5, b(1, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	468 ASSERT_DOUBLE_EQ(1, b(2, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	469 ASSERT_DOUBLE_EQ(0, b(2, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	470 ASSERT_DOUBLE_EQ(-1, b(2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	471 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	472
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	473
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	474 TEST(Matrix, Inverse4)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	475 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	476 OrthancStone::Matrix a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	477 a.resize(4, 4);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	478 a(0, 0) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	479 a(0, 1) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	480 a(0, 2) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	481 a(0, 3) = -3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	482 a(1, 0) = -2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	483 a(1, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	484 a(1, 2) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	485 a(1, 3) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	486 a(2, 0) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	487 a(2, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	488 a(2, 2) = -3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	489 a(2, 3) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	490 a(3, 0) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	491 a(3, 1) = -2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	492 a(3, 2) = -3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	493 a(3, 3) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	494
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	495 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	496 ASSERT_EQ(4u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	497 ASSERT_EQ(4u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	498
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	499 b *= 134.0; // This is the determinant
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	500
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	501 ASSERT_DOUBLE_EQ(8, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	502 ASSERT_DOUBLE_EQ(-44, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	503 ASSERT_DOUBLE_EQ(30, b(0, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	504 ASSERT_DOUBLE_EQ(-10, b(0, 3));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	505 ASSERT_DOUBLE_EQ(2, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	506 ASSERT_DOUBLE_EQ(-11, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	507 ASSERT_DOUBLE_EQ(41, b(1, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	508 ASSERT_DOUBLE_EQ(-36, b(1, 3));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	509 ASSERT_DOUBLE_EQ(16, b(2, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	510 ASSERT_DOUBLE_EQ(-21, b(2, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	511 ASSERT_DOUBLE_EQ(-7, b(2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	512 ASSERT_DOUBLE_EQ(-20, b(2, 3));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	513 ASSERT_DOUBLE_EQ(-28, b(3, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	514 ASSERT_DOUBLE_EQ(-47, b(3, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	515 ASSERT_DOUBLE_EQ(29, b(3, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	516 ASSERT_DOUBLE_EQ(-32, b(3, 3));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	517 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	518
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	519
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	520 TEST(FiniteProjectiveCamera, Calibration)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	521 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	522 unsigned int volumeWidth = 512;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	523 unsigned int volumeHeight = 512;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	524 unsigned int volumeDepth = 110;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	525
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	526 OrthancStone::Vector camera = OrthancStone::LinearAlgebra::CreateVector
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	527 (-1000, -5000, -static_cast<double>(volumeDepth) * 32);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	528
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	529 OrthancStone::Vector principalPoint = OrthancStone::LinearAlgebra::CreateVector
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	530 (volumeWidth/2, volumeHeight/2, volumeDepth * 2);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	531
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	532 OrthancStone::FiniteProjectiveCamera c(camera, principalPoint, 0, 512, 512, 1, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	533
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	534 double swapv[9] = { 1, 0, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	535 0, -1, 512,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	536 0, 0, 1 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	537 OrthancStone::Matrix swap;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	538 OrthancStone::LinearAlgebra::FillMatrix(swap, 3, 3, swapv);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	539
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	540 OrthancStone::Matrix p = OrthancStone::LinearAlgebra::Product(swap, c.GetMatrix());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	541 p /= p(2,3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	542
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	543 ASSERT_NEAR( 1.04437, p(0,0), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	544 ASSERT_NEAR(-0.0703111, p(0,1), 0.00000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	545 ASSERT_NEAR(-0.179283, p(0,2), 0.000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	546 ASSERT_NEAR( 61.7431, p(0,3), 0.0001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	547 ASSERT_NEAR( 0.11127, p(1,0), 0.000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	548 ASSERT_NEAR(-0.595541, p(1,1), 0.000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	549 ASSERT_NEAR( 0.872211, p(1,2), 0.000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	550 ASSERT_NEAR( 203.748, p(1,3), 0.001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	551 ASSERT_NEAR( 3.08593e-05, p(2,0), 0.0000000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	552 ASSERT_NEAR( 0.000129138, p(2,1), 0.000000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	553 ASSERT_NEAR( 9.18901e-05, p(2,2), 0.0000000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	554 ASSERT_NEAR( 1, p(2,3), 0.0000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	555 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	556
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	557
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	558 static bool IsEqualRotationVector(OrthancStone::Vector a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	559 OrthancStone::Vector b)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	560 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	561 if (a.size() != b.size() \|\|
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	562 a.size() != 3)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	563 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	564 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	565 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	566 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	567 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	568 OrthancStone::LinearAlgebra::NormalizeVector(a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	569 OrthancStone::LinearAlgebra::NormalizeVector(b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	570 return OrthancStone::LinearAlgebra::IsCloseToZero(boost::numeric::ublas::norm_2(a - b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	571 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	572 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	573
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	574
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	575 TEST(GeometryToolbox, AlignVectorsWithRotation)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	576 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	577 OrthancStone::Vector a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	578 OrthancStone::Matrix r;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	579
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	580 OrthancStone::LinearAlgebra::AssignVector(a, -200, 200, -846.63);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	581 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	582
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	583 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	584 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	585 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, a), b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	586
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	587 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	588 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	589 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, b), a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	590
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	591 OrthancStone::LinearAlgebra::AssignVector(a, 1, 0, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	592 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	593 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	594 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	595 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, a), b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	596
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	597 OrthancStone::LinearAlgebra::AssignVector(a, 0, 1, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	598 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	599 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	600 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	601 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, a), b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	602
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	603 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	604 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	605 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	606 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	607 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, a), b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	608
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	609 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	610 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	611 ASSERT_THROW(OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b), Orthanc::OrthancException);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	612
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	613 // TODO: Deal with opposite vectors
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	614
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	615 /*
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	616 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, -1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	617 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	618 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	619 OrthancStone::LinearAlgebra::Print(r);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	620 OrthancStone::LinearAlgebra::Print(boost::numeric::ublas::prod(r, a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	621 */
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	622 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	623
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	624
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	625 static bool IsEqualVectorL1(OrthancStone::Vector a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	626 OrthancStone::Vector b)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	627 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	628 if (a.size() != b.size())
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	629 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	630 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	631 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	632 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	633 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	634 for (size_t i = 0; i < a.size(); i++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	635 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	636 if (!OrthancStone::LinearAlgebra::IsNear(a[i], b[i], 0.0001))
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	637 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	638 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	639 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	640 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	641
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	642 return true;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	643 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	644 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	645
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	646
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	647 TEST(VolumeImageGeometry, Basic)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	648 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	649 using namespace OrthancStone;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	650
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	651 VolumeImageGeometry g;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	652 g.SetSizeInVoxels(10, 20, 30);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	653 g.SetVoxelDimensions(1, 2, 3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	654
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	655 Vector p = g.GetCoordinates(0, 0, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	656 ASSERT_EQ(3u, p.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	657 ASSERT_DOUBLE_EQ(-1.0 / 2.0, p[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	658 ASSERT_DOUBLE_EQ(-2.0 / 2.0, p[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	659 ASSERT_DOUBLE_EQ(-3.0 / 2.0, p[2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	660
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	661 p = g.GetCoordinates(1, 1, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	662 ASSERT_DOUBLE_EQ(-1.0 / 2.0 + 10.0 * 1.0, p[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	663 ASSERT_DOUBLE_EQ(-2.0 / 2.0 + 20.0 * 2.0, p[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	664 ASSERT_DOUBLE_EQ(-3.0 / 2.0 + 30.0 * 3.0, p[2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	665
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	666 VolumeProjection proj;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	667 ASSERT_TRUE(g.DetectProjection(proj, g.GetAxialGeometry().GetNormal()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	668 ASSERT_EQ(VolumeProjection_Axial, proj);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	669 ASSERT_TRUE(g.DetectProjection(proj, g.GetCoronalGeometry().GetNormal()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	670 ASSERT_EQ(VolumeProjection_Coronal, proj);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	671 ASSERT_TRUE(g.DetectProjection(proj, g.GetSagittalGeometry().GetNormal()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	672 ASSERT_EQ(VolumeProjection_Sagittal, proj);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	673
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	674 ASSERT_EQ(10u, g.GetProjectionWidth(VolumeProjection_Axial));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	675 ASSERT_EQ(20u, g.GetProjectionHeight(VolumeProjection_Axial));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	676 ASSERT_EQ(30u, g.GetProjectionDepth(VolumeProjection_Axial));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	677 ASSERT_EQ(10u, g.GetProjectionWidth(VolumeProjection_Coronal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	678 ASSERT_EQ(30u, g.GetProjectionHeight(VolumeProjection_Coronal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	679 ASSERT_EQ(20u, g.GetProjectionDepth(VolumeProjection_Coronal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	680 ASSERT_EQ(20u, g.GetProjectionWidth(VolumeProjection_Sagittal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	681 ASSERT_EQ(30u, g.GetProjectionHeight(VolumeProjection_Sagittal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	682 ASSERT_EQ(10u, g.GetProjectionDepth(VolumeProjection_Sagittal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	683
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	684 p = g.GetVoxelDimensions(VolumeProjection_Axial);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	685 ASSERT_EQ(3u, p.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	686 ASSERT_DOUBLE_EQ(1, p[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	687 ASSERT_DOUBLE_EQ(2, p[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	688 ASSERT_DOUBLE_EQ(3, p[2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	689 p = g.GetVoxelDimensions(VolumeProjection_Coronal);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	690 ASSERT_EQ(3u, p.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	691 ASSERT_DOUBLE_EQ(1, p[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	692 ASSERT_DOUBLE_EQ(3, p[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	693 ASSERT_DOUBLE_EQ(2, p[2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	694 p = g.GetVoxelDimensions(VolumeProjection_Sagittal);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	695 ASSERT_EQ(3u, p.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	696 ASSERT_DOUBLE_EQ(2, p[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	697 ASSERT_DOUBLE_EQ(3, p[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	698 ASSERT_DOUBLE_EQ(1, p[2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	699
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	700 // Loop over all the voxels of the volume
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	701 for (unsigned int z = 0; z < g.GetDepth(); z++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	702 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	703 const float zz = (0.5f + static_cast<float>(z)) / static_cast<float>(g.GetDepth()); // Z-center of the voxel
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	704
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	705 for (unsigned int y = 0; y < g.GetHeight(); y++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	706 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	707 const float yy = (0.5f + static_cast<float>(y)) / static_cast<float>(g.GetHeight()); // Y-center of the voxel
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	708
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	709 for (unsigned int x = 0; x < g.GetWidth(); x++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	710 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	711 const float xx = (0.5f + static_cast<float>(x)) / static_cast<float>(g.GetWidth()); // X-center of the voxel
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	712
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	713 const float sx = 1.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	714 const float sy = 2.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	715 const float sz = 3.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	716
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	717 Vector p = g.GetCoordinates(xx, yy, zz);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	718
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	719 Vector q = (g.GetAxialGeometry().MapSliceToWorldCoordinates(
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	720 static_cast<double>(x) * sx,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	721 static_cast<double>(y) * sy) +
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	722 z * sz * g.GetAxialGeometry().GetNormal());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	723 ASSERT_TRUE(IsEqualVectorL1(p, q));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	724
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	725 q = (g.GetCoronalGeometry().MapSliceToWorldCoordinates(
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	726 static_cast<double>(x) * sx,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	727 static_cast<double>(g.GetDepth() - 1 - z) * sz) +
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	728 y * sy * g.GetCoronalGeometry().GetNormal());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	729 ASSERT_TRUE(IsEqualVectorL1(p, q));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	730
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	731 /**
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	732 * WARNING: In sagittal geometry, the normal points to
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	733 * REDUCING X-axis in the 3D world. This is necessary to keep
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	734 * the right-hand coordinate system. Hence the "-".
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	735 **/
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	736 q = (g.GetSagittalGeometry().MapSliceToWorldCoordinates(
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	737 static_cast<double>(y) * sy,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	738 static_cast<double>(g.GetDepth() - 1 - z) * sz) +
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	739 x * sx * (-g.GetSagittalGeometry().GetNormal()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	740 ASSERT_TRUE(IsEqualVectorL1(p, q));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	741 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	742 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	743 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	744
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	745 ASSERT_EQ(0, (int) VolumeProjection_Axial);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	746 ASSERT_EQ(1, (int) VolumeProjection_Coronal);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	747 ASSERT_EQ(2, (int) VolumeProjection_Sagittal);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	748
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	749 for (int p = 0; p < 3; p++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	750 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	751 VolumeProjection projection = (VolumeProjection) p;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	752 const CoordinateSystem3D& s = g.GetProjectionGeometry(projection);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	753
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	754 ASSERT_THROW(g.GetProjectionSlice(projection, g.GetProjectionDepth(projection)), Orthanc::OrthancException);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	755
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	756 for (unsigned int i = 0; i < g.GetProjectionDepth(projection); i++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	757 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	758 CoordinateSystem3D plane = g.GetProjectionSlice(projection, i);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	759
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	760 if (projection == VolumeProjection_Sagittal)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	761 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	762 ASSERT_TRUE(IsEqualVectorL1(plane.GetOrigin(), s.GetOrigin() + static_cast<double>(i) *
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	763 (-s.GetNormal()) * g.GetVoxelDimensions(projection)[2]));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	764 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	765 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	766 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	767 ASSERT_TRUE(IsEqualVectorL1(plane.GetOrigin(), s.GetOrigin() + static_cast<double>(i) *
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	768 s.GetNormal() * g.GetVoxelDimensions(projection)[2]));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	769 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	770
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	771 ASSERT_TRUE(IsEqualVectorL1(plane.GetAxisX(), s.GetAxisX()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	772 ASSERT_TRUE(IsEqualVectorL1(plane.GetAxisY(), s.GetAxisY()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	773
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	774 unsigned int slice;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	775 VolumeProjection q;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	776 ASSERT_TRUE(g.DetectSlice(q, slice, plane));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	777 ASSERT_EQ(projection, q);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	778 ASSERT_EQ(i, slice);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	779 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	780 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	781 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	782
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	783
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	784 TEST(LinearAlgebra, ParseVectorLocale)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	785 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	786 OrthancStone::Vector v;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	787
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	788 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "1.2"));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	789 ASSERT_EQ(1u, v.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	790 ASSERT_DOUBLE_EQ(1.2, v[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	791
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	792 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "-1.2e+2"));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	793 ASSERT_EQ(1u, v.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	794 ASSERT_DOUBLE_EQ(-120.0, v[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	795
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	796 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "-1e-2\\2"));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	797 ASSERT_EQ(2u, v.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	798 ASSERT_DOUBLE_EQ(-0.01, v[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	799 ASSERT_DOUBLE_EQ(2.0, v[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	800
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	801 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "1.3671875\\1.3671875"));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	802 ASSERT_EQ(2u, v.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	803 ASSERT_DOUBLE_EQ(1.3671875, v[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	804 ASSERT_DOUBLE_EQ(1.3671875, v[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	805 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	806
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	807 TEST(GenericToolbox, NormalizeUuid)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	808 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	809 std::string ref("44ca5051-14ef-4d2f-8bd7-db20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	810
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	811 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	812 std::string u("44ca5051-14ef-4d2f-8bd7-db20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	813 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	814 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	815 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	816
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	817 // space left
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	818 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	819 std::string u(" 44ca5051-14ef-4d2f-8bd7-db20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	820 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	821 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	822 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	823
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	824 // space right
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	825 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	826 std::string u("44ca5051-14ef-4d2f-8bd7-db20bfb61fbb ");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	827 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	828 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	829 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	830
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	831 // space l & r
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	832 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	833 std::string u(" 44ca5051-14ef-4d2f-8bd7-db20bfb61fbb ");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	834 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	835 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	836 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	837
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	838 // space left + case
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	839 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	840 std::string u(" 44CA5051-14ef-4d2f-8bd7-dB20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	841 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	842 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	843 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	844
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	845 // space right + case
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	846 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	847 std::string u("44ca5051-14EF-4D2f-8bd7-db20bfb61fbB ");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	848 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	849 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	850 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	851
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	852 // space l & r + case
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	853 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	854 std::string u(" 44cA5051-14Ef-4d2f-8bD7-db20bfb61fbb ");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	855 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	856 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	857 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	858
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	859 // no
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	860 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	861 std::string u(" 44ca5051-14ef-4d2f-8bd7- db20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	862 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	863 ASSERT_NE(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	864 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	865
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	866 // no
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	867 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	868 std::string u("44ca5051-14ef-4d2f-8bd7-db20bfb61fb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	869 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	870 ASSERT_NE(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	871 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	872 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	873
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	874
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	875 TEST(CoordinateSystem3D, Basic)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	876 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	877 using namespace OrthancStone;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	878
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	879 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	880 CoordinateSystem3D c;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	881 ASSERT_FALSE(c.IsValid());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	882 ASSERT_DOUBLE_EQ(c.GetNormal()[0], 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	883 ASSERT_DOUBLE_EQ(c.GetNormal()[1], 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	884 ASSERT_DOUBLE_EQ(c.GetNormal()[2], 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	885
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	886 ASSERT_DOUBLE_EQ(0, c.ComputeDistance(LinearAlgebra::CreateVector(0, 0, 0)));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	887 ASSERT_DOUBLE_EQ(0, c.ComputeDistance(LinearAlgebra::CreateVector(5, 0, 0)));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	888 ASSERT_DOUBLE_EQ(0, c.ComputeDistance(LinearAlgebra::CreateVector(0, 5, 0)));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	889 ASSERT_DOUBLE_EQ(5, c.ComputeDistance(LinearAlgebra::CreateVector(0, 0, 5)));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	890 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	891
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	892 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	893 CoordinateSystem3D c("nope1", "nope2");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	894 ASSERT_FALSE(c.IsValid());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	895 ASSERT_DOUBLE_EQ(c.GetNormal()[0], 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	896 ASSERT_DOUBLE_EQ(c.GetNormal()[1], 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	897 ASSERT_DOUBLE_EQ(c.GetNormal()[2], 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	898 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	899
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	900 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	901 // https://www.vedantu.com/maths/perpendicular-distance-of-a-point-from-a-plane
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	902 CoordinateSystem3D c = CoordinateSystem3D::CreateFromPlaneGeneralForm(2, 4, -4, -6);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	903 ASSERT_DOUBLE_EQ(3, c.ComputeDistance(LinearAlgebra::CreateVector(0, 3, 6)));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	904 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	905
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	906 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	907 // https://mathinsight.org/distance_point_plane_examples
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	908 CoordinateSystem3D c = CoordinateSystem3D::CreateFromPlaneGeneralForm(2, -2, 5, 8);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	909 ASSERT_DOUBLE_EQ(39.0 / sqrt(33.0), c.ComputeDistance(LinearAlgebra::CreateVector(4, -4, 3)));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	910 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	911
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	912 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	913 // https://www.ck12.org/calculus/distance-between-a-point-and-a-plane/lesson/Distance-Between-a-Point-and-a-Plane-MAT-ALY/
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	914 const Vector a = LinearAlgebra::CreateVector(3, 6, 9);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	915 const Vector b = LinearAlgebra::CreateVector(9, 6, 3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	916 const Vector c = LinearAlgebra::CreateVector(6, -9, 9);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	917 CoordinateSystem3D d = CoordinateSystem3D::CreateFromThreePoints(a, b, c);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	918 ASSERT_DOUBLE_EQ(0, d.ComputeDistance(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	919 ASSERT_DOUBLE_EQ(0, d.ComputeDistance(b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	920 ASSERT_DOUBLE_EQ(0, d.ComputeDistance(c));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	921 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	922
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	923 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	924 // https://tutorial.math.lamar.edu/classes/calcii/eqnsofplanes.aspx
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	925 const Vector a = LinearAlgebra::CreateVector(1, -2, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	926 const Vector b = LinearAlgebra::CreateVector(3, 1, 4);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	927 const Vector c = LinearAlgebra::CreateVector(0, -1, 2);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	928 CoordinateSystem3D d = CoordinateSystem3D::CreateFromThreePoints(a, b, c);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	929 double r = d.GetNormal() [0] / 2.0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	930 ASSERT_DOUBLE_EQ(-8 * r, d.GetNormal() [1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	931 ASSERT_DOUBLE_EQ(5 * r, d.GetNormal() [2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	932 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	933 }

Mercurial > hg > orthanc-stone

annotate UnitTestsSources/GeometryToolboxTests.cpp @ 1762:604fc328dc10