kv

-

kernel handle returned by StartMaple

rt

-

Maple rtable object

i

-

ith indexing function in an rtable

Public Sub MyScalarMultiply(ByVal kv As Long)

    Dim rts As RTableSettings

    Dim rt As Long

    Dim scalar As Double

    Dim i, j, m, n As Long

    Dim index(0 To 1) As Long

    ' original rtable

    rt = EvalMapleStatement(kv, _

       "Matrix(4,4,(i,j)->i,datatype=float[8],shape=diagonal);")

    ' scalar to element-wise multiply

    scalar = 2

    RTableGetSettings kv, rts, rt

    m = RTableUpperBound(kv, rt, 1)

    n = RTableUpperBound(kv, rt, 2)

    If rts.data_type = RTABLE_FLOAT64 _

       And MapleNumArgs(kv, rts.index_functions) = 1 _

       And RTableIndFn(kv, rt, 1) = RTABLE_INDEX_DIAGONAL _

    Then

        ' handle diagonal shape more efficiently

        If m < n Then

            n = m

        End If

        For i = 1 To n

            index(0) = i

            index(1) = i

            RTableAssignFloat64 kv, rt, index(0), _

                RTableSelectFloat64(kv, rt, index(0)) * scalar

        Next i

    ElseIf rts.data_type = RTABLE_FLOAT64 Then

        For i = 1 To m

            index(0) = i

            For j = 1 To n

                index(1) = j

                RTableAssignFloat64 kv, rt, index(0), _

                    RTableSelectFloat64(kv, rt, index(0)) * scalar

            Next j

        Next i

    End If

    MapleALGEB_Printf1 kv, "rtable = %a", rt

End Sub