Menu ▾ ▴

GSL: Crashes sometimes under VC 7.1

Help
Ulrich
2006-04-01
2012-07-26

  • Ulrich

    Ulrich - 2006-04-01

    Hi,

    I am trying to run GSL 1.6 under VC 7.1. Most of the examples from the GSL manual work just fine, but some (naturally the ones I need) crash.

    Below is the example from Section 30.4 on Wavelet Transforms that crashes with "test.exe has encountered a problem and needs to close. We are sorry for the inconvenience."

    I compile and link as follows (there are no errors while compiling or linking):

    cl /c /DGSL_DLL /I "C:\Program Files\GnuWin32\include" test.c
    link /LIBPATH:"C:\Program Files\GnuWin32\lib" test.obj libgsl.lib

    Here is the error signature:

    AppName: test.exe AppVer: 0.0.0.0 ModName: libgsl.dll
    ModVer: 0.0.0.0 Offset: 000a7ece

    Any suggestions would be more than welcome.

    Thanks!!!
    Ulrich

    test.c contains:

    include <stdio.h>

     #include &lt;math.h&gt;
 #include &lt;gsl/gsl_sort.h&gt;
 #include &lt;gsl/gsl_wavelet.h&gt;

 int
 main (int argc, char **argv)
 {
   int i, n = 256, nc = 20;
   double *data = malloc (n * sizeof (double));
   double *abscoeff = malloc (n * sizeof (double));
   size_t *p = malloc (n * sizeof (size_t));

   FILE *f = fopen (argv[1], &quot;r&quot;);
   for (i = 0; i &lt; n; i++)
     {
       fscanf (f, &quot;%lg&quot;, &amp;data[i]);
     }
   fclose (f);

   {
     gsl_wavelet *w = gsl_wavelet_alloc (gsl_wavelet_daubechies, 4);
     gsl_wavelet_workspace *work = gsl_wavelet_workspace_alloc (n);

     gsl_wavelet_transform_forward (w, data, 1, n, work);

     for (i = 0; i &lt; n; i++)
       {
         abscoeff[i] = fabs (data[i]);
       }

     gsl_sort_index (p, abscoeff, 1, n);

     for (i = 0; (i + nc) &lt; n; i++)
       data[p[i]] = 0;

     gsl_wavelet_transform_inverse (w, data, 1, n, work);
   }

   for (i = 0; i &lt; n; i++)
     {
       printf (&quot;%g\n&quot;, data[i]);
     }
 }
     

    • Jerry S.

      Jerry S. - 2006-04-11

      The following works for me with the distributed files in GnuWin32 gsl-1.6 using VC7 or VC8 from the command line:

      1. use lib to create import library libgsl.lib from libgsl.def (in gsl-1.6-1-lib/lib):

        lib /DEF:libgsl.def

      2. compile the source with BOTH GSL_DLL and WIN32 defined:

        cl /c /I <path-to-gsl-headers> /DGSL_DLL /DWIN32 test.c

      3. link:

        link /LIBPATH:<path-to-libgsl.lib> test.obj libgsl.lib

      Proper definition of the import specification requires that both GSL_DLL and WIN32 be defined when gsl/gsl_types.h is parsed. This is a trap for VC users since WIN32 is not among the predefined macros. We'll look at modifying this in gl-1.7 so only GSL_DLL is required for VC also.

      If this doesn't solve your problem, please continue the conversation.

      Thanks again,
      Jerry

       

    • Jerry S.

      Jerry S. - 2006-04-04

      Ulrich;

      Thank you very much for the feedback. I have not seen this problem before, and am trying to reproduce it.

      Please provide more information about the gsl library you're using - did you compile the dll from scratch with msvc, or did you generate an import library and are using the dll in the binary distribution?

      Also, could you please provide the data file you're feeding to the wavelet transform? Either post it here or email to <jds2 at users dot sourceforge dot net>.

      Thanks,
      Jerry

       
MongoDB Logo MongoDB
Gen AI apps are built with MongoDB Atlas
Atlas offers built-in vector search and global availability across 125+ regions. Start building AI apps faster, all in one place.
Try Free →