rb-gsl 1.16.0.4 → 1.16.0.5

data/ext/gsl_native/const.c

@@ -0,0 +1,331 @@
+/*
+  const.c
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2006 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY.
+*/
+
+#include "include/rb_gsl_const.h"
+
+static void rb_gsl_const_mks(VALUE module);
+static void rb_gsl_const_cgs(VALUE module);
+static void rb_gsl_const_num(VALUE module);
+void Init_gsl_const_additional(VALUE mmks, VALUE mcgs, VALUE mnum);
+
+static void rb_gsl_const_mks(VALUE module)
+{
+  rb_define_const(module, "SPEED_OF_LIGHT",
+                  rb_float_new(GSL_CONST_MKSA_SPEED_OF_LIGHT));
+  rb_define_const(module, "GRAVITATIONAL_CONSTANT",
+                  rb_float_new(GSL_CONST_MKSA_GRAVITATIONAL_CONSTANT));
+  rb_define_const(module, "PLANCKS_CONSTANT_H",
+                  rb_float_new(GSL_CONST_MKSA_PLANCKS_CONSTANT_H));
+  rb_define_const(module, "PLANCKS_CONSTANT_HBAR",
+                  rb_float_new(GSL_CONST_MKSA_PLANCKS_CONSTANT_HBAR));
+  rb_define_const(module, "VACUUM_PERMEABILITY",
+                  rb_float_new(GSL_CONST_MKSA_VACUUM_PERMEABILITY));
+  rb_define_const(module, "VACUUM_PERMITTIVITY",
+                  rb_float_new(GSL_CONST_MKSA_VACUUM_PERMITTIVITY));
+  rb_define_const(module, "ASTRONOMICAL_UNIT",
+                  rb_float_new(GSL_CONST_MKSA_ASTRONOMICAL_UNIT));
+  rb_define_const(module, "LIGHT_YEAR", rb_float_new(GSL_CONST_MKSA_LIGHT_YEAR));
+  rb_define_const(module, "PARSEC", rb_float_new(GSL_CONST_MKSA_PARSEC));
+  rb_define_const(module, "GRAV_ACCEL", rb_float_new(GSL_CONST_MKSA_GRAV_ACCEL));
+  rb_define_const(module, "ELECTRON_VOLT",
+                  rb_float_new(GSL_CONST_MKSA_ELECTRON_VOLT));
+  rb_define_const(module, "MASS_ELECTRON",
+                  rb_float_new(GSL_CONST_MKSA_MASS_ELECTRON));
+  rb_define_const(module, "MASS_MUON", rb_float_new(GSL_CONST_MKSA_MASS_MUON));
+  rb_define_const(module, "MASS_PROTON", rb_float_new(GSL_CONST_MKSA_MASS_PROTON));
+  rb_define_const(module, "MASS_NEUTRON", rb_float_new(GSL_CONST_MKSA_MASS_NEUTRON));
+  rb_define_const(module, "RYDBERG", rb_float_new(GSL_CONST_MKSA_RYDBERG));
+  rb_define_const(module, "BOHR_MAGNETON",
+                  rb_float_new(GSL_CONST_MKSA_BOHR_MAGNETON));
+  rb_define_const(module, "NUCLEAR_MAGNETON",
+                  rb_float_new(GSL_CONST_MKSA_NUCLEAR_MAGNETON));
+  rb_define_const(module, "ELECTRON_MAGNETIC_MOMENT",
+                  rb_float_new(GSL_CONST_MKSA_ELECTRON_MAGNETIC_MOMENT));
+  rb_define_const(module, "PROTON_MAGNETIC_MOMENT",
+                  rb_float_new(GSL_CONST_MKSA_PROTON_MAGNETIC_MOMENT));
+  rb_define_const(module, "STANDARD_GAS_VOLUME",
+                  rb_float_new(GSL_CONST_MKSA_STANDARD_GAS_VOLUME));
+
+  rb_define_const(module, "MINUTE", rb_float_new(GSL_CONST_MKSA_MINUTE));
+  rb_define_const(module, "HOUR", rb_float_new(GSL_CONST_MKSA_HOUR));
+  rb_define_const(module, "DAY", rb_float_new(GSL_CONST_MKSA_DAY));
+  rb_define_const(module, "WEEK", rb_float_new(GSL_CONST_MKSA_WEEK));
+  rb_define_const(module, "INCH", rb_float_new(GSL_CONST_MKSA_INCH));
+  rb_define_const(module, "FOOT", rb_float_new(GSL_CONST_MKSA_FOOT));
+  rb_define_const(module, "YARD", rb_float_new(GSL_CONST_MKSA_YARD));
+  rb_define_const(module, "MILE", rb_float_new(GSL_CONST_MKSA_MILE));
+  rb_define_const(module, "NAUTICAL_MILE",
+                  rb_float_new(GSL_CONST_MKSA_NAUTICAL_MILE));
+  rb_define_const(module, "FATHOM", rb_float_new(GSL_CONST_MKSA_FATHOM));
+  rb_define_const(module, "MIL", rb_float_new(GSL_CONST_MKSA_MIL));
+  rb_define_const(module, "POINT", rb_float_new(GSL_CONST_MKSA_POINT));
+  rb_define_const(module, "TEXPOINT", rb_float_new(GSL_CONST_MKSA_TEXPOINT));
+  rb_define_const(module, "MICRON", rb_float_new(GSL_CONST_MKSA_MICRON));
+  rb_define_const(module, "ANGSTROM", rb_float_new(GSL_CONST_MKSA_ANGSTROM));
+  rb_define_const(module, "HECTARE", rb_float_new(GSL_CONST_MKSA_HECTARE));
+  rb_define_const(module, "ACRE", rb_float_new(GSL_CONST_MKSA_ACRE));
+  rb_define_const(module, "BARN", rb_float_new(GSL_CONST_MKSA_BARN));
+  rb_define_const(module, "LITER", rb_float_new(GSL_CONST_MKSA_LITER));
+  rb_define_const(module, "US_GALLON", rb_float_new(GSL_CONST_MKSA_US_GALLON));
+  rb_define_const(module, "QUART", rb_float_new(GSL_CONST_MKSA_QUART));
+  rb_define_const(module, "PINT", rb_float_new(GSL_CONST_MKSA_PINT));
+  rb_define_const(module, "CUP", rb_float_new(GSL_CONST_MKSA_CUP));
+  rb_define_const(module, "FLUID_OUNCE", rb_float_new(GSL_CONST_MKSA_FLUID_OUNCE));
+  rb_define_const(module, "TABLESPOON", rb_float_new(GSL_CONST_MKSA_TABLESPOON));
+  rb_define_const(module, "CANADIAN_GALLON",
+                  rb_float_new(GSL_CONST_MKSA_CANADIAN_GALLON));
+
+  rb_define_const(module, "UK_GALLON", rb_float_new(GSL_CONST_MKSA_UK_GALLON));
+  rb_define_const(module, "KILOMETERS_PER_HOUR",
+                  rb_float_new(GSL_CONST_MKSA_MILES_PER_HOUR));
+  rb_define_const(module, "MILES_PER_HOUR",
+                  rb_float_new(GSL_CONST_MKSA_KILOMETERS_PER_HOUR));
+  rb_define_const(module, "KNOT", rb_float_new(GSL_CONST_MKSA_KNOT));
+  rb_define_const(module, "POUND_MASS", rb_float_new(GSL_CONST_MKSA_POUND_MASS));
+  rb_define_const(module, "POUND_OUNCE", rb_float_new(GSL_CONST_MKSA_OUNCE_MASS));
+  rb_define_const(module, "POUND_TON", rb_float_new(GSL_CONST_MKSA_TON));
+  rb_define_const(module, "POUND_METRIC_TON",
+                  rb_float_new(GSL_CONST_MKSA_METRIC_TON));
+  rb_define_const(module, "POUND_UK_TON", rb_float_new(GSL_CONST_MKSA_UK_TON));
+  rb_define_const(module, "POUND_TROY_OUNCE",
+                  rb_float_new(GSL_CONST_MKSA_TROY_OUNCE));
+  rb_define_const(module, "CARAT", rb_float_new(GSL_CONST_MKSA_CARAT));
+  rb_define_const(module, "UNIFIED_ATOMIC_MASS",
+                  rb_float_new(GSL_CONST_MKSA_UNIFIED_ATOMIC_MASS));
+  rb_define_const(module, "GRAM_FORCE", rb_float_new(GSL_CONST_MKSA_GRAM_FORCE));
+  rb_define_const(module, "POUND_FORCE", rb_float_new(GSL_CONST_MKSA_POUND_FORCE));
+  rb_define_const(module, "KILOPOUND_FORCE",
+                  rb_float_new(GSL_CONST_MKSA_KILOPOUND_FORCE));
+  rb_define_const(module, "POUNDAL", rb_float_new(GSL_CONST_MKSA_POUNDAL));
+  rb_define_const(module, "CALORIE", rb_float_new(GSL_CONST_MKSA_CALORIE));
+  rb_define_const(module, "BTU", rb_float_new(GSL_CONST_MKSA_BTU));
+  rb_define_const(module, "THERM", rb_float_new(GSL_CONST_MKSA_THERM));
+  rb_define_const(module, "HORSEPOWER", rb_float_new(GSL_CONST_MKSA_HORSEPOWER));
+  rb_define_const(module, "BAR", rb_float_new(GSL_CONST_MKSA_BAR));
+  rb_define_const(module, "STD_ATMOSPHERE",
+                  rb_float_new(GSL_CONST_MKSA_STD_ATMOSPHERE));
+  rb_define_const(module, "TORR", rb_float_new(GSL_CONST_MKSA_TORR));
+  rb_define_const(module, "METER_OF_MERCURY",
+                  rb_float_new(GSL_CONST_MKSA_METER_OF_MERCURY));
+  rb_define_const(module, "INCH_OF_MERCURY",
+                  rb_float_new(GSL_CONST_MKSA_INCH_OF_MERCURY));
+  rb_define_const(module, "INCH_OF_WATER",
+                  rb_float_new(GSL_CONST_MKSA_INCH_OF_WATER));
+  rb_define_const(module, "PSI", rb_float_new(GSL_CONST_MKSA_PSI));
+  rb_define_const(module, "POISE", rb_float_new(GSL_CONST_MKSA_POISE));
+  rb_define_const(module, "STOKES", rb_float_new(GSL_CONST_MKSA_STOKES));
+  rb_define_const(module, "FARADAY", rb_float_new(GSL_CONST_MKSA_FARADAY));
+  rb_define_const(module, "ELECTRON_CHARGE",
+                  rb_float_new(GSL_CONST_MKSA_ELECTRON_CHARGE));
+  rb_define_const(module, "GAUSS", rb_float_new(GSL_CONST_MKSA_GAUSS));
+  rb_define_const(module, "STILB", rb_float_new(GSL_CONST_MKSA_STILB));
+  rb_define_const(module, "LUMEN", rb_float_new(GSL_CONST_MKSA_LUMEN));
+  rb_define_const(module, "LUX", rb_float_new(GSL_CONST_MKSA_LUX));
+  rb_define_const(module, "PHOT", rb_float_new(GSL_CONST_MKSA_PHOT));
+  rb_define_const(module, "FOOTCANDLE", rb_float_new(GSL_CONST_MKSA_FOOTCANDLE));
+  rb_define_const(module, "LAMBERT", rb_float_new(GSL_CONST_MKSA_LAMBERT));
+  rb_define_const(module, "CURIE", rb_float_new(GSL_CONST_MKSA_CURIE));
+  rb_define_const(module, "ROENTGEN", rb_float_new(GSL_CONST_MKSA_ROENTGEN));
+  rb_define_const(module, "RAD", rb_float_new(GSL_CONST_MKSA_RAD));
+  rb_define_const(module, "SOLAR_MASS", rb_float_new(GSL_CONST_MKSA_SOLAR_MASS));
+
+  rb_define_const(module, "BOLTZMANN", rb_float_new(GSL_CONST_MKSA_BOLTZMANN));
+  rb_define_const(module, "MOLAR_GAS", rb_float_new(GSL_CONST_MKSA_MOLAR_GAS));
+
+  rb_define_const(module, "BOHR_RADIUS", rb_float_new(GSL_CONST_MKSA_BOHR_RADIUS));
+  rb_define_const(module, "NEWTON", rb_float_new(GSL_CONST_MKSA_NEWTON));
+  rb_define_const(module, "DYNE", rb_float_new(GSL_CONST_MKSA_DYNE));
+  rb_define_const(module, "JOULE", rb_float_new(GSL_CONST_MKSA_JOULE));
+  rb_define_const(module, "ERG", rb_float_new(GSL_CONST_MKSA_ERG));
+
+  rb_define_const(module, "STEFAN_BOLTZMANN_CONSTANT",
+                  rb_float_new(GSL_CONST_MKSA_STEFAN_BOLTZMANN_CONSTANT));
+  rb_define_const(module, "THOMSON_CROSS_SECTION",
+                  rb_float_new(GSL_CONST_MKSA_THOMSON_CROSS_SECTION));
+  rb_define_const(module, "DEBYE", rb_float_new(GSL_CONST_MKSA_DEBYE));
+}
+
+
+static void rb_gsl_const_cgs(VALUE module)
+{
+  rb_define_const(module, "SPEED_OF_LIGHT",
+                  rb_float_new(GSL_CONST_CGSM_SPEED_OF_LIGHT));
+  rb_define_const(module, "GRAVITATIONAL_CONSTANT",
+                  rb_float_new(GSL_CONST_CGSM_GRAVITATIONAL_CONSTANT));
+  rb_define_const(module, "PLANCKS_CONSTANT_H",
+                  rb_float_new(GSL_CONST_CGSM_PLANCKS_CONSTANT_H));
+  rb_define_const(module, "PLANCKS_CONSTANT_HBAR",
+                  rb_float_new(GSL_CONST_CGSM_PLANCKS_CONSTANT_HBAR));
+  rb_define_const(module, "ASTRONOMICAL_UNIT",
+                  rb_float_new(GSL_CONST_CGSM_ASTRONOMICAL_UNIT));
+  rb_define_const(module, "LIGHT_YEAR",
+                  rb_float_new(GSL_CONST_CGSM_LIGHT_YEAR));
+  rb_define_const(module, "PARSEC",
+                  rb_float_new(GSL_CONST_CGSM_PARSEC));
+  rb_define_const(module, "GRAV_ACCEL",
+                  rb_float_new(GSL_CONST_CGSM_GRAV_ACCEL));
+  rb_define_const(module, "ELECTRON_VOLT",
+                  rb_float_new(GSL_CONST_CGSM_ELECTRON_VOLT));
+  rb_define_const(module, "MASS_ELECTRON",
+                  rb_float_new(GSL_CONST_CGSM_MASS_ELECTRON));
+  rb_define_const(module, "MASS_MUON", rb_float_new(GSL_CONST_CGSM_MASS_MUON));
+  rb_define_const(module, "MASS_PROTON", rb_float_new(GSL_CONST_CGSM_MASS_PROTON));
+  rb_define_const(module, "MASS_NEUTRON", rb_float_new(GSL_CONST_CGSM_MASS_NEUTRON));
+  rb_define_const(module, "RYDBERG", rb_float_new(GSL_CONST_CGSM_RYDBERG));
+
+  rb_define_const(module, "BOHR_MAGNETON",
+                  rb_float_new(GSL_CONST_CGSM_BOHR_MAGNETON));
+
+  rb_define_const(module, "NUCLEAR_MAGNETON",
+                  rb_float_new(GSL_CONST_CGSM_NUCLEAR_MAGNETON));
+  rb_define_const(module, "ELECTRON_MAGNETIC_MOMENT",
+                  rb_float_new(GSL_CONST_CGSM_ELECTRON_MAGNETIC_MOMENT));
+  rb_define_const(module, "PROTON_MAGNETIC_MOMENT",
+                  rb_float_new(GSL_CONST_CGSM_PROTON_MAGNETIC_MOMENT));
+  rb_define_const(module, "STANDARD_GAS_VOLUME",
+                  rb_float_new(GSL_CONST_CGSM_STANDARD_GAS_VOLUME));
+
+  rb_define_const(module, "MINUTE", rb_float_new(GSL_CONST_CGSM_MINUTE));
+  rb_define_const(module, "HOUR", rb_float_new(GSL_CONST_CGSM_HOUR));
+  rb_define_const(module, "DAY", rb_float_new(GSL_CONST_CGSM_DAY));
+  rb_define_const(module, "WEEK", rb_float_new(GSL_CONST_CGSM_WEEK));
+  rb_define_const(module, "INCH", rb_float_new(GSL_CONST_CGSM_INCH));
+  rb_define_const(module, "FOOT", rb_float_new(GSL_CONST_CGSM_FOOT));
+  rb_define_const(module, "YARD", rb_float_new(GSL_CONST_CGSM_YARD));
+  rb_define_const(module, "MILE", rb_float_new(GSL_CONST_CGSM_MILE));
+  rb_define_const(module, "NAUTICAL_MILE",
+                  rb_float_new(GSL_CONST_CGSM_NAUTICAL_MILE));
+  rb_define_const(module, "FATHOM", rb_float_new(GSL_CONST_CGSM_FATHOM));
+  rb_define_const(module, "MIL", rb_float_new(GSL_CONST_CGSM_MIL));
+  rb_define_const(module, "POINT", rb_float_new(GSL_CONST_CGSM_POINT));
+  rb_define_const(module, "TEXPOINT", rb_float_new(GSL_CONST_CGSM_TEXPOINT));
+  rb_define_const(module, "MICRON", rb_float_new(GSL_CONST_CGSM_MICRON));
+  rb_define_const(module, "ANGSTROM", rb_float_new(GSL_CONST_CGSM_ANGSTROM));
+  rb_define_const(module, "HECTARE", rb_float_new(GSL_CONST_CGSM_HECTARE));
+  rb_define_const(module, "ACRE", rb_float_new(GSL_CONST_CGSM_ACRE));
+  rb_define_const(module, "BARN", rb_float_new(GSL_CONST_CGSM_BARN));
+  rb_define_const(module, "LITER", rb_float_new(GSL_CONST_CGSM_LITER));
+  rb_define_const(module, "US_GALLON", rb_float_new(GSL_CONST_CGSM_US_GALLON));
+  rb_define_const(module, "QUART", rb_float_new(GSL_CONST_CGSM_QUART));
+  rb_define_const(module, "PINT", rb_float_new(GSL_CONST_CGSM_PINT));
+  rb_define_const(module, "CUP", rb_float_new(GSL_CONST_CGSM_CUP));
+  rb_define_const(module, "FLUID_OUNCE", rb_float_new(GSL_CONST_CGSM_FLUID_OUNCE));
+  rb_define_const(module, "TABLESPOON", rb_float_new(GSL_CONST_CGSM_TABLESPOON));
+  rb_define_const(module, "CANADIAN_GALLON",
+                  rb_float_new(GSL_CONST_CGSM_CANADIAN_GALLON));
+
+  rb_define_const(module, "UK_GALLON", rb_float_new(GSL_CONST_CGSM_UK_GALLON));
+  rb_define_const(module, "KILOMETERS_PER_HOUR",
+                  rb_float_new(GSL_CONST_CGSM_MILES_PER_HOUR));
+  rb_define_const(module, "MILES_PER_HOUR",
+                  rb_float_new(GSL_CONST_CGSM_KILOMETERS_PER_HOUR));
+  rb_define_const(module, "KNOT", rb_float_new(GSL_CONST_CGSM_KNOT));
+  rb_define_const(module, "POUND_MASS", rb_float_new(GSL_CONST_CGSM_POUND_MASS));
+  rb_define_const(module, "POUND_OUNCE", rb_float_new(GSL_CONST_CGSM_OUNCE_MASS));
+  rb_define_const(module, "POUND_TON", rb_float_new(GSL_CONST_CGSM_TON));
+  rb_define_const(module, "POUND_METRIC_TON",
+                  rb_float_new(GSL_CONST_CGSM_METRIC_TON));
+  rb_define_const(module, "POUND_UK_TON", rb_float_new(GSL_CONST_CGSM_UK_TON));
+  rb_define_const(module, "POUND_TROY_OUNCE",
+                  rb_float_new(GSL_CONST_CGSM_TROY_OUNCE));
+  rb_define_const(module, "CARAT", rb_float_new(GSL_CONST_CGSM_CARAT));
+  rb_define_const(module, "UNIFIED_ATOMIC_MASS",
+                  rb_float_new(GSL_CONST_CGSM_UNIFIED_ATOMIC_MASS));
+  rb_define_const(module, "GRAM_FORCE", rb_float_new(GSL_CONST_CGSM_GRAM_FORCE));
+  rb_define_const(module, "POUND_FORCE", rb_float_new(GSL_CONST_CGSM_POUND_FORCE));
+  rb_define_const(module, "KILOPOUND_FORCE",
+                  rb_float_new(GSL_CONST_CGSM_KILOPOUND_FORCE));
+  rb_define_const(module, "POUNDAL", rb_float_new(GSL_CONST_CGSM_POUNDAL));
+  rb_define_const(module, "CALORIE", rb_float_new(GSL_CONST_CGSM_CALORIE));
+  rb_define_const(module, "BTU", rb_float_new(GSL_CONST_CGSM_BTU));
+  rb_define_const(module, "THERM", rb_float_new(GSL_CONST_CGSM_THERM));
+  rb_define_const(module, "HORSEPOWER", rb_float_new(GSL_CONST_CGSM_HORSEPOWER));
+  rb_define_const(module, "BAR", rb_float_new(GSL_CONST_CGSM_BAR));
+  rb_define_const(module, "STD_ATMOSPHERE",
+                  rb_float_new(GSL_CONST_CGSM_STD_ATMOSPHERE));
+  rb_define_const(module, "TORR", rb_float_new(GSL_CONST_CGSM_TORR));
+  rb_define_const(module, "METER_OF_MERCURY",
+                  rb_float_new(GSL_CONST_CGSM_METER_OF_MERCURY));
+  rb_define_const(module, "INCH_OF_MERCURY",
+                  rb_float_new(GSL_CONST_CGSM_INCH_OF_MERCURY));
+  rb_define_const(module, "INCH_OF_WATER",
+                  rb_float_new(GSL_CONST_CGSM_INCH_OF_WATER));
+  rb_define_const(module, "PSI", rb_float_new(GSL_CONST_CGSM_PSI));
+  rb_define_const(module, "POISE", rb_float_new(GSL_CONST_CGSM_POISE));
+  rb_define_const(module, "STOKES", rb_float_new(GSL_CONST_CGSM_STOKES));
+  rb_define_const(module, "FARADAY", rb_float_new(GSL_CONST_CGSM_FARADAY));
+  rb_define_const(module, "ELECTRON_CHARGE",
+                  rb_float_new(GSL_CONST_CGSM_ELECTRON_CHARGE));
+  rb_define_const(module, "ELECTRON_CHARGE_ESU",
+                  rb_float_new(GSL_CONST_CGSM_ELECTRON_CHARGE*GSL_CONST_CGSM_SPEED_OF_LIGHT));
+  rb_define_const(module, "STILB", rb_float_new(GSL_CONST_CGSM_STILB));
+  rb_define_const(module, "LUMEN", rb_float_new(GSL_CONST_CGSM_LUMEN));
+  rb_define_const(module, "LUX", rb_float_new(GSL_CONST_CGSM_LUX));
+  rb_define_const(module, "PHOT", rb_float_new(GSL_CONST_CGSM_PHOT));
+  rb_define_const(module, "FOOTCANDLE", rb_float_new(GSL_CONST_CGSM_FOOTCANDLE));
+  rb_define_const(module, "LAMBERT", rb_float_new(GSL_CONST_CGSM_LAMBERT));
+  rb_define_const(module, "CURIE", rb_float_new(GSL_CONST_CGSM_CURIE));
+  rb_define_const(module, "ROENTGEN", rb_float_new(GSL_CONST_CGSM_ROENTGEN));
+  rb_define_const(module, "RAD", rb_float_new(GSL_CONST_CGSM_RAD));
+  rb_define_const(module, "SOLAR_MASS", rb_float_new(GSL_CONST_CGSM_SOLAR_MASS));
+
+  rb_define_const(module, "BOLTZMANN", rb_float_new(GSL_CONST_CGSM_BOLTZMANN));
+  rb_define_const(module, "MOLAR_GAS", rb_float_new(GSL_CONST_CGSM_MOLAR_GAS));
+
+  rb_define_const(module, "BOHR_RADIUS", rb_float_new(GSL_CONST_CGSM_BOHR_RADIUS));
+  rb_define_const(module, "NEWTON", rb_float_new(GSL_CONST_CGSM_NEWTON));
+  rb_define_const(module, "DYNE", rb_float_new(GSL_CONST_CGSM_DYNE));
+  rb_define_const(module, "JOULE", rb_float_new(GSL_CONST_CGSM_JOULE));
+  rb_define_const(module, "ERG", rb_float_new(GSL_CONST_CGSM_ERG));
+
+  rb_define_const(module, "STEFAN_BOLTZMANN_CONSTANT",
+                  rb_float_new(GSL_CONST_CGSM_STEFAN_BOLTZMANN_CONSTANT));
+  rb_define_const(module, "THOMSON_CROSS_SECTION",
+                  rb_float_new(GSL_CONST_CGSM_THOMSON_CROSS_SECTION));
+}
+
+static void rb_gsl_const_num(VALUE module)
+{
+  rb_define_const(module, "AVOGADRO", rb_float_new(GSL_CONST_NUM_AVOGADRO));
+  rb_define_const(module, "FINE_STRUCTURE",
+                  rb_float_new(GSL_CONST_NUM_FINE_STRUCTURE));
+  rb_define_const(module, "YOTTA", rb_float_new(GSL_CONST_NUM_YOTTA));
+  rb_define_const(module, "ZETTA", rb_float_new(GSL_CONST_NUM_ZETTA));
+  rb_define_const(module, "EXA", rb_float_new(GSL_CONST_NUM_EXA));
+  rb_define_const(module, "PETA", rb_float_new(GSL_CONST_NUM_PETA));
+  rb_define_const(module, "TERA", rb_float_new(GSL_CONST_NUM_TERA));
+  rb_define_const(module, "GIGA", rb_float_new(GSL_CONST_NUM_GIGA));
+  rb_define_const(module, "MEGA", rb_float_new(GSL_CONST_NUM_MEGA));
+  rb_define_const(module, "KILO", rb_float_new(GSL_CONST_NUM_KILO));
+  rb_define_const(module, "MILLI", rb_float_new(GSL_CONST_NUM_MILLI));
+  rb_define_const(module, "MICRO", rb_float_new(GSL_CONST_NUM_MICRO));
+  rb_define_const(module, "NANO", rb_float_new(GSL_CONST_NUM_NANO));
+  rb_define_const(module, "PICO", rb_float_new(GSL_CONST_NUM_PICO));
+  rb_define_const(module, "FEMTO", rb_float_new(GSL_CONST_NUM_FEMTO));
+  rb_define_const(module, "ATTO", rb_float_new(GSL_CONST_NUM_ATTO));
+  rb_define_const(module, "ZEPTO", rb_float_new(GSL_CONST_NUM_ZEPTO));
+  rb_define_const(module, "YOCTO", rb_float_new(GSL_CONST_NUM_YOCTO));
+}
+
+void Init_gsl_const(VALUE module)
+{
+  VALUE mgsl_const;
+  VALUE mgsl_const_mks, mgsl_const_cgs, mgsl_const_num;
+
+  mgsl_const = rb_define_module_under(module, "CONST");
+  mgsl_const_mks = rb_define_module_under(mgsl_const, "MKSA");
+  rb_gsl_const_mks(mgsl_const_mks);
+  mgsl_const_cgs = rb_define_module_under(mgsl_const, "CGSM");
+  rb_gsl_const_cgs(mgsl_const_cgs);
+  mgsl_const_num = rb_define_module_under(mgsl_const, "NUM");
+  rb_gsl_const_num(mgsl_const_num);
+  Init_gsl_const_additional(mgsl_const_mks, mgsl_const_cgs, mgsl_const_num);
+}

data/ext/{gsl → gsl_native}/const_additional.c

@@ -11,13 +11,6 @@
 
 #include "include/rb_gsl_const.h"
 
-#ifndef GSL_1_4_9_LATER
-#define MKS_STEFAN_BOLTZMANN_CONSTANT (5.6703993443638e-08)
-#define MKS_THOMSON_CROSS_SECTION (6.65245852869236e-29)
-#define CGS_STEFAN_BOLTZMANN_CONSTANT (5.6703993443638e-05)
-#define CGS_THOMSON_CROSS_SECTION (6.65245852869236e-25)
-#endif
-
 #define MKS_CLASSICAL_ELECTRON_RADIUS (2.81794028310825e-15)
 #define MKS_RADIATION_DENSITY_CONSTANT (7.56576650685962e-16)
 #define MKS_RADIATION_NUMBER_DENSITY_CONSTANT (20.2869161066108e-6)
@@ -43,47 +36,33 @@ static void rb_gsl_const_num(VALUE module);
 static void rb_gsl_const_mks(VALUE module)
 {
   rb_define_const(module, "RADIATION_DENSITY_CONSTANT",
-		  rb_float_new(MKS_RADIATION_DENSITY_CONSTANT));
-  rb_define_const(module, "RADIATION_NUMBER_DENSITY_CONSTANT", 
-		  rb_float_new(MKS_RADIATION_NUMBER_DENSITY_CONSTANT));
+                  rb_float_new(MKS_RADIATION_DENSITY_CONSTANT));
+  rb_define_const(module, "RADIATION_NUMBER_DENSITY_CONSTANT",
+                  rb_float_new(MKS_RADIATION_NUMBER_DENSITY_CONSTANT));
   rb_define_const(module, "CLASSICAL_ELECTRON_RADIUS",
-		  rb_float_new(MKS_CLASSICAL_ELECTRON_RADIUS));
+                  rb_float_new(MKS_CLASSICAL_ELECTRON_RADIUS));
   rb_define_const(module, "SOLAR_TIME", rb_float_new(MKS_SOLAR_TIME));
   rb_define_const(module, "SOLAR_GM", rb_float_new(MKS_SOLAR_GM));
 
   rb_define_const(module, "PLANCK_MASS", rb_float_new(MKS_PLANCK_MASS));
   rb_define_const(module, "PLANCK_LENGTH", rb_float_new(MKS_PLANCK_LENGTH));
   rb_define_const(module, "PLANCK_TIME", rb_float_new(MKS_PLANCK_TIME));
