- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
I have modified the numerical recipes ODE integrator. When I ran the original NR code, it worked fine. I have now changed something and I now get odd problems with my arrays.
In the odeint.f90, there is a statement (see http://tero.fis.uncor.edu/NR/f90/odeint.f90):
y(:)=ystart(:)
Now, when the array is passed to save_a_step (internal subroutine), I get an address error.
If I change the above line to:
y = ystart
the error disappears. Sadly, I think this is the cause of further errors, so I would like to know what is happening. What is the difference between:
y and y(:)
????
Thanks.
Link Copied
1 Reply
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
There should not be a difference. If you are usimg the latest version of the Intel compiler and can come up with an example that demonstrates the problem, please send it to Intel Premier Support.
I will note that the (:) syntax is not meaningful in this context and, I feel, is poor coding style. But it is legal Fortran.
I will note that the (:) syntax is not meaningful in this context and, I feel, is poor coding style. But it is legal Fortran.

Reply
Topic Options
- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page