00001 // This file is part of ff3d - http://www.freefem.org/ff3d 00002 // Copyright (C) 2001, 2002, 2003 Stéphane Del Pino 00003 00004 // This program is free software; you can redistribute it and/or modify 00005 // it under the terms of the GNU General Public License as published by 00006 // the Free Software Foundation; either version 2, or (at your option) 00007 // any later version. 00008 00009 // This program is distributed in the hope that it will be useful, 00010 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00011 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00012 // GNU General Public License for more details. 00013 00014 // You should have received a copy of the GNU General Public License 00015 // along with this program; if not, write to the Free Software Foundation, 00016 // Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. 00017 00018 // $Id: InfiniteCylinder.cpp,v 1.3 2006/10/01 23:25:58 delpinux Exp $ 00019 00020 #include <InfiniteCylinder.hpp> 00021 #include <Cylinder.hpp> 00022 00023 InfiniteCylinder:: 00024 InfiniteCylinder(const Cylinder& c) 00025 : Shape(InfiniteCylinder::infiniteCylinder), 00026 __center(c.__center), 00027 __unaryVector(c.__unaryVector), 00028 __radius2(c.__radius2) 00029 { 00030 this->setTransformationsList(c.transformationsList()); 00031 } 00032 00033 InfiniteCylinder:: 00034 InfiniteCylinder(const InfiniteCylinder& c) 00035 : Shape(c), 00036 __center(c.__center), 00037 __unaryVector(c.__unaryVector), 00038 __radius2(c.__radius2) 00039 { 00040 ; 00041 } 00042 00043 bool InfiniteCylinder:: 00044 __inShape(const TinyVector<3, real_t>& v) const 00045 { 00046 TinyVector<3, real_t> x = v-__center; 00047 00048 real_t ux = __unaryVector * x; 00049 TinyVector<3, real_t> normal (x - (ux * __unaryVector)); 00050 return (normal * normal < __radius2); 00051 } 00052 00053 ReferenceCounting<Shape> InfiniteCylinder:: 00054 __getCopy() const 00055 { 00056 return new InfiniteCylinder(*this); 00057 }
1.5.6