在R中大型稀疏矩阵上实现余弦相似性的快速方法

Question

我有一个大的稀疏矩阵（300,000 X 60,000）并试图计算用户之间的成对相似性。最终的相似性矩阵将是（300,000×300,000）。我尝试使用来自"qlcMatrix"的{​​{1}}和来自cosine的{​​{1}}，但似乎没有一个在时间方面有效运作。有没有办法实现或任何已经实现相同的包或者我应该转移到其他编程语言？