Inkscape::DrawingItem *SPClipPath::show(Inkscape::Drawing &drawing, unsigned int key) {
Inkscape::DrawingGroup *ai = new Inkscape::DrawingGroup(drawing);
display = sp_clippath_view_new_prepend(display, key, ai);
for ( SPObject *child = firstChild() ; child ; child = child->getNext() ) {
if (SP_IS_ITEM(child)) {
Inkscape::DrawingItem *ac = SP_ITEM(child)->invoke_show(drawing, key, SP_ITEM_REFERENCE_FLAGS);
if (ac) {
/* The order is not important in clippath */
ai->appendChild(ac);
}
}
}
if (clipPathUnits == SP_CONTENT_UNITS_OBJECTBOUNDINGBOX && display->bbox) {
Geom::Affine t = Geom::Scale(display->bbox->dimensions());
t.setTranslation(display->bbox->min());
ai->setChildTransform(t);
}
ai->setStyle(this->style);
return ai;
}
void SPClipPath::update(SPCtx* ctx, unsigned int flags) {
if (flags & SP_OBJECT_MODIFIED_FLAG) {
flags |= SP_OBJECT_PARENT_MODIFIED_FLAG;
}
flags &= SP_OBJECT_MODIFIED_CASCADE;
GSList *l = NULL;
for ( SPObject *child = this->firstChild(); child; child = child->getNext()) {
sp_object_ref(child);
l = g_slist_prepend(l, child);
}
l = g_slist_reverse(l);
while (l) {
SPObject *child = SP_OBJECT(l->data);
l = g_slist_remove(l, child);
if (flags || (child->uflags & (SP_OBJECT_MODIFIED_FLAG | SP_OBJECT_CHILD_MODIFIED_FLAG))) {
child->updateDisplay(ctx, flags);
}
sp_object_unref(child);
}
for (SPClipPathView *v = this->display; v != NULL; v = v->next) {
Inkscape::DrawingGroup *g = dynamic_cast<Inkscape::DrawingGroup *>(v->arenaitem);
if (this->clipPathUnits == SP_CONTENT_UNITS_OBJECTBOUNDINGBOX && v->bbox) {
Geom::Affine t = Geom::Scale(v->bbox->dimensions());
t.setTranslation(v->bbox->min());
g->setChildTransform(t);
} else {
g->setChildTransform(Geom::identity());
}
}
}
void draw(cairo_t *cr,
std::ostringstream *notify,
int width, int height, bool save, std::ostringstream *timer_stream) {
if (first_time)
{
first_time = false;
sliders[0].geometry(Point(50, 50), 100);
}
size_t const num_points = static_cast<size_t>(sliders[0].value());
D2<SBasis> B1 = b_handle.asBezier();
Piecewise<D2<SBasis> >B;
B.concat(Piecewise<D2<SBasis> >(B1));
// testing fuse_nearby_ends
std::vector< Piecewise<D2<SBasis> > > pieces;
pieces = fuse_nearby_ends(split_at_discontinuities(B),9);
Piecewise<D2<SBasis> > C;
for (unsigned i=0; i<pieces.size(); i++){
C.concat(pieces[i]);
}
// testing fuse_nearby_ends
cairo_set_line_width (cr, 2.);
cairo_set_source_rgba (cr, 0., 0.5, 0., 1);
//cairo_d2_sb(cr, B1);
//cairo_pw_d2_sb(cr, C);
//cairo_pw_d2_sb(cr, B);
cairo_stroke(cr);
Timer tm;
Timer::Time als_time = tm.lap();
cairo_set_source_rgba (cr, 0., 0., 0.9, 1);
//dot_plot(cr,uniform_B);
cairo_stroke(cr);
std::cout << B[0] << std::endl;
Geom::Affine translation;
Geom::Path original_path;
//original_bezier.append(B[0]);
//original_bezier.appendNew<CubicBezier> (B[0]);
CubicBezier original_bezier(b_handle.pts);
original_path.append(original_bezier);
std::vector<double> initial_t;
std::vector<Geom::Point> curve_points;
if (randomize_times) {
std::uniform_real_distribution<double> dist_t(0,1);
for (size_t ii = 0; ii < num_points; ++ii) {
double const t = dist_t(generator);
initial_t.push_back(t);
}
std::sort(initial_t.begin(), initial_t.end());
double const min = initial_t.front();
double const max = initial_t.back();
for (auto& t : initial_t) {
t = (t-min)/(max-min);
}
for (auto const t : initial_t) {
curve_points.push_back(original_bezier.pointAt(t));
}
}
else {
for (size_t ii = 0; ii < num_points; ++ii) {
double const t = static_cast<double>(ii) / (num_points-1);
Geom::Point const p = original_bezier.pointAt(t);
initial_t.push_back(t);
curve_points.push_back(p);
}
}
cairo_set_source_rgba (cr, 0., 0., .9, 1);
cairo_path(cr, original_path);
draw_text(cr, original_path.initialPoint(), "original curve and old fit");
Geom::CubicBezier fitted_new;
Geom::CubicBezier fitted_new_a;
Geom::Point very_old_version_raw[4];
bezier_fit_cubic(very_old_version_raw, curve_points.data(), curve_points.size(), 0.);
Geom::CubicBezier very_old_bezier(
very_old_version_raw[0],
very_old_version_raw[1],
very_old_version_raw[2],
very_old_version_raw[3]
);
Geom::Path very_old_version_path;
very_old_version_path.append(very_old_bezier);
cairo_set_source_rgba (cr, .7, .7, 0., 1);
cairo_stroke(cr);
cairo_path(cr, very_old_version_path);
cairo_set_source_rgba (cr, 0., 0., .9, 1);
cairo_stroke(cr);
cross_plot(cr, curve_points);
if(1) {
Geom::CubicBezier combination(very_old_bezier);
tm.ask_for_timeslice();
tm.start();
auto new_result_ig_a = experiment::fit_bezier(combination, curve_points);
als_time = tm.lap();
*notify << "Bezier fit a, old algorithm as initial guess, time = " << als_time << std::endl
<< "Worst residual: " << new_result_ig_a.first << " at t=" << new_result_ig_a.second << std::endl;
Geom::Path combination_path;
translation.setTranslation(Geom::Point(300,300));
combination_path.append(combination.transformed(translation));
cairo_set_source_rgba (cr, .0, .0, .9, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, combination_path);
draw_text(cr, combination_path.initialPoint(), "old fit as i.g.");
}
{
tm.ask_for_timeslice();
tm.start();
auto new_result = fit_bezier(fitted_new, curve_points);
als_time = tm.lap();
*notify << "Bezier fit, time = " << als_time << std::endl
<< "Worst residual: " << new_result.first << " at t=" << new_result.second << std::endl;
Geom::Path fitted_new_path;
translation.setTranslation(Geom::Point(300,0));
fitted_new_path.append(fitted_new.transformed(translation));
cairo_set_source_rgba (cr, .0, .9, .0, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, fitted_new_path);
draw_text(cr, fitted_new_path.initialPoint(), "new fit");
}
{
tm.ask_for_timeslice();
tm.start();
auto new_result_a = experiment::fit_bezier(fitted_new_a, curve_points);
als_time = tm.lap();
*notify << "Bezier fit a, time = " << als_time << std::endl
<< "Worst residual: " << new_result_a.first << " at t=" << new_result_a.second << std::endl;
Geom::Path fitted_new_a_path;
translation.setTranslation(Geom::Point(0,300));
fitted_new_a_path.append(fitted_new_a.transformed(translation));
cairo_set_source_rgba (cr, .9, .0, .0, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, fitted_new_a_path);
draw_text(cr, fitted_new_a_path.initialPoint(), "new fit (a)");
}
Geom::CubicBezier fixed_times_bezier;
{
tm.ask_for_timeslice();
tm.start();
auto fixed_times_result = experiment::fit_bezier_fixed_times(fixed_times_bezier, curve_points);
als_time = tm.lap();
*notify << "Bezier fit a (fixed times), time = " << als_time << std::endl
<< "Worst residual: " << fixed_times_result.first << " at t=" << fixed_times_result.second << std::endl;
Geom::Path fixed_times_path;
translation.setTranslation(Geom::Point(600,300));
fixed_times_path.append(fixed_times_bezier.transformed(translation));
cairo_set_source_rgba (cr, .9, .0, .0, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, fixed_times_path);
draw_text(cr, fixed_times_path.initialPoint(), "fixed t fit (a)");
}
Geom::CubicBezier fixed_times_ig_bezier = fixed_times_bezier;
{
tm.ask_for_timeslice();
tm.start();
auto fixed_times_ig_result = experiment::fit_bezier(fixed_times_ig_bezier, curve_points);
als_time = tm.lap();
*notify << "Bezier fit a (with fixed times as i.g.), time = " << als_time << std::endl
<< "Worst residual: " << fixed_times_ig_result.first << " at t=" << fixed_times_ig_result.second << std::endl;
Geom::Path fixed_times_path;
translation.setTranslation(Geom::Point(900,300));
fixed_times_path.append(fixed_times_ig_bezier.transformed(translation));
cairo_set_source_rgba (cr, .9, .0, .0, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, fixed_times_path);
draw_text(cr, fixed_times_path.initialPoint(), "new (a) with fixed t as i.g.");
}
Geom::CubicBezier icp_bezier;
{
tm.ask_for_timeslice();
tm.start();
auto icp_result = experiment::fit_bezier_icp(icp_bezier, curve_points);
als_time = tm.lap();
*notify << "Bezier fit icp, time = " << als_time << std::endl
<< "Worst residual: " << icp_result.first << " at t=" << icp_result.second << std::endl;
Geom::Path icp_path;
translation.setTranslation(Geom::Point(600,600));
icp_path.append(icp_bezier.transformed(translation));
cairo_set_source_rgba (cr, .9, .0, .0, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, icp_path);
draw_text(cr, icp_path.initialPoint(), "icp fit");
}
Geom::CubicBezier icp_ig_bezier(icp_bezier);
{
tm.ask_for_timeslice();
tm.start();
auto icp_ig_result = experiment::fit_bezier(icp_ig_bezier, curve_points);
als_time = tm.lap();
*notify << "Bezier fit with icp i.g., time = " << als_time << std::endl
<< "Worst residual: " << icp_ig_result.first << " at t=" << icp_ig_result.second << std::endl;
Geom::Path icp_ig_path;
translation.setTranslation(Geom::Point(900,600));
icp_ig_path.append(icp_ig_bezier.transformed(translation));
cairo_set_source_rgba (cr, .9, .0, .0, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, icp_ig_path);
draw_text(cr, icp_ig_path.initialPoint(), "bezier fit with icp as i.g.");
}
std::cout << "original: " << write_svg_path(original_path) << std::endl;
Geom::CubicBezier initial_guess(
curve_points.front(), curve_points.front(),
curve_points.back(), curve_points.back()
);
{
experiment::get_initial_guess(initial_guess, curve_points);
Geom::Path initial_guess_path;
translation.setTranslation(Geom::Point(600,0));
initial_guess_path.append(initial_guess.transformed(translation));
cairo_set_source_rgba (cr, .8, .0, .8, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, initial_guess_path);
draw_text(cr, initial_guess_path.initialPoint(), "initial guess");
}
cairo_stroke(cr);
Toy::draw(cr, notify, width, height, save,timer_stream);
}