extracblk.c

Go to the documentation of this file.

/**
 *
 * @file extracblk.c
 *
 * PaStiX analyse headers for extra symbolic structure functions.
 *
 * @copyright 1998-2025 Bordeaux INP, CNRS (LaBRI UMR 5800), Inria,
 *                      Univ. Bordeaux. All rights reserved.
 *
 * @version 6.4.0
 * @author Mathieu Faverge
 * @author Gregoire Pichon
 * @date 2024-07-05
 *
 * @addtogroup blend_dev_split
 * @{
 *
 **/
#include "common.h"
#include "symbol/symbol.h"
#include "cost.h"
#include "elimintree.h"
#include "cand.h"
#include "extracblk.h"
 
/**
 *******************************************************************************
 *
 * @brief Allocate the extracblk structure.
 *
 *******************************************************************************
 *
 * @param[inout] extracblk
 *          Allocate memory space to store extra cblk when they are created. The
 *          initial number is the one given at extraCblkInit().
 *
 *******************************************************************************/
static inline void
extraCblkAlloc( ExtraCblk_t *extracblk )
{
    pastix_int_t i;
 
    MALLOC_INTERN( extracblk->sptcblk, extracblk->cblknbr, pastix_int_t );
    MALLOC_INTERN( extracblk->sptcbnb, extracblk->cblknbr, pastix_int_t );
 
    /*
     * Unsplitted cblk will keep sptcblk to -1 and sptcbnb to 1
     * Splitted cblk will have sptcbnb to the number of cblk replacing the original one,
     * and sptcblk will contain the cblktab index of the first generated cblk
     */
    for(i=0; i<extracblk->cblknbr;i++)
    {
        extracblk->sptcblk[i] = -1;
        extracblk->sptcbnb[i] =  1;
    }
 
    /* We choose an arbitrary size for initial allocation of bloktab and cblktab (5%) */
    extracblk->sizcblk = ( extracblk->cblknbr + 20 ) / 20;
    MALLOC_INTERN( extracblk->cblktab, extracblk->cblknbr, symbol_cblk_t );
 
    return;
}
 
/**
 *******************************************************************************
 *
 * @brief Increment the number of extra cblk that can be stored.
 *
 * If there is not enough space, the structure is reallocated to a larger space.
 *
 *******************************************************************************
 *
 * @param[inout] extracblk
 *          Allocate memory space to store extra cblk when they are created. The
 *          initial number is the one given at extraCblkInit().
 *
 *******************************************************************************
 *
 * @return the current index for the next cblk added to the structure.
 *
 *******************************************************************************/
static inline pastix_int_t
extraCblkInc( ExtraCblk_t *extracblk )
{
    /* First cblk added */
    if ( extracblk->sizcblk == 0 ) {
        extraCblkAlloc( extracblk );
    }
 
    extracblk->curcblk++;
 
    /* Check that we have enough space and make it bigger if required */
    if( extracblk->curcblk >= extracblk->sizcblk )
    {
        /* Add 5% of original cblknbr to the cblktab */
        pastix_int_t extrasize = (extracblk->cblknbr + 20 ) / 20;
        symbol_cblk_t *tmp;
 
        assert( extracblk->curcblk == extracblk->sizcblk);
        tmp = extracblk->cblktab;
        extracblk->sizcblk += extrasize;
 
        MALLOC_INTERN( extracblk->cblktab, extracblk->sizcblk, symbol_cblk_t );
        memcpy(extracblk->cblktab, tmp, sizeof(symbol_cblk_t)*extracblk->curcblk);
 
        memFree_null(tmp);
    }
 
    return extracblk->curcblk;
}
 
/**
 *******************************************************************************
 *
 * @brief Initialize the extracblk structure.
 *
 *******************************************************************************
 *
 * @param[in]    cblknbr
 *          The starting number of cblk that the structure can hold.
 *
 * @param[inout] extracblk
 *          Pointer to the allocated extracblk structure to initialize.
 *
 *******************************************************************************/
void
extraCblkInit( pastix_int_t  cblknbr,
               ExtraCblk_t  *extracblk )
{
    extracblk->cblknbr = cblknbr;
    extracblk->addcblk = 0;
    extracblk->addblok = 0;
    extracblk->addblof = 0;
    extracblk->sptcblk = NULL;
    extracblk->sptcbnb = NULL;
    extracblk->curcblk = -1;
    extracblk->sizcblk = 0;
    extracblk->cblktab = NULL;
    return;
}
 
/**
 *******************************************************************************
 *
 * @brief Free the extracblk structure.
 *
 *******************************************************************************
 *
 * @param[inout] extracblk
 *          Pointer to the allocated extracblk structure to free.
 *
 *******************************************************************************/
