zhagcm.c

Go to the documentation of this file.

#include "morse.h"
#include "control/common.h"
#include "control/hicma_common.h"

int HICMA_zhagcm_Tile(
        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_context_t *morse;
    MORSE_sequence_t *sequence = NULL;
    MORSE_request_t request = MORSE_REQUEST_INITIALIZER;
    int status;

    morse = morse_context_self();
    if (morse == NULL) {
        morse_fatal_error("MORSE_zhagcm_Tile", "morse not initialized");
        return MORSE_ERR_NOT_INITIALIZED;
    }
    morse_sequence_create(morse, &sequence);
    HICMA_zhagcm_Tile_Async(
            uplo,
            AUV,
            Ark,
            numrows_matrix,
            numcols_matrix,
            numrows_block,
            numcols_block,
            maxrank, tol,
            sequence, &request );
    morse_sequence_wait(morse, sequence);
    //RUNTIME_desc_getoncpu(AD);
    status = sequence->status;
    morse_sequence_destroy(morse, sequence);
    return status;
}

int HICMA_zhagcm_Tile_Async(
        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 = morse_context_self();
    if (morse == NULL) {
        morse_fatal_error("MORSE_zhagcm_Tile", "morse not initialized");
        return MORSE_ERR_NOT_INITIALIZED;
    }
    if (sequence == NULL) {
        morse_fatal_error("MORSE_zhagcm_Tile", "NULL sequence");
        return MORSE_ERR_UNALLOCATED;
    }
    if (request == NULL) {
        morse_fatal_error("MORSE_zhagcm_Tile", "NULL request");
        return MORSE_ERR_UNALLOCATED;
    }
    /* Check sequence status */
    if (sequence->status == MORSE_SUCCESS)
        request->status = MORSE_SUCCESS;
    else
        return morse_request_fail(sequence, request, MORSE_ERR_SEQUENCE_FLUSHED);

    /* Check descriptors for correctness */
    if (morse_desc_check(A) != MORSE_SUCCESS) {
        morse_error("MORSE_zhagcm_Tile", "invalid descriptor");
        return morse_request_fail(sequence, request, MORSE_ERR_ILLEGAL_VALUE);
    }
    /* Check input arguments */
    /*if (A->nb != A->mb) {
        morse_error("MORSE_zhagcm_Tile", "only square tiles supported");
        return morse_request_fail(sequence, request, MORSE_ERR_ILLEGAL_VALUE);
    }*/

    /* Quick return */
    if (chameleon_min( A->m, A->n ) == 0)
        return MORSE_SUCCESS;

    hicma_pzhagcm(
            uplo,
            AUV,
            Ark,
            numrows_matrix,
            numcols_matrix,
            numrows_block,
            numcols_block,
            maxrank, tol,
            sequence, request);

    return MORSE_SUCCESS;
}