def test2(self): from pele.storage import Database from pele.systems import LJCluster np.random.seed(0) natoms = 13 system = LJCluster(natoms) pot = system.get_potential() mindist = system.get_mindist(niter=1) db = Database() db.addMinimum(pot.getEnergy(_x1), _x1) db.addMinimum(pot.getEnergy(_x2), _x2) m1, m2 = db.minima() connect = DoubleEndedConnect(m1, m2, pot, mindist, db, verbosity=10) connect.connect() self.assertTrue(connect.success()) path = connect.returnPath()
local_connect_params["tsSearchParams"] = tsSearchParams local_connect_params["NEBparams"] = NEBparams myconnect = DoubleEndedConnect( min1, min2, pot, mindist, database, local_connect_params=local_connect_params, ) if (myconnect.graph.areConnected(min1, min2)): print "ALERT: The minima are already connected in the database file", dbfile print " Delete it for a fresh run." myconnect.connect() print "" print "found a path!" #the path is now saved in the database. Lets retrieve it and #print it in a more visual format #now retrieve the path for printing print "" mints, S, energies = myconnect.returnPath() nmin = (len(mints)-1)/2 + 1 nts = nmin-1 print "the path has %d minima and %d transition states" % (nmin, nts) eofs = "path.EofS" print "saving energies to", eofs
NEBparams["image_density"] = 15 NEBparams["NEBquenchParams"] = NEBquenchParams local_connect_params = dict() local_connect_params["tsSearchParams"] = tsSearchParams local_connect_params["NEBparams"] = NEBparams myconnect = DoubleEndedConnect(min1, min2, pot, mindist, database, local_connect_params=local_connect_params, ) if (myconnect.graph.areConnected(min1, min2)): print "ALERT: The minima are already connected in the database file", dbfile print " Delete it for a fresh run." myconnect.connect() print "" print "found a path!" # the path is now saved in the database. Lets retrieve it and # print it in a more visual format # now retrieve the path for printing print "" mints, S, energies = myconnect.returnPath() nmin = (len(mints) - 1) / 2 + 1 nts = nmin - 1 print "the path has %d minima and %d transition states" % (nmin, nts) eofs = "path.EofS" print "saving energies to", eofs