text: // Copyright (c) 2000  
  // Utrecht University (The Netherlands),
  // ETH Zurich (Switzerland),
  // INRIA Sophia-Antipolis (France),
  // Max-Planck-Institute Saarbruecken (Germany),
  // and Tel-Aviv University (Israel).  All rights reserved. 
  //
  // This file is part of CGAL (www.cgal.org)
  //
  // $URL: https://github.com/CGAL/cgal/blob/releases/CGAL-5.0.2/Cartesian_kernel/include/CGAL/Cartesian/Scaling_rep_3.h $
  // $Id: Scaling_rep_3.h 52164b1 2019-10-19T15:34:59+02:00 Sébastien Loriot
  // SPDX-License-Identifier: LGPL-3.0-or-later OR LicenseRef-Commercial
  // 
  //
  // Author(s)     : Herve Bronnimann
  
  #ifndef CGAL_CARTESIAN_SCALING_REP_3_H
  #define CGAL_CARTESIAN_SCALING_REP_3_H
  
  namespace CGAL {
  
  template < class R >
  class Scaling_repC3 : public Aff_transformation_rep_baseC3<R>
  {
    friend class Aff_transformation_repC3<R>;
    friend class Translation_repC3<R>;
  
  public:
    typedef typename R::FT                                FT;
    typedef Aff_transformation_rep_baseC3<R>              Transformation_base_3;
    typedef Aff_transformation_repC3<R>                   Transformation_3;
    typedef Translation_repC3<R>                          Translation_3;
    typedef Scaling_repC3<R>                              Scaling_3;
    typedef typename Transformation_base_3::Point_3       Point_3;
    typedef typename Transformation_base_3::Vector_3      Vector_3;
    typedef typename Transformation_base_3::Direction_3   Direction_3;
    typedef typename Transformation_base_3::Aff_transformation_3
  	                                                Aff_transformation_3;
  
    Scaling_repC3() {}
    Scaling_repC3(const FT &s) : scalefactor_(s) {}
    virtual ~Scaling_repC3() {}
  
    virtual Point_3      transform(const Point_3 &p) const
    {
      return Point_3(scalefactor_ * p.x(),
                     scalefactor_ * p.y(),
                     scalefactor_ * p.z());
    }
  
    virtual Vector_3     transform(const Vector_3 &v) const
    {
      return Vector_3(scalefactor_ * v.x(), scalefactor_ * v.y(),
                      scalefactor_ * v.z());
    }
  
    virtual Direction_3  transform(const Direction_3 &d) const
    {
      return d;
    }
  
    virtual Aff_transformation_3 operator*(const Transformation_base_3 &t) const
    {
      return t.compose(*this);
    }
  
    virtual Aff_transformation_3 compose(const Transformation_3 &t) const
    {
      return Aff_transformation_3(scalefactor_ * t.t11,
                                  scalefactor_ * t.t12,
                                  scalefactor_ * t.t13,
                                  t.t14,
  				
                                  scalefactor_ * t.t21,
                                  scalefactor_ * t.t22,
                                  scalefactor_ * t.t23,
                                  t.t24,
  				
                                  scalefactor_ * t.t31,
                                  scalefactor_ * t.t32,
                                  scalefactor_ * t.t33,
                                  t.t34);
    }
  
    virtual Aff_transformation_3 compose(const Translation_3 &t) const
    {
      FT ft0(0);
      return Aff_transformation_3(scalefactor_,
                                  ft0,
                                  ft0,
  				t.translationvector_.x(),
  				
                                  ft0,
  				scalefactor_,
                                  ft0,
  				t.translationvector_.y(),
  				
                                  ft0,
                                  ft0,
  				scalefactor_,
  				t.translationvector_.z());
    }
  
    virtual Aff_transformation_3 compose(const Scaling_3 &t) const
    {
      return Aff_transformation_3(SCALING, scalefactor_*t.scalefactor_);
    }
  
    virtual Aff_transformation_3 inverse() const
    {
      return Aff_transformation_3(SCALING, FT(1)/scalefactor_);
    }
  
    virtual Aff_transformation_3 transpose() const
    {
      return Aff_transformation_3(SCALING, scalefactor_);
    }
  
    virtual bool is_even() const
    {
      return true;
    }
  
    virtual FT cartesian(int i, int j) const
    {
      if (i!=j) return FT(0);
      if (i==3) return FT(1);
      return scalefactor_;
    }
  
    virtual std::ostream &print(std::ostream &os) const
    {
      os << "Aff_transformationC3(" << scalefactor_ << ")";
      return os;
    }
  
  private:
    FT   scalefactor_;
  };
  
  } //namespace CGAL
  
  #endif // CGAL_CARTESIAN_SCALING_REP_3_H
  