whats good guys.
problem: i am giving an N x N sparse matrix as input to FEAST and I am trying to solve for d-dimensional manifolds.
Note that this dimensionality, d, is independent of N.
Now when d isn't too large; say <3000, FEAST works completely fine.
However, if I "up" the value of d to say, 9000, i get errors such as:
Intel MKL Extended Eigensolvers: Size subspace 9001
#Loop | #Eig | Trace | Error-Trace | Max-Residual
OMP: Error #34: System unable to allocate necessary resources for OMP thread:
OMP: System error #35: Resource temporarily unavailable
OMP: Hint: Try decreasing the value of OMP_NUM_THREADS.
Abort trap: 6
I have tried setting export OMP_STACKSIZE=1024m, and I have also tried using the -stack_size argument in clang++ to specify a stacksize that's pretty large. Neither of these solutions worked. I have also looked at how much memory is being consumed when I run FEAST as above with either settings and the limit was roughly the same. It suggests to me that I'm not setting OMP_stacksize right, or maybe OSX does this differently?
Using clang++ with c++11 threading, and also the MKL for all the math. any assistance on this situation would be great.