Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Function exp is not supported #139

Closed
kikiggg opened this issue Jan 10, 2023 · 4 comments
Closed

Function exp is not supported #139

kikiggg opened this issue Jan 10, 2023 · 4 comments

Comments

@kikiggg
Copy link

kikiggg commented Jan 10, 2023

Hi

I want to check StructuralIdentifiability some differential equation which include "exp" but it didnt work with message "Function exp is not supported"

Is there any solutiong for that? or internally not supported?

Thank you
Park

@iliailmer
Copy link
Collaborator

Hi @kikiggg , StructuralIdentifiability.jl supports ODE models with rational (polynomial) functions. You can try to modify your ODE so that exp is "implicit", for example

# original:
x' = A * exp(B * t) * x

# modified
x' = A * g * x
g' = B * g

@ChrisRackauckas
Copy link
Member

We could make that a ModelingToolkit pass of sorts, worth an issue.

@pogudingleb
Copy link
Collaborator

Small addition to the comment of @iliailmer
It is better to set g = A * exp(B * t), this way the resulting system will be

x' = g * x
g' = B * g

So the number degrees of freedom (num of states plus num of parameters will be the same), so the identifiability results can be transferred between the original system and the new one.

@ChrisRackauckas Here is a subtlety: one should ideally try to make the transformation carefully, keeping the dimension of the problem. Not impossible, definitely worth an issue, I will create one.

@pogudingleb
Copy link
Collaborator

I created a separate issue for handling non-rational models containing my thoughts on the matter. It subsumes the present issue, so I close it.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

4 participants