@@ -394,9 +394,9 @@ def spinorb_from_spatial(one_body_integrals, two_body_integrals):
394394 two_body_coefficients [2 * p , 2 * q , 2 * r , 2 * s ] = two_body_integrals [
395395 p , q , r , s
396396 ]
397- two_body_coefficients [
398- 2 * p + 1 , 2 * q + 1 , 2 * r + 1 , 2 * s + 1
399- ] = two_body_integrals [ p , q , r , s ]
397+ two_body_coefficients [2 * p + 1 , 2 * q + 1 , 2 * r + 1 , 2 * s + 1 ] = (
398+ two_body_integrals [ p , q , r , s ]
399+ )
400400
401401 # Truncate.
402402 one_body_coefficients [numpy .absolute (one_body_coefficients ) < EQ_TOLERANCE ] = 0.0
@@ -1025,11 +1025,9 @@ def get_molecular_hamiltonian(self, occupied_indices=None, active_indices=None):
10251025 one_body_integrals , two_body_integrals = self .get_integrals ()
10261026 constant = self .nuclear_repulsion
10271027 else :
1028- (
1029- core_adjustment ,
1030- one_body_integrals ,
1031- two_body_integrals ,
1032- ) = self .get_active_space_integrals (occupied_indices , active_indices )
1028+ (core_adjustment , one_body_integrals , two_body_integrals ) = (
1029+ self .get_active_space_integrals (occupied_indices , active_indices )
1030+ )
10331031 constant = self .nuclear_repulsion + core_adjustment
10341032
10351033 one_body_coefficients , two_body_coefficients = spinorb_from_spatial (
0 commit comments