1 files changed, 77 insertions, 0 deletions
diff --git a/kpovmodeler/pmsorsegment.h b/kpovmodeler/pmsorsegment.h
new file mode 100644
index 00000000..f0c5a2ea
--- /dev/null
+++ b/kpovmodeler/pmsorsegment.h
@@ -0,0 +1,77 @@
+//-*-C++-*-
+/*
+**************************************************************************
+                                 description
+                             --------------------
+    copyright            : (C) 2002 by Andreas Zehender
+    email                : [email protected]
+**************************************************************************
+
+**************************************************************************
+*                                                                        *
+*  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; either version 2 of the License, or     *
+*  (at your option) any later version.                                   *
+*                                                                        *
+**************************************************************************/
+
+
+#ifndef PMSORSEGMENT_H
+#define PMSORSEGMENT_H
+
+#ifdef HAVE_CONFIG_H
+#include <config.h>
+#endif
+
+#include "pmvector.h"
+
+/**
+ * Helper class for sors
+ *
+ * Each instance of this class represents one sor segment. A point
+ * on the segment is given by the equation
+ *
+ * fi(t) = A[i] * t^3 + B[i] * t^2 + C[i] * t + D[i]
+ *
+ * with t ranging from 0 to 1.
+ */
+class PMSorSegment
+{
+public:
+   /**
+    * Standard constructor
+    */
+   PMSorSegment( )
+   {
+      m_a = m_b = m_c = m_d = 0.0;
+   }
+   /**
+    * Constructor that calculates the segment
+    */
+   PMSorSegment( const PMVector& p0, const PMVector& p1,
+                 const PMVector& p2, const PMVector& p3 )
+   {
+      calculateSor( p0, p1, p2, p3 );
+   }
+   
+   /**
+    * Returns a 2D vector with the point on the sor segment
+    */
+   PMVector point( double t ) const;
+   /**
+    * Returns the gradient on the sor
+    */
+   //PMVector gradient( double t ) const;
+   
+   /**
+    * Calculates the sor parameters
+    */
+   void calculateSor( const PMVector& p0, const PMVector& p1,
+                      const PMVector& p2, const PMVector& p3 );
+private:
+   double m_a, m_b, m_c, m_d;
+   double m_t, m_s;
+};
+
+#endif