-
-#ifndef GSL_1_4_9_LATER
-  rb_define_const(module, "STEFAN_BOLTZMANN_CONSTANT", 
-		  rb_float_new(MKS_STEFAN_BOLTZMANN_CONSTANT));
-  rb_define_const(module, "THOMSON_CROSS_SECTION", 
-		  rb_float_new(MKS_THOMSON_CROSS_SECTION));
-#endif
 }
 
 static void rb_gsl_const_cgs(VALUE module)
 {
-  rb_define_const(module, "RADIATION_DENSITY_CONSTANT", 
-		  rb_float_new(CGS_RADIATION_DENSITY_CONSTANT));
-  rb_define_const(module, "RADIATION_NUMBER_DENSITY_CONSTANT", 
-		  rb_float_new(CGS_RADIATION_NUMBER_DENSITY_CONSTANT));
-  rb_define_const(module, "CLASSICAL_ELECTRON_RADIUS", 
-		  rb_float_new(CGS_CLASSICAL_ELECTRON_RADIUS));
+  rb_define_const(module, "RADIATION_DENSITY_CONSTANT",
+                  rb_float_new(CGS_RADIATION_DENSITY_CONSTANT));
+  rb_define_const(module, "RADIATION_NUMBER_DENSITY_CONSTANT",
+                  rb_float_new(CGS_RADIATION_NUMBER_DENSITY_CONSTANT));
+  rb_define_const(module, "CLASSICAL_ELECTRON_RADIUS",
+                  rb_float_new(CGS_CLASSICAL_ELECTRON_RADIUS));
   rb_define_const(module, "SOLAR_TIME", rb_float_new(CGS_SOLAR_TIME));
   rb_define_const(module, "SOLAR_GM", rb_float_new(CGS_SOLAR_GM));
 
   rb_define_const(module, "PLANCK_MASS", rb_float_new(CGS_PLANCK_MASS));
   rb_define_const(module, "PLANCK_LENGTH", rb_float_new(CGS_PLANCK_LENGTH));
   rb_define_const(module, "PLANCK_TIME", rb_float_new(CGS_PLANCK_TIME));
-
-#ifndef GSL_1_4_9_LATER
-  rb_define_const(module, "STEFAN_BOLTZMANN_CONSTANT", 
-		  rb_float_new(CGS_STEFAN_BOLTZMANN_CONSTANT));
-  rb_define_const(module, "THOMSON_CROSS_SECTION", 
-		  rb_float_new(CGS_THOMSON_CROSS_SECTION));
-#endif
 }
 
 static void rb_gsl_const_num(VALUE module)
