forked from evangelistalab/scf_orbital_overlap
-
Notifications
You must be signed in to change notification settings - Fork 0
/
pymodule.py
61 lines (51 loc) · 1.76 KB
/
pymodule.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
#
#@BEGIN LICENSE
#
# scf_orbital_overlap by Psi4 Developer, a plugin to:
#
# PSI4: an ab initio quantum chemistry software package
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License along
# with this program; if not, write to the Free Software Foundation, Inc.,
# 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
#
#@END LICENSE
#
import psi4
import re
import os
import inputparser
import math
import warnings
from driver import *
from wrappers import *
from molutil import *
import p4util
from p4xcpt import *
def run_scf_orbital_overlap(name, **kwargs):
r"""Function encoding sequence of PSI module and plugin calls so that
scf_orbital_overlap can be called via :py:func:`~driver.energy`. For post-scf plugins.
>>> energy('scf_orbital_overlap')
"""
lowername = name.lower()
kwargs = p4util.kwargs_lower(kwargs)
# Your plugin's psi4 run sequence goes here
#psi4.set_global_option('BASIS', 'sto-3g')
psi4.set_local_option('SCF_ORBITAL_OVERLAP', 'PRINT', 1)
returnvalue = psi4.plugin('scf_orbital_overlap.so')
psi4.set_variable('CURRENT ENERGY', returnvalue)
# Integration with driver routines
procedures['energy']['scf_orbital_overlap'] = run_scf_orbital_overlap
def exampleFN():
# Your Python code goes here
pass