Alexey Cherkaev writes:
> Found the problem: "double precision" on my machine (presumably, any 64-bit
> machine) with "-fdefault-real-8" flag will be "real*16" (something like
> non-existent "Float128" in Julia).
In Fortran, if you want your code to be portable and avoid these type of
issues it
Found the problem: "double precision" on my machine (presumably, any 64-bit
machine) with "-fdefault-real-8" flag will be "real*16" (something like
non-existent "Float128" in Julia).
On Wednesday, November 9, 2016 at 4:05:59 PM UTC+2, Alexey Cherkaev wrote:
>
> On a smaller example I can see tha
On a smaller example I can see that it is indeed is because floats are not
double floats.
It's seem to be the statement
implicit double precision (a-h,o-z)
in Fortran code that throws things upside down. I'm not so familiar with
the Fortran standard and especially with gfortran extensions, but
Hi,
Alexey Cherkaev writes:
> To cut long story short, I have a problem of passing floating-point numbers
> from
> Julia to Fortran subroutine. Documentation suggests, that each parameter to
> Fortran's program needs "Ref{T}" wrapper as they are passed by reference (I've
> seen that older code a
I have encountered a rather strange error as I'm trying to create a wrapper
for Fortran's TWPBVP (ODE BVP) solver.
Full code is available at https://github.com/mobius-eng/TWPBVP.jl
After downloading (to julia's package directory), run
Pkg.build("TWPBVP")
and
Pkg.test("TWPBVP")
The issue (rig