IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 424707ff authored by Maximilian Reininghaus's avatar Maximilian Reininghaus :vulcan:
Browse files

geometry example usage

parent 853ddcee
No related branches found
No related tags found
No related merge requests found
/*
from within ngc/build/ directory:
g++ -O3 -I ../third_party/ -I<path to eigen> -I ../ -Wall -pedantic \
-std=c++14 ../Main/geometry_example.cc \
../Framework/Geometry/CoordinateSystem.cc -o geometry_example
*/
#include <Framework/Geometry/Vector.h>
#include <Framework/Geometry/Sphere.h>
#include <Framework/Geometry/Point.h>
#include <Framework/Geometry/CoordinateSystem.h>
#include <phys/units/quantity.hpp>
#include <phys/units/io.hpp>
#include <iostream>
#include <cstdlib>
int main()
{
using namespace phys::units;
using namespace phys::units::io;
using namespace phys::units::literals;
CoordinateSystem root;
Point const p1(root, {0_m, 0_m, 0_m});
CoordinateSystem cs2 = root.translate({0_m, 0_m, 1_m});
Point const p2(cs2, {0_m, 0_m, 0_m});
auto const diff = p2 - p1;
auto const norm = diff.squaredNorm();
std::cout << diff.getComponents() << std::endl;
std::cout << norm << std::endl;
Sphere s(p1, 10_m);
std::cout << s.isInside(p2) << std::endl;
Sphere s2(p1, 3_um);
std::cout << s2.isInside(p2) << std::endl;
return EXIT_SUCCESS;
}
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment