numerix 1.0.0

data/ext/numerix/matrix4x4.h

@@ -0,0 +1,90 @@
+#ifndef NUMERIX_MATRIX4X4_H
+#define NUMERIX_MATRIX4X4_H 1
+
+#include "common.h"
+#include "matrix_base.h"
+#include "plane.h"
+#include "quaternion.h"
+
+#define MATRIX4X4() \
+    Matrix4x4 *m;   \
+    Data_Get_Struct(self, Matrix4x4, m)
+
+typedef struct CanonicalBasis {
+    Vector3 row0;
+    Vector3 row1;
+    Vector3 row2;
+} CanonicalBasis;
+
+typedef struct VectorBasis {
+    Vector3 *element0;
+    Vector3 *element1;
+    Vector3 *element2;
+} VectorBasis;
+
+void Init_matrix4x4(VALUE outer);
+static VALUE rb_matrix4x4_allocate(VALUE klass);
+VALUE rb_matrix4x4_initialize(int argc, VALUE *argv, VALUE self);
+
+// Instance
+VALUE rb_matrix4x4_identity_p(VALUE self);
+VALUE rb_matrix4x4_translation(VALUE self);
+VALUE rb_matrix4x4_translation_set(VALUE self, VALUE value);
+VALUE rb_matrix4x4_determinant(VALUE self);
+VALUE rb_matrix4x4_add(VALUE self, VALUE other);
+VALUE rb_matrix4x4_subtract(VALUE self, VALUE other);
+VALUE rb_matrix4x4_multiply(VALUE self, VALUE other);
+VALUE rb_matrix4x4_negate(VALUE self);
+VALUE rb_matrix4x4_equal(VALUE self, VALUE other);
+VALUE rb_matrix4x4_pow(VALUE self, VALUE exponent);
+VALUE rb_matrix4x4_row(VALUE self, VALUE row);
+VALUE rb_matrix4x4_column(VALUE self, VALUE column);
+VALUE rb_matrix4x4_each_row(VALUE self);
+VALUE rb_matrix4x4_each_column(VALUE self);
+VALUE rb_matrix4x4_aref(int argc, VALUE *argv, VALUE self);
+VALUE rb_matrix4x4_aset(int argc, VALUE *argv, VALUE self);
+VALUE rb_matrix4x4_invert(VALUE self);
+VALUE rb_matrix4x4_transform(VALUE self, VALUE quaternion);
+VALUE rb_matrix4x4_transform_bang(VALUE self, VALUE quaternion);
+VALUE rb_matrix4x4_transpose(VALUE self);
+VALUE rb_matrix4x4_lerp(VALUE self, VALUE other, VALUE amount);
+VALUE rb_matrix4x4_lerp_bang(VALUE self, VALUE other, VALUE amount);
+VALUE rb_matrix4x4_map(VALUE self);
+VALUE rb_matrix4x4_map_bang(VALUE self);
+
+// Conversion
+VALUE rb_matrix4x4_to_s(VALUE self);
+VALUE rb_matrix4x4_to_a(VALUE self);
+VALUE rb_matrix4x4_to_h(VALUE self);
+// TODO to_mat3x2
+
+// // Class
+VALUE rb_matrix4x4_identity(VALUE klass);
+VALUE rb_matrix4x4_create_billboard(VALUE klass, VALUE position, VALUE camera, VALUE up, VALUE forward);
+VALUE rb_matrix4x4_create_constrained_billboard(VALUE klass, VALUE position, VALUE rotate, VALUE camera, VALUE up, VALUE forward);
+VALUE rb_matrix4x4_create_translation(int argc, VALUE *argv, VALUE klass);
+VALUE rb_matrix4x4_create_scale(int argc, VALUE *argv, VALUE klass);
+VALUE rb_matrix4x4_create_rotation_x(int argc, VALUE *argv, VALUE klass);
+VALUE rb_matrix4x4_create_rotation_y(int argc, VALUE *argv, VALUE klass);
+VALUE rb_matrix4x4_create_rotation_z(int argc, VALUE *argv, VALUE klass);
+VALUE rb_matrix4x4_from_axis_angle(VALUE klass, VALUE axis, VALUE angle);
+VALUE rb_matrix4x4_create_perspective_fov(VALUE klass, VALUE fieldOfView, VALUE aspectRatio, VALUE nearDist, VALUE farDist);
+VALUE rb_matrix4x4_create_perspective(VALUE klass, VALUE width, VALUE height, VALUE nearDist, VALUE farDist);
+VALUE rb_matrix4x4_create_perspective_off_center(VALUE klass, VALUE left, VALUE right, VALUE bottom, VALUE top, VALUE nearDist, VALUE farDist);
+VALUE rb_matrix4x4_create_orthographic(VALUE klass, VALUE width, VALUE height, VALUE zNearPlane, VALUE zFarPlane);
+VALUE rb_matrix4x4_create_orthographic_off_center(VALUE klass, VALUE left, VALUE right, VALUE bottom, VALUE top, VALUE zNearPlane, VALUE zFarPlane);
+VALUE rb_matrix4x4_create_look_at(VALUE klass, VALUE camPos, VALUE camTarget, VALUE camUp);
+VALUE rb_matrix4x4_create_world(VALUE klass, VALUE position, VALUE forward, VALUE up);
+VALUE rb_matrix4x4_from_quaternion(VALUE klass, VALUE quaternion);
+VALUE rb_matrix4x4_from_yaw_pitch_roll(VALUE klass, VALUE yaw, VALUE pitch, VALUE roll);
+VALUE rb_matrix4x4_create_shadow(VALUE klass, VALUE light, VALUE plane);
+VALUE rb_matrix4x4_create_reflection(VALUE klass, VALUE plane);
+VALUE rb_matrix4x4_lerp_s(VALUE klass, VALUE matrix1, VALUE matrix2, VALUE amount);
+
+// Internal
+static inline void numerix_vec3_normalize(Vector3 *v);
+static inline void numerix_vec3_cross(Vector3 *v1, Vector3 *v2, Vector3 *result);
+static inline float numerix_vec3_dot(Vector3 *v1, Vector3 *v2);
+static inline void numerix_cross_norm(Vector3 *v1, Vector3 *v2, Vector3 *result);
+
+#endif /* NUMERIX_MATRIX4X4_H */

