line.h

/******************************************************************************
 *   Copyright (C) 2006-2024 by the GIMLi development team                    *
 *   Carsten Rücker carsten@resistivity.net                                   *
 *                                                                            *
 *   Licensed under the Apache License, Version 2.0 (the "License");          *
 *   you may not use this file except in compliance with the License.         *
 *   You may obtain a copy of the License at                                  *
 *                                                                            *
 *       http://www.apache.org/licenses/LICENSE-2.0                           *
 *                                                                            *
 *   Unless required by applicable law or agreed to in writing, software      *
 *   distributed under the License is distributed on an "AS IS" BASIS,        *
 *   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. *
 *   See the License for the specific language governing permissions and      *
 *   limitations under the License.                                           *
 *                                                                            *
 ******************************************************************************/
 
#ifndef GIMLI_LINE__H
#define GIMLI_LINE__H
 
#include "gimli.h"
#include "pos.h"
 
namespace GIMLI{


 
class DLLEXPORT Line {
public:
    Line();

    Line(const RVector3 & p0);

    Line(const RVector3 & p0, const RVector3 & p1);

    Line(const Line & line);

    ~Line();

    Line & operator = (const Line & line);

    inline const RVector3 operator()(double t) const { return this->at(t); }

    inline bool operator == (const Line & line) const { return this->compare(line); }

    inline bool operator != (const Line & line) const { return !(*this == line); }

    bool checkValidity(double tol=TOLERANCE);

    inline bool valid() const { return valid_; }

    bool compare(const Line & line, double tol=TOLERANCE) const;

    const RVector3 & p0() const { return _p0; }

    const RVector3 & p1() const { return _p1; }

    bool intersectRay(const RVector3 & start, const RVector3 & dir,
                      RVector3 & pos, double tol=TOLERANCE) const;

    RVector3 intersect(const RVector3 & start, const RVector3 & dir, double tol=TOLERANCE) const;

    RVector3 intersect(const Line & line, double tol=TOLERANCE) const;

    double nearest(const RVector3 & p) const;

    inline RVector3 at(double t) const { return _p0 + (_p1 - _p0) * t; }

    double distance(const RVector3 & pos) const;

    double t(const RVector3 & pos, double tol=TOLERANCE) const ;

    bool touch1(const RVector3 & pos, int & pFunIdx, double tol=TOLERANCE) const ;

    bool touch(const RVector3 & pos, double tol=TOLERANCE) const ;

    double length() const { return _p0.dist(_p1); }
 
protected:

    void copy_(const Line & line);
 
    RVector3 _p0;
    RVector3 _p1;
    bool valid_;
};
 
std::ostream & operator << (std::ostream & str, const Line & l);
 
} // namespace GIMLI;
 
#endif // GIMLI_LINE__H