Why Segmentation fault is happening in this openmp code?

Question

main program:

program main                                                                                                                                            


        

Answer 1

Segmentation fault is due to stack memory limit when using OpenMP. Using the solutions from the previous answer did not solve the problem for me on my Windows OS. Using memory allocation into heap rather than stack memory seems to work:

integer, parameter :: nmax = 202000  
real(dp), dimension(:), allocatable :: e_in
integer i

allocate(e_in(nmax))

e_in = 0

! rest of code

deallocate(e_in)

Plus this would not involve changing any default environment parameters.

Acknowledgement to and refer to ohm314's solution here: large array using heap memory allocation