"x0 is infeasible" error for lc_geometry (two gaussian) and ebai (mlp) fit

[aydee218]

I'm getting the error "x0 is infeasible." while attempting to fit some light curve using lc_geometry (two gaussian) and ebai (mlp method).
Usually it has to do with bounds and initial values but i have no idea how to fix it using phoebe.
is it even fixable?
Note : my data can be fitted using the polyfit for lc_geometry and knn for ebai.

[kecnry]

I have not seen this happen before and the error isn't directly from phoebe, but I'm guessing from one of the fitting algorithms. Can you post the full error message traceback(s) so we can see where from phoebe this is being triggered? And it would also be super helpful if you could save and attach the bundle (or you can email it personally if you aren't comfortable posting the data publicly) along with the solver labels that trigger the error.

[aydee218]

Here for LC geometry

2G model CE failed with error: x0 is infeasible.
2G model CGE failed with error: x0 is infeasible.
2G model CG12E1 failed with error: x0 is infeasible.
2G model CG12E2 failed with error: x0 is infeasible.

---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
Cell In[24], line 1
----> 1 b.run_solver('lcgeom_solver', solution='lcgeom_solution',overwrite=True)

File ~/.local/lib/python3.10/site-packages/phoebe/parameters/parameters.py:424, in send_if_client.<locals>._send_if_client(self, *args, **kwargs)
    421     return ret_
    423 else:
--> 424     return fctn(self, *args, **kwargs)

File ~/.local/lib/python3.10/site-packages/phoebe/frontend/bundle.py:13649, in Bundle.run_solver(self, solver, solution, detach, return_changes, **kwargs)
  13645     return job_param
  13648 solver_class = getattr(_solverbackends, '{}Backend'.format(solver_ps.kind.title()))
> 13649 params = solver_class().run(self, solver_ps.solver, compute, solution=solution, **{k:v for k,v in kwargs.items() if k not in ['compute']})
  13650 metawargs = {'context': 'solution',
  13651              'solver': solver_ps.solver,
  13652              'compute': compute,
  13653              'kind': solver_ps.kind,
  13654              'solution': solution}
  13657 comment_param = StringParameter(qualifier='comments', value=kwargs.get('comments', solver_ps.get_value(qualifier='comments', default='', **_skip_filter_checks)), description='User-provided comments for this solution.  Feel free to place any notes here.')

File ~/.local/lib/python3.10/site-packages/phoebe/solverbackends/solverbackends.py:555, in BaseSolverBackend.run(self, b, solver, compute, **kwargs)
    552     rpacketlists_per_worker = [self.run_worker(**packet)]
    554 else:
--> 555     rpacketlists_per_worker = [self.run_worker(**packet)]
    557 logger.debug("rank:{}/{} calling _fill_solution".format(mpi.myrank, mpi.nprocs))
    558 return self._fill_solution(solution_ps, rpacketlists_per_worker)

File ~/.local/lib/python3.10/site-packages/phoebe/solverbackends/solverbackends.py:641, in Lc_GeometryBackend.run_worker(self, b, solver, compute, **kwargs)
    637 model_classes = {'two-gaussian': ligeor.models.TwoGaussianModel,
    638                 'polyfit': ligeor.models.Polyfit}
    640 model = model_classes.get(analytical_model)(phases=phases, fluxes=fluxes, sigmas=sigmas)
--> 641 model.fit()
    642 eclipse_dict = model.compute_eclipse_params(interactive=interactive)
    644 fitted_params = [t0_supconj_param, ecc_param, per0_param, rsum_param, teffratio_param]

File ~/.local/lib/python3.10/site-packages/phoebe/dependencies/ligeor/models/twogaussian.py:102, in TwoGaussianModel.fit(self, fit_funcs, param_vals)
     99     self.fit_twoGaussian_models()
    100     # compute all model light curves
--> 102 self.compute_twoGaussian_models()
    103 # compute corresponding BIC values
    104 self.compute_twoGaussian_models_BIC()

File ~/.local/lib/python3.10/site-packages/phoebe/dependencies/ligeor/models/twogaussian.py:621, in TwoGaussianModel.compute_twoGaussian_models(self)
    618 models = {}
    620 for fkey in self.fits.keys():
--> 621     models[fkey] = self.twogfuncs[fkey](self.phases, *self.fits[fkey])
    623 self.models = models

TypeError: TwoGaussianModel.ce() missing 2 required positional arguments: 'Aell' and 'phi0'

and for EBAI mlp

2G model CE failed with error: x0 is infeasible.
2G model CGE failed with error: x0 is infeasible.
2G model CG12E1 failed with error: x0 is infeasible.
2G model CG12E2 failed with error: x0 is infeasible.

---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
Cell In[32], line 1
----> 1 b.run_solver('ebai_mlp', solution='ebai_sol')

File ~/.local/lib/python3.10/site-packages/phoebe/parameters/parameters.py:424, in send_if_client.<locals>._send_if_client(self, *args, **kwargs)
    421     return ret_
    423 else:
--> 424     return fctn(self, *args, **kwargs)

File ~/.local/lib/python3.10/site-packages/phoebe/frontend/bundle.py:13649, in Bundle.run_solver(self, solver, solution, detach, return_changes, **kwargs)
  13645     return job_param
  13648 solver_class = getattr(_solverbackends, '{}Backend'.format(solver_ps.kind.title()))
> 13649 params = solver_class().run(self, solver_ps.solver, compute, solution=solution, **{k:v for k,v in kwargs.items() if k not in ['compute']})
  13650 metawargs = {'context': 'solution',
  13651              'solver': solver_ps.solver,
  13652              'compute': compute,
  13653              'kind': solver_ps.kind,
  13654              'solution': solution}
  13657 comment_param = StringParameter(qualifier='comments', value=kwargs.get('comments', solver_ps.get_value(qualifier='comments', default='', **_skip_filter_checks)), description='User-provided comments for this solution.  Feel free to place any notes here.')

File ~/.local/lib/python3.10/site-packages/phoebe/solverbackends/solverbackends.py:555, in BaseSolverBackend.run(self, b, solver, compute, **kwargs)
    552     rpacketlists_per_worker = [self.run_worker(**packet)]
    554 else:
--> 555     rpacketlists_per_worker = [self.run_worker(**packet)]
    557 logger.debug("rank:{}/{} calling _fill_solution".format(mpi.myrank, mpi.nprocs))
    558 return self._fill_solution(solution_ps, rpacketlists_per_worker)

File ~/.local/lib/python3.10/site-packages/phoebe/solverbackends/solverbackends.py:1250, in EbaiBackend.run_worker(self, b, solver, compute, **kwargs)
   1247 else:
   1248     lcModel = ligeor.models.Polyfit(phases=phases_shifted[s], fluxes=fluxes[s], sigmas=sigmas[s])
-> 1250 lcModel.fit()
   1251 ebai_phases = np.linspace(-0.5,0.5,ebai_phase_bins)
   1252 ebai_fluxes = lcModel.compute_model(ebai_phases, best_fit=True)

File ~/.local/lib/python3.10/site-packages/phoebe/dependencies/ligeor/models/twogaussian.py:102, in TwoGaussianModel.fit(self, fit_funcs, param_vals)
     99     self.fit_twoGaussian_models()
    100     # compute all model light curves
--> 102 self.compute_twoGaussian_models()
    103 # compute corresponding BIC values
    104 self.compute_twoGaussian_models_BIC()

File ~/.local/lib/python3.10/site-packages/phoebe/dependencies/ligeor/models/twogaussian.py:621, in TwoGaussianModel.compute_twoGaussian_models(self)
    618 models = {}
    620 for fkey in self.fits.keys():
--> 621     models[fkey] = self.twogfuncs[fkey](self.phases, *self.fits[fkey])
    623 self.models = models

TypeError: TwoGaussianModel.ce() missing 2 required positional arguments: 'Aell' and 'phi0'

[kecnry]

Can you try disabling binning phase_bin = False in the solver options and see if that circumvents the error? This seems similar to #743 / #749, even though the error messages are slightly different.

If that doesn't seem to solve the issue, please email me the bundle (call b.save('reproduce.bundle') immediately before either of those calls to b.run_solver) and I can try to dig deeper. You should be able to find my email on my GitHub profile.

[aydee218]

I did disable the binning in my case. i'll check the 2 discussions referenced.
If the problem persist, I'll email you the bundle !