kv

-

kernel handle of type MKernelVector

rt

-

type ALGEB rtable object

dim

-

integer specifying dimension

    #include "maplec.h"

    ALGEB M_DECL MyReindexRTable( MKernelVector kv, ALGEB *args )

    {

    M_INT argc, i, num_dims, n, bounds[RTABLE_MAX_BOUNDS];

    ALGEB rt1, rt2, *data1, *data2;

    RTableSettings rts;

    argc = MapleNumArgs(kv,(ALGEB)args);

    if( argc != 1 ) {

        MapleRaiseError(kv,"one argument expected");

            return( NULL );

        }

    if( !IsMapleRTable(kv,args[1]) ) {

        MapleRaiseError(kv,"rtable expected for parameter 1");

            return( NULL );

        }

    rt1 = args[1];

    RTableGetSettings(kv,&rts,rt1);

    if( rts.data_type != RTABLE_DAG ) {

        MapleRaiseError(kv,"DAG rtable expected for parameter 1");

            return( NULL );

        }

    if( rts.storage == RTABLE_SPARSE ) {

        MapleRaiseError(kv,"dense rtable expected for parameter 1");

            return( NULL );

        }

    num_dims = RTableNumDimensions(kv,rt1);

    for( i=1; i<=num_dims; ++i ) {

        bounds[2*i-2] = 1;

        bounds[2*i-1] = RTableUpperBound(kv,rt1,i)

         - RTableLowerBound(kv,rt1,i) + 1;

    }

    rt2 = RTableCreate(kv,&rts,NULL,bounds);

    n = RTableNumElements(kv,rt1);

    data1 = (ALGEB*)RTableDataBlock(kv,rt1);

    data2 = (ALGEB*)RTableDataBlock(kv,rt2);

    for( i=0; i<n; ++i ) {

        data2[i] = data1[i];

    }

    return( rt2 );

    }