def METHOD_NAME():
    global mol, mo, nao
    mol = gto.Mole()
    mol.verbose = 0
    mol.output = None#'out_h2o'
    mol.atom.extend([
        ['O' , (0. , 0.     , 0.)],
        [1   , (0. , -0.757 , 0.587)],
        [1   , (0. , 0.757  , 0.587)] ])
    mol.basis = 'cc-pvdz'
    mol.build()
    nao = mol.nao_nr()
    numpy.random.seed(15)
    mo = numpy.random.random((nao,nao))
    mo = mo.copy(order='F')