Документация

using System;

public interface IMultiply
{
    // Scalar multiplication
    System.Double multiply(System.Double x, System.Double y);

    // Multiply vector by a scalar, return a vector
    System.Double[] multiply(System.Double[] x, System.Double y);

    // Matrix multiplication
    System.Double[,] multiply(System.Double[,] x, System.Double[,] y);
}

compiler.build.dotNETAssembly('multiply.m', ...
    'AssemblyName','Multiply', ...
    'ClassName','Arithmetic');

ntswrap('-c','Multiply.Arithmetic', ...
        '-a','IMultiply.dll', ...
        '-i','IMultiply');

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;

namespace MultiplyApp
{
    class Program
    {
        delegate String VectorToString<T>(IEnumerable<T> v);
        delegate IEnumerable<String> MatrixToString<T>(T[,] m);

        static void Main(string[] args)
        {
            Console.WriteLine("\nStarting application...\n");
            
            //Create an instance of the Type-safe API
            IMultiply m = new ArithmeticIMultiply();
            
            // Scalar multiplication
            double x = 17, y = 3.14159;
            double z = m.multiply(x, y);
            System.Console.WriteLine("{0} * {1} = {2}\n", x, y, z);

            // Vector times scalar
            double[] v = new double[] { 2.5, 81, 64 };
            double s = 11;
            double[] d = m.multiply(v, s);

            VectorToString<double> vec2str = (vec =>
              vec.Select(n => n.ToString()).Aggregate((str, next) => str + " " + next));

            System.Console.WriteLine("[ {0} ] * {1} = [ {2} ]\n",
                vec2str(v), s, vec2str(d));

            // Matrix multiplication
            double[,] magic = new double[,]{  // 3x3 magic square
                     { 8, 1, 6 },
                     { 3, 5, 7 },
                     { 4, 9, 2 } };
            double[,] squareSquared = m.multiply(magic, magic);

            MatrixToString<double> mat2str = mat =>
                mat.EnumerateRows<double>().Select(r => vec2str(r));

            PrintParallel(mat2str(magic), " * ".Select(c => c.ToString()),
                          mat2str(magic), " = ".Select(c => c.ToString()),
                          mat2str(squareSquared));

            Console.WriteLine("\nClosing application...");
        }

        public static void PrintParallel<T>(params IEnumerable<T>[] sources)
        {
            int max = sources.Select(s => s.Count()).Max();
            for (int i = 0; i < max; i++)
            {
                foreach (var src in sources)
                    System.Console.Write("{0} ", src.ElementAt(i));
                System.Console.WriteLine();
            }
        }
    }
    public static class ArrayExtensions
    {
        public static IEnumerable<IEnumerable<T>> EnumerateRows<T>(this Array a)
        {
            return Enumerable.Range(0, a.GetLength(1)).Select(row =>
                a.ToIEnumerable<T>().Skip(row * a.GetLength(0)).Take(a.GetLength(0)));
        }

        public static IEnumerable<T> ToIEnumerable<T>(this Array a)
        {
            foreach (var item in a)
                yield return (T)item;
        }

    }
}