def METHOD_NAME(self):
    profile = {"Mo": {"r": 0.6, "w": 1.0}}
    self.describer1 = BispectrumCoefficients(
        rcutfac=4.6, twojmax=6, element_profile=profile, quadratic=False, pot_fit=True
    )
    model1 = SKLModel(describer=self.describer1, model=LinearRegression())
    self.potential1 = SNAPotential(model=model1, name="test")
    self.describer2 = BispectrumCoefficients(
        rcutfac=4.6, twojmax=6, element_profile=profile, quadratic=True, pot_fit=True
    )
    model2 = SKLModel(describer=self.describer2, model=LinearRegression())
    self.potential2 = SNAPotential(model=model2, name="test")
    self.test_pool = test_datapool
    self.test_structures = []
    self.test_energies = []
    self.test_forces = []
    self.test_stresses = []
    for d in self.test_pool:
        self.test_structures.append(d["structure"])
        self.test_energies.append(d["outputs"]["energy"])
        self.test_forces.append(d["outputs"]["forces"])
        self.test_stresses.append(d["outputs"]["virial_stress"])
    self.test_struct = self.test_pool[-1]["structure"]