Dear Chao Cao, Thank you very much for your message. Your situation is the same of my. But I run the code on IBM PPC970 cluster. I can try to use other mathematical libraries, may be it helps. Best regards, Sergey