I'm using FGMRES to solve some large scale problems.And now I have encountered with some errors when I allocate the memory for the tmp array.
The MKL references say that,the FGMRES subroutine must have a parameter named tmp,which length is N* (2 * N+ 1) + (N* (N+ 9)) / 2 + 1,where N is the problem's dimension.
My question is,when my problem's dimension is about 4e11,I can do it well.But when I increase the number of unknows and when it comes to about 2e12,I can't allocate the memory for array tmp.I catch the error "bad alloc".
My computer has 256G RAM,and I also found that when I allocate the memory successfully,the program didn't occupy the memory immediately.So I think it might be some other errors.
And the last question,if I cannot do the FGMRES with large problems,does it mean the FGMRES is not suitable for that?
Thanks in advance,
Solo