intento llamar una funcion desde Fortran hacia C++,
cuando defino n_points=10^6 no genera problema, pero con n_points= 10^8 aparece Violación de segmento (`core' generado)
.
Aqui el codigo fortran pi_rand.f90
module monte_carlo_mod
use, intrinsic :: iso_c_binding
implicit none
contains
function monte_carlo_pi(n_points) result(pi_estimate) bind(c, name="monte_carlo_pi")
integer :: n_points ! = 100000000
integer :: i, inside_circle
real :: x, y, pi_estimate, radius
real :: rand
integer :: seed
seed = 12345
inside_circle = 0
call random_seed(seed)
do i = 1, n_points
call random_number(x)
call random_number(y)
x = 2.d0 * x - 1.d0
y = 2.d0 * y - 1.d0
if (x**2 + y**2 <= 1.d0) then
inside_circle = inside_circle + 1
end if
end do
pi_estimate = 4.d0 * real(inside_circle) / real(n_points)
print *, 'Estimated value of Pi:', pi_estimate
end function monte_carlo_pi
end module
Aqui codigo en c++ main.cpp
// main.cpp
#include <iostream>
// Declarar la función Fortran que se va a utilizar
extern "C" {
double monte_carlo_pi(int n_points);
}
int main() {
int n_points = 100000000;
double pi_estimate;
// Llamar a la función Fortran
pi_estimate = monte_carlo_pi(n_points);
// Mostrar el resultado
// std::cout << "Estimated value of Pi: " << pi_estimate << std::endl;
return 0;
}
Codigos para compilar que use :
gfortran -c pi_rand.f90
g++ -c main.cpp
g++ -o salida main.o pi_rand.o -lfortran
time ./salida
.
Me da error
Violación de segmento (`core' generado)
. Que hice mal ?