kv

-

kernel handle of type MKernelVector

rts

-

pointer to an RTableSettings structure

pdata

-

pointer to array data (optional)

bounds

-

array of lower and upper bounds

    #include "maplec.h"

    ALGEB M_DECL MyIdentity( MKernelVector kv, ALGEB *args )

    {

    M_INT argc, n, i;

    RTableSettings rts;

    M_INT bounds[4];

    ALGEB rt;

    INTEGER32 *data;

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

    if( argc != 1 ) {

        MapleRaiseError(kv,"one argument expected");

            return( NULL );

        }

    n = MapleToM_INT(kv,args[1]);

    RTableGetDefaults(kv,&rts);

    rts.num_dimensions = 2;

    rts.subtype = RTABLE_MATRIX;

    rts.data_type = RTABLE_INTEGER32;

    bounds[0] = 1;

    bounds[1] = n;

    bounds[2] = 1;

    bounds[3] = n;

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

    data = (INTEGER32*)RTableDataBlock(kv,rt);

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

        data[MATRIX_OFFSET_FORTRAN_RECT(i,i,n,n)] = 1;

    }

    return( rt );

    }