Saving MULTI-Dimensonal Structure Data: from Matalb to Mex Fortran INPUTS-

Gang, I can pass and save a single structure of data in my Mex fortran file. My problem is how do I save multidimensional structure data.... Right now it is only saving the LAST dimension of data...

This is what I have in my Fortran file. The problem with this is only THE LAST array of inputs gets saved in ALL indexes.

 % code
 Do II = 1,10
   Do JJ = 1,5
      tmp = mxGetField(prhs(1),JJ,'Time')
      call mxCopyPtrToReal4(mxGetData(tmp),Data_Block(II,JJ)%Time,1)
   EndDo
 EndDo

So what happens is 'Time' is ALL the same value for ALL indexes of the Data_Block Strucure..

The II or JJ indexes doesn't seem to increment or something.... Any Help ????

Answer 1

[Geoff Hayes]

Rick - your code to get the time data from the structure depends only on JJ

 tmp = mxGetField(prhs(1),JJ,'Time')

Why are you not using II as well?

What does your input data look like? Are all time fields distinct?

[Rick]

Geoff, I guess the main thing is that I really don't know how to use the II indexing as well in the "tmp".

    % code
The input data structure from matlab Looks like this:
Block_data(10).GROUP1(5).Time

    % code
The actual call to this MEX file looks like this:
Call_Wrapper(Block_Data(10).GROUP1)

[Geoff Hayes]

Rick - since you are passing in only one instance of the GROUP1 i.e. that from the tenth element of Block_data, then all you are passing in is a struct array of five elements. To verify this, before calling your Call_Wrapper, type

 size(Block_Data(10).GROUP1)

in the Command Window to see the dimensions of the object that you are going to pass.

[Rick]

Thanx Geoff...... Yeah that is correct. I changed it to: Call_Wrapper(Block_Data(10).GROUP1) So now I am sending all 10....... :-)

So now how do I capture all 10 instances along with the 5 GROUP1 instances.... of "Time" back in he Fortran wrapper. My original statements

 Do II = 1,10
   Do JJ = 1,5
      tmp = mxGetField(prhs(1),JJ,'Time')
      call mxCopyPtrToReal4(mxGetData(tmp),Data_Block(II,JJ)%Time,1)
   EndDo
 EndDo

still do not work..........

[Geoff Hayes]

Rick - I don't see the change, or are you passing

 Call_Wrapper(Block_Data)

instead?

[Rick]

Yes, I am sorry. I meant what you just wrote I changed it to: Call_Wrapper(Block_Data)

[Geoff Hayes]

Rick - I created some dummy test data as

 Block_Data = [];
 for k=1:10
      Block_Data(k).GROUP1 = []; 
      for j=1:5
           Block_Data(k).GROUP1(j).time = double(100*k+j);
      end
 end

Now in the main body of the Fortran code, I did the following

 C     get the size of the input array.
       mrows = mxGetM(prhs(1))
       ncols = mxGetN(prhs(1))

 C     iterate over each column (GROUP1)
       Do n = 1,ncols

 C         get the field GROUP1
           gmx = mxGetField(prhs(1),n,"GROUP1")

 C         get the dimensions of the GROUP1
           gmrows = mxGetM(gmx)
           gncols = mxGetN(gmx)

 C         iterate over each column (time)
           Do gn = 1,gncols

 C             get the time data
               tmx = mxGetPr(mxGetField(gmx,gn,"time"))
               call mxCopyPtrToReal8(tmx,tmx_input,1)

 C             write the time data to the console
               write(line,*) 'time = ',tmx_input
               k=mexPrintf(line//achar(13)) 

          EndDo

       EndDo

Once I compiled the code (with the above body), I ran the function as

 Call_Wrapper(Block_Data)

and I observed all the correct times being written to the console. Try incorporating the above code/logic into your program, substituting the "write" commands with that which updates your array.

[Rick]

Thank you Thank you......... I will try this later......... and get back with you..........Thank you for all your time I have taken up...