void
extraCblkExit( ExtraCblk_t *extracblk )
{
    if ( extracblk->sizcblk > 0 ) {
        memFree_null( extracblk->sptcblk );
        memFree_null( extracblk->sptcbnb );
        memFree_null( extracblk->cblktab );
    }
    extracblk->curcblk = -1;
    extracblk->sizcblk = 0;
    return;
}
 
/**
 *******************************************************************************
 *
 * @brief Add a new additional cblk defined by its first and last columns.
 *
 *******************************************************************************
 *
 * @param[inout] extracblk
 *          Pointer to the extracblk structure to add the cblk.
 *
 * @param[in] fcolnum
 *          Index of the first column in the new cblk.
 *
 * @param[in] lcolnum
 *          Index of the last column in the new cblk.
 *
 * @param[in] selevtx
 *          TODO
 *
 *******************************************************************************/
void
extraCblkAdd( ExtraCblk_t *extracblk,
              pastix_int_t fcolnum,
              pastix_int_t lcolnum,
              int8_t       selevtx )
{
    pastix_int_t curcblk = extraCblkInc( extracblk );
 
    extracblk->cblktab[ curcblk ].fcolnum = fcolnum;
    extracblk->cblktab[ curcblk ].lcolnum = lcolnum;
    extracblk->cblktab[ curcblk ].bloknum = -1;
    extracblk->cblktab[ curcblk ].selevtx = selevtx;
}
 
/**
 *******************************************************************************
 *
 * @brief Merge the existing symbol structure with the additional information
 * from the extracblk structure.
 *
 *******************************************************************************
 *
 * @param[in] extracblk
 *          Pointer to the extracblk structure that contains information about
 *          splited cblk.
 *
 * @param[inout] newsymb
 *          Symbol matrix to update. On exit, the symbol matrix structure is
 *          extended by the splited cblk described in extracblk structure.
 *
 * @param[inout] candtab
 *          On entry, the candtab aray associated to the input symbol matrix.
 *          On exit, the updated candtab array with the extended number of cblk.
 *
 *******************************************************************************/
