topic Re: How to create Fortran array referencing SafeArrayAccessData in Intel® Fortran Compiler

How to create Fortran array referencing SafeArrayAccessData

MWind2 — Fri, 31 Dec 2021 06:29:13 GMT

VBArray INOUT length known
integer(C_INTPTR_T) :: dptr
...
iRes=SafeArrayAccessData(VBArray,dptr)
dptr address is starting point of one dimensional c array of real{8}

I want to put it in a Fortran array so I can see it in debugger for one reason, but I want the Fortran array to reference the underlying data and not be a copy. Intention would be changes to data in Fortran array would be preserved in dptr array.

Re: How to create Fortran array referencing SafeArrayAccessData

Steve_Lionel — Fri, 31 Dec 2021 17:30:17 GMT

Something like this:

use, intrinsic :: iso_c_binding ... real(c_double), pointer :: array_data(:) ... call c_f_pointer(transfer(dptr,c_null_ptr),array_data,[array_size])

where array_size is the number of elements in the array (can be a variable or a constant). After this, array_data is your one-dimensional array, no copy done.

Re: How to create Fortran array referencing SafeArrayAccessData

MWind2 — Fri, 31 Dec 2021 21:00:02 GMT

I seem to be getting an Exception thrown at 0x00007FF88E912D83 (df.dll) in EXCEL.EXE: 0xC0000005: Access violation writing location 0x0000000000000000 because of the "c_null_ptr" in the transfer. Data is "adptr"(changed from dptr as dptr used as another pointer).

Re: How to create Fortran array referencing SafeArrayAccessData

Steve_Lionel — Fri, 31 Dec 2021 21:08:51 GMT

C_NULL_PTR as the second argument to TRANSFER is just the "mold" - it provides the type to interpret the source (first argument) as. It is not accessed as data. If you got an access violation, then adptr was zero. I suggest adding some logging to your code (write to a file perhaps.)

Re: How to create Fortran array referencing SafeArrayAccessData

MWind2 — Fri, 31 Dec 2021 23:06:34 GMT

I attached a vs2022 folder with all inside. Excel is 64 bit and projects.

Re: How to create Fortran array referencing SafeArrayAccessData

Steve_Lionel — Sat, 01 Jan 2022 02:10:43 GMT

Two errors. First, you omitted the size argument I specified to C_F_POINTER. Second, you treated the array as if it were zero-origin, but the array you get from C_F_POINTER is one-origin. (This could be finessed with pointer assignment redeclaring the bounds, but that seems unnecessary.)

I attach my modified df.f90. I didn't change anything else.

Re: How to create Fortran array referencing SafeArrayAccessData

MWind2 — Sat, 01 Jan 2022 03:02:32 GMT

Thank you! I spent all day trying for magic incantations. Intention here was to show that a Excel vba array of doubles used in calculations could be easily sent to a fortran dll and the results sent back in a position(s) in the array for Re: Reading Data from an Excel Spreadsheet? - Page 2 - Intel Communities. One change was made in formula

array_data(4)=PI*array_data(1)*array_data(1)*array_data(2)*array_data(3)/4.0 for weight of cylinder agree with others I had Option 0 on in vba, but array in fortran was [1,4]. Without mucking it up, how to set array_data to be [0,3]?

Re: How to create Fortran array referencing SafeArrayAccessData

Steve_Lionel — Sat, 01 Jan 2022 14:28:58 GMT

Add a new array pointer declaration as follows:

real(8), pointer :: array_data(:), array_data1(:)

then replace the C_F_POINTER call with these two lines:

call c_f_pointer(cptr,array_data1,[4]) array_data(0:) => array_data1

I had thought I could do:

array_data(0:) => array_data

but that didn't work, for reasons still unclear to me (it caused array_data to become disassociated - this MIGHT be a bug, and I will investigate further.)

Edit: I can't reproduce the problem in a simple test. Shrug....

Re: How to create Fortran array referencing SafeArrayAccessData

MWind2 — Sat, 01 Jan 2022 20:22:13 GMT

Thank you! Sometimes I'm "cset" in my ways if permitted, and I would have still been looking for Rosetta stone.

While trying to find things that applied in books, I was thinking that most write to the overview, not that there is anything wrong with that, but a compendium of use questions and answers would be more useful at times. It seems forums nearly do the job, but there are problems with that like too much archiving or deleting of still relevant material, besides signal to noise.

