Java-реализация разложения по сингулярным значениям для больших разреженных матриц

Мне просто интересно, знает ли кто-нибудь о Java-реализации разложения по сингулярным значениям (SVD) для больших разреженных матриц? Мне нужна эта реализация для скрытого семантического анализа (LSA).

Я пробовал пакеты из UJMP и JAMA, но они задыхаются, когда количество row> = 1000 и col> = 500. Если кто-нибудь может указать мне на псевдокод или что-то еще там, это было бы очень признательно.