data/ext/numerix/matrix_base.c

@@ -0,0 +1,307 @@
+
+#include "matrix_base.h"
+
+VALUE rb_cMatrixBase;
+VALUE rb_cMatrix3x2;
+VALUE rb_cMatrix4x4;
+
+void Init_matrix_base(VALUE outer) {
+    rb_cMatrixBase = rb_define_class_under(outer, "MatrixBase", rb_cNumerixStruct);
+    rb_cMatrix3x2 = rb_define_class_under(outer, "Matrix3x2", rb_cMatrixBase);
+    rb_cMatrix4x4 = rb_define_class_under(outer, "Matrix4x4", rb_cMatrixBase);
+
+    rb_define_method(rb_cMatrixBase, "initialize", rb_numerix_abstract_initialize, 0);
+
+    rb_define_method(rb_cMatrix3x2, "m11", rb_matrix3x2_m11, 0);
+    rb_define_method(rb_cMatrix3x2, "m12", rb_matrix3x2_m12, 0);
+    rb_define_method(rb_cMatrix3x2, "m21", rb_matrix3x2_m21, 0);
+    rb_define_method(rb_cMatrix3x2, "m22", rb_matrix3x2_m22, 0);
+    rb_define_method(rb_cMatrix3x2, "m31", rb_matrix3x2_m31, 0);
+    rb_define_method(rb_cMatrix3x2, "m32", rb_matrix3x2_m32, 0);
+
+    rb_define_method(rb_cMatrix3x2, "m11=", rb_matrix3x2_m11_set, 1);
+    rb_define_method(rb_cMatrix3x2, "m12=", rb_matrix3x2_m12_set, 1);
+    rb_define_method(rb_cMatrix3x2, "m21=", rb_matrix3x2_m21_set, 1);
+    rb_define_method(rb_cMatrix3x2, "m22=", rb_matrix3x2_m22_set, 1);
+    rb_define_method(rb_cMatrix3x2, "m31=", rb_matrix3x2_m31_set, 1);
+    rb_define_method(rb_cMatrix3x2, "m32=", rb_matrix3x2_m32_set, 1);
+
+    rb_define_method(rb_cMatrix4x4, "m11", rb_matrix4x4_m11, 0);
+    rb_define_method(rb_cMatrix4x4, "m12", rb_matrix4x4_m12, 0);
+    rb_define_method(rb_cMatrix4x4, "m13", rb_matrix4x4_m13, 0);
+    rb_define_method(rb_cMatrix4x4, "m14", rb_matrix4x4_m14, 0);
+    rb_define_method(rb_cMatrix4x4, "m21", rb_matrix4x4_m21, 0);
+    rb_define_method(rb_cMatrix4x4, "m22", rb_matrix4x4_m22, 0);
+    rb_define_method(rb_cMatrix4x4, "m23", rb_matrix4x4_m23, 0);
+    rb_define_method(rb_cMatrix4x4, "m24", rb_matrix4x4_m24, 0);
+    rb_define_method(rb_cMatrix4x4, "m31", rb_matrix4x4_m31, 0);
+    rb_define_method(rb_cMatrix4x4, "m32", rb_matrix4x4_m32, 0);
+    rb_define_method(rb_cMatrix4x4, "m33", rb_matrix4x4_m33, 0);
+    rb_define_method(rb_cMatrix4x4, "m34", rb_matrix4x4_m34, 0);
+    rb_define_method(rb_cMatrix4x4, "m41", rb_matrix4x4_m41, 0);
+    rb_define_method(rb_cMatrix4x4, "m42", rb_matrix4x4_m42, 0);
+    rb_define_method(rb_cMatrix4x4, "m43", rb_matrix4x4_m43, 0);
+    rb_define_method(rb_cMatrix4x4, "m44", rb_matrix4x4_m44, 0);
+
+    rb_define_method(rb_cMatrix4x4, "m11=", rb_matrix4x4_m11_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m12=", rb_matrix4x4_m12_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m13=", rb_matrix4x4_m13_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m14=", rb_matrix4x4_m14_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m21=", rb_matrix4x4_m21_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m22=", rb_matrix4x4_m22_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m23=", rb_matrix4x4_m23_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m24=", rb_matrix4x4_m24_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m31=", rb_matrix4x4_m31_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m32=", rb_matrix4x4_m32_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m33=", rb_matrix4x4_m33_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m34=", rb_matrix4x4_m34_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m41=", rb_matrix4x4_m41_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m42=", rb_matrix4x4_m42_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m43=", rb_matrix4x4_m43_set, 1);
+    rb_define_method(rb_cMatrix4x4, "m44=", rb_matrix4x4_m44_set, 1);
+
+    PRIVATE_CLASS_METHOD(rb_cMatrixBase, "new");
+    rb_include_module(rb_cMatrixBase, rb_mEnumerable);
+}
+
+VALUE rb_matrix3x2_m11(VALUE self) {
+    MATRIX3X2();
+    return DBL2NUM(m->m11);
+}
+
+VALUE rb_matrix3x2_m12(VALUE self) {
+    MATRIX3X2();
+    return DBL2NUM(m->m12);
+}
+
+VALUE rb_matrix3x2_m21(VALUE self) {
+    MATRIX3X2();
+    return DBL2NUM(m->m21);
+}
+
+VALUE rb_matrix3x2_m22(VALUE self) {
+    MATRIX3X2();
+    return DBL2NUM(m->m22);
+}
+
+VALUE rb_matrix3x2_m31(VALUE self) {
+    MATRIX3X2();
+    return DBL2NUM(m->m31);
+}
+
+VALUE rb_matrix3x2_m32(VALUE self) {
+    MATRIX3X2();
+    return DBL2NUM(m->m32);
+}
+
+VALUE rb_matrix4x4_m11(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m11);
+}
+
+VALUE rb_matrix4x4_m12(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m12);
+}
+
+VALUE rb_matrix4x4_m13(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m13);
+}
+
+VALUE rb_matrix4x4_m14(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m14);
+}
+
+VALUE rb_matrix4x4_m21(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m21);
+}
+
+VALUE rb_matrix4x4_m22(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m22);
+}
+
+VALUE rb_matrix4x4_m23(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m23);
+}
+
+VALUE rb_matrix4x4_m24(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m24);
+}
+
+VALUE rb_matrix4x4_m31(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m31);
+}
+
+VALUE rb_matrix4x4_m32(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m32);
+}
+
+VALUE rb_matrix4x4_m33(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m33);
+}
+
+VALUE rb_matrix4x4_m34(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m34);
+}
+
+VALUE rb_matrix4x4_m41(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m41);
+}
+
+VALUE rb_matrix4x4_m42(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m42);
+}
+
+VALUE rb_matrix4x4_m43(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m43);
+}
+
+VALUE rb_matrix4x4_m44(VALUE self) {
+    MATRIX4X4();
+    return DBL2NUM(m->m44);
+}
+
+VALUE rb_matrix3x2_m11_set(VALUE self, VALUE value) {
+    MATRIX3X2();
+    m->m11 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix3x2_m12_set(VALUE self, VALUE value) {
+    MATRIX3X2();
+    m->m12 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix3x2_m21_set(VALUE self, VALUE value) {
+    MATRIX3X2();
+    m->m21 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix3x2_m22_set(VALUE self, VALUE value) {
+    MATRIX3X2();
+    m->m22 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix3x2_m31_set(VALUE self, VALUE value) {
+    MATRIX3X2();
+    m->m31 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix3x2_m32_set(VALUE self, VALUE value) {
+    MATRIX3X2();
+    m->m32 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m11_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m11 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m12_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m12 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m13_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m13 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m14_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m14 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m21_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m21 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m22_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m22 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m23_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m23 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m24_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m24 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m31_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m31 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m32_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m32 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m33_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m33 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m34_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m34 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m41_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m41 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m42_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m42 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m43_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m43 = NUM2FLT(value);
+    return value;
+}
+
+VALUE rb_matrix4x4_m44_set(VALUE self, VALUE value) {
+    MATRIX4X4();
+    m->m44 = NUM2FLT(value);
+    return value;
+}

data/ext/numerix/matrix_base.h

@@ -0,0 +1,70 @@
+#ifndef NUMERIX_MATRIX_BASE_H
+#define NUMERIX_MATRIX_BASE_H 1
+
+#include "common.h"
+
+#define MATRIX3X2() \
+    Matrix3x2 *m;   \
+    Data_Get_Struct(self, Matrix3x2, m)
+
+#define MATRIX4X4() \
+    Matrix4x4 *m;   \
+    Data_Get_Struct(self, Matrix4x4, m)
+
+void Init_matrix_base(VALUE outer);
+
+VALUE rb_matrix3x2_m11(VALUE self);
+VALUE rb_matrix3x2_m12(VALUE self);
+VALUE rb_matrix3x2_m21(VALUE self);
+VALUE rb_matrix3x2_m22(VALUE self);
+VALUE rb_matrix3x2_m31(VALUE self);
+VALUE rb_matrix3x2_m32(VALUE self);
+
+VALUE rb_matrix4x4_m11(VALUE self);
+VALUE rb_matrix4x4_m12(VALUE self);
+VALUE rb_matrix4x4_m13(VALUE self);
+VALUE rb_matrix4x4_m14(VALUE self);
+
+VALUE rb_matrix4x4_m21(VALUE self);
+VALUE rb_matrix4x4_m22(VALUE self);
+VALUE rb_matrix4x4_m23(VALUE self);
+VALUE rb_matrix4x4_m24(VALUE self);
+
+VALUE rb_matrix4x4_m31(VALUE self);
+VALUE rb_matrix4x4_m32(VALUE self);
+VALUE rb_matrix4x4_m33(VALUE self);
+VALUE rb_matrix4x4_m34(VALUE self);
+
+VALUE rb_matrix4x4_m41(VALUE self);
+VALUE rb_matrix4x4_m42(VALUE self);
+VALUE rb_matrix4x4_m43(VALUE self);
+VALUE rb_matrix4x4_m44(VALUE self);
+
+VALUE rb_matrix3x2_m11_set(VALUE self, VALUE value);
+VALUE rb_matrix3x2_m12_set(VALUE self, VALUE value);
+VALUE rb_matrix3x2_m21_set(VALUE self, VALUE value);
+VALUE rb_matrix3x2_m22_set(VALUE self, VALUE value);
+VALUE rb_matrix3x2_m31_set(VALUE self, VALUE value);
+VALUE rb_matrix3x2_m32_set(VALUE self, VALUE value);
+
+VALUE rb_matrix4x4_m11_set(VALUE self, VALUE value);
+VALUE rb_matrix4x4_m12_set(VALUE self, VALUE value);
+VALUE rb_matrix4x4_m13_set(VALUE self, VALUE value);
+VALUE rb_matrix4x4_m14_set(VALUE self, VALUE value);
+
+VALUE rb_matrix4x4_m21_set(VALUE self, VALUE value);
+VALUE rb_matrix4x4_m22_set(VALUE self, VALUE value);
+VALUE rb_matrix4x4_m23_set(VALUE self, VALUE value);
+VALUE rb_matrix4x4_m24_set(VALUE self, VALUE value);
+
+VALUE rb_matrix4x4_m31_set(VALUE self, VALUE value);
+VALUE rb_matrix4x4_m32_set(VALUE self, VALUE value);
+VALUE rb_matrix4x4_m33_set(VALUE self, VALUE value);
+VALUE rb_matrix4x4_m34_set(VALUE self, VALUE value);
+
+VALUE rb_matrix4x4_m41_set(VALUE self, VALUE value);
+VALUE rb_matrix4x4_m42_set(VALUE self, VALUE value);
+VALUE rb_matrix4x4_m43_set(VALUE self, VALUE value);
+VALUE rb_matrix4x4_m44_set(VALUE self, VALUE value);
+
+#endif /* NUMERIX_MATRIX_BASE_H */

data/ext/numerix/numerix.c

@@ -0,0 +1,33 @@
+#include "numerix.h"
+
+VALUE rb_mNumerix;
+VALUE rb_eNumerixError;
+
+void Init_numerix(void) {
+    rb_mNumerix = rb_define_module("Numerix");
+    rb_eNumerixError = rb_define_class_under(rb_mNumerix, "NumerixError", rb_eStandardError);
+
+    Init_basic_vector(rb_mNumerix);
+    Init_numerix_structure(rb_mNumerix);
+    Init_vector_base(rb_mNumerix);
+    Init_matrix_base(rb_mNumerix);
+    Init_vector2(rb_mNumerix);
+    Init_vector3(rb_mNumerix);
+    Init_vector4(rb_mNumerix);
+    Init_matrix3x2(rb_mNumerix);
+    Init_matrix4x4(rb_mNumerix);
+    Init_quaternion(rb_mNumerix);
+    Init_plane(rb_mNumerix);
+
+    rb_define_const(rb_mNumerix, "SIZEOF_VECTOR2", INT2NUM(sizeof(Vector2)));
+    rb_define_const(rb_mNumerix, "SIZEOF_VECTOR3", INT2NUM(sizeof(Vector3)));
+    rb_define_const(rb_mNumerix, "SIZEOF_VECTOR4", INT2NUM(sizeof(Vector4)));
+    rb_define_const(rb_mNumerix, "SIZEOF_MATRIX3X2", INT2NUM(sizeof(Matrix3x2)));
+    rb_define_const(rb_mNumerix, "SIZEOF_MATRiX4X4", INT2NUM(sizeof(Matrix4x4)));
+    rb_define_const(rb_mNumerix, "SIZEOF_QUATERNION", INT2NUM(sizeof(Quaternion)));
+    rb_define_const(rb_mNumerix, "SIZEOF_PLANE", INT2NUM(sizeof(Plane)));
+}
+
+VALUE rb_numerix_abstract_initialize(VALUE self) {
+    rb_raise(rb_eNumerixError, "cannot create instance of abstract class %s", rb_class2name(CLASS_OF(self)));
+}

data/ext/numerix/numerix.h

@@ -0,0 +1,19 @@
+#ifndef NUMERIX_H
+#define NUMERIX_H 1
+
+#include "common.h"
+#include "matrix3x2.h"
+#include "matrix4x4.h"
+#include "matrix_base.h"
+#include "plane.h"
+#include "quaternion.h"
+#include "structure.h"
+#include "vector.h"
+#include "vector2.h"
+#include "vector3.h"
+#include "vector4.h"
+#include "vector_base.h"
+
+VALUE rb_numerix_abstract_initialize(VALUE self);
+
+#endif /* NUMERIX_H */