Changeset 40

Show
Ignore:
Timestamp:
08/28/07 18:01:02 (1 year ago)
Author:
pernet
Message:

(hopefully) last bug fix for ftrsm ftrmm new algorithm, + linbox related tweaks.

Location:
include
Files:
6 modified

Legend:

Unmodified
Added
Removed
  • include/config-blas.h

    r18 r40  
    2222 
    2323 
    24 #ifndef __LINBOX_CONFIG_BLAS 
    25 #define __LINBOX_CONFIG_BLAS 
     24#ifndef _LINBOX_CONFIG_BLAS 
     25#define _LINBOX_CONFIG_BLAS 
    2626 
    2727// #ifndef __LINBOX_CONFIGURATION 
  • include/fflas-ffpack/fflas.h

    r38 r40  
    2121#endif 
    2222 
    23 #ifdef __LINBOX_CONFIG_H  
     23#ifdef _LINBOX_CONFIG_H  
    2424#include "linbox/config-blas.h" 
    2525#include "linbox/field/unparametric.h" 
  • include/fflas-ffpack/fflas_bounds.inl

    r36 r40  
    99 */ 
    1010 
    11 #ifdef __LINBOX_CONFIG_H 
     11#ifdef _LINBOX_CONFIG_H 
    1212#define FFLAS_INT_TYPE Integer 
    1313#else 
     
    4242                        //FFLAS_INT_TYPE c = (p-1)*(ex)/2; //bound for a centered representation 
    4343                        long long c; 
    44 #ifndef __LINBOX_CONFIG_H 
     44#ifndef _LINBOX_CONFIG_H 
    4545                        if (F.balanced) 
    4646                                c = (p-1)*(ex)/2; // balanced representation 
     
    176176                F.characteristic(pi); 
    177177                static FFLAS_INT_TYPE p=pi; 
    178 #ifdef __LINBOX_CONFIG_H 
     178#ifdef _LINBOX_CONFIG_H 
    179179                static size_t nmax = bound_compute_double(pi); 
    180180#else 
     
    184184                        return nmax; 
    185185                else 
    186 #ifdef __LINBOX_CONFIG_H 
     186#ifdef _LINBOX_CONFIG_H 
    187187                        return nmax= bound_compute_double (p=pi); //(F.balanced) ? bound_compute_balanced(p=pi) : bound_compute(p=pi); 
    188188#else 
     
    200200                F.characteristic(pi); 
    201201                static FFLAS_INT_TYPE p=pi; 
    202 #ifdef __LINBOX_CONFIG_H 
     202#ifdef _LINBOX_CONFIG_H 
    203203                static size_t nmax = bound_compute_float(pi); 
    204204#else 
     
    208208                        return nmax; 
    209209                else 
    210 #ifdef __LINBOX_CONFIG_H 
     210#ifdef _LINBOX_CONFIG_H 
    211211                        return nmax= bound_compute_float (p=pi); //(F.balanced) ? bound_compute_balanced(p=pi) : bound_compute(p=pi); 
    212212#else 
  • include/fflas-ffpack/fflas_ftrmm_src.inl

    r38 r40  
    155155        if (!M || !N ) return; 
    156156         
    157         static typename Field::Element one, Mone; 
     157        static typename Field::Element one; 
    158158        F.init(one, 1.0); 
    159         F.neg(Mone, one); 
    160159         
    161160        size_t nsplit = DotProdBound (F, 0, one, 
     
    206205                 typename Field::Element * B, const size_t ldb) { 
    207206         
    208         static typename Field::Element Mone; 
    209207        static typename Field::Element one; 
    210208        F.init(one, 1.0); 
    211         F.neg(Mone,one); 
    212209        if (__FFLAS__Na == 1) 
    213210#ifdef __FFLAS__NONUNIT 
  • include/fflas-ffpack/fflas_ftrsm_src.inl

    r38 r40  
    221221#ifdef __FFLAS__RIGHT 
    222222                fgemm (D, FflasNoTrans, Mjoin (Fflas, __FFLAS__TRANS), __FFLAS__Mb2, __FFLAS__Nb2, nsplit, 
    223                        Mone, __FFLAS__B1, ldb, __FFLAS__A2, lda, one, __FFLAS__B2, ldb); 
     223                       -1.0, __FFLAS__B1, ldb, __FFLAS__A2, lda, one, __FFLAS__B2, ldb); 
    224224#else 
    225225                fgemm (D, Mjoin (Fflas, __FFLAS__TRANS), FflasNoTrans, __FFLAS__Mb2, __FFLAS__Nb2, nsplit, 
    226                        Mone, __FFLAS__A2, lda, __FFLAS__B1, ldb, one, __FFLAS__B2, ldb); 
     226                       -1.0, __FFLAS__A2, lda, __FFLAS__B1, ldb, one, __FFLAS__B2, ldb); 
    227227#endif 
    228228 
  • include/fflas-ffpack/ffpack.h

    r34 r40  
    1212#define __FFPACK_H 
    1313 
    14 #ifdef __LINBOX_CONFIG_H 
     14#ifdef _LINBOX_CONFIG_H 
    1515#include "linbox/fflas/fflas.h" 
    1616#else 
     
    2121#include <vector> 
    2222 
    23 #ifdef __LINBOX_CONFIG_H 
     23#ifdef _LINBOX_CONFIG_H 
    2424namespace LinBox{ 
    2525#endif 
     
    10791079#include "ffpack_krylovelim.inl" 
    10801080#include "ffpack_frobenius.inl" 
    1081 #ifdef __LINBOX_CONFIG_H 
     1081#ifdef _LINBOX_CONFIG_H 
    10821082} 
    10831083#endif