toeplitz_utils.c

Go to the documentation of this file.

 
#include "toeplitz.h"
#include <time.h>
 
extern int PRINT_RANK;
 
// set of utilitaries routines - fd@apc
 
 
int defineTpltz(Tpltz *Nm1, int64_t nrow, int m_cw, int m_rw,
                Block *tpltzblocks, int nb_blocks_loc, int nb_blocks_tot,
                int64_t idp, int local_V_size, Flag flag_stgy, MPI_Comm comm) {
 
    Nm1->nrow = nrow; // glob //recup du fichier params apres (en variables
                      // globales)
    Nm1->m_cw          = m_cw;          // glob
    Nm1->m_rw          = m_rw;          // glob
    Nm1->tpltzblocks   = tpltzblocks;   // toep
    Nm1->nb_blocks_loc = nb_blocks_loc; // toep
    Nm1->nb_blocks_tot = nb_blocks_tot; // toep
    Nm1->idp           = idp;           // comput
    Nm1->local_V_size  = local_V_size;  // comput
    Nm1->flag_stgy     = flag_stgy;     // param
    Nm1->comm          = comm;          // param
 
 
    return 0;
}
 
 
int defineBlocks_avg(Block *tpltzblocks, double *T, int nb_blocks_loc,
                     int n_block_avg, int lambda_block_avg, int64_t id0) {
 
    int i;
 
 
    for (i = 0; i < nb_blocks_loc; i++) tpltzblocks[i].n = n_block_avg;
 
    for (i = 0; i < nb_blocks_loc; i++)
        tpltzblocks[i].lambda = lambda_block_avg;
 
    tpltzblocks[0].idv = (int64_t) (id0 / n_block_avg) * n_block_avg;
    for (i = 1; i < nb_blocks_loc; i++)
        tpltzblocks[i].idv =
                (int64_t) tpltzblocks[i - 1].idv + tpltzblocks[i - 1].n;
 
    for (i = 0; i < nb_blocks_loc; i++) { tpltzblocks[i].T_block = (T); }
 
 
    return 0;
}
 
 
//=============================================
 
 
int createRandomT(double *T, int Tsize) {
 
    int i;
    srand(time(NULL)); // init seed
 
    // input matrix definition of T
    for (i = 0; i < Tsize; i++) T[i] = rand() / ((double) RAND_MAX);
 
    return 0;
}
 
 
int createTbasic1(double *T, int Tsize) {
 
    int i;
    srand(Tsize);
 
    // input matrix definition of T
    for (i = 0; i < Tsize; i++) T[i] = 1.0 + rand() / ((double) RAND_MAX);
 
    return 0;
}
 
 
int createTbasic2(double *T, int Tsize) {
 
    int i;
    srand(Tsize);
 
    // input matrix definition of T
    for (i = 0; i < Tsize; i++) {
        if (i == 0) {
            T[i] = 10.;
        } else if (i == 1) {
            T[i] = 2.;
        } else if (i == 2) {
            T[i] = 3.;
        } else {
            T[i] = rand() / ((double) RAND_MAX);
        }
    }
 
    return 0;
}
 
 
int createTbasic3(double *T, int Tsize) {
 
    int i;
    srand(Tsize);
 
    // input matrix definition of T
    for (i = 0; i < Tsize; i++) {
        if (i == 0) {
            T[i] = 2.;
        } else if (i == 1) {
            T[i] = -1.;
        } else if (i == 2) {
            T[i] = 0.;
        } else {
            T[i] = 0.; // rand()/((double) RAND_MAX);
        }
    }
 
 
    return 0;
}
 
 
int createTfrominvtt(double *T, int Tsize) {
 
    int i;
 
    // #include "invtt_params.h"
 
    double *invtt;
 
    T = invtt;
    //  createinvtt(invtt);
 
 
    return 0;
}