> Can you make a patch for the level-2 cblas functions only, I'm
> reasonably happy with those but not sure about the level-3.
Ok, but I think it will take a while. Now I don't have much time.
> I think it would be better to move the pos argument inside the macro, as it
> shouldn't be needed outside it (example below).
I prefer maintain the pos argument outside (I think this was
discussed). If pos is outside we could use the same macro in the cblas
functions and in the GSL wreppers of cblas to check the errors and
return a GSL error code. The "problem" with cblas is that an incorrect
input argument STOPS the execution of the program (in GSL this will
cause memory leaks, for example).
> #define CBLAS_ERROR_GEMV(order,TransA,M,N,alpha,A,lda,X,incX,beta,Y,incY) \
> { int pos = 0; \
> CHECK_ORDER(pos,1,order); \
> CHECK_TRANSPOSE(pos,2,TransA); \
> CHECK_DIM(pos,3,M); \
> CHECK_DIM(pos,4,N); \
> if((order)==CblasRowMajor) { \
> if((lda)<CBLAS_MAX(1,N)) { \
> pos = 7; \
> } \
> } else if((order)==CblasColMajor) { \
> if((lda)<CBLAS_MAX(1,M)) { \
> pos = 7; \
> } \
> } \
> CHECK_STRIDE(pos,9,incX); \
> CHECK_STRIDE(pos,12,incY); \
> REPORT_ERROR(pos) ; }
>
--
*****************************************
José Luis García Pallero
[email protected]
(o<
/ / \
V_/_
Use Debian GNU/Linux and enjoy!
*****************************************
_______________________________________________
Help-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/help-gsl
