def test_init_input_density_matrix(self):
        # Setup
        Job = InitJob("test_data/JobDef_input_density_matrix.json")
        input_rho_file = Job.Def['input_rho']
        with open(input_rho_file, 'r') as file_handle:
            input_rho = np.matrix(json.load(file_handle))

        # Action
        electronic = Electronic(Job)

        # Result
        assert np.array_equal(electronic.rho, input_rho)
        assert np.array_equal(electronic.rhotot, input_rho)