solver_copy.c

Go to the documentation of this file.

/**
 *
 * @file solver_copy.c
 *
 * PaStiX solver matrix copy and reallocation functions.
 *
 * @copyright 2004-2023 Bordeaux INP, CNRS (LaBRI UMR 5800), Inria,
 *                      Univ. Bordeaux. All rights reserved.
 *
 * @version 6.3.2
 * @author Pascal Henon
 * @author Mathieu Faverge
 * @author Tony Delarue
 * @date 2023-07-21
 *
 **/
#include "common.h"
#include "queue.h"
#include "blend/solver.h"
 
/**
 *******************************************************************************
 *
 * @ingroup blend_dev_solver_null
 *
 * @brief Copy the solver matrix data structure from solvin to solvout.
 *
 * Every data is copied, event the coefficient if they are allocated and
 * initialized.
 * It is also used to reallocate the data in a contiguous way after the
 * initialization that allocates all internal arrays in multiple step which
 * might results in fragmentation.
 * @warning This function is not able to copy a solver matrix with low rank
 * blocks yet.
 *
 *******************************************************************************
 *
 * @param[in] solvin
 *          The solver matrix structure to duplicate.
 *
 * @param[out] solvout
 *          The allocated pointer to the solver matrix structure that will
 *          contain the copy.
 *
 * @param[in] flttype
 *          The floating point arithmetic used in the input solver matrix to
 *          know the size of the memory space to duplicate for the coefficients.
 *
 *******************************************************************************/
static inline void
solver_copy( const SolverMatrix *solvin,
             SolverMatrix       *solvout,
             pastix_coeftype_t   flttype )
{
    SolverCblk *solvcblk;
    SolverBlok *solvblok;
    pastix_int_t i;
 
    /* Copy tasktab */
    MALLOC_INTERN(solvout->tasktab, solvout->tasknbr, Task);
    memcpy(solvout->tasktab, solvin->tasktab, solvout->tasknbr*sizeof(Task));
 
    /* Copy cblktab and bloktab */
    MALLOC_INTERN(solvout->cblktab, solvout->cblknbr+1, SolverCblk);
    memcpy(solvout->cblktab, solvin->cblktab,
           (solvout->cblknbr+1)*sizeof(SolverCblk));
 
    MALLOC_INTERN(solvout->bloktab, solvout->bloknbr+1, SolverBlok);
    memcpy(solvout->bloktab, solvin->bloktab,
           (solvout->bloknbr+1)*sizeof(SolverBlok));
 
    MALLOC_INTERN(solvout->browtab, solvout->brownbr, pastix_int_t);
    memcpy(solvout->browtab, solvin->browtab,
           solvout->brownbr*sizeof(pastix_int_t));
 
    if ( solvin->gcbl2loc ) {
        MALLOC_INTERN(solvout->gcbl2loc, solvout->gcblknbr, pastix_int_t);
        memcpy(solvout->gcbl2loc, solvin->gcbl2loc,
               solvout->gcblknbr*sizeof(pastix_int_t));
    }
    else {
        solvout->gcbl2loc = NULL;
    }
 
    solvblok = solvout->bloktab;
    for (solvcblk = solvout->cblktab; solvcblk  < solvout->cblktab + solvout->cblknbr; solvcblk++) {
        pastix_int_t bloknbr = (solvcblk+1)->fblokptr - solvcblk->fblokptr;
        solvcblk->fblokptr = solvblok;
        solvblok += bloknbr;
 
        if ( flttype == PastixPattern ) {
            solvcblk->lcoeftab = NULL;
            solvcblk->ucoeftab = NULL;
        }
        else {
            if ( solvcblk->cblktype & CBLK_COMPRESSED ) {
                /* Not handled for now */
            }
            else {
                void *lcoeftab = solvcblk->lcoeftab;
                void *ucoeftab = solvcblk->ucoeftab;
                size_t size = cblk_colnbr( solvcblk ) * solvcblk->stride
                            * pastix_size_of( flttype );
 
                if ( ucoeftab ) {
                    MALLOC_INTERN( solvcblk->lcoeftab, 2 * size, char );
                    solvcblk->ucoeftab = (char*)lcoeftab + size;
                    memcpy(solvcblk->lcoeftab, lcoeftab, size );
                    memcpy(solvcblk->ucoeftab, ucoeftab, size );
                }
                else {
                    MALLOC_INTERN( solvcblk->lcoeftab, size, char );
                    memcpy(solvcblk->lcoeftab, lcoeftab, size );
                    solvcblk->ucoeftab = NULL;
                }
            }
        }
    }
    solvcblk->fblokptr = solvblok;
 
    /* Copy ttsktab & ttsknbr */
    if (solvout->bublnbr>0)
    {
        MALLOC_INTERN(solvout->ttsknbr, solvout->bublnbr, pastix_int_t);
        memcpy(solvout->ttsknbr, solvin->ttsknbr, solvout->bublnbr*sizeof(pastix_int_t));
        MALLOC_INTERN(solvout->ttsktab, solvout->bublnbr, pastix_int_t*);
 
        for (i=0;i<solvout->bublnbr;i++)
        {
            solvout->ttsktab[i] = NULL;
            MALLOC_INTERN(solvout->ttsktab[i], solvout->ttsknbr[i], pastix_int_t);
            memcpy(solvout->ttsktab[i], solvin->ttsktab[i],
                   solvout->ttsknbr[i]*sizeof(pastix_int_t));
        }
    }
    else
    {
        solvout->ttsknbr = NULL;
        solvout->ttsktab = NULL;
    }
}
 
