summaryrefslogtreecommitdiff
path: root/graphics/epix/sphere.cc
diff options
context:
space:
mode:
authorNorbert Preining <norbert@preining.info>2019-09-02 13:46:59 +0900
committerNorbert Preining <norbert@preining.info>2019-09-02 13:46:59 +0900
commite0c6872cf40896c7be36b11dcc744620f10adf1d (patch)
tree60335e10d2f4354b0674ec22d7b53f0f8abee672 /graphics/epix/sphere.cc
Initial commit
Diffstat (limited to 'graphics/epix/sphere.cc')
-rw-r--r--graphics/epix/sphere.cc131
1 files changed, 131 insertions, 0 deletions
diff --git a/graphics/epix/sphere.cc b/graphics/epix/sphere.cc
new file mode 100644
index 0000000000..5f172e0fb0
--- /dev/null
+++ b/graphics/epix/sphere.cc
@@ -0,0 +1,131 @@
+/*
+ * sphere.cc -- ePiX::Sphere class and mathematical operators
+ *
+ * This file is part of ePiX, a C++ library for creating high-quality
+ * figures in LaTeX
+ *
+ * Version 1.1.9
+ * Last Change: July 30, 2007
+ */
+
+/*
+ * Copyright (C) 2001, 2002, 2003, 2004, 2005, 2006, 2007
+ * Andrew D. Hwang <rot 13 nujnat at zngupf dot ubylpebff dot rqh>
+ * Department of Mathematics and Computer Science
+ * College of the Holy Cross
+ * Worcester, MA, 01610-2395, USA
+ */
+
+/*
+ * ePiX 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; either version 2 of the License, or
+ * (at your option) any later version.
+ *
+ * ePiX 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 ePiX; if not, write to the Free Software Foundation, Inc.,
+ * 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+ */
+
+#include <cmath>
+
+#include "errors.h"
+
+#include "triples.h"
+
+#include "functions.h"
+
+#include "camera.h"
+#include "circle.h"
+
+#include "sphere.h"
+
+namespace ePiX {
+
+ Sphere::Sphere(const P& ctr, double rad)
+ : m_ctr(ctr), m_rad(rad), m_malformed(false) { }
+
+ Sphere::Sphere(const P& ctr, const P& pt)
+ : m_ctr(ctr), m_rad(norm(pt-ctr)), m_malformed(false) { }
+
+ P Sphere::center() const
+ {
+ return m_ctr;
+ }
+
+ double Sphere::radius() const
+ {
+ return m_rad;
+ }
+
+ bool Sphere::malformed() const
+ {
+ return m_malformed;
+ }
+
+ // translation
+ Sphere& Sphere::shift(const P& arg)
+ {
+ if (!m_malformed)
+ m_ctr += arg;
+ return *this;
+ }
+
+ Sphere& Sphere::move_to(const P& arg)
+ {
+ if (!m_malformed)
+ m_ctr = arg;
+ return *this;
+ }
+
+ Sphere& Sphere::scale(double& arg)
+ {
+ if (!m_malformed)
+ m_rad *= arg;
+ return *this;
+ }
+
+
+ void Sphere::draw() const
+ {
+ if (m_malformed)
+ return;
+
+ // else
+ P location(cam().viewpt());
+
+ P dir(m_ctr - location);
+ double dist(norm(dir));
+
+ if (dist <= m_rad)
+ epix_warning("Cannot draw sphere enclosing camera");
+
+ else
+ {
+ dir *= 1/dist;
+ double x(m_rad*m_rad/dist);
+
+ Circle temp(m_ctr - x*dir, sqrt((m_rad-x)*(m_rad+x)), dir);
+ temp.draw();
+ }
+ } // end of Sphere::photo()
+
+
+ //// global functions ////
+ void sphere(const P& ctr, double rad)
+ {
+ Sphere S(ctr, rad);
+ S.draw();
+ }
+
+ void sphere(const P& ctr, const P& pt)
+ {
+ Sphere S(ctr, pt);
+ S.draw();
+ }
+} // end of namespace