def METHOD_NAME(self):
    part = self.chol.part
    ngs = self.chol.ngs
    Xaoik = numpy.zeros((part.nkk,ngs,part.nij),
                         dtype=numpy.complex128)
    Xaolj = numpy.zeros((part.nkk,ngs,part.nij),
                        dtype=numpy.complex128)
    self.chol.generate_orbital_products(1, self.X, Xaoik, Xaolj)
    self.assertAlmostEqual(numpy.max(numpy.abs(Xaoik)),
                           0.01239256218954048, places=8)
    self.assertAlmostEqual(numpy.max(numpy.abs(Xaolj)),
                           23.93386691184637, places=8)