void
extraCblkMerge( const ExtraCblk_t *extracblk,
                symbol_matrix_t   *newsymb,
                Cand             **candtab )
{
    pastix_int_t  i, j, k, l;
    pastix_int_t  curbloknum, curcblknum;
    pastix_int_t  lastcblksplit;
    pastix_int_t  addblok = 0;
    pastix_int_t *newnum      = NULL;
    pastix_int_t *extranewnum = NULL;
 
    symbol_matrix_t *oldsymb;
 
#if defined(PASTIX_BLEND_PROPMAP_2STEPS)
    Cand *oldcand = *candtab;
    Cand *newcand;
#else
    (void)candtab;
#endif
 
    symbol_cblk_t *curcblk;
    symbol_blok_t *curblok;
 
    /* No splitted cblk: partition remains the same */
    if( extracblk->addcblk == 0 ) {
        return;
    }
 
    /* Backup the old symbol */
    MALLOC_INTERN(oldsymb, 1, symbol_matrix_t);
    memcpy( oldsymb, newsymb, sizeof(symbol_matrix_t) );
 
    /* Allocate new cblktab */
    newsymb->cblknbr = oldsymb->cblknbr + extracblk->addcblk;
    MALLOC_INTERN(newsymb->cblktab, newsymb->cblknbr+1, symbol_cblk_t);
 
    newsymb->browtab = NULL;
 
#if defined(PASTIX_BLEND_PROPMAP_2STEPS)
    /* Allocate new candtab */
    newcand = candInit( newsymb->cblknbr );
 
    /* Copy the root cand */
    newcand[-1] = oldcand[-1];
#endif
 
    /*
     * We use the sptcbnb array to get the new numbering of the former cblk
     * in the new symbolic matrix
     * newnum[i+1] becomes the new number of the first cblk generated from the
     * split of former cblk number i.
     */
    MALLOC_INTERN(newnum,      oldsymb->cblknbr+1,   pastix_int_t);
    MALLOC_INTERN(extranewnum, extracblk->curcblk+1, pastix_int_t);
    newnum[0] = 0;
    memcpy(newnum+1, extracblk->sptcbnb, (oldsymb->cblknbr) * sizeof(pastix_int_t));
#ifndef NDEBUG
    memset(extranewnum, 0xff, (extracblk->curcblk+1) * sizeof(pastix_int_t));
#endif
 
    /* Compute number of blocks that will be generated,
     * and copy main information of cblktab and candtab */
    lastcblksplit = -1;
    for(i=0; i<oldsymb->cblknbr; i++)
    {
        pastix_int_t fbloknum = oldsymb->cblktab[i  ].bloknum;
        pastix_int_t lbloknum = oldsymb->cblktab[i+1].bloknum;
        pastix_int_t sptcbnbw = extracblk->sptcbnb[i];
 
        /*
         * First we compute the number of extra blocks that will be generated
         */
 
        /* Diagonal block */
        addblok += (((sptcbnbw+1) * sptcbnbw) / 2) - 1;
        for(j=fbloknum+1; j<lbloknum; j++)
        {
            pastix_int_t fcblknum = oldsymb->bloktab[j].fcblknm;
            pastix_int_t sptfcbnb = extracblk->sptcbnb[fcblknum];
            pastix_int_t sptcbnbh = 0;
 
            /* If facing cblk is splitted */
            if ( sptfcbnb > 1 )
            {
                symbol_cblk_t  *newfcblk =  &(extracblk->cblktab[ extracblk->sptcblk[fcblknum] ]);
                pastix_int_t frownum  = oldsymb->bloktab[j].frownum;
                pastix_int_t lrownum  = oldsymb->bloktab[j].lrownum;
 
                /* Compute how many times the block is splitted horizontally */
                for(k = 0; k < sptfcbnb; k++, newfcblk++)
                {
                    /* This block doesn't face this new cblk */
                    if ( frownum > newfcblk->lcolnum ) {
                        continue;
                    }
 
                    /* No more facing cblk will be found */
                    if ( lrownum < newfcblk->fcolnum ) {
                        break;
                    }
 
                    assert( frownum <= lrownum );
                    sptcbnbh++;
                    frownum = newfcblk->lcolnum+1;
                }
            }
            else
                sptcbnbh = 1;
 
            /*
             * The number of extra blocks is the number of times the block
             * is psplitted horizontally times the number of time the cblk
             * is splitted vertically minu itself
             */
            addblok += sptcbnbw * sptcbnbh - 1;
        }
 
        /*
         * Second, we create newnum/extranewnum arrays and copy information into
         * cblktab and candtab
         */
        {
            /* This cblk is splitted, we generate new cblktab from extra */
            pastix_int_t newcblknum = newnum[i];
            if (sptcbnbw > 1) {
                pastix_int_t nbcblk2copy = (i - lastcblksplit - 1);
                pastix_int_t sptcblk = extracblk->sptcblk[i];
 
                /* Copy the previous unchanged cblks from oldsymb */
                if ( nbcblk2copy > 0 ) {
                    lastcblksplit++;
                    memcpy( newsymb->cblktab + newnum[ lastcblksplit ],
                            oldsymb->cblktab + lastcblksplit,
                            nbcblk2copy * sizeof(symbol_cblk_t) );
 
#if defined(PASTIX_BLEND_PROPMAP_2STEPS)
                    memcpy( newcand + newnum[ lastcblksplit ],
                            oldcand + lastcblksplit,
                            nbcblk2copy * sizeof(Cand) );
#endif
                }
 
                /* Copy the new cblk from extracblk */
                assert( (sptcblk >= 0) && (sptcblk <= extracblk->curcblk) );
                memcpy( newsymb->cblktab   + newcblknum,
                        extracblk->cblktab + sptcblk,
                        sptcbnbw * sizeof(symbol_cblk_t) );
 
                /* Initialize extranewnum and duplicate the cand for each new cblk */
                for(j=0; j<sptcbnbw; j++, sptcblk++) {
                    extranewnum[sptcblk] = newcblknum+j;
 
                    assert( (extranewnum[sptcblk] >= 0) &&
                            (extranewnum[sptcblk] <  newsymb->cblknbr) );
 
#if defined(PASTIX_BLEND_PROPMAP_2STEPS)
                    memcpy( newcand + extranewnum[ sptcblk ],
                            oldcand + i, sizeof(Cand) );
#endif
                }
 
                lastcblksplit = i;
            }
            /* Update newnum of following cblk (newnum is allocated with one extra space) */
            newnum[i+1] += newcblknum;
        }
    }
 
    /* Copy last unsplitted block */
    if ( lastcblksplit < (oldsymb->cblknbr-1) )
    {
        pastix_int_t nbcblk2copy = oldsymb->cblknbr - lastcblksplit - 1;
        lastcblksplit++;
        memcpy( newsymb->cblktab + newnum[ lastcblksplit ],
                oldsymb->cblktab + lastcblksplit,
                nbcblk2copy * sizeof(symbol_cblk_t) );
 
#if defined(PASTIX_BLEND_PROPMAP_2STEPS)
        memcpy( newcand + newnum[ lastcblksplit ],
                oldcand + lastcblksplit,
                nbcblk2copy * sizeof(Cand) );
#endif
    }
 
#if defined(PASTIX_BLEND_PROPMAP_2STEPS)
    candExit(oldcand);
    *candtab = newcand;
#endif
 
    /* Allocate new bloktab */
    newsymb->bloknbr = oldsymb->bloknbr + addblok;
    MALLOC_INTERN(newsymb->bloktab, newsymb->bloknbr,   symbol_blok_t);
 
    /* Fill in the new symbolic matrix resulting from the splitting of the former one */
    curbloknum = 0;
    curcblknum = 0;
    curcblk = newsymb->cblktab;
    curblok = newsymb->bloktab;
 
#if defined(PASTIX_SYMBOL_DUMP_SYMBMTX)
    symbol_cblk_t *cblk = newsymb->cblktab;
#endif
    for(i=0; i<oldsymb->cblknbr; i++)
    {
        pastix_int_t fbloknum = oldsymb->cblktab[i  ].bloknum;
        pastix_int_t lbloknum = oldsymb->cblktab[i+1].bloknum;
        pastix_int_t sptcbnbw = extracblk->sptcbnb[i];
 
        /* For each new column created by the cblk split */
        for(j=0; j<sptcbnbw; j++, curcblknum++, curcblk++)
        {
            /* Store diagonal bloknum */
            curcblk->bloknum = curbloknum;
 
#if defined(PASTIX_SYMBOL_DUMP_SYMBMTX)
            cblk->split_cblk = sptcbnbw-j-1;
            cblk++;
#endif
 
            /* Create odb due to the splitting of the diag blok */
            for(k=0; k<(sptcbnbw-j); k++, curbloknum++, curblok++)
            {
                curblok->frownum = curcblk[k].fcolnum;
                curblok->lrownum = curcblk[k].lcolnum;
                curblok->lcblknm = curcblknum;
                curblok->fcblknm = curcblknum + k;
            }
            /* Next cblk will have one block less on the diagonal */
 
            /* Create other off diagonal blocks */
            for(k=fbloknum+1; k<lbloknum; k++)
            {
                pastix_int_t frownum  = oldsymb->bloktab[k].frownum;
                pastix_int_t lrownum  = oldsymb->bloktab[k].lrownum;
                pastix_int_t fcblknum = oldsymb->bloktab[k].fcblknm;
                pastix_int_t sptfcblk = extracblk->sptcblk[fcblknum];
                pastix_int_t sptfcbnb = extracblk->sptcbnb[fcblknum];
 
                /* If facing cblk is splitted */
                if ( sptfcbnb > 1 )
                {
                    pastix_int_t   newfcblknum = extranewnum[ sptfcblk ];
                    symbol_cblk_t *newfcblk    = &(extracblk->cblktab[ sptfcblk ]);
 
                    assert( newfcblknum != -1 );
                    assert( newfcblk != NULL );
 
                    /* Create new blocks facing this cblk */
                    for(l=0; l<sptfcbnb; l++, newfcblk++)
                    {
                        /* This block doesn't face this new cblk */
                        if ( frownum > newfcblk->lcolnum ) {
                            continue;
                        }
 
                        /* No more facing cblk will be found */
                        if ( lrownum < newfcblk->fcolnum ) {
                            break;
                        }
 
                        assert( frownum <= lrownum );
                        assert( frownum >= newfcblk->fcolnum );
 
                        curblok->frownum = frownum;
                        curblok->lrownum = pastix_imin( lrownum, newfcblk->lcolnum );
                        curblok->lcblknm = curcblknum;
                        curblok->fcblknm = newfcblknum + l;
                        curblok++; curbloknum++;
 
                        frownum = newfcblk->lcolnum+1;
                    }
                }
                else
                {
                    curblok->frownum = frownum;
                    curblok->lrownum = lrownum;
                    curblok->lcblknm = curcblknum;
                    curblok->fcblknm = newnum[fcblknum];
                    curblok++; curbloknum++;
                }
            }
        }
    }
 
    assert(curcblknum == newsymb->cblknbr);
    assert(curbloknum == newsymb->bloknbr);
    assert((curcblk - newsymb->cblktab) == newsymb->cblknbr);
    assert((curblok - newsymb->bloktab) == newsymb->bloknbr);
 
    /* Free old versions and temporary buffer */
    pastixSymbolExit(oldsymb);
    memFree_null(oldsymb);
    memFree_null(newnum);
    memFree_null(extranewnum);
 
    /* Virtual cblk to avoid side effect in the loops on cblk bloks */
    curcblk[0].fcolnum = curcblk[-1].lcolnum + 1;
    curcblk[0].lcolnum = curcblk[-1].lcolnum + 1;
    curcblk[0].bloknum = curbloknum;
 
    pastixSymbolBuildRowtab( newsymb );
 
    return;
}
 
/**
 *@}
 */