cosmos 4.0.3-java

data/ext/cosmos/ext/cosmos_io/extconf.rb

@@ -0,0 +1,13 @@
+require 'mkmf'
+
+unless $CFLAGS.gsub!(/ -O[\dsz]?/, ' -O3')
+  $CFLAGS << ' -O3'
+end
+if CONFIG['CC'] =~ /gcc/
+  $CFLAGS << ' -Wall'
+  if $DEBUG && !$CFLAGS.gsub!(/ -O[\dsz]?/, ' -O0 -ggdb')
+    $CFLAGS << ' -O0 -ggdb'
+  end
+end
+
+create_makefile 'cosmos/ext/cosmos_io'

data/ext/cosmos/ext/crc/crc.c

@@ -0,0 +1,341 @@
+/*
+# Copyright 2014 Ball Aerospace & Technologies Corp.
+# All Rights Reserved.
+#
+# This program is free software; you can modify and/or redistribute it
+# under the terms of the GNU General Public License
+# as published by the Free Software Foundation; version 3 with
+# attribution addendums as found in the LICENSE.txt
+*/
+
+#include "ruby.h"
+#include "stdio.h"
+
+VALUE mCosmos;
+VALUE cCrc;
+VALUE cCrc16;
+VALUE cCrc32;
+VALUE cCrc64;
+
+static ID id_ivar_seed = 0;
+static ID id_ivar_xor = 0;
+static ID id_ivar_reflect = 0;
+static ID id_ivar_table = 0;
+
+static const unsigned char BIT_REVERSE_TABLE[] =
+{
+  0x00, 0x80, 0x40, 0xC0, 0x20, 0xA0, 0x60, 0xE0, 0x10, 0x90, 0x50, 0xD0, 0x30, 0xB0, 0x70, 0xF0,
+  0x08, 0x88, 0x48, 0xC8, 0x28, 0xA8, 0x68, 0xE8, 0x18, 0x98, 0x58, 0xD8, 0x38, 0xB8, 0x78, 0xF8,
+  0x04, 0x84, 0x44, 0xC4, 0x24, 0xA4, 0x64, 0xE4, 0x14, 0x94, 0x54, 0xD4, 0x34, 0xB4, 0x74, 0xF4,
+  0x0C, 0x8C, 0x4C, 0xCC, 0x2C, 0xAC, 0x6C, 0xEC, 0x1C, 0x9C, 0x5C, 0xDC, 0x3C, 0xBC, 0x7C, 0xFC,
+  0x02, 0x82, 0x42, 0xC2, 0x22, 0xA2, 0x62, 0xE2, 0x12, 0x92, 0x52, 0xD2, 0x32, 0xB2, 0x72, 0xF2,
+  0x0A, 0x8A, 0x4A, 0xCA, 0x2A, 0xAA, 0x6A, 0xEA, 0x1A, 0x9A, 0x5A, 0xDA, 0x3A, 0xBA, 0x7A, 0xFA,
+  0x06, 0x86, 0x46, 0xC6, 0x26, 0xA6, 0x66, 0xE6, 0x16, 0x96, 0x56, 0xD6, 0x36, 0xB6, 0x76, 0xF6,
+  0x0E, 0x8E, 0x4E, 0xCE, 0x2E, 0xAE, 0x6E, 0xEE, 0x1E, 0x9E, 0x5E, 0xDE, 0x3E, 0xBE, 0x7E, 0xFE,
+  0x01, 0x81, 0x41, 0xC1, 0x21, 0xA1, 0x61, 0xE1, 0x11, 0x91, 0x51, 0xD1, 0x31, 0xB1, 0x71, 0xF1,
+  0x09, 0x89, 0x49, 0xC9, 0x29, 0xA9, 0x69, 0xE9, 0x19, 0x99, 0x59, 0xD9, 0x39, 0xB9, 0x79, 0xF9,
+  0x05, 0x85, 0x45, 0xC5, 0x25, 0xA5, 0x65, 0xE5, 0x15, 0x95, 0x55, 0xD5, 0x35, 0xB5, 0x75, 0xF5,
+  0x0D, 0x8D, 0x4D, 0xCD, 0x2D, 0xAD, 0x6D, 0xED, 0x1D, 0x9D, 0x5D, 0xDD, 0x3D, 0xBD, 0x7D, 0xFD,
+  0x03, 0x83, 0x43, 0xC3, 0x23, 0xA3, 0x63, 0xE3, 0x13, 0x93, 0x53, 0xD3, 0x33, 0xB3, 0x73, 0xF3,
+  0x0B, 0x8B, 0x4B, 0xCB, 0x2B, 0xAB, 0x6B, 0xEB, 0x1B, 0x9B, 0x5B, 0xDB, 0x3B, 0xBB, 0x7B, 0xFB,
+  0x07, 0x87, 0x47, 0xC7, 0x27, 0xA7, 0x67, 0xE7, 0x17, 0x97, 0x57, 0xD7, 0x37, 0xB7, 0x77, 0xF7,
+  0x0F, 0x8F, 0x4F, 0xCF, 0x2F, 0xAF, 0x6F, 0xEF, 0x1F, 0x9F, 0x5F, 0xDF, 0x3F, 0xBF, 0x7F, 0xFF
+};
+
+/*
+ * Bit Reverse an unsigned char
+ */
+static unsigned char bit_reverse_8(unsigned char value)
+{
+  return BIT_REVERSE_TABLE[value];
+}
+
+/*
+ * Bit Reverse an unsigned short
+ */
+static unsigned short bit_reverse_16(unsigned short value)
+{
+  return (BIT_REVERSE_TABLE[value & 0xff] << 8) |
+    (BIT_REVERSE_TABLE[(value >> 8) & 0xff]);
+}
+
+/*
+ * Bit Reverse an unsigned int
+ */
+static unsigned int bit_reverse_32(unsigned int value)
+{
+  return (BIT_REVERSE_TABLE[value & 0xff] << 24) |
+    (BIT_REVERSE_TABLE[(value >> 8) & 0xff] << 16) |
+    (BIT_REVERSE_TABLE[(value >> 16) & 0xff] << 8) |
+    (BIT_REVERSE_TABLE[(value >> 24) & 0xff]);
+}
+
+/*
+ * Bit Reverse an unsigned long long
+ */
+static unsigned long long bit_reverse_64(unsigned long long value)
+{
+  return ((unsigned long long)BIT_REVERSE_TABLE[value & 0x00000000000000ffULL] << 56) |
+    ((unsigned long long)BIT_REVERSE_TABLE[(value >> 8) & 0x00000000000000ffULL] << 48) |
+    ((unsigned long long)BIT_REVERSE_TABLE[(value >> 16) & 0x00000000000000ffULL] << 40) |
+    ((unsigned long long)BIT_REVERSE_TABLE[(value >> 24) & 0x00000000000000ffULL] << 32) |
+    ((unsigned long long)BIT_REVERSE_TABLE[(value >> 32) & 0x00000000000000ffULL] << 24) |
+    ((unsigned long long)BIT_REVERSE_TABLE[(value >> 40) & 0x00000000000000ffULL] << 16) |
+    ((unsigned long long)BIT_REVERSE_TABLE[(value >> 48) & 0x00000000000000ffULL] << 8) |
+    ((unsigned long long)BIT_REVERSE_TABLE[(value >> 56) & 0x00000000000000ffULL]);
+}
+
+/*
+ * Calculate a 16-bit CRC
+ */
+static VALUE crc16_calculate(int argc, VALUE* argv, VALUE self)
+{
+  volatile VALUE param_data = Qnil;
+  volatile VALUE param_seed = Qnil;
+  unsigned char* data = NULL;
+  unsigned short* table = NULL;
+  int i = 0;
+  long length = 0;
+  unsigned short crc = 0;
+
+  switch (argc)
+  {
+    case 1:
+      Check_Type(argv[0], T_STRING);
+      param_data = argv[0];
+      param_seed = rb_ivar_get(self, id_ivar_seed);
+      break;
+    case 2:
+      Check_Type(argv[0], T_STRING);
+      param_data = argv[0];
+      if (argv[1] == Qnil)
+      {
+        param_seed = rb_ivar_get(self, id_ivar_seed);
+      }
+      else
+      {
+        param_seed = argv[1];
+      }
+      break;
+    default:
+      /* Invalid number of arguments given */
+      rb_raise(rb_eArgError, "wrong number of arguments (%d for 1..2)", argc);
+      break;
+  };
+
+  crc = NUM2UINT(param_seed);
+  data = (unsigned char*) RSTRING_PTR(param_data);
+  length = RSTRING_LEN(param_data);
+  table = (unsigned short*) RSTRING_PTR(rb_ivar_get(self, id_ivar_table));
+
+  if (RTEST(rb_ivar_get(self, id_ivar_reflect)))
+  {
+    for (i = 0; i < length; i++)
+    {
+      crc = (crc << 8) ^ table[(crc >> 8) ^ bit_reverse_8(data[i])];
+    }
+
+    if (RTEST(rb_ivar_get(self, id_ivar_xor)))
+    {
+      return UINT2NUM(bit_reverse_16(crc ^ 0xFFFF));
+    }
+    else
+    {
+      return UINT2NUM(bit_reverse_16(crc));
+    }
+  }
+  else
+  {
+    for (i = 0; i < length; i++)
+    {
+      crc = (crc << 8) ^ table[(crc >> 8) ^ data[i]];
+    }
+
+    if (RTEST(rb_ivar_get(self, id_ivar_xor)))
+    {
+      return UINT2NUM(crc ^ 0xFFFF);
+    }
+    else
+    {
+      return UINT2NUM(crc);
+    }
+  }
+}
+
+/*
+ * Calculate a 32-bit CRC
+ */
+static VALUE crc32_calculate(int argc, VALUE* argv, VALUE self)
+{
+  volatile VALUE param_data = Qnil;
+  volatile VALUE param_seed = Qnil;
+  unsigned char* data = NULL;
+  unsigned int* table = NULL;
+  int i = 0;
+  long length = 0;
+  unsigned int crc = 0;
+
+  switch (argc)
+  {
+    case 1:
+      Check_Type(argv[0], T_STRING);
+      param_data = argv[0];
+      param_seed = rb_ivar_get(self, id_ivar_seed);
+      break;
+    case 2:
+      Check_Type(argv[0], T_STRING);
+      param_data = argv[0];
+      if (argv[1] == Qnil)
+      {
+        param_seed = rb_ivar_get(self, id_ivar_seed);
+      }
+      else
+      {
+        param_seed = argv[1];
+      }
+      break;
+    default:
+      /* Invalid number of arguments given */
+      rb_raise(rb_eArgError, "wrong number of arguments (%d for 1..2)", argc);
+      break;
+  };
+
+  crc = NUM2UINT(param_seed);
+  data = (unsigned char*) RSTRING_PTR(param_data);
+  length = RSTRING_LEN(param_data);
+  table = (unsigned int*) RSTRING_PTR(rb_ivar_get(self, id_ivar_table));
+
+  if (RTEST(rb_ivar_get(self, id_ivar_reflect)))
+  {
+    for (i = 0; i < length; i++)
+    {
+      crc = (crc << 8) ^ table[((crc >> 24) ^ bit_reverse_8(data[i])) & 0x000000FF];
+    }
+
+    if (RTEST(rb_ivar_get(self, id_ivar_xor)))
+    {
+      return UINT2NUM(bit_reverse_32(crc ^ 0xFFFFFFFF));
+    }
+    else
+    {
+      return UINT2NUM(bit_reverse_32(crc));
+    }
+  }
+  else
+  {
+    for (i = 0; i < length; i++)
+    {
+      crc = (crc << 8) ^ table[((crc >> 24) ^ data[i]) & 0x000000FF];
+    }
+
+    if (RTEST(rb_ivar_get(self, id_ivar_xor)))
+    {
+      return UINT2NUM(crc ^ 0xFFFFFFFF);
+    }
+    else
+    {
+      return UINT2NUM(crc);
+    }
+  }
+}
+
+/*
+ * Calculate a 64-bit CRC
+ */
+static VALUE crc64_calculate(int argc, VALUE* argv, VALUE self)
+{
+  volatile VALUE param_data = Qnil;
+  volatile VALUE param_seed = Qnil;
+  unsigned char* data = NULL;
+  unsigned long long* table = NULL;
+  int i = 0;
+  long length = 0;
+  unsigned long long crc = 0;
+
+  switch (argc)
+  {
+    case 1:
+      Check_Type(argv[0], T_STRING);
+      param_data = argv[0];
+      param_seed = rb_ivar_get(self, id_ivar_seed);
+      break;
+    case 2:
+      Check_Type(argv[0], T_STRING);
+      param_data = argv[0];
+      if (argv[1] == Qnil)
+      {
+        param_seed = rb_ivar_get(self, id_ivar_seed);
+      }
+      else
+      {
+        param_seed = argv[1];
+      }
+      break;
+    default:
+      /* Invalid number of arguments given */
+      rb_raise(rb_eArgError, "wrong number of arguments (%d for 1..2)", argc);
+      break;
+  };
+
+  crc = NUM2ULL(param_seed);
+  data = (unsigned char*) RSTRING_PTR(param_data);
+  length = RSTRING_LEN(param_data);
+  table = (unsigned long long*) RSTRING_PTR(rb_ivar_get(self, id_ivar_table));
+
+  if (RTEST(rb_ivar_get(self, id_ivar_reflect)))
+  {
+    for (i = 0; i < length; i++)
+    {
+      crc = (crc << 8) ^ table[((crc >> 56) ^ bit_reverse_8(data[i])) & 0x00000000000000FFULL];
+    }
+
+    if (RTEST(rb_ivar_get(self, id_ivar_xor)))
+    {
+      return ULL2NUM(bit_reverse_64(crc ^ 0xFFFFFFFFFFFFFFFFULL));
+    }
+    else
+    {
+      return ULL2NUM(bit_reverse_64(crc));
+    }
+  }
+  else
+  {
+    for (i = 0; i < length; i++)
+    {
+      crc = (crc << 8) ^ table[((crc >> 56) ^ data[i]) & 0x00000000000000FFULL];
+    }
+
+    if (RTEST(rb_ivar_get(self, id_ivar_xor)))
+    {
+      return ULL2NUM(crc ^ 0xFFFFFFFFFFFFFFFFULL);
+    }
+    else
+    {
+      return ULL2NUM(crc);
+    }
+  }
+}
+
+/*
+ * Initialize methods for Crc
+ */
+void Init_crc()
+{
+  id_ivar_seed = rb_intern("@seed");
+  id_ivar_xor = rb_intern("@xor");
+  id_ivar_reflect = rb_intern("@reflect");
+  id_ivar_table = rb_intern("@table");
+
+  mCosmos = rb_define_module("Cosmos");
+
+  cCrc = rb_define_class_under(mCosmos, "Crc", rb_cObject);
+
+  cCrc16 = rb_define_class_under(mCosmos, "Crc16", cCrc);
+  rb_define_method(cCrc16, "calc", crc16_calculate,  -1);
+
+  cCrc32 = rb_define_class_under(mCosmos, "Crc32", cCrc);
+  rb_define_method(cCrc32, "calc", crc32_calculate,  -1);
+
+  cCrc64 = rb_define_class_under(mCosmos, "Crc64", cCrc);
+  rb_define_method(cCrc64, "calc", crc64_calculate,  -1);
+}

