AlgObj3DInt.h

//===========================================================================
// GoTools - SINTEF Geometry Tools
//
// GoTools module: Intersections, version 1.0
//
// Copyright (C) 2000-2007 SINTEF ICT, Applied Mathematics, Norway.
//
// This program is free software; you can redistribute it and/or          
// modify it under the terms of the GNU General Public License            
// as published by the Free Software Foundation version 2 of the License. 
//
// This program is distributed in the hope that it will be useful,        
// but WITHOUT ANY WARRANTY; without even the implied warranty of         
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the          
// GNU General Public License for more details.                           
//
// You should have received a copy of the GNU General Public License      
// along with this program; if not, write to the Free Software            
// Foundation, Inc.,                                                      
// 59 Temple Place - Suite 330,                                           
// Boston, MA  02111-1307, USA.                                           
//
// Contact information: E-mail: tor.dokken@sintef.no                      
// SINTEF ICT, Department of Applied Mathematics,                         
// P.O. Box 124 Blindern,                                                 
// 0314 Oslo, Norway.                                                     
//
// Other licenses are also available for this software, notably licenses
// for:
// - Building commercial software.                                        
// - Building software whose source code you wish to keep private.        
//===========================================================================
#ifndef _ALGOBJ3DINT_H
#define _ALGOBJ3DINT_H


#include "AlgObjectInt.h"
#include "Array.h"
#include "BernsteinTetrahedralPoly.h"
#include "BaryCoordSystem.h"
#include <vector>


namespace Go {




struct Alg3DElem {

    double factor_;

    Array<int, 3> degrees_; // (deg_x, deg_y, deg_z).

    Alg3DElem(double factor, int degree_x, int degree_y, int degree_z)
        : factor_(factor), 
          degrees_(Array<int, 3>(degree_x, degree_y, degree_z))
    {}
};



class AlgObj3DInt : public AlgObjectInt {
public:
    AlgObj3DInt(int degree);

    AlgObj3DInt(const std::vector<Alg3DElem>& terms);

    AlgObj3DInt(const BernsteinTetrahedralPoly& implicit,
                const BaryCoordSystem3D& bc);

    virtual ~AlgObj3DInt();

    int numTerms()
    { return terms_.size(); }

    int degree()
    { return degree_; }

    Alg3DElem term(int index);

    void getImplicit(BernsteinTetrahedralPoly& impl, BaryCoordSystem3D& bc)
    {
        impl = implicit_;
        bc = bc_;
    }

    bool usingPowerBasis()
    { return power_basis_; }


protected:
    int degree_; // Total degree, i.e. the highest sum of the degrees
                 // of a term in the polynomial.

    std::vector<Alg3DElem> terms_;
    bool power_basis_; // The other option is Bernstein basis &
                       // barycentric coordinates.

    // @@sbr Not sure which to use!
    // Should the approximative impclicitization take place inside
    // this function?  Would assume that it was to be calculated on
    // the outside and only stored here.  Could then use class in
    // Implicitization, which depends on that structure not to change.
    // But tu implement easy support for converting from power to
    // Bernstein basis this approach makes sense.
    BernsteinTetrahedralPoly implicit_; // deg &
                                        // (deg+1)*(deg+2)*(deg+3)/6
                                        // basis elements.
    BaryCoordSystem3D bc_; // 4 3D-corners. Requires boundingbox for
                           // alg obj.

private:

//     // We need a utility function to change from power basis to
//     // Bernstein, using the format in implicitization.
//     BernsteinTetrahedralPoly powerToBernstein();

//     // We convert x to the Bernstein basis.
//     SplineSurface unitSf();

};


} // namespace Go


#endif // _ALGOBJ3DINT_H

---