On Saturday 17 September 2005 20:25, Sergey Plis wrote: > Thank you for the help. I am not sure if the change is necessary as > long as what is implemented now works everything should be fine. > Especially, taking into account that the definition happens only once > and the usage semantics does not change and is very convenient. > > However, the problem described in my next posted topic still > persists. Also now I get the following error when calling ]]* > > a[[ b[[ ]]* bigforth: symbol lookup error: /usr/lib/libgsl.so.0: > undefined symbol: cblas_dgemm
Strange. What you do is correct. It works here. The attached test stuff prints out a multiplied matrix, with the correct values. Maybe you need to change libgslblas.so with libgslblas.so.0, but both works here (libgsl from SuSE 9.3). I get the error only if I omit the depends libgslblas part. -- Bernd Paysan "If you want it done right, you have to do it yourself" http://www.jwdt.com/~paysan/
import float float also
Module GSL
also DOS
library libgslcblas libgslcblas.so
library libgsl libgsl.so.0 depends libgslcblas
legacy off
libgsl gsl_matrix_alloc int int (int) gsl_matrix_alloc
( n m -- *gsl_matrix )
libgsl gsl_matrix_calloc int int (int) gsl_matrix_calloc
( n m -- *gsl_matrix )
libgsl gsl_matrix_free int (void) gsl_matrix_free
( *gsl_matrix -- )
libgsl gsl_matrix_get int int int (fp) gsl_matrix_get
( *gsl_matrix i j -- df )
\ ******* this is the function that gives trouble **********
libgsl gsl_matrix_set df int int int (void/fp) gsl_matrix_set
( df *gsl_matrix i j -- )
libgsl gsl_matrix_set_all df int (void/fp) gsl_matrix_set_all
( df *gsl_matrix i j -- )
libgsl gsl_blas_dgemm int df int int df int int (int/fp) gsl_blas_dgemm
( alpha beta int int *gsl_matrix *gsl_matrix *gsl_matrix -- n )
struct{ cell size1 cell size2 } gsl_matrix
111 Constant CblasNoTrans
: ]]* ( *gsl_matrix *gsl_matrix -- *gsl_matrix )
!1 !0 CblasNoTrans CblasNoTrans 2swap 2dup
gsl_matrix size2 @ swap gsl_matrix size1 @ swap
gsl_matrix_alloc dup >r
gsl_blas_dgemm drop r> ;
legacy on
previous
Module;
gsl also
3 3 gsl_matrix_alloc constant a[[
3 3 gsl_matrix_alloc constant b[[
3e a[[ gsl_matrix_set_all
2e b[[ gsl_matrix_set_all
a[[ b[[ ]]* constant c[[
9 0 [dO] [i] 3 mod 0= [if] cr [thEN] c[[ [i] 3 /mod gsl_matrix_get f. [LOOP]
pgpDyY7k0ETqK.pgp
Description: PGP signature
