Assigning matrices is terribly inefficient (unneeded expand)

[asodja]

The following model calles function EulerAngles with same arguments for every element of the matrix:

model TransformMatrix
  function EulerAngles 
    "Computes transformation matrix for given Euler angles"
    input Real phi;
    input Real theta;
    input Real psi;
    output Real A[4, 4];
  algorithm 
    A := diagonal(vector(ones(1, 4)));
    A[1, 1] := cos(psi)*cos(phi) - cos(theta)*sin(phi)*sin(psi);
    A[1, 2] := cos(psi)*sin(phi) + cos(theta)*cos(phi)*sin(psi);
    A[1, 3] := sin(psi)*sin(theta);
    A[2, 1] := -sin(psi)*cos(phi) - cos(theta)*sin(phi)*cos(psi);
    A[2, 2] := -sin(psi)*sin(phi) + cos(theta)*cos(phi)*cos(psi);
    A[2, 3] := cos(psi)*sin(theta);
    A[3, 1] := sin(theta)*sin(phi);
    A[3, 2] := -sin(theta)*cos(phi);
    A[3, 3] := cos(theta);
  end EulerAngles;
  
  parameter Real phi = 0;
  parameter Real theta = 0;
  parameter Real psi = 0;
  parameter Real T[4, 4]=EulerAngles(phi,theta,psi);  
end TransformMatrix;