#include <stdlib.h>
#include <stdio.h>
+#include <string.h>
#include <memory.h>
#include <math.h>
+#include <ctype.h>
#include "drawer.h"
static char* getToken(const char**p)
{
const char*start;
char*result;
- while(**p && strchr(" ,\t\n\r", **p)) {
+ while(**p && strchr(" ,()\t\n\r", **p)) {
(*p)++;
}
start = *p;
- while(**p && !strchr(" ,\t\n\r", **p)) {
+
+ /*
+ SVF pathdata can exclude whitespace after L and M commands.
+ Ref: http://www.w3.org/TR/SVG11/paths.html#PathDataGeneralInformation
+ This allows us to use svg files output from gnuplot.
+ Also checks for relative MoveTo and LineTo (m and l).
+ 051106 Magnus Lundin, lundin@mlu.mine.nu
+ */
+ if (strchr("LMlm", **p) && (isdigit(*(*p+1))||strchr("+-", *(*p+1)))) {
+ (*p)++;
+ }
+ else while(**p && !strchr(" ,()\t\n\r", **p)) {
(*p)++;
}
result = malloc((*p)-start+1);
draw->pos = *to;
}
+/* convenience routine */
+static void draw_conicTo2(drawer_t*draw, double x1, double y1, double x2, double y2)
+{
+ FPOINT c1,c2;
+ c1.x = x1;
+ c1.y = y1;
+ c2.x = x2;
+ c2.y = y2;
+ draw_conicTo(draw, &c1, &c2);
+}
+/* convenience routine */
+static void draw_moveTo2(drawer_t*draw, double x, double y)
+{
+ FPOINT c;
+ c.x = x; c.y = y;
+ draw->moveTo(draw, &c);
+}
+/* convenience routine */
+static void draw_lineTo2(drawer_t*draw, double x, double y)
+{
+ FPOINT c;
+ c.x = x; c.y = y;
+ draw->lineTo(draw, &c);
+}
+
+
void draw_string(drawer_t*draw, const char*string)
{
const char*p = string;
char*token = getToken(&p);
if(!token || !*token)
break;
- if(!strncmp(token, "moveTo", 6)) {
+ if(!strncmp(token, "moveTo", 6) ||
+ !strncmp(token, "M", 1) //svg
+ ) {
FPOINT to;
- to.x = atoi(getToken(&p));
- to.y = atoi(getToken(&p));
+ to.x = atof(getToken(&p));
+ to.y = atof(getToken(&p));
draw->moveTo(draw, &to);
}
- else if(!strncmp(token, "lineTo", 6)) {
+ else if(!strncmp(token, "lineTo", 6) ||
+ !strncmp(token, "L", 1) //svg
+ ) {
FPOINT to;
- to.x = atoi(getToken(&p));
- to.y = atoi(getToken(&p));
+ to.x = atof(getToken(&p));
+ to.y = atof(getToken(&p));
draw->lineTo(draw, &to);
}
else if(!strncmp(token, "curveTo", 7) || !strncmp(token, "splineTo", 8)) {
FPOINT mid,to;
- mid.x = atoi(getToken(&p));
- mid.y = atoi(getToken(&p));
- to.x = atoi(getToken(&p));
- to.y = atoi(getToken(&p));
+ mid.x = atof(getToken(&p));
+ mid.y = atof(getToken(&p));
+ to.x = atof(getToken(&p));
+ to.y = atof(getToken(&p));
draw->splineTo(draw, &mid, &to);
}
- else if(!strncmp(token, "cubicTo", 5)) {
+ else if(!strncmp(token, "conicTo", 5)) {
+ FPOINT mid,to;
+ mid.x = atof(getToken(&p));
+ mid.y = atof(getToken(&p));
+ to.x = atof(getToken(&p));
+ to.y = atof(getToken(&p));
+ draw_conicTo(draw, &mid, &to);
+ }
+ else if(!strncmp(token, "circle", 6)) {
+ int mx,my,r;
+ double r2;
+ mx = atof(getToken(&p));
+ my = atof(getToken(&p));
+ r = atof(getToken(&p));
+ r2 = 0.70710678118654757*r;
+ draw_moveTo2(draw, mx, my-r);
+ draw_conicTo2(draw, mx+r2, my-r2, mx+r, my);
+ draw_conicTo2(draw, mx+r2, my+r2, mx, my+r);
+ draw_conicTo2(draw, mx-r2, my+r2, mx-r, my);
+ draw_conicTo2(draw, mx-r2, my-r2, mx, my-r);
+ }
+ else if(!strncmp(token, "box", 3)) {
+ int x1,y1,x2,y2;
+ x1 = atof(getToken(&p));
+ y1 = atof(getToken(&p));
+ x2 = atof(getToken(&p));
+ y2 = atof(getToken(&p));
+ draw_moveTo2(draw, x1, y1);
+ draw_lineTo2(draw, x1, y2);
+ draw_lineTo2(draw, x2, y2);
+ draw_lineTo2(draw, x2, y1);
+ draw_lineTo2(draw, x1, y1);
+ }
+ else if(!strncmp(token, "cubicTo", 5) ||
+ !strncmp(token, "C", 1) //svg
+ ) {
FPOINT mid1,mid2,to;
- mid1.x = atoi(getToken(&p));
- mid1.y = atoi(getToken(&p));
- mid2.x = atoi(getToken(&p));
- mid2.y = atoi(getToken(&p));
- to.x = atoi(getToken(&p));
- to.y = atoi(getToken(&p));
+ mid1.x = atof(getToken(&p));
+ mid1.y = atof(getToken(&p));
+ mid2.x = atof(getToken(&p));
+ mid2.y = atof(getToken(&p));
+ to.x = atof(getToken(&p));
+ to.y = atof(getToken(&p));
draw_cubicTo(draw, &mid1, &mid2, &to);
}
- else fprintf(stderr, "drawer: Warning: unknown primitive '%s'", token);
+ else if(!strncmp(token, "z", 1) //svg
+ ) {
+ // ignore
+ }
+ else
+ fprintf(stderr, "drawer: Warning: unknown primitive '%s'\n", token);
free(token);
}
char left = 0,recurse=0;
int t;
int probes = 15;
+ double dx,dy;
/* create simple approximation: a qspline which run's through the
qspline point at 0.5 */
test.control.y += test.end.y;
}
+#define PROBES
+#ifdef PROBES
/* measure the spline's accurancy, by taking a number of probes */
-
for(t=0;t<probes;t++) {
struct SPLINEPOINT qr1,qr2,cr1,cr2;
double pos = 0.5/(probes*2)*(t*2+1);
recurse=1;break;
}
}
+#else // quadratic error: *much* faster!
+
+ /* convert control point representation to
+ d*x^3 + c*x^2 + b*x + a */
+ dx= s->end.x - s->control2.x*3 + s->control1.x*3 - s->start.x;
+ dy= s->end.y - s->control2.y*3 + s->control1.y*3 - s->start.y;
+
+ /* we need to do this for the subspline between [start,end], not [0,1]
+ as a transformation of t->a*t+b does nothing to highest coefficient
+ of the spline except multiply it with a^3, we just need to modify
+ d here. */
+ {double m = end-start;
+ dx*=m*m*m;
+ dy*=m*m*m;
+ }
+
+ /* use the integral over (f(x)-g(x))^2 between 0 and 1
+ to measure the approximation quality.
+ (it boils down to const*d^2)
+ */
+ recurse = (dx*dx + dy*dy > quality2);
+#endif
if(recurse && istep>1 && size-level > num) {
istep >>= 1;
{
struct qspline q[128];
struct cspline c;
- double quality = 80;
- double maxerror = (500-(quality*5)>1?500-(quality*5):1)/20.0;
-
- int num = approximate3(&c, q, 128, maxerror*maxerror);
- int t;
+ //double quality = 80;
+ double maxerror = 1;//(500-(quality*5)>1?500-(quality*5):1)/20.0;
+ int t,num;
c.start.x = draw->pos.x;
c.start.y = draw->pos.y;
c.control2.y = control2->y;
c.end.x = to->x;
c.end.y = to->y;
+
+ num = approximate3(&c, q, 128, maxerror*maxerror);
for(t=0;t<num;t++) {
FPOINT mid;