data/ext/cosmos/ext/crc/extconf.rb

@@ -0,0 +1,13 @@
+require 'mkmf'
+
+unless $CFLAGS.gsub!(/ -O[\dsz]?/, ' -O3')
+  $CFLAGS << ' -O3'
+end
+if CONFIG['CC'] =~ /gcc/
+  $CFLAGS << ' -Wall'
+  if $DEBUG && !$CFLAGS.gsub!(/ -O[\dsz]?/, ' -O0 -ggdb')
+    $CFLAGS << ' -O0 -ggdb'
+  end
+end
+
+create_makefile 'cosmos/ext/crc'

data/ext/cosmos/ext/line_graph/extconf.rb

@@ -0,0 +1,13 @@
+require 'mkmf'
+
+unless $CFLAGS.gsub!(/ -O[\dsz]?/, ' -O3')
+  $CFLAGS << ' -O3'
+end
+if CONFIG['CC'] =~ /gcc/
+  $CFLAGS << ' -Wall'
+  if $DEBUG && !$CFLAGS.gsub!(/ -O[\dsz]?/, ' -O0 -ggdb')
+    $CFLAGS << ' -O0 -ggdb'
+  end
+end
+
+create_makefile 'cosmos/ext/line_graph'

data/ext/cosmos/ext/line_graph/line_graph.c

