Dear Pierre,

As I said earlier the dimension check as it is implemented is not
meaningful. So I disabled it and added the solver_doc.ic to the patch
file.

Some tests still fail. Also find attached results of the tests. Test
that failed are
block_test.py
interpolation_test.py
statistics_test.py

I am currently interested in multimininization. Does any of these test
failures affect multimin routines?

cheers,
Krishna.

On 5/31/06, Krishna Mohan Gundu <[EMAIL PROTECTED]> wrote:
Dear Pierre,

I went through the code and identified the problem. I don't know the
solution but it should give you a hint.

src/init/block_helpers.c : 94
          if(dimensions[i] == -1 ){

this statement fails eventhough dimensions[i] is equal to -1 although
both are of the same type, dont ask me why!! (NOTE: in error messages
we have the number 4294967295 = ulong32(-1)). Even if this succeeds we
get an error saying "No one cares about its first dimension!" That
leads to the following code.

typemaps/block_conversion_functions.h: 37-38
        ((pyvector =  PyGSL_PyArray_PREPARE_gsl_vector_view(input,
basetype ## _py_array_type, \

PyGSL_NON_CONTIGUOUS | flag, -1, argnum, NULL)) != NULL)

I do not understand why you are passing "-1" as the size of the vector
to be compared. Can you explain?

cheers,
Krishna.

--- /dev/null	2004-09-18 16:50:22.000000000 -0700
+++ testing/src/solvers/solver_doc.ic	2006-06-01 00:02:08.000000000 -0700
@@ -0,0 +1,13 @@
+/* -*- C -*- */
+#ifndef SOLVER_DOC_IC
+#define SOLVER_DOC_IC 1
+static const char PyGSL_solver_module_doc[] = "XXX Missing";
+static const char PyGSL_solver_type_doc[] = "XXX Missing";
+
+
+
+
+
+
+
+#endif /* SOLVER_DOC_IC */
--- src/init/.block_helpers.c~	2006-02-23 05:03:18.000000000 -0700
+++ src/init/block_helpers.c	2006-06-01 00:09:31.000000000 -0700
@@ -89,7 +89,7 @@
      }
      
 
-     
+#if 0
      for(i=0; i<nd; ++i){
 	  if(dimensions[i] == -1 ){
 	       switch(i){
@@ -114,6 +114,7 @@
 	       goto fail;
 	  }
      }
+#endif
 
      if(  ((PyArrayObject *) (a_array))->data == NULL){
 	  gsl_error("Got an array object were the data was NULL!", filename, __LINE__, GSL_ESANITY);

Attachment: pygsl-tests-2006-07-01.res
Description: Binary data

Reply via email to