Concatonation of Matrices in octave API HOW?

[John Weatherwax]

Hello,
 I have several matrices, each one represents the solutions to a system of 
PDE's at several timesteps.  I currently get each one by a call to a fortran 
subroutine that integrates the PDE from time t_{n} to time t_{n+1}.  I would 
like to be able to "concatonate" these matrices into a NDArray (or similar) 
structure to then pass out of the C++ gateway routine I have written that 
interfaces with the fortran code but I cannot seem to be able to figure out 
how to do this.  What I need is something like:

NDArray uAll;
for( int i=1; i<=nTimesteps; i++){
 Matrix uOneStep( numberOfPDEs, numberOfGridPts );
 double* uOneStepP = uOneStep.fortran_vec();
 F77_FUNC(pdecode,PDECODE)(uOneStepP,... other args...)
 uAll.insert( uOneStep )
}
return uAll;

Thus uAll would have its "page" dimension each timeslice.  This may not be 
the best way to do this.  Any "other" method to get this multidemisional 
information out would be helpful too.

TIA,

Wax

_________________________________________________________________
Is your PC infected? Get a FREE online computer virus scan from McAfee® 
Security. http://clinic.mcafee.com/clinic/ibuy/campaign.asp?cid=3963



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------