alglib 0.1.1 → 0.1.2

data/Rakefile

@@ -1,19 +1,28 @@
 #!/usr/bin/ruby
 # -*- ruby -*-
-
+MAKEFILE="ext/Makefile.#{RUBY_VERSION}"
+BASEDIR=File.dirname(__FILE__)
 require 'rubygems'
 require 'hoe'
-
+require 'fileutils'
 task :default => [:compile]
 task :compile => ["ext/alglib_ext.so"]
+task :clean do
+    FileList.new("ext/**/*.o","ext/**/*.so","ext/Makefile.*").each {|file|
+    FileUtils.rm(file)
+    }
+end
 
-file "ext/Makefile" => "ext/extconf.rb" do |t|
-    system %(cd ext && ruby extconf.rb)
+file MAKEFILE => "ext/extconf.rb" do |t|
+  Dir.chdir("#{BASEDIR}/ext")
+  ruby %{extconf.rb}
+  Dir.chdir("#{BASEDIR}")
 end
 
-file "ext/alglib_ext.so" => Dir.glob("ext/*.i")+["ext/Makefile", "ext/alglib/alglib_util.cpp", "ext/alglib/alglib_util.h"] do |t|
-    system %(cd ext && rake)
-    
+file "ext/alglib_ext.so" => Dir.glob("ext/*.i")+[MAKEFILE, "ext/alglib/alglib_util.cpp", "ext/alglib/alglib_util.h"] do |t|
+    Dir.chdir "#{BASEDIR}/ext"
+    ruby %{-S rake}
+    Dir.chdir BASEDIR
 end
 require File.dirname(__FILE__)+'/lib/alglib_version.rb'
 Hoe.spec 'alglib' do

data/ext/Rakefile

@@ -1,10 +1,9 @@
 require 'rake/testtask'
 require 'rake/rdoctask'
-
-
+MAKEFILE="Makefile.#{RUBY_VERSION}"
 task :default => [:compile]
 task :test => ["test.rb"] do |t|
-    system %(ruby test.rb)
+    ruby %(test.rb)
 end
 task :compile => ["alglib.so"]
 
@@ -15,10 +14,10 @@ end
 
 
 cpp=Dir.glob("alglib/*.cpp")+Dir.glob("alglib/*.h")
