Integrator.hpp

//
// Integrator Tool for SAL3D.
// Copyright 2010 AQSENSE, S.L.
//

#ifndef AQSENSE_SAL3DPP_INTEGRATOR_HPP
#define AQSENSE_SAL3DPP_INTEGRATOR_HPP

#include <sal3d/integrator.h>
#include <sal3dpp/Mesh.hpp>
#include <sal3dpp/COP.hpp>
#include <stdio.h>

//REVIEW Public API, doxygen

namespace sal3d
{
    class Integrator
    {
    public:

        Integrator(const sal3d::COP& aCop)
        {
            _integrator.p = 0;
            sal3d_cop cop = aCop.c_cop();

            Error e;
            sal3d_integrator_new_from_cop(&_integrator, cop, &e.e);
            if (e.e.value < 0)
                throw e;
        };

        Integrator(const std::vector<sal3d::COP> &aCops)
        {
            _integrator.p = 0;
            std::vector<sal3d_cop> ccops;
            for(unsigned int i = 0; i < aCops.size(); i++)
            {
                ccops.push_back(aCops[i].c_cop());
            }

            Error e;
            sal3d_integrator_new_from_cops(&_integrator,
                    &ccops[0], static_cast<int>(ccops.size()), &e.e);
            if (e.e.value < 0)
                throw e;
        };

        Integrator(const Integrator &aIntegrator):
            _integrator(aIntegrator._integrator)
        {
            Error e;
            sal3d_integrator_share(_integrator, &e.e);
            if (e.e.value < 0)
                throw e;
        }

        explicit Integrator(sal3d_integrator aInteg):
            _integrator(_integrator)
        {
            Error e;
            sal3d_integrator_is_valid(aInteg, &e.e);
            if (e.e.value < 0)
                throw e;
        }

        ~Integrator()
        {
            sal3d_integrator_release(_integrator);
        }

        Mesh
        Reconstruct(unsigned int aDepth, unsigned int aSPN, int aVerbose = 0) const
        {
            sal3d_integrator_params params = sal3d_integrator_params_default;
            params.iDepth = aDepth;
            params.iVerbose = aVerbose;
            params.iSamplesPerNode = aSPN;

            Error e;
            sal3d_integrator_reconstruct_surface(_integrator, &params, &e.e);
            if (e.e.value < 0)
                throw e;

            sal3d_mesh cmesh;
            sal3d_integrator_get_mesh(_integrator, &cmesh, &e.e);
            if (e.e.value < 0)
                throw e;

            Mesh mesh(cmesh); 

            return mesh;
        }

        sal3d_integrator
        c_integrator() const
        {
            return _integrator;
        }

    private:
        sal3d_integrator _integrator;
    };
}

#endif