Hi

The following function declaration doesn't work. It takes object flowfield and returns object matrix_square.

Code:
matrix_square eigenvalue(flow_field, flow_field);


matrix_square eigenvalue_plus(flow_field &ob, flow_field &ob1)
{

  int i;
  matrix_square temp;

    for (i = 0; i < Grids; i++)
      temp.ptr[i] = ob.ptr[i];							//Element 11
    for (i = 4*Grids; i < 5*Grids; i++)
      temp.ptr[i] = ob.ptr[i - 4*Grids] + ob1.ptr[i - 4*Grids];		//Element 22
    for (i = 8*Grids; i < 9*Grids; i++)
      temp.ptr[i] = 0;								//Element 33

return temp;
}