Matrix exponential

Dear folks,

I have a system of ODE of the form y_dot = A(y) * y,
where y_dot is time derivative of vector y and A(y) a tensor (n by n Matrix form).
For numerical approximation, I need to calculate exponential of matrix A, which is exp(del_t * A)
where del_t is time step.
Do you have any idea how to implement this using the PDE module? or is there any other method to do this?

Thank you.
Wonseok Yoon