linear_algebra.cpp

Go to the documentation of this file.

#include <linear_algebra.h>
 
#include <cmath>            // sqrt
 
double dist(Point3D& p1, Point3D& p2) {
    return norm((p1-p2)*(p1-p2));
}
 
double norm2(const Point3D& p) {
    return (p.x*p.x+p.y*p.y+p.z*p.z);
}
 
double norm(const Point3D& p) {
    return sqrt(norm2(p));
}