Re: How to create Fortran array referencing SafeArrayAccessData

MWind2 — Mon, 03 Jan 2022 22:28:17 GMT

Is it possible to do this for the two dimensional array? So far I have had no success?

The bounds part works.

Re: How to create Fortran array referencing SafeArrayAccessData

Steve_Lionel — Mon, 03 Jan 2022 23:35:01 GMT

Look again at the VB.NET-Safearrays sample, from which you borrowed some code. It illustrates how to access arrays of any dimension.

Re: How to create Fortran array referencing SafeArrayAccessData

MWind2 — Tue, 04 Jan 2022 01:41:00 GMT

The looping into each one array goes fine, but combining or assigning those to a Fortran 2d array is where I am not getting.

Re: How to create Fortran array referencing SafeArrayAccessData

Steve_Lionel — Tue, 04 Jan 2022 19:34:59 GMT

Please zip and attach what you have now, and I'll take a look.

Re: How to create Fortran array referencing SafeArrayAccessData

MWind2 — Tue, 04 Jan 2022 21:02:59 GMT

I've been using x1.xlsm, dfad and dcpp1, no progress with fortran.

Re: How to create Fortran array referencing SafeArrayAccessData

Steve_Lionel — Tue, 04 Jan 2022 21:41:38 GMT

You have a lot of code commented out, but I see this:

cptr=transfer(adptr,c_null_ptr) call c_f_pointer(cptr,array_data,[12]) array_data(0:) => array_data1

That should be array_data1 in the second line.

Where is the code where you're "assigning to a 2D Fortran array"?

Re: How to create Fortran array referencing SafeArrayAccessData

MWind2 — Tue, 04 Jan 2022 22:21:10 GMT

I was getting nowhere but locking exceptions, so this was what I was doing . Trying to do in fortran what was done in dcpp1 and so far getting nowhere.

Re: How to create Fortran array referencing SafeArrayAccessData

Steve_Lionel — Tue, 04 Jan 2022 23:53:38 GMT

I have not tried to edit your sources, but 1) You still have that error on the call to c_f_pointer, and 2) unlike the C++ code, you unlock the SafeArray before accessing it.

Re: How to create Fortran array referencing SafeArrayAccessData

MWind2 — Wed, 05 Jan 2022 02:22:54 GMT

When I use

readloop: do ires=SafeArrayGetElemsize(VBArray) ires = SafeArrayGetElement (VBArray, indexes(1), loc(adptr)) ! cptr=transfer(adptr,c_null_ptr) call c_f_pointer(cptr,array_data,[ilen]) call c_f_pointer(cptr,array_data1,[ilen]) array_data(0:) => array_data1 ! ! array_data1(3)=PI*array_data1(0)*array_data1(0)*array_data1(1)*array_data1(2) do i = nbounds, 1, -1 indexes(i) = indexes(i) + 1 if (indexes(i) <= bounds(i)%ub) exit indexes(i) = bounds(i)%lb if (i == 1) exit readloop end do end do readloop iRes=SafeArrayUnlock(VBArray)

I get an element size of 8, which does not seem to be an array size. Maybe I'm reading it wrong but GetElement may be just a double.

Re: How to create Fortran array referencing SafeArrayAccessData

MWind2 — Wed, 05 Jan 2022 03:00:35 GMT

Got a fortran array with such but maybe order wrong:

iRes=SafeArrayLock(VBArray) indexes = bounds%lb ! Initialize to all lower bounds ilen=bounds(2)%ub-bounds(2)%lb ires=SafeArrayGetElemsize(VBArray) !ires = SafeArrayPtrOfIndex(VBArray, indx, cptr0) ires=SafeArrayAccessData(VBArray,adptr) !ires = SafeArrayGetElement (VBArray, indexes(1), loc(adptr)) ! cptr=transfer(adptr,c_null_ptr) call c_f_pointer(cptr,array_data,[12]) !call c_f_pointer(cptr,array_data1,[12]) !array_data(0:) => array_data1 ! ! array_data1(3)=PI*array_data1(0)*array_data1(0)*array_data1(1)*array_data1(2) iRes=SafeArrayUnlock(VBArray)

Re: How to create Fortran array referencing SafeArrayAccessData

MWind2 — Wed, 05 Jan 2022 04:34:58 GMT

Not pretty but curious.