@@ -102,8 +81,8 @@ void Init_gsl_const_additional(VALUE mmks, VALUE mcgs, VALUE mnum)
 #undef MKS_STEFAN_BOLTZMANN_CONSTANT
 #undef MKS_RADIATION_DENSITY_CONSTANT
 #undef MKS_RADIATION_NUMBER_DENSITY_CONSTANT
-#undef CGS_CLASSICAL_ELECTRON_RADIUS 
-#undef CGS_STEFAN_BOLTZMANN_CONSTANT 
+#undef CGS_CLASSICAL_ELECTRON_RADIUS
+#undef CGS_STEFAN_BOLTZMANN_CONSTANT
 #undef CGS_RADIATION_DENSITY_CONSTANT
 #undef CGS_RADIATION_NUMBER_DENSITY_CONSTANT
 #undef CGS_THOMSON_CROSS_SECTION

data/ext/gsl_native/cqp.c

@@ -0,0 +1,283 @@
+#ifdef HAVE_GSL_GSL_CQP_H
+
+#include "include/rb_gsl.h"
+#include "gsl/gsl_cqp.h"
+
+static const gsl_cqpminimizer_type* type_by_string(VALUE t);
+static const gsl_cqpminimizer_type* get_type(VALUE t)
+{
+
+  switch (TYPE(t)) {
+  case T_STRING:
+    return type_by_string(t);
+    break;
+  default:
+    rb_raise(rb_eTypeError, "Wrong argument type %s.", rb_class2name(CLASS_OF(t)));
+  }
+}
+
+static const gsl_cqpminimizer_type* type_by_string(VALUE t)
+{
+  char *name;
+  name = STR2CSTR(t);
+  if (strcmp(name, "mg_pdip") == 0) {
+    return gsl_cqpminimizer_mg_pdip;
+  } else {
+    rb_raise(rb_eRuntimeError, "Unknown minimizer type %s.", name);
+  }
+  return NULL; /* never reach here */
+}
+
+static VALUE rb_cqpminimizer_alloc(VALUE klass,  VALUE t, VALUE n, VALUE me, VALUE mi)
+{
+  gsl_cqpminimizer *m;
+  m = gsl_cqpminimizer_alloc(get_type(t), (size_t) FIX2INT(n), (size_t) FIX2INT(me), (size_t) FIX2INT(mi));
+  return Data_Wrap_Struct(klass, 0, gsl_cqpminimizer_free, m);
+}
+
+static VALUE rb_cqpminimizer_set(VALUE obj, VALUE data)
+{
+  gsl_cqpminimizer *m;
+  gsl_cqp_data *d;
+  Data_Get_Struct(obj, gsl_cqpminimizer, m);
+  Data_Get_Struct(data, gsl_cqp_data, d);
+  gsl_cqpminimizer_set(m, d);
+  return Qtrue;
+}
+
+static VALUE rb_cqpminimizer_name(VALUE obj)
+{
+  gsl_cqpminimizer *m;
+  Data_Get_Struct(obj, gsl_cqpminimizer, m);
+  return rb_str_new2(gsl_cqpminimizer_name(m));
+}
+
+static VALUE rb_cqpminimizer_iterate(VALUE obj)
+{
+  gsl_cqpminimizer *m;
+  Data_Get_Struct(obj, gsl_cqpminimizer, m);
+  return INT2FIX(gsl_cqpminimizer_iterate(m));
+}
+
+static VALUE rb_cqpminimizer_x(VALUE obj)
+{
+  gsl_cqpminimizer *m;
+  Data_Get_Struct(obj, gsl_cqpminimizer, m);
+  return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, gsl_cqpminimizer_x(m));
+}
+
+static VALUE rb_cqpminimizer_lm_eq(VALUE obj)
+{
+  gsl_cqpminimizer *m;
+  Data_Get_Struct(obj, gsl_cqpminimizer, m);
+  return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, gsl_cqpminimizer_lm_eq(m));
+}
+static VALUE rb_cqpminimizer_lm_ineq(VALUE obj)
+{
+  gsl_cqpminimizer *m;
+  Data_Get_Struct(obj, gsl_cqpminimizer, m);
+  return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, gsl_cqpminimizer_lm_ineq(m));
+}
+static VALUE rb_cqpminimizer_f(VALUE obj)
+{
+  gsl_cqpminimizer *m;
+  Data_Get_Struct(obj, gsl_cqpminimizer, m);
+  return rb_float_new(gsl_cqpminimizer_f(m));
+}
+static VALUE rb_cqpminimizer_gap(VALUE obj)
+{
+  gsl_cqpminimizer *m;
+  Data_Get_Struct(obj, gsl_cqpminimizer, m);
+  return rb_float_new(gsl_cqpminimizer_gap(m));
+}
+static VALUE rb_cqpminimizer_residuals_norm(VALUE obj)
+{
+  gsl_cqpminimizer *m;
+  Data_Get_Struct(obj, gsl_cqpminimizer, m);
+  return rb_float_new(gsl_cqpminimizer_residuals_norm(m));
+}
+/*
+static VALUE rb_cqpminimizer_minimum(VALUE obj)
+{
+  gsl_cqpminimizer *m;
+  Data_Get_Struct(obj, gsl_cqpminimizer, m);
+  return rb_float_new(gsl_cqpminimizer_minimum(m));
+}
+*/
+static VALUE rb_cqpminimizer_test_convergence(VALUE obj, VALUE g, VALUE r)
+{
+  gsl_cqpminimizer *m;
+  Data_Get_Struct(obj, gsl_cqpminimizer, m);
+  return INT2FIX(gsl_cqpminimizer_test_convergence(m, NUM2DBL(g), NUM2DBL(r)));
+}
+static VALUE rb_cqpminimizer_test_infeasibility(VALUE obj, VALUE e)
+{
+  gsl_cqpminimizer *m;
+  Data_Get_Struct(obj, gsl_cqpminimizer, m);
+  return INT2FIX(gsl_cqp_minimizer_test_infeasibility(m, NUM2DBL(e)));
+}
+
+static VALUE rb_cqp_data_alloc(VALUE klass)
+{
+  gsl_cqp_data *d;
+  d = (gsl_cqp_data*) malloc(sizeof(gsl_cqp_data));
+  return Data_Wrap_Struct(klass, 0, free, d);
+}
+
+static VALUE rb_cqp_data_Q(VALUE obj)
+{
+  gsl_cqp_data *d;
+  Data_Get_Struct(obj, gsl_cqp_data, d);
+  return Data_Wrap_Struct(cgsl_matrix_view, 0, NULL,   d->Q);
+}
+
+static VALUE rb_cqp_data_q(VALUE obj)
+{
+  gsl_cqp_data *d;
+  Data_Get_Struct(obj, gsl_cqp_data, d);
+  return Data_Wrap_Struct(cgsl_vector_view, 0, NULL,   d->q);
+}
+
+static VALUE rb_cqp_data_A(VALUE obj)
+{
+  gsl_cqp_data *d;
+  Data_Get_Struct(obj, gsl_cqp_data, d);
+  return Data_Wrap_Struct(cgsl_matrix_view, 0, NULL,   d->A);
+}
+
+static VALUE rb_cqp_data_b(VALUE obj)
+{
+  gsl_cqp_data *d;
+  Data_Get_Struct(obj, gsl_cqp_data, d);
+  return Data_Wrap_Struct(cgsl_vector_view, 0, NULL,   d->b);
+}
+
+static VALUE rb_cqp_data_C(VALUE obj)
+{
+  gsl_cqp_data *d;
+  Data_Get_Struct(obj, gsl_cqp_data, d);
+  return Data_Wrap_Struct(cgsl_matrix_view, 0, NULL,   d->C);
+}
+
+static VALUE rb_cqp_data_d(VALUE obj)
+{
+  gsl_cqp_data *d;
+  Data_Get_Struct(obj, gsl_cqp_data, d);
+  return Data_Wrap_Struct(cgsl_vector_view, 0, NULL,   d->d);
+}
+
+static VALUE rb_cqp_data_set_Q(VALUE obj, VALUE mm)
+{
+  gsl_cqp_data *d;
+  gsl_matrix *m;
+  Data_Get_Struct(obj, gsl_cqp_data, d);
+  CHECK_MATRIX(mm);
+  Data_Get_Struct(mm, gsl_matrix, m);
+  d->Q = m;
+  return Qtrue;
+}
+
+static VALUE rb_cqp_data_set_q(VALUE obj, VALUE vv)
+{
+  gsl_cqp_data *d;
+  gsl_vector *v;
+  Data_Get_Struct(obj, gsl_cqp_data, d);
+  CHECK_VECTOR(vv);
+  Data_Get_Struct(vv, gsl_vector, v);
+  d->q = v;
+  return Qtrue;
+}
+
+static VALUE rb_cqp_data_set_A(VALUE obj, VALUE mm)
+{
+  gsl_cqp_data *d;
+  gsl_matrix *m;
+  Data_Get_Struct(obj, gsl_cqp_data, d);
+  CHECK_MATRIX(mm);
+  Data_Get_Struct(mm, gsl_matrix, m);
+  d->A = m;
+  return Qtrue;
+}
+
+static VALUE rb_cqp_data_set_b(VALUE obj, VALUE vv)
+{
+  gsl_cqp_data *d;
+  gsl_vector *v;
+  Data_Get_Struct(obj, gsl_cqp_data, d);
+  CHECK_VECTOR(vv);
+  Data_Get_Struct(vv, gsl_vector, v);
+  d->b = v;
+  return Qtrue;
+}
+
+static VALUE rb_cqp_data_set_C(VALUE obj, VALUE mm)
+{
+  gsl_cqp_data *d;
+  gsl_matrix *m;
+  Data_Get_Struct(obj, gsl_cqp_data, d);
+  CHECK_MATRIX(mm);
+  Data_Get_Struct(mm, gsl_matrix, m);
+  d->C = m;
+  return Qtrue;
+}
+
+static VALUE rb_cqp_data_set_d(VALUE obj, VALUE vv)
+{
+  gsl_cqp_data *d;
+  gsl_vector *v;
+  Data_Get_Struct(obj, gsl_cqp_data, d);
+  CHECK_VECTOR(vv);
+  Data_Get_Struct(vv, gsl_vector, v);
+  d->d = v;
+  return Qtrue;
+}
+
+void Init_cqp(VALUE module)
+{
+  VALUE mCQP, cMinimizer, cData;
+
+  mCQP = rb_define_module_under(module, "CQP");
+  cMinimizer = rb_define_class_under(mCQP, "Minimizer", cGSL_Object);
+  cData = rb_define_class_under(mCQP, "Data", cGSL_Object);
+
+  rb_define_singleton_method(cMinimizer, "alloc", rb_cqpminimizer_alloc, 4);
+
+  rb_define_method(cMinimizer, "set", rb_cqpminimizer_set, 1);
+  rb_define_method(cMinimizer, "name", rb_cqpminimizer_name, 0);
+  rb_define_method(cMinimizer, "iterate", rb_cqpminimizer_iterate, 0);
+  rb_define_method(cMinimizer, "x", rb_cqpminimizer_x, 0);
+  rb_define_method(cMinimizer, "lm_eq", rb_cqpminimizer_lm_eq, 0);
+  rb_define_method(cMinimizer, "lm_ineq", rb_cqpminimizer_lm_ineq, 0);
+  rb_define_method(cMinimizer, "f", rb_cqpminimizer_f, 0);
+  rb_define_method(cMinimizer, "gap", rb_cqpminimizer_gap, 0);
+  rb_define_method(cMinimizer, "residuals_norm", rb_cqpminimizer_residuals_norm, 0);
+/*  rb_define_method(cMinimizer, "minimum", rb_cqpminimizer_minimum, 0);  */
+  rb_define_method(cMinimizer, "test_convergence", rb_cqpminimizer_test_convergence, 2);
+  rb_define_method(cMinimizer, "test_infeasibility", rb_cqpminimizer_test_infeasibility, 1);
+
+  /*****/
+  rb_define_singleton_method(cData, "alloc", rb_cqp_data_alloc, 0);
+
+  rb_define_method(cData, "Q", rb_cqp_data_Q, 0);
+  rb_define_method(cData, "q", rb_cqp_data_q, 0);
+  rb_define_method(cData, "A", rb_cqp_data_A, 0);
+  rb_define_method(cData, "b", rb_cqp_data_b, 0);
+  rb_define_method(cData, "C", rb_cqp_data_C, 0);
+  rb_define_method(cData, "d", rb_cqp_data_d, 0);
+
+  rb_define_method(cData, "set_Q", rb_cqp_data_set_Q, 1);
+  rb_define_method(cData, "set_q", rb_cqp_data_set_q, 1);
+  rb_define_method(cData, "set_A", rb_cqp_data_set_A, 1);
+  rb_define_method(cData, "set_b", rb_cqp_data_set_b, 1);
+  rb_define_method(cData, "set_C", rb_cqp_data_set_C, 1);
+  rb_define_method(cData, "set_d", rb_cqp_data_set_d, 1);
+  rb_define_alias(cData, "Q=", "set_Q");
+  rb_define_alias(cData, "q=", "set_q");
+  rb_define_alias(cData, "A=", "set_A");
+  rb_define_alias(cData, "b=", "set_b");
+  rb_define_alias(cData, "C=", "set_C");
+  rb_define_alias(cData, "d=", "set_d");
+}
+
+#endif
+