search

PALab / RUS

Resonant Ultrasound Spectroscopy
12 stars 7 forks source link

Shape issue in forward and inverse problem #18

Open dangdx-code opened 8 months ago

dangdx-code commented 8 months ago

Dear Dr. Shepherd, Here is Dang. I am a Ph.D student in Paul Scherrer Institut (PSI). Currently, I am in charging to build up a RUS in my lab. Your computing is great to resolve the forward and inverse problems. However, I am still confusing the shape number when I try to plug to forward and inverse. Rectangular shape is 0 in forward and 2 in inverse, Right? At some points, the reverse program does not working well. I attach the struggle, which I got. I also attach the reference parameter, which I used for running test. I am looking forward your response soon. Thank you so much in advance.

P/s: reference source for runing test, table S1: https://static-content.springer.com/esm/art%3A10.1038%2Fs43246-020-00083-1/MediaObjects/43246_2020_83_MOESM1_ESM.pdf

Reverse running test

kaspervanwijk commented 8 months ago

To my knowledge, the shape parameter is consistent between the forward and inverse codes.

But in your example, there is an error about dimensions. I suggest reading the fortran version manual and starting with the example in there. There is an art to the inversion process, due to the fact that any experimental data misses some resonant peaks.

Cheers, Kasper van Wijk

From: dangdx-code @.> Sent: Tuesday, February 20, 2024 2:29:53 AM To: PALab/RUS @.> Cc: Subscribed @.***> Subject: [PALab/RUS] Shape issue in forward and inverse problem (Issue #18)

Dear Dr. Shepherd, Here is Dang. I am a Ph.D student in Paul Scherrer Institut (PSI). Currently, I am in charging to build up a RUS in my lab. Your computing is great to resolve the forward and inverse problems. However, I am still confusing the shape number when I try to plug to forward and inverse. Rectangular shape is 0 in forward and 2 in inverse, Right? At some points, the reverse program does not working well. I attach the struggle, which I got. I also attach the reference parameter, which I used for running test. I am looking forward your response soon. Thank you so much in advance.

P/s: reference source for runing test, table S1: https://static-content.springer.com/esm/art%3A10.1038%2Fs43246-020-00083-1/MediaObjects/43246_2020_83_MOESM1_ESM.pdfhttps://static-content.springer.com/esm/art%3A10.1038%2Fs43246-020-00083-1/MediaObjects/43246_2020_83_MOESM1_ESM.pdf

Reverse.running.test.png (view on web)https://github.com/PALab/RUS/assets/125368137/8b2c0201-a397-4e7e-a956-33be95b2715a

— Reply to this email directly, view it on GitHubhttps://github.com/PALab/RUS/issues/18, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ABG3TXQEY3JQ5V4HZ6KVQ2LYUNHVDAVCNFSM6AAAAABDPPA67SVHI2DSMVQWIX3LMV43ASLTON2WKOZSGE2DEMZXHE3TCMY. You are receiving this because you are subscribed to this thread.Message ID: @.***>

dangdx-code commented 8 months ago

Hi Kasper, Thank you so much for your response and guidance. when I run the help in both cases. In the forward, shape = 0 is rectangular but it is spherical in the case of inverse. and the unit of dimension is centimeter, I try out but it did not work at all. Best, Dang

kaspervanwijk commented 8 months ago

If the definition of the rectangle is inconsistent, then that is ugly but not a real bug. About cm/m: if you run the example on the aluminium, you know the units. From memory, they are in m.

Cheers, Kasper

From: dangdx-code @.> Sent: Wednesday, February 21, 2024 3:22:42 AM To: PALab/RUS @.> Cc: Kasper van Wijk @.>; Comment @.> Subject: Re: [PALab/RUS] Shape issue in forward and inverse problem (Issue #18)

Hi Kasper, Thank you so much for your response and guidance. when I run the help in both cases. In the forward, shape = 0 is rectangular but it is spherical in the case of inverse. and the unit of dimension is centimeter, I try out but it did not work at all. Best, Dang

— Reply to this email directly, view it on GitHubhttps://github.com/PALab/RUS/issues/18#issuecomment-1954321445, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ABG3TXWV4DWJ3OEDTS74KPDYUSWTFAVCNFSM6AAAAABDPPA67SVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTSNJUGMZDCNBUGU. You are receiving this because you commented.Message ID: @.***>

dangdx-code commented 8 months ago

Thank you so much for your great guidance. Could you please guide me to solve issue of list index out of range, which has showed like as image? And, I am still confusing about the example when I run the forward problem, I figured out that it did not need the elastic tensor. Could you please help me explain why (maybe, is it my stupid question)? Thank you very much for valuable discussion. Regards, Dang