foo.c

#include "SDL.h"
#include <cairo.h>

#include <math.h>

int fps = 30;
int w = 320;
int h = 240;

typedef float vec[4];
typedef float *vec_p;
typedef float mat[4 * 4];
typedef float *mat_p;

void mat_identity(mat m) {
  int i;

  for (i = 0; i < 16; i = i + 1) {
    m[i] = 0;
  }
  m[ 0] = 1; m[ 5] = 1; m[10] = 1; m[15] = 1;
}

void mat_load(mat m, mat n) {
  int i;

  for (i = 0; i < 16; i = i + 1) {
    m[i] = n[i];
  }
}

void mat_mul(mat m, mat n) {
  mat o;
  int i, j;

  for (j = 0; j < 4; j = j + 1) {
    for (i = 0; i < 4; i = i + 1) {
      o[4*j+i] = m[4*j+0] * n[4*0+i] + m[4*j+1] * n[4*1+i]
        + m[4*j+2] * n[4*2+i] + m[4*j+3] * n[4*3+i];
    }
  }
  for (i = 0; i < 16; i = i + 1) {
    m[i] = o[i];
  }
}

void vec_load(vec v, vec w) {
  int i;

  for (i = 0; i < 4; i = i + 1) {
    v[i] = w[i];
  }
}

void vec_mat_mul(vec v, mat m) {
  float x, y, z, w;
  x = v[0] * m[ 0] + v[1] * m[ 1] + v[2] * m[ 2] + v[3] * m[ 3];   
  y = v[0] * m[ 4] + v[1] * m[ 5] + v[2] * m[ 6] + v[3] * m[ 7];
  z = v[0] * m[ 8] + v[1] * m[ 9] + v[2] * m[10] + v[3] * m[11];
  w = v[0] * m[12] + v[1] * m[13] + v[2] * m[14] + v[3] * m[15];
  v[0] = x; v[1] = y; v[2] = z; v[3] = w;
}

int main(int argc, char **argv) {
  SDL_Surface *sdl_surface;
  cairo_t *cr;

  Uint8 *keystate;
  Uint32 next_frame;

  mat matrix;

  vec verts[8] = {
    -1/2.0, -1/2.0, -1/2.0, 1,
     1/2.0, -1/2.0, -1/2.0, 1,
    -1/2.0,  1/2.0, -1/2.0, 1,
     1/2.0,  1/2.0, -1/2.0, 1,
    -1/2.0, -1/2.0,  1/2.0, 1,
     1/2.0, -1/2.0,  1/2.0, 1,
    -1/2.0,  1/2.0,  1/2.0, 1,
     1/2.0,  1/2.0,  1/2.0, 1
    };

  vec_p lines[] = {
    verts[0], verts[1],
    verts[1], verts[3],
    verts[3], verts[2],
    verts[2], verts[0],

    verts[4], verts[5],
    verts[5], verts[7],
    verts[7], verts[6],
    verts[6], verts[4],

    verts[0], verts[4],
    verts[1], verts[5],
    verts[2], verts[6],
    verts[3], verts[7],
    0};

  int running;
  float angle;

  SDL_Init(SDL_INIT_VIDEO | SDL_INIT_TIMER);
  SDL_ShowCursor(0);
  SDL_SetVideoMode(w, h, 32, 0);
  sdl_surface = SDL_GetVideoSurface();

  { /* Initialize Canvas */
    cairo_surface_t *surface;
    surface = cairo_image_surface_create_for_data(
      sdl_surface->pixels,
      CAIRO_FORMAT_RGB24,
      sdl_surface->w,
      sdl_surface->h,
      sdl_surface->pitch
      );
    cr = cairo_create(surface);
    cairo_surface_destroy(surface);
  }

  // Cartesian
  cairo_translate(cr, w/2.0, h/2.0);
  cairo_scale(cr, 1, -1);

  // fixed scale
  cairo_scale(cr, h, h);

  /* Initialize Delay */
  next_frame = 1024.0 / fps;

  /* Game Logic */
  running = 1;
  {
    mat_p m = matrix;
    float a = M_PI / 4.0;
    mat_identity(m);
    m[11] = 4;
    m[ 5] =  cos(a); m[ 6] = -sin(a);
    m[ 9] =  sin(a); m[10] =  cos(a);
  }
  angle = 0;

  SDL_LockSurface(sdl_surface);
  while (running) {

    /* Render Frame */
    cairo_set_operator(cr, CAIRO_OPERATOR_CLEAR);
    cairo_paint(cr);
    cairo_set_operator(cr, CAIRO_OPERATOR_OVER);

    {
      mat m, n; vec v;
      float x, y, a = angle;
      int i;

      mat_identity(n);
      n[0] =  cos(a); n[1] = -sin(a);
      n[4] =  sin(a); n[5] =  cos(a);

      mat_load(m, matrix);
      mat_mul(m, n);

      i = 0;
      while (lines[i]) {
        vec_load(v, lines[i+0]);
        vec_mat_mul(v, m);
        x = v[0] / v[2]; y = v[1] / v[2];
        cairo_move_to(cr, x, y);

        vec_load(v, lines[i+1]);
        vec_mat_mul(v, m);
        x = v[0] / v[2]; y = v[1] / v[2];
        cairo_line_to(cr, x, y);

        i = i + 2;
      }

      cairo_set_source_rgb(cr, 1.0, 1.0, 1.0);
      cairo_set_line_width(cr, 1/64.0);
      cairo_set_line_cap(cr, CAIRO_LINE_CAP_ROUND);
      cairo_stroke(cr);
    }

    /* Update Display */
    SDL_UnlockSurface(sdl_surface);
    SDL_Flip(sdl_surface);
    SDL_LockSurface(sdl_surface);

    /* Delay */
    {
      Uint32 now;
      now = SDL_GetTicks();
      if (now < next_frame) {
        SDL_Delay(next_frame - now);
      }
      next_frame = next_frame + 1024.0 / fps;
    }

    /* Game Logic */
    SDL_PumpEvents();
    keystate = SDL_GetKeyState(NULL);
    if (keystate[SDLK_q]) {
      running = 0;
    }

    angle = angle + M_PI / 2.0 / fps;

  }
  SDL_UnlockSurface(sdl_surface);

  cairo_destroy(cr);
  SDL_Quit();

  return 0;
}