diff options
Diffstat (limited to 'src/3rdparty/2geom/src/toys/bounds-test.cpp')
| -rw-r--r-- | src/3rdparty/2geom/src/toys/bounds-test.cpp | 171 |
1 files changed, 171 insertions, 0 deletions
diff --git a/src/3rdparty/2geom/src/toys/bounds-test.cpp b/src/3rdparty/2geom/src/toys/bounds-test.cpp new file mode 100644 index 0000000..f600162 --- /dev/null +++ b/src/3rdparty/2geom/src/toys/bounds-test.cpp @@ -0,0 +1,171 @@ +#include <2geom/d2.h> +#include <2geom/sbasis.h> +#include <2geom/bezier-to-sbasis.h> + +#include <toys/path-cairo.h> +#include <toys/toy-framework-2.h> + +#include <time.h> +#include <vector> +using std::vector; +using namespace Geom; +using namespace std; + +static void plot(cairo_t* cr, SBasis const &B,double vscale=1,double a=0,double b=1){ + D2<SBasis> plot; + plot[0]=SBasis(Linear(150+a*300,150+b*300)); + plot[1]=B*(-vscale); + plot[1]+=300; + cairo_d2_sb(cr, plot); + cairo_stroke(cr); +} +static void plot_bar(cairo_t* cr, double height, double vscale=1,double a=0,double b=1){ + cairo_move_to(cr, Geom::Point(150+300*a,-height*vscale+300)); + cairo_line_to(cr, Geom::Point(150+300*b,-height*vscale+300)); + cairo_stroke(cr); +} + +class BoundsTester: public Toy { + unsigned size; + PointSetHandle hand; + void draw(cairo_t *cr, std::ostringstream *notify, int width, int height, bool save, std::ostringstream *timer_stream) override { + + for (unsigned i=0;i<size;i++){ + hand.pts[i ][0]=150+15*(i-size); + hand.pts[i+size][0]=450+15*(i+1); + cairo_move_to(cr, Geom::Point(hand.pts[i ][0],150)); + cairo_line_to(cr, Geom::Point(hand.pts[i ][0],450)); + cairo_move_to(cr, Geom::Point(hand.pts[i+size][0],150)); + cairo_line_to(cr, Geom::Point(hand.pts[i+size][0],450)); + } + cairo_move_to(cr, Geom::Point(0,300)); + cairo_line_to(cr, Geom::Point(600,300)); + + cairo_set_line_width (cr, 1); + cairo_set_source_rgba (cr, 0.2, 0.2, 0.2, 1); + cairo_stroke(cr); + + SBasis B(size, Linear()); + for (unsigned i=0;i<size;i++){ + B[i] = Linear(-(hand.pts[i ][1]-300)*std::pow(4.,(int)i), + -(hand.pts[i+size][1]-300)*std::pow(4.,(int)i) ); + } + B.normalize(); + plot(cr,B,1); + cairo_set_source_rgba (cr, 0., 0., 0.8, 1); + cairo_stroke(cr); + + Interval bnds = *bounds_local(B,Interval(0.,.5)); + plot_bar(cr,bnds.min(),1,.0,.5); + plot_bar(cr,bnds.max(),1,.0,.5); + cairo_set_source_rgba (cr, 0.4, 0., 0., 1); + cairo_stroke(cr); + bnds = *bounds_exact(B); + plot_bar(cr,bnds.min()); + plot_bar(cr,bnds.max()); + cairo_set_source_rgba (cr, 0.9, 0., 0., 1); + cairo_stroke(cr); + +/* +This is a multi-root test... +*/ + hand.pts[2*size ][0]=150; + hand.pts[2*size+1][0]=150; + hand.pts[2*size+2][0]=150; + hand.pts[2*size ][1]=std::max(hand.pts[2*size ][1],hand.pts[2*size+1][1]); + hand.pts[2*size+1][1]=std::max(hand.pts[2*size+1][1],hand.pts[2*size+2][1]); + vector<double> levels; + levels.push_back((300-hand.pts[2*size ][1])); + levels.push_back((300-hand.pts[2*size+1][1])); + levels.push_back((300-hand.pts[2*size+2][1])); + for (double level : levels) plot_bar(cr,level); + + cairo_set_source_rgba (cr, 0., 0.5, 0., 1); + + *notify<<"Use hand.pts to set the coefficients of the s-basis."<<std::endl; + + vector<double>my_roots; + +// cairo_set_source_rgba (cr, 0.9, 0., 0.8, 1); +// for (unsigned i=0;i<levels.size();i++){ +// my_roots.clear(); +// my_roots=roots(B-Linear(levels[i])); +// for(unsigned j=0;j<my_roots.size();j++){ +// draw_cross(cr,Point(150+300*my_roots[j],300-levels[i])); +// } +// } + +// cairo_set_source_rgba (cr, 0.9, 0., 0.8, 1); + + vector<vector<double> > sols=multi_roots(B,levels,.001,.001); + //map<double,unsigned> sols=multi_roots(B,levels); + //for(map<double,unsigned>::iterator sol=sols.begin();sol!=sols.end();sol++){ + // draw_handle(cr,Point(150+300*(*sol).first,300-levels[(*sol).second])); + //} + + for (unsigned i=0;i<sols.size();i++){ + for (unsigned j=0;j<sols[i].size();j++){ + draw_handle(cr,Point(150+300*sols[i][j],300-levels[i])); + } + } + cairo_set_source_rgba (cr, 0.9, 0., 0.8, 1); + +/* + + clock_t end_t; + unsigned iterations = 0; + + my_roots.clear(); + end_t = clock()+clock_t(0.1*CLOCKS_PER_SEC); + iterations = 0; + while(end_t > clock()) { + my_roots.clear(); + for (unsigned i=0;i<levels.size();i++){ + my_roots=roots(B-Linear(levels[i])); + } + iterations++; + } + *notify << 1000*0.1/iterations <<" ms = roots time"<< std::endl; + + sols.clear(); + end_t = clock()+clock_t(0.1*CLOCKS_PER_SEC); + iterations = 0; + while(end_t > clock()) { + sols.clear(); + sols=multi_roots(B,levels); + iterations++; + } + *notify << 1000*0.1/iterations <<" ms = multi roots time"<< std::endl; +*/ + Toy::draw(cr, notify, width, height, save,timer_stream); + } + +public: + BoundsTester(){ + size=5; + if(hand.pts.empty()) { + for(unsigned i = 0; i < 2*size; i++) + hand.pts.emplace_back(0,150+150+uniform()*300*0); + } + hand.pts.emplace_back(150,300+ 50+uniform()*100); + hand.pts.emplace_back(150,300- 50+uniform()*100); + hand.pts.emplace_back(150,300-150+uniform()*100); + handles.push_back(&hand); + } +}; + +int main(int argc, char **argv) { + init(argc, argv, new BoundsTester); + return 0; +} + +/* + Local Variables: + mode:c++ + c-file-style:"stroustrup" + c-file-offsets:((innamespace . 0)(inline-open . 0)(case-label . +)) + indent-tabs-mode:nil + fill-column:99 + End: +*/ +// vim: filetype = cpp:expandtab:shiftwidth = 4:tabstop = 8:softtabstop = 4:encoding = utf-8:textwidth = 99 : |