diff options
| author | Timothy Pearson <kb9vqf@pearsoncomputing.net> | 2012-07-11 14:15:27 -0500 |
|---|---|---|
| committer | Timothy Pearson <kb9vqf@pearsoncomputing.net> | 2012-07-11 14:15:27 -0500 |
| commit | b85a292ce06475d560bfa1195b63a8bfe211f22d (patch) | |
| tree | 463d71be55ff807513139f1de106aef6bdd7b4db /lib/tqwtplot3d/examples/mesh2/src/functions.h | |
| parent | ce039289815e2802fdeca8d384126c807ca9cb58 (diff) | |
| download | ulab-b85a292ce06475d560bfa1195b63a8bfe211f22d.tar.gz ulab-b85a292ce06475d560bfa1195b63a8bfe211f22d.zip | |
Add 0.2.7 release of qwtplot3d for future TQt3 conversion and use
Diffstat (limited to 'lib/tqwtplot3d/examples/mesh2/src/functions.h')
| -rw-r--r-- | lib/tqwtplot3d/examples/mesh2/src/functions.h | 214 |
1 files changed, 214 insertions, 0 deletions
diff --git a/lib/tqwtplot3d/examples/mesh2/src/functions.h b/lib/tqwtplot3d/examples/mesh2/src/functions.h new file mode 100644 index 0000000..1a86e0d --- /dev/null +++ b/lib/tqwtplot3d/examples/mesh2/src/functions.h @@ -0,0 +1,214 @@ +#ifndef __EXAMPLE_H__
+#define __EXAMPLE_H__
+
+#include <math.h>
+#include "../../../include/qwt3d_parametricsurface.h"
+#include "../../../include/qwt3d_function.h"
+
+using namespace Qwt3D;
+
+
+class Rosenbrock : public Function
+{
+public:
+
+ Rosenbrock(SurfacePlot& pw)
+ :Function(pw)
+ {
+ }
+
+ double operator()(double x, double y)
+ {
+ return 0.7 * log10((1-x)*(1-x) + 10 * (y - x*x)*(y - x*x));
+ }
+// QString name() const {return "Almost {\\it Rosenbrock}\\,:\\quad$\\frac{\\ln((1-x)^2 + 100(y-x^2)^2)}{8}$";}
+};
+
+class Hat : public Function
+{
+public:
+
+ Hat(SurfacePlot& pw)
+ :Function(pw)
+ {
+ //setMinZ(0.3);
+ setDomain(0,10,0,10);
+ }
+
+ double operator()(double x, double y)
+ {
+ return 1.0 / (x*x+y*y+0.5);
+ //return x*x*y/(x*x*x*x+y*y);
+ }
+};
+
+
+class Ripple : public Function
+{
+public:
+
+ Ripple(SurfacePlot& pw)
+ :Function(pw)
+ {
+ double l = 12;
+ setDomain(-l,l,-l,l);
+ }
+
+ double operator()(double x, double y)
+ {
+ return (cos(sqrt(x*x+y*y) + cos(sqrt(((x+.913*2*Qwt3D::PI)*(x+.913*2*Qwt3D::PI))+y*y))
+ + cos(sqrt(((x-.913*2*Qwt3D::PI)*(x-.913*2*Qwt3D::PI))+(y*y))))*4);
+ }
+};
+
+class Saddle : public Function
+{
+public:
+
+ Saddle()
+ :Function()
+ {
+ // setMaxZ(0.8);
+ }
+
+ double operator()(double x, double y)
+ {
+ return x*x - y*y;
+ }
+// QString name() const {return "$x^2-y^2$";}
+};
+
+class Mex : public Function
+{
+public:
+
+ Mex()
+ :Function()
+ {
+ // setMaxZ(0.8);
+ }
+
+ double operator()(double x, double y)
+ {
+ double n = sqrt(x*x+y*y);
+
+ if (n < DBL_MIN)
+ return 20;
+
+ return 20 * sin(sqrt(x*x+y*y)) / n;
+ }
+// QString name() const {return "$\\frac{20\\sin\\sqrt{x^2+y^2}}{\\sqrt{x^2+y^2}}$";}
+};
+
+class Torus : public ParametricSurface
+{
+public:
+
+ Torus(SurfacePlot& pw)
+ :ParametricSurface(pw)
+ {
+ setMesh(41,31);
+ setDomain(-2*Qwt3D::PI, 0,-2*Qwt3D::PI,0);
+ setPeriodic(true,true);
+ }
+
+
+ Triple operator()(double u, double v)
+ {
+ double x,y,z;
+ double c = 1.9;
+ x = (c + cos(v)) * cos(u);
+ y = (c + cos(v)) * sin(u);
+ z = sin(v) + cos(v);
+ return Triple(x,y,z);
+ }
+};
+
+class Seashell : public ParametricSurface
+{
+public:
+
+ Seashell(SurfacePlot& pw)
+ :ParametricSurface(pw)
+ {
+ setMesh(41,131);
+ setDomain(0,2*Qwt3D::PI,0,2*Qwt3D::PI);
+ setPeriodic(true,true);
+ }
+
+
+ Triple operator()(double u, double v)
+ {
+ double x,y,z;
+ double a = 1;
+ double b = 6;
+ double c = 0.5;
+ int n = 3;
+
+ double f = v/(2*Qwt3D::PI);
+
+ x = a*(1-f)*cos(n*v)*(1+cos(u)) + c*cos(n*v) ;
+ y = a*(1-f)*sin(n*v)*(1+cos(u)) + c*sin(n*v) ;
+ z = b*f + a*(1-f)*sin(u);
+ return Triple(x,y,z);
+ }
+};
+
+class Boy : public ParametricSurface
+{
+public:
+
+ Boy(SurfacePlot& pw)
+ :ParametricSurface(pw)
+ {
+ setMesh(141,131);
+ setDomain(0,Qwt3D::PI,0,Qwt3D::PI);
+ setPeriodic(true,true);
+ }
+
+
+ Triple operator()(double u, double v)
+ {
+ double x,y,z;
+ double a = 2/3.;
+ double b = sqrt(2.);
+
+ x = a*(cos(u)*cos(2*v)+b*sin(u)*cos(v))*cos(u) / (b-sin(2*u)*sin(3*v));
+ y = a*(cos(u)*sin(2*v)-b*sin(u)*sin(v))*cos(u) / (b-sin(2*u)*sin(3*v));
+ z = b*cos(u)*cos(u) / (b-sin(2*u)*sin(2*v));
+
+ return Triple(x,y,z);
+ }
+};
+
+class Dini : public ParametricSurface
+{
+public:
+
+ Dini(SurfacePlot& pw)
+ :ParametricSurface(pw)
+ {
+ setMesh(141,35);
+ setDomain(0,5*Qwt3D::PI,0.001, 2);
+ setPeriodic(true,true);
+ }
+
+
+ Triple operator()(double u, double v)
+ {
+ double x,y,z;
+ double a = 5;
+ double b = 1;
+
+
+ x=a*cos(u)*sin(v);
+ y=a*sin(u)*sin(v);
+ z=a*(cos(v)+log(tan(v/2)))+b*u;
+
+ return Triple(x,y,z);
+ }
+};
+
+void createCone(Qwt3D::TripleField& conepos, Qwt3D::CellField& conecell);
+
+#endif
|