@@ -122,34 +122,28 @@ def test_binary_solution_convergence(self):
122122
123123 bdg = burnman .Solution (
124124 "bdg" ,
125- solution_model = burnman .classes .solutionmodel .IdealSolution (bdg_endmembers )
125+ solution_model = burnman .classes .solutionmodel .IdealSolution (bdg_endmembers ),
126126 )
127127 ppv = burnman .Solution (
128128 "ppv" ,
129- solution_model = burnman .classes .solutionmodel .IdealSolution (ppv_endmembers )
129+ solution_model = burnman .classes .solutionmodel .IdealSolution (ppv_endmembers ),
130130 )
131131
132132 composition = {"Mg" : 0.9 , "Fe" : 0.1 , "Si" : 1.0 , "O" : 3.0 }
133133 assemblage = burnman .Composite (
134- phases = [bdg , ppv ],
135- fractions = [0.9 , 0.1 ],
136- name = "MgSiO3-pv-ppv-assemblage"
134+ phases = [bdg , ppv ], fractions = [0.9 , 0.1 ], name = "MgSiO3-pv-ppv-assemblage"
137135 )
138136 bdg .set_composition ([0.9 , 0.1 ])
139137 ppv .set_composition ([0.9 , 0.1 ])
140138
141139 temperatures = np .linspace (1000 , 4000 , 4 )
142- assemblage .set_state (120.e9 , temperatures [0 ])
140+ assemblage .set_state (120.0e9 , temperatures [0 ])
143141
144- equality_constraints = [
145- ("T" , temperatures ),
146- ("phase_fraction" , (ppv , 0.0 ))
147- ]
142+ equality_constraints = [("T" , temperatures ), ("phase_fraction" , (ppv , 0.0 ))]
148143
149144 sol , prm = equilibrate (composition , assemblage , equality_constraints )
150145 self .assertFalse (any (not s .success for s in sol ))
151146
152147
153-
154148if __name__ == "__main__" :
155149 unittest .main ()
0 commit comments