#include "common.h"

void* alloc(size_t size) {

  void *p = calloc(1, size);

  if (!p) DIE("cannot allocate %llu bytes of memory", (unsigned long long)size);

  return p;

}

void* alloc2d(size_t cnt, size_t size) {

  if (cnt < 0) return NULL;

  size_t ps = cnt*sizeof(void*);

  char **p = (char**)alloc(ps + cnt*size);

  char *d = ((char*)p) + ps;

  if (size) for(size_t i = 0; i < cnt; ++i) p[i] = d + i*size;

  return p;

}

char* vstrprintf(const char *fmt, va_list args) {

  va_list copy;

  va_copy(copy, args);

  int size = vsnprintf(NULL, 0, fmt, copy);

  va_end(copy);

  if (size < 0) size = 0;

  ++size;

  char *buf = alloc(size+1);

  vsnprintf(buf, size, fmt, args);

  return buf;

}

char* strprintf(const char *format, ...) {

  va_list args;

  va_start(args,format);

  char* s = vstrprintf(format, args);

  va_end(args);

  return s;

}

char* unquote(char *s, char prefix, char suffix) {

  char *s0 = NULL, *s1 = s;

  for(char *c = s; *c; ++c) {

    if (isspace(*c)) continue;

    if (!s0) s0 = c;

    s1 = c+1;

  }

  if (!s0) s0 = s;

  if (s0 < s1 && *s0 == prefix) ++s0;

  if (s0 < s1 && *(s1-1) == suffix) --s1;

  memmove(s, s0, s1 - s0);

  s[s1-s0] = 0;

  return s;

}

char* readline(FILE *f) {

  size_t len = 0, blen = 1024;

  char *buf = alloc(blen+1);

  while(1) {

    int ch = fgetc(f);

    if (ch <= 0 || ch == '\n')

      break;

    if (len+1 >= blen) {

      size_t bl = blen + blen/8;

      buf = realloc(buf, bl);

      if (!buf) DIE("cannot re-allocate %llu bytes of memory", (unsigned long long)bl);

      memset(buf + blen, 0, bl - blen);

      blen = bl;

    }

    buf[len++] = ch;

  }

  buf[len++] = 0;

  buf = realloc(buf, len);

  if (!buf) DIE("cannot re-allocate %llu bytes of memory", (unsigned long long)len);

  return buf;

}

void cubicto(double x0, double y0, double x1, double y1, double x2, double y2, double x3, double y3, int levels) {

  if (--levels <= 0) { lineTo(x3, y3); return; }

  double xx0 = (x0 + x1)/2, yy0 = (y0 + y1)/2;

  double xx1 = (x1 + x2)/2, yy1 = (y1 + y2)/2;

  double xx2 = (x2 + x3)/2, yy2 = (y2 + y3)/2;

  double xxx0 = (xx0 + xx1)/2, yyy0 = (yy0 + yy1)/2;

  double xxx1 = (xx1 + xx2)/2, yyy1 = (yy1 + yy2)/2;

  double xxxx = (xxx0 + xxx1)/2, yyyy = (yyy0 + yyy1)/2;

  cubicto(x0, y0, xx0, yy0, xxx0, yyy0, xxxx, yyyy, levels);

  cubicto(xxxx, yyyy, xxx1, yyy1, xx2, yy2, x3, y3, levels);

}