/**
 *******************************************************************************
 *
 * @ingroup blend_dev_solver
 *
 * @brief Generate a copy of a solver matrix structure.
 *
 * Every data is copied, event the coefficient if they are allocated and
 * initialized.
 * @warning This function is not able to copy a solver matrix with low rank
 * blocks yet.
 *
 *******************************************************************************
 *
 * @param[in] solvin
 *          The solver matrix structure to duplicate.
 *
 * @param[in] flttype
 *          The floating point arithmetic used in the input solver matrix to
 *          know the size of the memory space to duplicate for the coefficients.
 *
 *******************************************************************************
 *
 * @return The pointer to the solver matrix internally allocated and that is a
 *         copy of the input solver. This pointer is NULL if the copy failed.
 *
 *******************************************************************************/
SolverMatrix *
solverCopy( const SolverMatrix *solvin,
            pastix_coeftype_t   flttype )
{
    SolverMatrix *solvout;
 
    MALLOC_INTERN(solvout, 1, SolverMatrix);
    memcpy(solvout, solvin, sizeof(SolverMatrix));
 
    solver_copy( solvin, solvout, flttype );
 
    return solvout;
}
 
/**
 *******************************************************************************
 *
 * @ingroup blend_dev_solver
 *
 * @brief Realloc in a contiguous way a given solver structure.
 *
 * All internal data of the solver structure are reallocated in a contiguous
 * manner to avoid the possible fragmentation from the initialization at
 * runtime.
 * @warning This function is not able to copy a solver matrix with low rank
 * blocks yet.
 *
 *******************************************************************************
 *
 * @param[inout] solvmtx
 *          On entry, the solver matrix to reallocate.
 *          On exit, the solver matrix with all internal data reallocated.
 *
 *******************************************************************************/
void
solverRealloc( SolverMatrix *solvmtx )
{
    SolverMatrix *tmp;
 
    MALLOC_INTERN(tmp, 1, SolverMatrix);
    /** copy general info **/
    memcpy(tmp, solvmtx, sizeof(SolverMatrix));
 
    solver_copy( tmp, solvmtx, PastixPattern );
 
    /** Free the former solver matrix **/
    solverExit(tmp);
    memFree_null(tmp);
}