search

pypest / pyemu

python modules for model-independent uncertainty analyses, data-worth analyses, and interfacing with PEST(++)
BSD 3-Clause "New" or "Revised" License
169 stars 94 forks source link

issue with zero_order_tikhonov #410

Closed hjia1005 closed 1 year ago

hjia1005 commented 1 year ago

I used pyemu.helpers.zero_order_tikhonov(pst,parbounds=True) to set up the Tikhonov regularisation. All parameters, observations and prior information in pest control file are in lower letter. image

However, when I start to run pestpp-glm, I received this error: image

I opened 008.rec and found this: image

and this:

image

Looks like there is a conflict between upper letter and lower letter. I have no idea why those lower letter in the pest control file are converted to upper letters in the rec file.

Any way to fix it? Thanks in advance.

hjia1005 commented 1 year ago

If I remove the below option from the pst file, then pestpp-glm runs normally without the precalculated jacobian matrix. ++base_jacobian(007_reuse.jco)

But the same thing shows up in the rec file.

image

This is what the jacobian matrix looks like - all lower letters: image

It has exactly the same parameters and observations as the pest file.

Why the jacobian matrix is not compatible with the prior information?

jtwhite79 commented 1 year ago

Im not sure its a case issue - pest++ tries to work entirely in upper case internally (for historical reasons) while pyemu is lower case - its a heaps of fun...

what version of pestpp-glm are you using? I just tested this use case on 5.2 and its seems to find the pi equations in the jcb file...

hjia1005 commented 1 year ago

@jtwhite79 , thanks for your reply.

Yes I was using version 5.1.21. Just downloaded the latest version 5.2.1, still have the same issue. Very strange.

jtwhite79 commented 1 year ago

can you post your pst and jcb files so I can have a look?

hjia1005 commented 1 year ago

Can I send them to you via email?

On Tue, 28 Feb 2023 at 10:47, J Dub @.***> wrote:

can you post your pst and jcb files so I can have a look?

— Reply to this email directly, view it on GitHub https://github.com/pypest/pyemu/issues/410#issuecomment-1447359630, or unsubscribe https://github.com/notifications/unsubscribe-auth/AH7GLMQVRUMQLCVNKUBLPXDWZVDI3ANCNFSM6AAAAAAVI633N4 . You are receiving this because you authored the thread.Message ID: @.***>

hjia1005 commented 1 year ago

Solved, thank you very much! @jtwhite79