Changeset 27

Show
Ignore:
Timestamp:
07/09/07 10:20:31 (1 year ago)
Author:
pernet
Message:

update LUdivine calling parameters

Location:
include/fflas-ffpack
Files:
8 modified

Legend:

Unmodified
Added
Removed

  • include/fflas-ffpack/fflas.h

    r22 r27  
    2121#endif 
    2222 
    23 #ifdef _LINBOX_CONFIG_H  
    24 #include "linbox/config-blas.h" 
     23#ifdef _LINBOX_LINBOX_CONFIG_H  
     24#include "linbox/linbox/config-blas.h" 
    2525#include "linbox/field/unparametric.h" 
    2626 
    … …  
    708708#include "fflas_fcopy.inl" 
    709709 
    710 #ifdef _LINBOX_CONFIG_H 
     710#ifdef _LINBOX_LINBOX_CONFIG_H 
    711711} 
    712712#endif 

  • include/fflas-ffpack/fflas_bounds.inl

    r22 r27  
    99 */ 
    1010 
    11 #ifdef _LINBOX_CONFIG_H 
     11#ifdef _LINBOX_LINBOX_CONFIG_H 
    1212#define FFLAS_INT_TYPE Integer 
    1313#else 

  • include/fflas-ffpack/ffpack.h

    r26 r27  
    1212#define __FFPACK_H 
    1313 
    14 #ifdef _LINBOX_CONFIG_H 
     14#ifdef _LINBOX_LINBOX_CONFIG_H 
    1515#include "linbox/fflas/fflas.h" 
    1616#else 
    … …  
    2121#include <vector> 
    2222 
    23 #ifdef _LINBOX_CONFIG_H 
     23#ifdef _LINBOX_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_LINBOX_CONFIG_H 
    10821082} 
    10831083#endif 

  • include/fflas-ffpack/ffpack_charpoly_kgfast.inl

    r1 r27  
    5555                        size_t * Q = new size_t[mc]; 
    5656 
    57                         if ( (r = LUdivine( F, FflasNonUnit, mc, mc,  
     57                        if ( (r = LUdivine( F, FflasNonUnit, FflasNoTrans, mc, mc,  
    5858                                            LUP, mc, P, Q, FfpackLQUP)) < mc ){ 
    5959                                * kg_mc = mc; 

  • include/fflas-ffpack/ffpack_charpoly_kgfastgeneralized.inl

    r1 r27  
    142142                                Q[i]=0; 
    143143 
    144                         size_t r = LUdivine (F, FflasNonUnit, lambda + me, ncols, LUP, ncols,  
     144                        size_t r = LUdivine (F, FflasNonUnit, FflasNoTrans, lambda + me, ncols, LUP, ncols,  
    145145                                           P, Q, FfpackLQUP); 
    146146#if DEBUG 

  • include/fflas-ffpack/ffpack_charpoly_kglu.inl

    r1 r27  
    128128        for ( i=0;i<2*N;++i)  
    129129                Q[i]=0; 
    130         LUdivine( F, FflasNonUnit, 2*N, N, X, N, P, Q, FfpackLQUP); 
     130        LUdivine( F, FflasNonUnit, FflasNoTrans, 2*N, N, X, N, P, Q, FfpackLQUP); 
    131131         
    132132        k = newD( F,d, KeepOn, l, N, X, Q, m); 
    … …  
    219219                for ( i=0;i<2*N;++i)  
    220220                        Q[i]=0; 
    221                 LUdivine( F, FflasNonUnit, nrowX, N, X, N, P, Q, FfpackLQUP); 
     221                LUdivine( F, FflasNonUnit, FflasNoTrans, nrowX, N, X, N, P, Q, FfpackLQUP); 
    222222                 
    223223                // Recompute the degrees of the list factors 

  • include/fflas-ffpack/ffpack_frobenius.inl

    r20 r27  
    6868        for (size_t i=0; i<N; ++i) 
    6969                Pk[i] = 0; 
    70         size_t R = LUdivine(F, FflasNonUnit, noc*c, N, K, ldk, Pk, Qk, FfpackLQUP); 
     70        size_t R = LUdivine(F, FflasNonUnit, FflasNoTrans, noc*c, N, K, ldk, Pk, Qk, FfpackLQUP); 
    7171 
    7272        size_t row_idx = 0; 
    … …  
    310310                size_t *P=new size_t[Mk]; 
    311311                size_t *Q=new size_t[Mk]; 
    312                 if (LUdivine (F, FflasNonUnit, Mk, Mk , K2 + (Ncurr-Mk)*ldk, ldk, P, Q, FfpackLQUP) < Mk){ 
     312                if (LUdivine (F, FflasNonUnit, FflasNoTrans, Mk, Mk , K2 + (Ncurr-Mk)*ldk, ldk, P, Q, FfpackLQUP) < Mk){ 
    313313                        // should never happen (not a LAS VEGAS check) 
    314314                        std::cerr<<"FAIL R2 < MK"<<std::endl; 

  • include/fflas-ffpack/ffpack_ludivine.inl

    r25 r27  
    730730//      tim.clear(); 
    731731//      tim.start(); 
    732         q1 = LUdivine( F, FflasNonUnit, mloc, no2, NW, ld1, P1, Q1, FfpackLQUP, cutoff); 
     732        q1 = LUdivine( F, FflasNonUnit, FflasNoTrans, mloc, no2, NW, ld1, P1, Q1, FfpackLQUP, cutoff); 
    733733         
    734734//      tim.stop(); 
    … …  
    804804        mloc = M-mo2; 
    805805        nloc = N-no2; 
    806         q2 = LUdivine( F, FflasNonUnit, mloc, nloc, SE, ld4, P2, Q2, FfpackLQUP, cutoff); 
     806        q2 = LUdivine( F, FflasNonUnit, FflasNoTrans, mloc, nloc, SE, ld4, P2, Q2, FfpackLQUP, cutoff); 
    807807#if DEBUG 
    808808        std::cerr<<"  E1 = L2.Q2.U2.P2"<<std::endl; 
    … …  
    866866        mloc = M-mo2-q2; 
    867867        nloc = no2-q1; 
    868         q3 = LUdivine( F, FflasNonUnit, mloc, nloc, SW+q2*ld3+q1, ld3, P1+q1, Q2+q2, FfpackLQUP, cutoff); 
     868        q3 = LUdivine( F, FflasNonUnit, FflasNoTrans, mloc, nloc, SW+q2*ld3+q1, ld3, P1+q1, Q2+q2, FfpackLQUP, cutoff); 
    869869         
    870870        // Updating P1,Q2 
    … …  
    878878        nloc = N-no2-q2; 
    879879 
    880         q3b = LUdivine( F, FflasNonUnit, mloc, nloc, NE+q1*ld2+q2, ld2, P2+q2, Q1+q1, FfpackLQUP, cutoff); 
     880        q3b = LUdivine( F, FflasNonUnit, FflasNoTrans, mloc, nloc, NE+q1*ld2+q2, ld2, P2+q2, Q1+q1, FfpackLQUP, cutoff); 
    881881         
    882882        // Updating P2, Q1 
    … …  
    955955//              for (size_t j=0;j<mo2-q1;++j) 
    956956//                      rP4[j]=0; 
    957                 q4 = LUdivine( F, FflasNonUnit, mloc, nloc, NW+(q1+q3b)*ld1+q1+q3, ld1, P1+q1+q3, Q1+q1+q3b, FfpackLQUP, cutoff); 
     957                q4 = LUdivine( F, FflasNonUnit, FflasNoTrans, mloc, nloc, NW+(q1+q3b)*ld1+q1+q3, ld1, P1+q1+q3, Q1+q1+q3b, FfpackLQUP, cutoff); 
    958958 
    959959                // Updating P