Matrix/Array multiplication

General Discussion
1

Hello,

I am using a Gaussian quadrature method for computing statistical moments.
I am running into a surprising issue, the following bit of code does not work :

real[int] A = [1, 2, 3];
real[int] B = [4, 5, 6];
real[int, int] C = [[1, 2, 3], [4, 5, 6], [7, 8, 9]];
real D = A'*C*B;

While this one does :

real[int] A = [1, 2, 3];
real[int] B = [4, 5, 6];
real[int, int] C = [[1, 2, 3], [4, 5, 6], [7, 8, 9]];
real[int] E = C*B;
real D = A'*E;

What is the issue with the first bit of code, and what has significantly changed in the second code that makes it work ?
Thank you in advance,
Hugo

2

FreeFEM does not know how to compute a series of operations grouped a single statement, as done in your first code snippet, so you need to break them down into multiple statements, as done in your second code snippet.