-file "alglib.so"=>cpp + ["Makefile"] do |t|
-	system %(make)
+file "alglib.so"=>cpp + [MAKEFILE] do |t|
+	system %(make -f #{MAKEFILE})
 end
 
-file "Makefile" => ["extconf.rb"] do |t|
-	system %(ruby extconf.rb)
+file MAKEFILE => ["extconf.rb"] do |t|
+	ruby %( extconf.rb)
 end

data/ext/alglib.i

@@ -2,8 +2,8 @@
 %{
 #include "ap.h"
 #include "alglib_util.h"
-#include "iostream"
-
+#include <iostream>
+#include <cstdio>
 %}
 
 %include "typemaps.i"
@@ -20,5 +20,6 @@
 %include "correlation.i"
 %include "correlationtests.i"
 %include "linreg.i"
+%include "mannwhitneyu.i"
 %include "logit.i"
-%include "pca.i"
+%include "pca.i"

data/ext/alglib/alglib.cpp

@@ -1,6 +1,6 @@
 /* ----------------------------------------------------------------------------
  * This file was automatically generated by SWIG (http://www.swig.org).
- * Version 1.3.35
+ * Version 1.3.36
  * 
  * This file is not intended to be easily readable and contains a number of 
  * coding conventions designed to improve portability and efficiency. Do not make
@@ -72,6 +72,12 @@ template <typename T> T SwigValueInit() {
 # endif
 #endif
 
+#ifndef SWIG_MSC_UNSUPPRESS_4505
+# if defined(_MSC_VER)
+#   pragma warning(disable : 4505) /* unreferenced local function has been removed */
+# endif 
+#endif
+
 #ifndef SWIGUNUSEDPARM
 # ifdef __cplusplus
 #   define SWIGUNUSEDPARM(p)
@@ -175,6 +181,12 @@ template <typename T> T SwigValueInit() {
 # endif
 #endif
 
+#ifndef SWIG_MSC_UNSUPPRESS_4505
+# if defined(_MSC_VER)
+#   pragma warning(disable : 4505) /* unreferenced local function has been removed */
+# endif 
+#endif
+
 #ifndef SWIGUNUSEDPARM
 # ifdef __cplusplus
 #   define SWIGUNUSEDPARM(p)
@@ -1804,7 +1816,7 @@ static VALUE mAlglib_ext;
 #define SWIG_RUBY_THREAD_END_BLOCK
 
 
-#define SWIGVERSION 0x010335 
+#define SWIGVERSION 0x010336 
 #define SWIG_VERSION SWIGVERSION
 
 
@@ -1817,8 +1829,8 @@ static VALUE mAlglib_ext;
 
 #include "ap.h"
 #include "alglib_util.h"
-#include "iostream"
-
+#include <iostream>
+#include <cstdio>
 
 
 #ifdef __cplusplus
@@ -2059,6 +2071,9 @@ VALUE spearmanrankcorrelationsignificance(double r,
 #include "linreg.h"
 
 
+#include "mannwhitneyu.h"
+
+
 #include "logit.h"
 
 
@@ -2083,8 +2098,8 @@ _wrap_new_Integer1dArray__SWIG_0(int argc, VALUE *argv, VALUE self) {
   if ((argc < 0) || (argc > 0)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
   }
-  result = (ap::template_1d_array< int > *)new ap::template_1d_array< int >();DATA_PTR(self) = result;
-  
+  result = (ap::template_1d_array< int > *)new ap::template_1d_array< int >();
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -2116,9 +2131,9 @@ _wrap_Integer1dArray_allocate(VALUE self) {
 SWIGINTERN VALUE
 _wrap_new_Integer1dArray__SWIG_1(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< int > *arg1 = 0 ;
-  ap::template_1d_array< int > *result = 0 ;
   void *argp1 ;
   int res1 = 0 ;
+  ap::template_1d_array< int > *result = 0 ;
   
   if ((argc < 1) || (argc > 1)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
@@ -2131,8 +2146,8 @@ _wrap_new_Integer1dArray__SWIG_1(int argc, VALUE *argv, VALUE self) {
     SWIG_exception_fail(SWIG_ValueError, Ruby_Format_TypeError("invalid null reference ", "ap::template_1d_array< int > const &","ap::template_1d_array<(int)>", 1, argv[0])); 
   }
   arg1 = reinterpret_cast< ap::template_1d_array< int > * >(argp1);
-  result = (ap::template_1d_array< int > *)new ap::template_1d_array< int >((ap::template_1d_array< int > const &)*arg1);DATA_PTR(self) = result;
-  
+  result = (ap::template_1d_array< int > *)new ap::template_1d_array< int >((ap::template_1d_array< int > const &)*arg1);
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -2184,11 +2199,11 @@ SWIGINTERN VALUE
 _wrap_Integer1dArray___getitem__(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< int > *arg1 = (ap::template_1d_array< int > *) 0 ;
   int arg2 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -2226,13 +2241,13 @@ _wrap_Integer1dArray___setitem__(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< int > *arg1 = (ap::template_1d_array< int > *) 0 ;
   int arg2 ;
   int arg3 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
   int val3 ;
   int ecode3 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 2) || (argc > 2)) {
@@ -2273,9 +2288,9 @@ Size or Length of the Integer1dArray.
 SWIGINTERN VALUE
 _wrap_Integer1dArray_size(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< int > *arg1 = (ap::template_1d_array< int > *) 0 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -2408,9 +2423,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_Integer1dArray_getcontent(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< int > *arg1 = (ap::template_1d_array< int > *) 0 ;
-  int *result = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int *result = 0 ;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -2433,11 +2448,11 @@ SWIGINTERN VALUE
 _wrap_Integer1dArray_getlowbound__SWIG_0(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< int > *arg1 = (ap::template_1d_array< int > *) 0 ;
   int arg2 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -2464,9 +2479,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_Integer1dArray_getlowbound__SWIG_1(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< int > *arg1 = (ap::template_1d_array< int > *) 0 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -2534,11 +2549,11 @@ SWIGINTERN VALUE
 _wrap_Integer1dArray_gethighbound__SWIG_0(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< int > *arg1 = (ap::template_1d_array< int > *) 0 ;
   int arg2 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -2565,9 +2580,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_Integer1dArray_gethighbound__SWIG_1(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< int > *arg1 = (ap::template_1d_array< int > *) 0 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -2640,8 +2655,8 @@ _wrap_new_Real1dArray__SWIG_0(int argc, VALUE *argv, VALUE self) {
   if ((argc < 0) || (argc > 0)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
   }
-  result = (ap::template_1d_array< double,true > *)new ap::template_1d_array< double,true >();DATA_PTR(self) = result;
-  
+  result = (ap::template_1d_array< double,true > *)new ap::template_1d_array< double,true >();
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -2673,9 +2688,9 @@ _wrap_Real1dArray_allocate(VALUE self) {
 SWIGINTERN VALUE
 _wrap_new_Real1dArray__SWIG_1(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< double,true > *arg1 = 0 ;
-  ap::template_1d_array< double,true > *result = 0 ;
   void *argp1 ;
   int res1 = 0 ;
+  ap::template_1d_array< double,true > *result = 0 ;
   
   if ((argc < 1) || (argc > 1)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
@@ -2688,8 +2703,8 @@ _wrap_new_Real1dArray__SWIG_1(int argc, VALUE *argv, VALUE self) {
     SWIG_exception_fail(SWIG_ValueError, Ruby_Format_TypeError("invalid null reference ", "ap::template_1d_array< double,true > const &","ap::template_1d_array<(double,true)>", 1, argv[0])); 
   }
   arg1 = reinterpret_cast< ap::template_1d_array< double,true > * >(argp1);
-  result = (ap::template_1d_array< double,true > *)new ap::template_1d_array< double,true >((ap::template_1d_array< double,true > const &)*arg1);DATA_PTR(self) = result;
-  
+  result = (ap::template_1d_array< double,true > *)new ap::template_1d_array< double,true >((ap::template_1d_array< double,true > const &)*arg1);
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -2741,11 +2756,11 @@ SWIGINTERN VALUE
 _wrap_Real1dArray___getitem__(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< double,true > *arg1 = (ap::template_1d_array< double,true > *) 0 ;
   int arg2 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -2783,13 +2798,13 @@ _wrap_Real1dArray___setitem__(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< double,true > *arg1 = (ap::template_1d_array< double,true > *) 0 ;
   int arg2 ;
   double arg3 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
   double val3 ;
   int ecode3 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 2) || (argc > 2)) {
@@ -2830,9 +2845,9 @@ Size or Length of the Real1dArray.
 SWIGINTERN VALUE
 _wrap_Real1dArray_size(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< double,true > *arg1 = (ap::template_1d_array< double,true > *) 0 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -2965,9 +2980,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_Real1dArray_getcontent(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< double,true > *arg1 = (ap::template_1d_array< double,true > *) 0 ;
-  double *result = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double *result = 0 ;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -2990,11 +3005,11 @@ SWIGINTERN VALUE
 _wrap_Real1dArray_getlowbound__SWIG_0(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< double,true > *arg1 = (ap::template_1d_array< double,true > *) 0 ;
   int arg2 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -3021,9 +3036,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_Real1dArray_getlowbound__SWIG_1(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< double,true > *arg1 = (ap::template_1d_array< double,true > *) 0 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -3091,11 +3106,11 @@ SWIGINTERN VALUE
 _wrap_Real1dArray_gethighbound__SWIG_0(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< double,true > *arg1 = (ap::template_1d_array< double,true > *) 0 ;
   int arg2 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -3122,9 +3137,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_Real1dArray_gethighbound__SWIG_1(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< double,true > *arg1 = (ap::template_1d_array< double,true > *) 0 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -3197,8 +3212,8 @@ _wrap_new_Boolean1dArray__SWIG_0(int argc, VALUE *argv, VALUE self) {
   if ((argc < 0) || (argc > 0)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
   }
-  result = (ap::template_1d_array< bool > *)new ap::template_1d_array< bool >();DATA_PTR(self) = result;
-  
+  result = (ap::template_1d_array< bool > *)new ap::template_1d_array< bool >();
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -3230,9 +3245,9 @@ _wrap_Boolean1dArray_allocate(VALUE self) {
 SWIGINTERN VALUE
 _wrap_new_Boolean1dArray__SWIG_1(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< bool > *arg1 = 0 ;
-  ap::template_1d_array< bool > *result = 0 ;
   void *argp1 ;
   int res1 = 0 ;
+  ap::template_1d_array< bool > *result = 0 ;
   
   if ((argc < 1) || (argc > 1)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
@@ -3245,8 +3260,8 @@ _wrap_new_Boolean1dArray__SWIG_1(int argc, VALUE *argv, VALUE self) {
     SWIG_exception_fail(SWIG_ValueError, Ruby_Format_TypeError("invalid null reference ", "ap::template_1d_array< bool > const &","ap::template_1d_array<(bool)>", 1, argv[0])); 
   }
   arg1 = reinterpret_cast< ap::template_1d_array< bool > * >(argp1);
-  result = (ap::template_1d_array< bool > *)new ap::template_1d_array< bool >((ap::template_1d_array< bool > const &)*arg1);DATA_PTR(self) = result;
-  
+  result = (ap::template_1d_array< bool > *)new ap::template_1d_array< bool >((ap::template_1d_array< bool > const &)*arg1);
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -3298,11 +3313,11 @@ SWIGINTERN VALUE
 _wrap_Boolean1dArray___getitem__(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< bool > *arg1 = (ap::template_1d_array< bool > *) 0 ;
   int arg2 ;
-  bool result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  bool result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -3340,13 +3355,13 @@ _wrap_Boolean1dArray___setitem__(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< bool > *arg1 = (ap::template_1d_array< bool > *) 0 ;
   int arg2 ;
   bool arg3 ;
-  bool result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
   bool val3 ;
   int ecode3 = 0 ;
+  bool result;
   VALUE vresult = Qnil;
   
   if ((argc < 2) || (argc > 2)) {
@@ -3387,9 +3402,9 @@ Size or Length of the Boolean1dArray.
 SWIGINTERN VALUE
 _wrap_Boolean1dArray_size(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< bool > *arg1 = (ap::template_1d_array< bool > *) 0 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -3522,9 +3537,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_Boolean1dArray_getcontent(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< bool > *arg1 = (ap::template_1d_array< bool > *) 0 ;
-  bool *result = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  bool *result = 0 ;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -3547,11 +3562,11 @@ SWIGINTERN VALUE
 _wrap_Boolean1dArray_getlowbound__SWIG_0(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< bool > *arg1 = (ap::template_1d_array< bool > *) 0 ;
   int arg2 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -3578,9 +3593,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_Boolean1dArray_getlowbound__SWIG_1(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< bool > *arg1 = (ap::template_1d_array< bool > *) 0 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -3648,11 +3663,11 @@ SWIGINTERN VALUE
 _wrap_Boolean1dArray_gethighbound__SWIG_0(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< bool > *arg1 = (ap::template_1d_array< bool > *) 0 ;
   int arg2 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -3679,9 +3694,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_Boolean1dArray_gethighbound__SWIG_1(int argc, VALUE *argv, VALUE self) {
   ap::template_1d_array< bool > *arg1 = (ap::template_1d_array< bool > *) 0 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -3754,8 +3769,8 @@ _wrap_new_Integer2dArray__SWIG_0(int argc, VALUE *argv, VALUE self) {
   if ((argc < 0) || (argc > 0)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
   }
-  result = (ap::template_2d_array< int > *)new ap::template_2d_array< int >();DATA_PTR(self) = result;
-  
+  result = (ap::template_2d_array< int > *)new ap::template_2d_array< int >();
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -3787,9 +3802,9 @@ _wrap_Integer2dArray_allocate(VALUE self) {
 SWIGINTERN VALUE
 _wrap_new_Integer2dArray__SWIG_1(int argc, VALUE *argv, VALUE self) {
   ap::template_2d_array< int > *arg1 = 0 ;
-  ap::template_2d_array< int > *result = 0 ;
   void *argp1 ;
   int res1 = 0 ;
+  ap::template_2d_array< int > *result = 0 ;
   
   if ((argc < 1) || (argc > 1)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
@@ -3802,8 +3817,8 @@ _wrap_new_Integer2dArray__SWIG_1(int argc, VALUE *argv, VALUE self) {
     SWIG_exception_fail(SWIG_ValueError, Ruby_Format_TypeError("invalid null reference ", "ap::template_2d_array< int > const &","ap::template_2d_array<(int)>", 1, argv[0])); 
   }
   arg1 = reinterpret_cast< ap::template_2d_array< int > * >(argp1);
-  result = (ap::template_2d_array< int > *)new ap::template_2d_array< int >((ap::template_2d_array< int > const &)*arg1);DATA_PTR(self) = result;
-  
+  result = (ap::template_2d_array< int > *)new ap::template_2d_array< int >((ap::template_2d_array< int > const &)*arg1);
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -3847,13 +3862,13 @@ _wrap_Integer2dArray___call__(int argc, VALUE *argv, VALUE self) {
   ap::template_2d_array< int > *arg1 = (ap::template_2d_array< int > *) 0 ;
   int arg2 ;
   int arg3 ;
-  int *result = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
   int val3 ;
   int ecode3 = 0 ;
+  int *result = 0 ;
   VALUE vresult = Qnil;
   
   if ((argc < 2) || (argc > 2)) {
@@ -4040,11 +4055,11 @@ SWIGINTERN VALUE
 _wrap_Integer2dArray_getlowbound(int argc, VALUE *argv, VALUE self) {
   ap::template_2d_array< int > *arg1 = (ap::template_2d_array< int > *) 0 ;
   int arg2 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -4072,11 +4087,11 @@ SWIGINTERN VALUE
 _wrap_Integer2dArray_gethighbound(int argc, VALUE *argv, VALUE self) {
   ap::template_2d_array< int > *arg1 = (ap::template_2d_array< int > *) 0 ;
   int arg2 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -4114,13 +4129,13 @@ _wrap_Integer2dArray___getitem__(int argc, VALUE *argv, VALUE self) {
   ap::template_2d_array< int > *arg1 = (ap::template_2d_array< int > *) 0 ;
   int arg2 ;
   int arg3 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
   int val3 ;
   int ecode3 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 2) || (argc > 2)) {
@@ -4212,8 +4227,8 @@ _wrap_new_Real2dArray__SWIG_0(int argc, VALUE *argv, VALUE self) {
   if ((argc < 0) || (argc > 0)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
   }
-  result = (ap::template_2d_array< double,true > *)new ap::template_2d_array< double,true >();DATA_PTR(self) = result;
-  
+  result = (ap::template_2d_array< double,true > *)new ap::template_2d_array< double,true >();
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -4245,9 +4260,9 @@ _wrap_Real2dArray_allocate(VALUE self) {
 SWIGINTERN VALUE
 _wrap_new_Real2dArray__SWIG_1(int argc, VALUE *argv, VALUE self) {
   ap::template_2d_array< double,true > *arg1 = 0 ;
-  ap::template_2d_array< double,true > *result = 0 ;
   void *argp1 ;
   int res1 = 0 ;
+  ap::template_2d_array< double,true > *result = 0 ;
   
   if ((argc < 1) || (argc > 1)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
@@ -4260,8 +4275,8 @@ _wrap_new_Real2dArray__SWIG_1(int argc, VALUE *argv, VALUE self) {
     SWIG_exception_fail(SWIG_ValueError, Ruby_Format_TypeError("invalid null reference ", "ap::template_2d_array< double,true > const &","ap::template_2d_array<(double,true)>", 1, argv[0])); 
   }
   arg1 = reinterpret_cast< ap::template_2d_array< double,true > * >(argp1);
-  result = (ap::template_2d_array< double,true > *)new ap::template_2d_array< double,true >((ap::template_2d_array< double,true > const &)*arg1);DATA_PTR(self) = result;
-  
+  result = (ap::template_2d_array< double,true > *)new ap::template_2d_array< double,true >((ap::template_2d_array< double,true > const &)*arg1);
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -4305,13 +4320,13 @@ _wrap_Real2dArray___call__(int argc, VALUE *argv, VALUE self) {
   ap::template_2d_array< double,true > *arg1 = (ap::template_2d_array< double,true > *) 0 ;
   int arg2 ;
   int arg3 ;
-  double *result = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
   int val3 ;
   int ecode3 = 0 ;
+  double *result = 0 ;
   VALUE vresult = Qnil;
   
   if ((argc < 2) || (argc > 2)) {
@@ -4498,11 +4513,11 @@ SWIGINTERN VALUE
 _wrap_Real2dArray_getlowbound(int argc, VALUE *argv, VALUE self) {
   ap::template_2d_array< double,true > *arg1 = (ap::template_2d_array< double,true > *) 0 ;
   int arg2 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -4530,11 +4545,11 @@ SWIGINTERN VALUE
 _wrap_Real2dArray_gethighbound(int argc, VALUE *argv, VALUE self) {
   ap::template_2d_array< double,true > *arg1 = (ap::template_2d_array< double,true > *) 0 ;
   int arg2 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -4572,13 +4587,13 @@ _wrap_Real2dArray___getitem__(int argc, VALUE *argv, VALUE self) {
   ap::template_2d_array< double,true > *arg1 = (ap::template_2d_array< double,true > *) 0 ;
   int arg2 ;
   int arg3 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
   int val3 ;
   int ecode3 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 2) || (argc > 2)) {
@@ -4664,9 +4679,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_real1d_to_array(int argc, VALUE *argv, VALUE self) {
   ap::real_1d_array *arg1 = (ap::real_1d_array *) 0 ;
-  VALUE result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  VALUE result;
   VALUE vresult = Qnil;
   
   if ((argc < 1) || (argc > 1)) {
@@ -4713,9 +4728,9 @@ _wrap_pearsoncorrelation(int argc, VALUE *argv, VALUE self) {
   ap::real_1d_array *arg1 = 0 ;
   ap::real_1d_array *arg2 = 0 ;
   int arg3 ;
-  double result;
   int val3 ;
   int ecode3 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 3) || (argc > 3)) {
@@ -4753,9 +4768,9 @@ _wrap_spearmanrankcorrelation(int argc, VALUE *argv, VALUE self) {
   ap::real_1d_array *arg1 = 0 ;
   ap::real_1d_array *arg2 = 0 ;
   int arg3 ;
-  double result;
   int val3 ;
   int ecode3 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 3) || (argc > 3)) {
@@ -4792,11 +4807,11 @@ SWIGINTERN VALUE
 _wrap_pearsoncorrelationsignificance(int argc, VALUE *argv, VALUE self) {
   double arg1 ;
   int arg2 ;
-  VALUE result;
   double val1 ;
   int ecode1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  VALUE result;
   VALUE vresult = Qnil;
   
   if ((argc < 2) || (argc > 2)) {
@@ -4824,11 +4839,11 @@ SWIGINTERN VALUE
 _wrap_spearmanrankcorrelationsignificance(int argc, VALUE *argv, VALUE self) {
   double arg1 ;
   int arg2 ;
-  VALUE result;
   double val1 ;
   int ecode1 = 0 ;
   int val2 ;
   int ecode2 = 0 ;
+  VALUE result;
   VALUE vresult = Qnil;
   
   if ((argc < 2) || (argc > 2)) {
@@ -4877,7 +4892,6 @@ _wrap_LinearModel_w_set(int argc, VALUE *argv, VALUE self) {
   }
   arg2 = reinterpret_cast< ap::real_1d_array * >(argp2);
   if (arg1) (arg1)->w = *arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -4887,9 +4901,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LinearModel_w_get(int argc, VALUE *argv, VALUE self) {
   linearmodel *arg1 = (linearmodel *) 0 ;
-  ap::real_1d_array *result = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  ap::real_1d_array *result = 0 ;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -4932,8 +4946,8 @@ _wrap_new_LinearModel(int argc, VALUE *argv, VALUE self) {
   if ((argc < 0) || (argc > 0)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
   }
-  result = (linearmodel *)new linearmodel();DATA_PTR(self) = result;
-  
+  result = (linearmodel *)new linearmodel();
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -4970,7 +4984,6 @@ _wrap_LrReport_c_set(int argc, VALUE *argv, VALUE self) {
   }
   arg2 = reinterpret_cast< ap::real_2d_array * >(argp2);
   if (arg1) (arg1)->c = *arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -4980,9 +4993,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LrReport_c_get(int argc, VALUE *argv, VALUE self) {
   lrreport *arg1 = (lrreport *) 0 ;
-  ap::real_2d_array *result = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  ap::real_2d_array *result = 0 ;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -5024,7 +5037,6 @@ _wrap_LrReport_rmserror_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->rmserror = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -5034,9 +5046,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LrReport_rmserror_get(int argc, VALUE *argv, VALUE self) {
   lrreport *arg1 = (lrreport *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -5078,7 +5090,6 @@ _wrap_LrReport_avgerror_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->avgerror = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -5088,9 +5099,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LrReport_avgerror_get(int argc, VALUE *argv, VALUE self) {
   lrreport *arg1 = (lrreport *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -5132,7 +5143,6 @@ _wrap_LrReport_avgrelerror_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->avgrelerror = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -5142,9 +5152,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LrReport_avgrelerror_get(int argc, VALUE *argv, VALUE self) {
   lrreport *arg1 = (lrreport *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -5186,7 +5196,6 @@ _wrap_LrReport_cvrmserror_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->cvrmserror = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -5196,9 +5205,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LrReport_cvrmserror_get(int argc, VALUE *argv, VALUE self) {
   lrreport *arg1 = (lrreport *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -5240,7 +5249,6 @@ _wrap_LrReport_cvavgerror_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->cvavgerror = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -5250,9 +5258,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LrReport_cvavgerror_get(int argc, VALUE *argv, VALUE self) {
   lrreport *arg1 = (lrreport *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -5294,7 +5302,6 @@ _wrap_LrReport_cvavgrelerror_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->cvavgrelerror = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -5304,9 +5311,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LrReport_cvavgrelerror_get(int argc, VALUE *argv, VALUE self) {
   lrreport *arg1 = (lrreport *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -5348,7 +5355,6 @@ _wrap_LrReport_ncvdefects_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< int >(val2);
   if (arg1) (arg1)->ncvdefects = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -5358,9 +5364,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LrReport_ncvdefects_get(int argc, VALUE *argv, VALUE self) {
   lrreport *arg1 = (lrreport *) 0 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -5402,7 +5408,6 @@ _wrap_LrReport_cvdefects_set(int argc, VALUE *argv, VALUE self) {
   }
   arg2 = reinterpret_cast< ap::integer_1d_array * >(argp2);
   if (arg1) (arg1)->cvdefects = *arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -5412,9 +5417,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LrReport_cvdefects_get(int argc, VALUE *argv, VALUE self) {
   lrreport *arg1 = (lrreport *) 0 ;
-  ap::integer_1d_array *result = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  ap::integer_1d_array *result = 0 ;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -5457,8 +5462,8 @@ _wrap_new_LrReport(int argc, VALUE *argv, VALUE self) {
   if ((argc < 0) || (argc > 0)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
   }
-  result = (lrreport *)new lrreport();DATA_PTR(self) = result;
-  
+  result = (lrreport *)new lrreport();
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -5872,9 +5877,9 @@ SWIGINTERN VALUE
 _wrap_lrprocess(int argc, VALUE *argv, VALUE self) {
   linearmodel *arg1 = 0 ;
   ap::real_1d_array *arg2 = 0 ;
-  double result;
   void *argp1 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 2) || (argc > 2)) {
@@ -5903,6 +5908,96 @@ fail:
 }
 
 
+SWIGINTERN VALUE
+_wrap_mannwhitneyutest(int argc, VALUE *argv, VALUE self) {
+  ap::real_1d_array *arg1 = 0 ;
+  int arg2 ;
+  ap::real_1d_array *arg3 = 0 ;
+  int arg4 ;
+  double *arg5 = 0 ;
+  double *arg6 = 0 ;
+  double *arg7 = 0 ;
+  int val2 ;
+  int ecode2 = 0 ;
+  int val4 ;
+  int ecode4 = 0 ;
+  double temp5 ;
+  int res5 = 0 ;
+  double temp6 ;
+  int res6 = 0 ;
+  double temp7 ;
+  int res7 = 0 ;
+  
+  if ((argc < 7) || (argc > 7)) {
+    rb_raise(rb_eArgError, "wrong # of arguments(%d for 7)",argc); SWIG_fail;
+  }
+  {
+    arg1 = new ap::real_1d_array();
+    arg1->setlength(RARRAY_LEN(argv[0]));
+    for(long i=0;i<RARRAY_LEN(argv[0]);i++) {
+      arg1->operator()(i)=rb_num2dbl(RARRAY_PTR(argv[0])[i]);
+    }
+  }
+  ecode2 = SWIG_AsVal_int(argv[1], &val2);
+  if (!SWIG_IsOK(ecode2)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode2), Ruby_Format_TypeError( "", "int","mannwhitneyutest", 2, argv[1] ));
+  } 
+  arg2 = static_cast< int >(val2);
+  {
+    arg3 = new ap::real_1d_array();
+    arg3->setlength(RARRAY_LEN(argv[2]));
+    for(long i=0;i<RARRAY_LEN(argv[2]);i++) {
+      arg3->operator()(i)=rb_num2dbl(RARRAY_PTR(argv[2])[i]);
+    }
+  }
+  ecode4 = SWIG_AsVal_int(argv[3], &val4);
+  if (!SWIG_IsOK(ecode4)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode4), Ruby_Format_TypeError( "", "int","mannwhitneyutest", 4, argv[3] ));
+  } 
+  arg4 = static_cast< int >(val4);
+  if (!(SWIG_IsOK((res5 = SWIG_ConvertPtr(argv[4],SWIG_as_voidptrptr(&arg5),SWIGTYPE_p_double,0))))) {
+    double val; 
+    int ecode = SWIG_AsVal_double(argv[4], &val);
+    if (!SWIG_IsOK(ecode)) {
+      SWIG_exception_fail(SWIG_ArgError(ecode), Ruby_Format_TypeError( "", "double","mannwhitneyutest", 5, argv[4] ));
+    }
+    temp5 = static_cast< double >(val);
+    arg5 = &temp5;
+    res5 = SWIG_AddTmpMask(ecode);
+  }
+  if (!(SWIG_IsOK((res6 = SWIG_ConvertPtr(argv[5],SWIG_as_voidptrptr(&arg6),SWIGTYPE_p_double,0))))) {
+    double val; 
+    int ecode = SWIG_AsVal_double(argv[5], &val);
+    if (!SWIG_IsOK(ecode)) {
+      SWIG_exception_fail(SWIG_ArgError(ecode), Ruby_Format_TypeError( "", "double","mannwhitneyutest", 6, argv[5] ));
+    }
+    temp6 = static_cast< double >(val);
+    arg6 = &temp6;
+    res6 = SWIG_AddTmpMask(ecode);
+  }
+  if (!(SWIG_IsOK((res7 = SWIG_ConvertPtr(argv[6],SWIG_as_voidptrptr(&arg7),SWIGTYPE_p_double,0))))) {
+    double val; 
+    int ecode = SWIG_AsVal_double(argv[6], &val);
+    if (!SWIG_IsOK(ecode)) {
+      SWIG_exception_fail(SWIG_ArgError(ecode), Ruby_Format_TypeError( "", "double","mannwhitneyutest", 7, argv[6] ));
+    }
+    temp7 = static_cast< double >(val);
+    arg7 = &temp7;
+    res7 = SWIG_AddTmpMask(ecode);
+  }
+  mannwhitneyutest((ap::template_1d_array< double,true > const &)*arg1,arg2,(ap::template_1d_array< double,true > const &)*arg3,arg4,*arg5,*arg6,*arg7);
+  if (SWIG_IsNewObj(res5)) delete arg5;
+  if (SWIG_IsNewObj(res6)) delete arg6;
+  if (SWIG_IsNewObj(res7)) delete arg7;
+  return Qnil;
+fail:
+  if (SWIG_IsNewObj(res5)) delete arg5;
+  if (SWIG_IsNewObj(res6)) delete arg6;
+  if (SWIG_IsNewObj(res7)) delete arg7;
+  return Qnil;
+}
+
+
 swig_class cLogitModel;
 
 SWIGINTERN VALUE
@@ -5928,7 +6023,6 @@ _wrap_LogitModel_w_set(int argc, VALUE *argv, VALUE self) {
   }
   arg2 = reinterpret_cast< ap::real_1d_array * >(argp2);
   if (arg1) (arg1)->w = *arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -5938,9 +6032,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitModel_w_get(int argc, VALUE *argv, VALUE self) {
   logitmodel *arg1 = (logitmodel *) 0 ;
-  ap::real_1d_array *result = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  ap::real_1d_array *result = 0 ;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -5983,8 +6077,8 @@ _wrap_new_LogitModel(int argc, VALUE *argv, VALUE self) {
   if ((argc < 0) || (argc > 0)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
   }
-  result = (logitmodel *)new logitmodel();DATA_PTR(self) = result;
-  
+  result = (logitmodel *)new logitmodel();
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -6021,7 +6115,6 @@ _wrap_LogitMcState_brackt_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< bool >(val2);
   if (arg1) (arg1)->brackt = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6031,9 +6124,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_brackt_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  bool result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  bool result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6075,7 +6168,6 @@ _wrap_LogitMcState_stage1_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< bool >(val2);
   if (arg1) (arg1)->stage1 = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6085,9 +6177,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_stage1_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  bool result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  bool result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6129,7 +6221,6 @@ _wrap_LogitMcState_infoc_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< int >(val2);
   if (arg1) (arg1)->infoc = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6139,9 +6230,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_infoc_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6183,7 +6274,6 @@ _wrap_LogitMcState_dg_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->dg = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6193,9 +6283,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_dg_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6237,7 +6327,6 @@ _wrap_LogitMcState_dgm_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->dgm = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6247,9 +6336,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_dgm_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6291,7 +6380,6 @@ _wrap_LogitMcState_dginit_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->dginit = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6301,9 +6389,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_dginit_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6345,7 +6433,6 @@ _wrap_LogitMcState_dgtest_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->dgtest = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6355,9 +6442,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_dgtest_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6399,7 +6486,6 @@ _wrap_LogitMcState_dgx_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->dgx = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6409,9 +6495,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_dgx_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6453,7 +6539,6 @@ _wrap_LogitMcState_dgxm_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->dgxm = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6463,9 +6548,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_dgxm_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6507,7 +6592,6 @@ _wrap_LogitMcState_dgy_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->dgy = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6517,9 +6601,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_dgy_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6561,7 +6645,6 @@ _wrap_LogitMcState_dgym_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->dgym = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6571,9 +6654,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_dgym_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6615,7 +6698,6 @@ _wrap_LogitMcState_finit_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->finit = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6625,9 +6707,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_finit_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6669,7 +6751,6 @@ _wrap_LogitMcState_ftest1_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->ftest1 = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6679,9 +6760,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_ftest1_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6723,7 +6804,6 @@ _wrap_LogitMcState_fm_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->fm = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6733,9 +6813,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_fm_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6777,7 +6857,6 @@ _wrap_LogitMcState_fx_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->fx = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6787,9 +6866,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_fx_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6831,7 +6910,6 @@ _wrap_LogitMcState_fxm_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->fxm = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6841,9 +6919,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_fxm_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6885,7 +6963,6 @@ _wrap_LogitMcState_fy_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->fy = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6895,9 +6972,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_fy_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6939,7 +7016,6 @@ _wrap_LogitMcState_fym_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->fym = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -6949,9 +7025,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_fym_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -6993,7 +7069,6 @@ _wrap_LogitMcState_stx_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->stx = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -7003,9 +7078,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_stx_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -7047,7 +7122,6 @@ _wrap_LogitMcState_sty_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->sty = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -7057,9 +7131,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_sty_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -7101,7 +7175,6 @@ _wrap_LogitMcState_stmin_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->stmin = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -7111,9 +7184,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_stmin_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -7155,7 +7228,6 @@ _wrap_LogitMcState_stmax_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->stmax = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -7165,9 +7237,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_stmax_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -7209,7 +7281,6 @@ _wrap_LogitMcState_width_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->width = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -7219,9 +7290,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_width_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -7263,7 +7334,6 @@ _wrap_LogitMcState_width1_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->width1 = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -7273,9 +7343,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_width1_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -7317,7 +7387,6 @@ _wrap_LogitMcState_xtrapf_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< double >(val2);
   if (arg1) (arg1)->xtrapf = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -7327,9 +7396,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_LogitMcState_xtrapf_get(int argc, VALUE *argv, VALUE self) {
   logitmcstate *arg1 = (logitmcstate *) 0 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -7372,8 +7441,8 @@ _wrap_new_LogitMcState(int argc, VALUE *argv, VALUE self) {
   if ((argc < 0) || (argc > 0)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
   }
-  result = (logitmcstate *)new logitmcstate();DATA_PTR(self) = result;
-  
+  result = (logitmcstate *)new logitmcstate();
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -7410,7 +7479,6 @@ _wrap_MNLReport_ngrad_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< int >(val2);
   if (arg1) (arg1)->ngrad = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -7420,9 +7488,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_MNLReport_ngrad_get(int argc, VALUE *argv, VALUE self) {
   mnlreport *arg1 = (mnlreport *) 0 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -7464,7 +7532,6 @@ _wrap_MNLReport_nhess_set(int argc, VALUE *argv, VALUE self) {
   } 
   arg2 = static_cast< int >(val2);
   if (arg1) (arg1)->nhess = arg2;
-  
   return Qnil;
 fail:
   return Qnil;
@@ -7474,9 +7541,9 @@ fail:
 SWIGINTERN VALUE
 _wrap_MNLReport_nhess_get(int argc, VALUE *argv, VALUE self) {
   mnlreport *arg1 = (mnlreport *) 0 ;
-  int result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  int result;
   VALUE vresult = Qnil;
   
   if ((argc < 0) || (argc > 0)) {
@@ -7519,8 +7586,8 @@ _wrap_new_MNLReport(int argc, VALUE *argv, VALUE self) {
   if ((argc < 0) || (argc > 0)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
   }
-  result = (mnlreport *)new mnlreport();DATA_PTR(self) = result;
-  
+  result = (mnlreport *)new mnlreport();
+  DATA_PTR(self) = result;
   return self;
 fail:
   return Qnil;
@@ -7617,9 +7684,9 @@ SWIGINTERN VALUE
 _wrap_mnlprocess(int argc, VALUE *argv, VALUE self) {
   logitmodel *arg1 = 0 ;
   ap::real_1d_array *arg2 = 0 ;
-  VALUE result;
   void *argp1 = 0 ;
   int res1 = 0 ;
+  VALUE result;
   VALUE vresult = Qnil;
   
   if ((argc < 2) || (argc > 2)) {
@@ -7882,13 +7949,13 @@ _wrap_mnlavgce(int argc, VALUE *argv, VALUE self) {
   logitmodel *arg1 = 0 ;
   ap::real_2d_array *arg2 = 0 ;
   int arg3 ;
-  double result;
   void *argp1 = 0 ;
   int res1 = 0 ;
   void *argp2 ;
   int res2 = 0 ;
   int val3 ;
   int ecode3 = 0 ;
+  double result;
   VALUE vresult = Qnil;
   
   if ((argc < 3) || (argc > 3)) {
@@ -8465,6 +8532,7 @@ SWIGEXPORT void Init_alglib_ext(void) {
   rb_define_module_function(mAlglib_ext, "lrunpack", VALUEFUNC(_wrap_lrunpack), -1);
   rb_define_module_function(mAlglib_ext, "lrpack", VALUEFUNC(_wrap_lrpack), -1);
   rb_define_module_function(mAlglib_ext, "lrprocess", VALUEFUNC(_wrap_lrprocess), -1);
+  rb_define_module_function(mAlglib_ext, "mannwhitneyutest", VALUEFUNC(_wrap_mannwhitneyutest), -1);
   
   cLogitModel.klass = rb_define_class_under(mAlglib_ext, "LogitModel", rb_cObject);
   SWIG_TypeClientData(SWIGTYPE_p_logitmodel, (void *) &cLogitModel);