pzhagcm.c

Go to the documentation of this file.

/*
 *
 * file pzplrnt.c
 *
 * MORSE auxiliary routines
 * MORSE is a software package provided by Univ. of Tennessee,
 * Univ. of California Berkeley and Univ. of Colorado Denver
 *
 * version 2.5.0
 * comment This file has been automatically generated
 *         from Plasma 2.5.0 for MORSE 1.0.0
 * author Mathieu Faverge
 * author Emmanuel Agullo
 * author Cedric Castagnede
 * date 2010-11-15
 *
 **/
#include "morse.h"
#include "control/common.h"
#include "hicma_runtime_z.h"
extern int store_only_diagonal_tiles;

void hicma_pzhagcm(
        MORSE_enum uplo,
        MORSE_desc_t *AUV,
        MORSE_desc_t *Ark,
        int numrows_matrix,
        int numcols_matrix,
        int numrows_block,
        int numcols_block,
        int maxrank, double tol,
        MORSE_sequence_t *sequence, MORSE_request_t *request )
{
    MORSE_desc_t *A = AUV; // FIXME
    MORSE_context_t *morse;
    MORSE_option_t options;

    int m, n;
    int tempmm, tempnn;

    morse = morse_context_self();
    if (sequence->status != MORSE_SUCCESS)
        return;
    RUNTIME_options_init(&options, morse, sequence, request);

    for (m = 0; m < A->mt; m++) {
        tempmm = m == A->mt-1 ? numrows_matrix-m*numrows_block : numrows_block;
        int ldamUV = BLKLDD(AUV, m);

        //for (n = 0; n < A->mt; n++) {
        //    tempnn = n == A->mt-1 ? A->m-n*A->mb : A->mb;
        for (n = 0; n < A->nt; n++) { //I hope this change does not break anything
            tempnn = n == A->nt-1 ? numcols_matrix-n*numcols_block : numcols_block;

            // if(m<n)
            //     continue;
            if(uplo == MorseLower && m <= n)
                continue;
            else
            if(uplo == MorseUpper && m >= n)
                continue;
            //printf("Tile (%d,%d) ldamUV=%d A->m=%d A->n=%d A->mb=%d A->nb=%d tempmm=%d tempnn=%d (%dx%d) (%dx%d) (%dx%d)\n", m, n, ldamUV, A->m, A->n, A->mb, A->nb, tempmm, tempnn, numrows_matrix, numcols_matrix, numrows_block, numcols_block, A->mt, A->nt);

            HICMA_TASK_zhagcm(
                    &options,
                    tempmm, tempnn,
                    AUV,
                    Ark,
                    m, n, 
                    ldamUV, ldamUV,
                    maxrank, tol
                    );
        }
    }
    RUNTIME_options_finalize(&options, morse);
    //MORSE_TASK_dataflush_all(); removed in newer chameleon
}