How can I hook into MATLAB to evaluate symbolic derivatives from C code?

Question

How can I use MATLAB's Symbolic Toolkit to symbolically evaluate derivatives from C or C++?

Answer 1

Here's a working example, in C, along with with Makefile for gcc\Linux. If anybody can provide build instructions for Windows or Mac, please do.

First the C code:

#define _GNU_SOURCE
#include <stdio.h>
#include <ctype.h>
#include <stdlib.h>
#include <string.h>
#include <engine.h>

int main (int argc, char *argv[])
{
    Engine    *engine;
    mxArray   *f;
    mxChar    *output;
    char      *input;
    char      *command;
    size_t     input_size;
    int        i;

    engine = engOpen ("matlab -nojvm");
    engEvalString (engine, "syms x");

    input = NULL;
    printf (">> "); getline (&input, &input_size, stdin);
    while (!feof (stdin)) {
     // remove the \n character
     input[strlen (input) - 1] = '\0';
     command = (char *)malloc (strlen (input) + 19);
     sprintf (command, "f=simple(diff(%s,x))", input);

     // execute the `diff` command on the user input expression
     engEvalString (engine, command);

     // the MATLAB engine does not understand the Symbolic Toolbox
     // therefore you have to convert the output expression into
     // some textual form (ccode,fortran,matlabFunction,latex)
     engEvalString (engine, "fstr=ccode(f)");
     f = engGetVariable (engine, "fstr");
     output = mxGetChars (f);

     // the expression is prefixed with some left-hand line (eg, "t0=")
     // so, discard everything to the left of the =, and display just
     // the derivative
     for (i = 0; output[i] != '='; i++);
     i++; // skip the =

     // `mxChar` is a 16-bit character, and `printf ("%ls")` is giving
     // me errors, so go through the string one character at a time, until
     // the semicolon is found (it is C-code), ignoring non-printable
     // characters along the way.
     for (; output[i] != ';'; i++) {
      if (isgraph (output[i])) {
       putchar (output[i]);
      }
     }
     putchar ('\n');
     printf (">> "); getline (&input, &input_size, stdin);
    }

    printf ("exiting...\n");
    engClose (engine);
    return EXIT_SUCCESS;
}

If you need to use the expression, instead of merely echoing it, you will have to write your own parser. But since it's a simple C expression (with the math.h functions) it is fairly trivial.

And here's my Makefile, nothing fancy .. just remember to set your paths appropriately.

CC=gcc
CFLAGS=-O0 -g -Wall -I$(MATLAB_DIR)/extern/include
LDFLAGS=-L$(MATLAB_DIR)/bin/glnxa64 -leng -lmx -lmat -lut -licudata -licuuc -licui18n -licuio -lhdf5
MATLAB_DIR=/opt/matlab/2008a

%.o: %.c
    $(CC) $(CFLAGS) -c $< -o $@

diff_test: diff_test.o
    $(CC) $< -o $@ $(LDFLAGS)

And finally a quick demo:

# ./diff_test 
>> sin(x)
cos(x)
>> 2*x^2*tan(x)
2.0*x*(2.0*tan(x)+x+x*pow(tan(x),2.0))
>> log(x^2)
2.0/x
>> ^D
exiting...