plucker.cc

Go to the documentation of this file.

#include "plucker.hh"

using namespace std;

Plucker::Plucker()
: coordinates_({0,0,0,0,0,0}),
  scale_(0),
  computed_(false),
  points_{nullptr,nullptr}
{}


Plucker::Plucker(const Node* a, const Node* b)
{
    ASSERT_DBG(a);
    ASSERT_DBG(b);
    points_[0] = a;
    points_[1] = b;
    coordinates_(arma::span(0,2)) = point(1) - point(0);
    
    // Check empty
    ASSERT_DBG(arma::norm(coordinates_(arma::span(0,2)),2) > 0);

    scale_ = 0;
    scale_ = std::max(  scale_, std::fabs(coordinates_[0]));
    scale_ = std::max(  scale_, std::fabs(coordinates_[1]));
    scale_ = std::max(  scale_, std::fabs(coordinates_[2]));
    
    computed_ = false;
}

Plucker::Plucker(const Node* a, const Node* b, bool compute_pc)
: Plucker(a,b)
{
    if(compute_pc) compute();
}


double Plucker::operator*(const Plucker &b){
    return (coordinates_[0]*b[3]) + (coordinates_[1]*b[4]) + (coordinates_[2]*b[5]) + (coordinates_[3]*b[0]) + (coordinates_[4]*b[1]) +(coordinates_[5]*b[2]);
}


void Plucker::compute(){
    if(computed_) return;
    
    ASSERT_DBG(points_[0]);
    ASSERT_DBG(points_[1]);
    coordinates_[3] = coordinates_[1]*point(0)[2] - coordinates_[2]*point(0)[1];
    coordinates_[4] = coordinates_[2]*point(0)[0] - coordinates_[0]*point(0)[2];
    coordinates_[5] = coordinates_[0]*point(0)[1] - coordinates_[1]*point(0)[0];
    computed_ = true;
}

ostream& operator<<(ostream& os, const Plucker& p)
{
    if(p.computed_){
        os <<"(" << p.coordinates_[0] << "," << p.coordinates_[1] << "," << p.coordinates_[2] << "," 
           << p.coordinates_[3] << "," << p.coordinates_[4] << "," << p.coordinates_[5] << ")";
    }else{
        os << "NULL (Plucker coords have not been computed)";
    }
    return os;
}