zhagdm.c

Go to the documentation of this file.

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

int HICMA_zhagdm_Tile(
        MORSE_enum   uplo,
        MORSE_desc_t *Dense
        )
{
    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_zhagdm_Tile", "morse not initialized");
        return MORSE_ERR_NOT_INITIALIZED;
    }
    morse_sequence_create(morse, &sequence);
    HICMA_zhagdm_Tile_Async(
            uplo,
            Dense,
            sequence, &request );
    MORSE_Desc_Flush( Dense, sequence );
    morse_sequence_wait(morse, sequence);
    /*RUNTIME_desc_getoncpu(Dense);*/
    status = sequence->status;
    morse_sequence_destroy(morse, sequence);
    return status;
}

int HICMA_zhagdm_Tile_Async(
        MORSE_enum       uplo,
        MORSE_desc_t *Dense,
        MORSE_sequence_t *sequence,
        MORSE_request_t  *request)
{
    MORSE_context_t *morse;

    morse = morse_context_self();
    if (morse == NULL) {
        morse_fatal_error("MORSE_zhagdm_Tile", "morse not initialized");
        return MORSE_ERR_NOT_INITIALIZED;
    }
    if (sequence == NULL) {
        morse_fatal_error("MORSE_zhagdm_Tile", "NULL sequence");
        return MORSE_ERR_UNALLOCATED;
    }
    if (request == NULL) {
        morse_fatal_error("MORSE_zhagdm_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(Dense) != MORSE_SUCCESS) {
        morse_error("MORSE_zhagdm_Tile", "invalid descriptor");
        return morse_request_fail(sequence, request, MORSE_ERR_ILLEGAL_VALUE);
    }
    /* Check input arguments */
    /*if (A->nb != A->mb) {
        morse_error("MORSE_zhagdm_Tile", "only square tiles supported");
        return morse_request_fail(sequence, request, MORSE_ERR_ILLEGAL_VALUE);
    }*/

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

    hicma_pzhagdm(
            uplo,
            Dense,
            sequence, request);

    return MORSE_SUCCESS;
}

int HICMA_zhagdmdiag_Tile(
        MORSE_enum   uplo,
        MORSE_desc_t *Dense
        )
{
    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_zhagdm_Tile", "morse not initialized");
        return MORSE_ERR_NOT_INITIALIZED;
    }
    morse_sequence_create(morse, &sequence);
    HICMA_zhagdmdiag_Tile_Async(
            uplo,
            Dense,
            sequence, &request );
    MORSE_Desc_Flush( Dense, sequence );
    morse_sequence_wait(morse, sequence);
    /*RUNTIME_desc_getoncpu(Dense);*/
    status = sequence->status;
    morse_sequence_destroy(morse, sequence);
    return status;
}

int HICMA_zhagdmdiag_Tile_Async(
        MORSE_enum       uplo,
        MORSE_desc_t *Dense,
        MORSE_sequence_t *sequence,
        MORSE_request_t  *request)
{
    MORSE_context_t *morse;

    morse = morse_context_self();
    if (morse == NULL) {
        morse_fatal_error("MORSE_zhagdm_Tile", "morse not initialized");
        return MORSE_ERR_NOT_INITIALIZED;
    }
    if (sequence == NULL) {
        morse_fatal_error("MORSE_zhagdm_Tile", "NULL sequence");
        return MORSE_ERR_UNALLOCATED;
    }
    if (request == NULL) {
        morse_fatal_error("MORSE_zhagdm_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(Dense) != MORSE_SUCCESS) {
        morse_error("MORSE_zhagdm_Tile", "invalid descriptor");
        return morse_request_fail(sequence, request, MORSE_ERR_ILLEGAL_VALUE);
    }
    /* Check input arguments */
    /*if (A->nb != A->mb) {
        morse_error("MORSE_zhagdm_Tile", "only square tiles supported");
        return morse_request_fail(sequence, request, MORSE_ERR_ILLEGAL_VALUE);
    }*/

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

    hicma_pzhagdmdiag(
            uplo,
            Dense,
            sequence, request);

    return MORSE_SUCCESS;
}