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;
}