@@ -0,0 +1,501 @@
+/*
+# Copyright 2014 Ball Aerospace & Technologies Corp.
+# All Rights Reserved.
+#
+# This program is free software; you can modify and/or redistribute it
+# under the terms of the GNU General Public License
+# as published by the Free Software Foundation; version 3 with
+# attribution addendums as found in the LICENSE.txt
+*/
+
+#include "ruby.h"
+#include "stdio.h"
+
+#ifndef RFLOAT_VALUE
+  #define RFLOAT_VALUE(v) (RFLOAT(v)->value)
+#endif
+
+/* Cosmos module reference */
+VALUE mCosmos = Qnil;
+
+/* LineClip class reference */
+VALUE cLineClip = Qnil;
+
+/* Reference to LineGraph class */
+VALUE cLineGraph = Qnil;
+
+/* Reference to other needed classes/modules */
+static VALUE mQt = Qnil;
+static VALUE cQtBase = Qnil;
+static VALUE cQtWidget = Qnil;
+
+/* Method Ids */
+static ID id_method_to_f = 0;
+static ID id_method_left = 0;
+static ID id_method_right = 0;
+static ID id_method_addLineColor = 0;
+static ID id_method_addRectColorFill = 0;
+
+/* Instance Variable Ids */
+static ID id_ivar_x_max = 0;
+static ID id_ivar_x_min = 0;
+static ID id_ivar_x_scale = 0;
+static ID id_ivar_graph_left_x = 0;
+static ID id_ivar_left_y_max = 0;
+static ID id_ivar_left_y_min = 0;
+static ID id_ivar_left_y_scale = 0;
+static ID id_ivar_right_y_max = 0;
+static ID id_ivar_right_y_min = 0;
+static ID id_ivar_right_y_scale = 0;
+static ID id_ivar_graph_top_y = 0;
+static ID id_ivar_lines = 0;
+static ID id_ivar_show_lines = 0;
+static ID id_ivar_point_size = 0;
+
+/* Constant Ids */
+static ID id_LEFT = 0;
+static ID id_RIGHT = 0;
+
+/* Enumeration of line position codes */
+const long LINE_CLIP_OUTCODES_TOP = 0x1;
+const long LINE_CLIP_OUTCODES_BOTTOM = 0x2;
+const long LINE_CLIP_OUTCODES_RIGHT = 0x4;
+const long LINE_CLIP_OUTCODES_LEFT = 0x8;
+
+/*
+ * This is part of the line clipping algorithm.  This function returns
+ * a code that indicates where a point is in relation to the viewable area.
+ */
+static int cal_code (double x, double y, double xmin, double ymin, double xmax, double ymax) {
+  int code = 0;
+
+  if (y > ymax) {
+    code |= LINE_CLIP_OUTCODES_TOP;
+  } else if (y < ymin) {
+    code |= LINE_CLIP_OUTCODES_BOTTOM;
+  }
+
+  if (x > xmax) {
+    code |= LINE_CLIP_OUTCODES_RIGHT;
+  } else if (x < xmin) {
+    code |= LINE_CLIP_OUTCODES_LEFT;
+  }
+
+  return code;
+}
+
+/*
+ * Internal function to perform clipping
+ */
+static VALUE line_clip_internal(double x0, double y0, double x1, double y1, double xmin, double ymin, double xmax, double ymax, double* result_x0, double* result_y0, double* result_x1, double* result_y1, volatile VALUE* result_clipped0, volatile VALUE* result_clipped1) {
+  int code0 = 0;
+  int code1 = 0;
+  int codeout = 0;
+  volatile VALUE accept = Qfalse;
+  volatile VALUE done = Qfalse;
+  volatile VALUE clipped0 = Qfalse;
+  volatile VALUE clipped1 = Qfalse;
+  double x = 0.0;
+  double y = 0.0;
+
+  code0 = cal_code(x0, y0, xmin, ymin, xmax, ymax);
+  code1 = cal_code(x1, y1, xmin, ymin, xmax, ymax);
+
+  while (1) {
+    if ((code0 | code1) == 0) {
+      /* Both points are within the viewable area.  The entire line can be
+       * graphed. */
+      accept = Qtrue;
+      done = Qtrue;
+    } else if ((code0 & code1) != 0) {
+      /* The entire line is outside of the viewable area.  No part of the
+       * line can be graphed. */
+      accept = Qfalse;
+      clipped0 = Qtrue;
+      clipped1 = Qtrue;
+      done = Qtrue;
+    } else {
+      /* Part of the line is inside the viewable area.  Figure out which part
+       * of the line can be drawn. */
+      x = 0.0;
+      y = 0.0;
+
+      if (code0 != 0) {
+        codeout = code0;
+        clipped0 = Qtrue;
+      } else {
+        codeout = code1;
+        clipped1 = Qtrue;
+      }
+
+      if ((codeout & LINE_CLIP_OUTCODES_TOP) != 0) {
+        x = x0 + (x1 - x0) * (ymax - y0) / (y1 - y0);
+        y = ymax;
+      } else if ((codeout & LINE_CLIP_OUTCODES_BOTTOM) != 0) {
+        x = x0 + (x1 - x0) * (ymin - y0) / (y1 - y0);
+        y = ymin;
+      } else if ((codeout & LINE_CLIP_OUTCODES_RIGHT) != 0) {
+        y = y0 + (y1 - y0) * (xmax - x0) / (x1 - x0);
+        x = xmax;
+      } else {
+        y = y0 + (y1 - y0) * (xmin - x0) / (x1 - x0);
+        x = xmin;
+      }
+
+      if (codeout == code0) {
+        x0 = x;
+        y0 = y;
+        code0 = cal_code(x0, y0, xmin, ymin, xmax, ymax);
+      } else {
+        x1 = x;
+        y1 = y;
+        code1 = cal_code(x1, y1, xmin, ymin, xmax, ymax);
+      }
+    }
+
+    if (done == Qtrue) {
+      break;
+    }
+  }
+
+  *result_x0 = x0;
+  *result_y0 = y0;
+  *result_x1 = x1;
+  *result_y1 = y1;
+  *result_clipped0 = clipped0;
+  *result_clipped1 = clipped1;
+
+  return accept;
+}
+
+/*
+ * This is a line-clipping algorithm.  It takes two points and a viewable
+ * area.  It returns the part of the line that is within the viewable area.
+ * If no part of the line is viewable, it returns nil
+ */
+static VALUE line_clip(VALUE self, VALUE x0, VALUE y0, VALUE x1, VALUE y1, VALUE xmin, VALUE ymin, VALUE xmax, VALUE ymax) {
+  volatile VALUE result = Qnil;
+  volatile VALUE result_clipped0 = Qnil;
+  volatile VALUE result_clipped1 = Qnil;
+  volatile VALUE return_value = Qnil;
+  double double_x0 = 0.0;
+  double double_y0 = 0.0;
+  double double_x1 = 0.0;
+  double double_y1 = 0.0;
+  double double_xmin = 0.0;
+  double double_ymin = 0.0;
+  double double_xmax = 0.0;
+  double double_ymax = 0.0;
+  double result_x0 = 0.0;
+  double result_y0 = 0.0;
+  double result_x1 = 0.0;
+  double result_y1 = 0.0;
+
+  double_x0 = RFLOAT_VALUE(rb_funcall(x0, id_method_to_f, 0));
+  double_y0 = RFLOAT_VALUE(rb_funcall(y0, id_method_to_f, 0));
+  double_x1 = RFLOAT_VALUE(rb_funcall(x1, id_method_to_f, 0));
+  double_y1 = RFLOAT_VALUE(rb_funcall(y1, id_method_to_f, 0));
+  double_xmin = RFLOAT_VALUE(rb_funcall(xmin, id_method_to_f, 0));
+  double_ymin = RFLOAT_VALUE(rb_funcall(ymin, id_method_to_f, 0));
+  double_xmax = RFLOAT_VALUE(rb_funcall(xmax, id_method_to_f, 0));
+  double_ymax = RFLOAT_VALUE(rb_funcall(ymax, id_method_to_f, 0));
+
+  result = line_clip_internal(double_x0, double_y0, double_x1, double_y1, double_xmin, double_ymin, double_xmax, double_ymax, &result_x0, &result_y0, &result_x1, &result_y1, &result_clipped0, &result_clipped1);
+
+  if (result == Qtrue)
+  {
+    return_value = rb_ary_new2(6);
+    rb_ary_push(return_value, rb_float_new(result_x0));
+    rb_ary_push(return_value, rb_float_new(result_y0));
+    rb_ary_push(return_value, rb_float_new(result_x1));
+    rb_ary_push(return_value, rb_float_new(result_y1));
+    rb_ary_push(return_value, result_clipped0);
+    rb_ary_push(return_value, result_clipped1);
+  }
+
+  return return_value;
+}
+
+/*
+ * rounds a double to the nearest integer
+ */
+static long round_nearest (double value) {
+  return ((long) (value + 0.5));
+}
+
+/*
+ * scale_value_to_graph_y internal
+ */
+static long scale_value_to_graph_y_internal (double y, double y_max, double y_scale, long graph_top_y) {
+  return (round_nearest((y_max - y) * y_scale) + graph_top_y);
+}
+
+/*
+ * This function converts a y value to a y coordinate on the graph
+ */
+static VALUE scale_value_to_graph_y(int argc, VALUE* argv, VALUE self) {
+  volatile VALUE y = Qnil;
+  ID id_axis = 0;
+  long long_graph_top_y = 0;
+  double double_y = 0.0;
+  double double_y_max = 0.0;
+  double double_y_scale = 0.0;
+
+  switch (argc) {
+    case 1:
+      y = argv[0];
+      id_axis = id_LEFT;
+      break;
+    case 2:
+      y = argv[0];
+      id_axis = SYM2ID(argv[1]);
+      break;
+    default:
+      /* Invalid number of arguments given */
+      rb_raise(rb_eArgError, "wrong number of arguments (%d for 1..2)", argc);
+      break;
+  };
+
+  long_graph_top_y = FIX2INT(rb_ivar_get(self, id_ivar_graph_top_y));
+  double_y = RFLOAT_VALUE(rb_funcall(y, id_method_to_f, 0));
+
+  if (id_axis == id_LEFT) {
+    double_y_max = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_left_y_max),   id_method_to_f, 0));
+    double_y_scale = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_left_y_scale), id_method_to_f, 0));
+  } else /* id_axis == id_RIGHT */ {
+    double_y_max = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_right_y_max),   id_method_to_f, 0));
+    double_y_scale = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_right_y_scale), id_method_to_f, 0));
+  }
+
+  return INT2FIX(scale_value_to_graph_y_internal(double_y, double_y_max, double_y_scale, long_graph_top_y));
+}
+
+/*
+ * Internal function converts an x value to an x coordinate on the graph
+ */
+static long scale_value_to_graph_x_internal (double x, double x_min, double x_scale, long graph_left_x) {
+  return (round_nearest((x - x_min) * x_scale) + graph_left_x);
+}
+
+/*
+ * This function converts an x value to an x coordinate on the graph
+ */
+static VALUE scale_value_to_graph_x(VALUE self, VALUE x) {
+  long long_graph_left_x = 0;
+  double double_x = 0.0;
+  double double_x_min = 0.0;
+  double double_x_scale = 0.0;
+
+  long_graph_left_x = FIX2INT(rb_ivar_get(self, id_ivar_graph_left_x));
+  double_x = RFLOAT_VALUE(rb_funcall(x, id_method_to_f, 0));
+  double_x_min = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_x_min), id_method_to_f, 0));
+  double_x_scale = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_x_scale), id_method_to_f, 0));
+
+  return INT2FIX(scale_value_to_graph_x_internal(double_x, double_x_min, double_x_scale, long_graph_left_x));
+}
+
+/*
+ * Internal version to draw a line
+ */
+static void draw_line_internal(VALUE dc, double x1, double y1, double x2, double y2, double x_min, double y_min, double x_max, double y_max, double x_scale, double y_scale, long graph_left_x, long graph_top_y, ID id_axis, VALUE show_line, VALUE point_size, VALUE color) {
+  volatile VALUE result = Qnil;
+  volatile VALUE clipped1 = Qnil;
+  volatile VALUE clipped2 = Qnil;
+  long x1_scaled = 0;
+  long y1_scaled = 0;
+  long x2_scaled = 0;
+  long y2_scaled = 0;
+  double clipped_x1 = 0.0;
+  double clipped_y1 = 0.0;
+  double clipped_x2 = 0.0;
+  double clipped_y2 = 0.0;
+
+  /* Calculate potentially clipped version of line */
+  result = line_clip_internal(x1, y1, x2, y2, x_min, y_min, x_max, y_max, &clipped_x1, &clipped_y1, &clipped_x2, &clipped_y2, &clipped1, &clipped2);
+
+  if (result == Qtrue) /* Line is visible so draw it */ {
+    /* Scale to graph coordinates */
+    x1_scaled = scale_value_to_graph_x_internal(clipped_x1, x_min, x_scale, graph_left_x);
+    y1_scaled = scale_value_to_graph_y_internal(clipped_y1, y_max, y_scale, graph_top_y);
+    x2_scaled = scale_value_to_graph_x_internal(clipped_x2, x_min, x_scale, graph_left_x);
+    y2_scaled = scale_value_to_graph_y_internal(clipped_y2, y_max, y_scale, graph_top_y);
+
+    /* Draw the line */
+    if (RTEST(show_line)) {
+      rb_funcall(dc, id_method_addLineColor, 5, INT2FIX(x1_scaled), INT2FIX(y1_scaled), INT2FIX(x2_scaled), INT2FIX(y2_scaled), color);
+    }
+
+    /* Draw point at line */
+    if (FIX2INT(point_size) > 0) {
+      /* Only show point if second point wasn't clipped */
+      if (!RTEST(clipped2)) {
+        rb_funcall(dc, id_method_addRectColorFill, 5, INT2FIX(x2_scaled - 2), INT2FIX(y2_scaled - 2), point_size, point_size, color);
+      }
+    }
+  }
+}
+
+/*
+ * Draws a line between two points that is clipped to fit the visible graph if necessary
+ */
+static VALUE draw_line(VALUE self, VALUE dc, VALUE x1, VALUE y1, VALUE x2, VALUE y2, VALUE show_line, VALUE point_size, VALUE axis, VALUE color) {
+  long long_graph_left_x = 0;
+  long long_graph_top_y = 0;
+  ID id_axis = 0;
+  double double_x1 = 0.0;
+  double double_y1 = 0.0;
+  double double_x2 = 0.0;
+  double double_y2 = 0.0;
+  double double_x_min = 0.0;
+  double double_y_min = 0.0;
+  double double_x_max = 0.0;
+  double double_y_max = 0.0;
+  double double_x_scale = 0.0;
+  double double_y_scale = 0.0;
+
+  id_axis = SYM2ID(axis);
+  double_x_max = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_x_max), id_method_to_f, 0));
+  double_x_min = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_x_min), id_method_to_f, 0));
+  double_x1 = RFLOAT_VALUE(rb_funcall(x1, id_method_to_f, 0));
+  double_y1 = RFLOAT_VALUE(rb_funcall(y1, id_method_to_f, 0));
+  double_x2 = RFLOAT_VALUE(rb_funcall(x2, id_method_to_f, 0));
+  double_y2 = RFLOAT_VALUE(rb_funcall(y2, id_method_to_f, 0));
+  double_x_scale = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_x_scale), id_method_to_f, 0));
+  long_graph_left_x = FIX2INT(rb_ivar_get(self, id_ivar_graph_left_x));
+  long_graph_top_y = FIX2INT(rb_ivar_get(self, id_ivar_graph_top_y));
+
+  if (id_axis == id_LEFT) {
+    double_y_max = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_left_y_max), id_method_to_f, 0));
+    double_y_min = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_left_y_min), id_method_to_f, 0));
+    double_y_scale = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_left_y_scale), id_method_to_f, 0));
+  } else /* id_axis == id_RIGHT */ {
+    double_y_max = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_right_y_max), id_method_to_f, 0));
+    double_y_min = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_right_y_min), id_method_to_f, 0));
+    double_y_scale = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_right_y_scale), id_method_to_f, 0));
+  }
+
+  draw_line_internal(dc, double_x1, double_y1, double_x2, double_y2, double_x_min, double_y_min, double_x_max, double_y_max, double_x_scale, double_y_scale, long_graph_left_x, long_graph_top_y, id_axis, show_line, point_size, color);
+
+  return Qnil;
+}
+
+/*
+ * Draws all lines for the given axis
+ */
+static VALUE draw_lines (VALUE self, VALUE dc, VALUE axis) {
+  long long_graph_left_x = 0;
+  long long_graph_top_y = 0;
+  long num_lines = 0;
+  long line_index = 0;
+  long line_length = 0;
+  long point_index = 0;
+  ID id_axis = 0;
+  volatile VALUE lines = Qnil;
+  volatile VALUE line = Qnil;
+  volatile VALUE x_values = Qnil;
+  volatile VALUE y_values = Qnil;
+  volatile VALUE color = Qnil;
+  volatile VALUE show_lines = Qnil;
+  volatile VALUE point_size = Qnil;
+  double double_x1 = 0.0;
+  double double_y1 = 0.0;
+  double double_x2 = 0.0;
+  double double_y2 = 0.0;
+  double double_x_min = 0.0;
+  double double_y_min = 0.0;
+  double double_x_max = 0.0;
+  double double_y_max = 0.0;
+  double double_x_scale = 0.0;
+  double double_y_scale = 0.0;
+
+  id_axis = SYM2ID(axis);
+  double_x_max = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_x_max), id_method_to_f, 0));
+  double_x_min = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_x_min), id_method_to_f, 0));
+  double_x_scale = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_x_scale), id_method_to_f, 0));
+  long_graph_left_x = FIX2INT(rb_ivar_get(self, id_ivar_graph_left_x));
+  long_graph_top_y = FIX2INT(rb_ivar_get(self, id_ivar_graph_top_y));
+
+  if (id_axis == id_LEFT) {
+    lines = rb_funcall(rb_ivar_get(self, id_ivar_lines), id_method_left, 0);
+    double_y_max = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_left_y_max), id_method_to_f, 0));
+    double_y_min = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_left_y_min), id_method_to_f, 0));
+    double_y_scale = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_left_y_scale), id_method_to_f, 0));
+  } else {
+    lines = rb_funcall(rb_ivar_get(self, id_ivar_lines), id_method_right, 0);
+    double_y_max = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_right_y_max), id_method_to_f, 0));
+    double_y_min = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_right_y_min), id_method_to_f, 0));
+    double_y_scale = RFLOAT_VALUE(rb_funcall(rb_ivar_get(self, id_ivar_right_y_scale), id_method_to_f, 0));
+  }
+
+  show_lines = rb_ivar_get(self, id_ivar_show_lines);
+  point_size = rb_ivar_get(self, id_ivar_point_size);
+
+  num_lines = RARRAY_LEN(lines);
+  for (line_index = 0; line_index < num_lines; line_index++) {
+    line = rb_ary_entry(lines, line_index);
+    x_values = rb_ary_entry(line, 0);
+    y_values = rb_ary_entry(line, 1);
+    color = rb_ary_entry(line, 6);
+
+    /* Get the first point of the line */
+    double_x1 = RFLOAT_VALUE(rb_ary_entry(x_values, 0));
+    double_y1 = RFLOAT_VALUE(rb_ary_entry(y_values, 0));
+
+    /* Loop over each data point of the line */
+    line_length = RARRAY_LEN(x_values);
+    for (point_index = 0; point_index < line_length; point_index++) {
+      double_x2 = RFLOAT_VALUE(rb_ary_entry(x_values, point_index));
+      double_y2 = RFLOAT_VALUE(rb_ary_entry(y_values, point_index));
+
+      draw_line_internal(dc, double_x1, double_y1, double_x2, double_y2, double_x_min, double_y_min, double_x_max, double_y_max, double_x_scale, double_y_scale, long_graph_left_x, long_graph_top_y, id_axis, show_lines, point_size, color);
+
+      double_x1 = double_x2;
+      double_y1 = double_y2;
+    }
+  }
+
+  return Qnil;
+}
+
+/*
+ * Initialize methods for C LineGraph
+ */
+void Init_line_graph (void)
+{
+  rb_require("cosmos");
+  mQt = rb_define_module("Qt");
+  cQtBase = rb_define_class_under(mQt, "Base", rb_cObject);
+  cQtWidget = rb_define_class_under(mQt, "Widget", cQtBase);
+
+  id_method_left = rb_intern("left");
+  id_method_right = rb_intern("right");
+  id_method_to_f = rb_intern("to_f");
+  id_method_addLineColor = rb_intern("addLineColor");
+  id_method_addRectColorFill = rb_intern("addRectColorFill");
+  id_ivar_x_max = rb_intern("@x_max");
+  id_ivar_x_min = rb_intern("@x_min");
+  id_ivar_x_scale = rb_intern("@x_scale");
+  id_ivar_graph_left_x = rb_intern("@graph_left_x");
+  id_ivar_graph_top_y = rb_intern("@graph_top_y");
+  id_ivar_left_y_max = rb_intern("@left_y_max");
+  id_ivar_left_y_min = rb_intern("@left_y_min");
+  id_ivar_left_y_scale = rb_intern("@left_y_scale");
+  id_ivar_right_y_max = rb_intern("@right_y_max");
+  id_ivar_right_y_min = rb_intern("@right_y_min");
+  id_ivar_right_y_scale = rb_intern("@right_y_scale");
+  id_ivar_lines = rb_intern("@lines");
+  id_ivar_show_lines = rb_intern("@show_lines");
+  id_ivar_point_size = rb_intern("@point_size");
+  id_LEFT = rb_intern("LEFT");
+  id_RIGHT = rb_intern("RIGHT");
+
+  mCosmos = rb_define_module("Cosmos");
+  cLineClip = rb_define_class_under(mCosmos, "LineClip", rb_cObject);
+  rb_define_singleton_method(cLineClip, "line_clip", line_clip, 8);
+
+  cLineGraph = rb_define_class_under(mCosmos, "LineGraph", cQtWidget);
+  rb_define_method(cLineGraph, "scale_value_to_graph_x", scale_value_to_graph_x,  1);
+  rb_define_method(cLineGraph, "scale_value_to_graph_y", scale_value_to_graph_y, -1);
+  rb_define_method(cLineGraph, "draw_line", draw_line, 9);
+  rb_define_method(cLineGraph, "draw_lines", draw_lines, 2);
+}