AMDG Allen Chen wrote:

Six years ago, I design my own matrix library like matlab. I overload operators, it make it conveninent to do matrix manupulations. For example

matrix A, B, C, D; A = B*C +D;

It is easy to use and read.

Why does not ublas overload such operators as +, -, *, /, ^, [], (),...

It does. However, operator* only works for multiplying by a scalar. For matrix-matrix multiplication use prod. In Christ, Steven Watanabe