def METHOD_NAME(self):
    residual = self.residual_reductor.reduce()
    initial_residual = self.initial_residual_reductor.reduce()
    if self.fom.output_functional.linear:
        output_adjoint = self.fom.output_functional.H
        output_adjoint_riesz_range = estimate_image(vectors=(output_adjoint,), orthonormalize=True,
                                                    product=self.products['RB'], riesz_representatives=True)
        projected_output_adjoint = project(output_adjoint, output_adjoint_riesz_range, None)
    else:
        projected_output_adjoint = None
    return ParabolicRBEstimator(residual, self.residual_reductor.residual_range_dims, initial_residual,
                                self.initial_residual_reductor.residual_range_dims, self.coercivity_estimator,
                                projected_output_adjoint)