Rank: Advanced Member Groups: Registered
Joined: 25/10/2017(UTC) Posts: 109 Was thanked: 7 time(s) in 5 post(s)

Hi. I need help with a iterative calculation. This is my first time trying to build a iterative loop in SMath. I need to calculate a coefficient that can only be calculated when the difference between X.0 and X.c is below a limit. I am trying to solve this as a function of r and rho. Can someone help me get started on how to formulate the iteration? problem.sm (12kb) downloaded 10 time(s).Best regards,




Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,647 Was thanked: 915 time(s) in 575 post(s)



1 user thanked mkraska for this useful post.



Rank: Advanced Member Groups: Registered
Joined: 25/10/2017(UTC) Posts: 109 Was thanked: 7 time(s) in 5 post(s)

Thank you for your suggestion Martin!
I need the output X to be a function of r and rho. Is that possible?
Some random questions... 1. why do you use arctg instead of atan? 2. does the iterative process have to be wrapped into a for loop where you need to give it a range, and ultimately an upper boundry (n.max)?




Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,647 Was thanked: 915 time(s) in 575 post(s)

Originally Posted by: hanskl Thank you for your suggestion Martin!
I need the output X to be a function of r and rho. Is that possible?
Some random questions... 1. why do you use arctg instead of atan? 2. does the iterative process have to be wrapped into a for loop where you need to give it a range, and ultimately an upper boundry (n.max)? 0. I think you would need to wrap the full procedure into a function of r and rho with X as return value. Alternatively, you might come up with a recursive definition of the procedure. 1. It is just the setting world versus european style for functions. It is the same function, just displayed according to the GUI settings. 2. the for loop just came in handy for generating the indices for storage of iteration history. You could use while 1 as well. I also thought it was a good idea to have an upper limit for the iterations in case the algorithm might diverge. 




Rank: Advanced Member Groups: Registered
Joined: 25/10/2017(UTC) Posts: 109 Was thanked: 7 time(s) in 5 post(s)

Originally Posted by: mkraska Originally Posted by: hanskl Thank you for your suggestion Martin!
I need the output X to be a function of r and rho. Is that possible?
Some random questions... 1. why do you use arctg instead of atan? 2. does the iterative process have to be wrapped into a for loop where you need to give it a range, and ultimately an upper boundry (n.max)? 0. I think you would need to wrap the full procedure into a function of r and rho with X as return value. Alternatively, you might come up with a recursive definition of the procedure. 1. It is just the setting world versus european style for functions. It is the same function, just displayed according to the GUI settings. 2. the for loop just came in handy for generating the indices for storage of iteration history. You could use while 1 as well. I also thought it was a good idea to have an upper limit for the iterations in case the algorithm might diverge. Thanks again for the replies, on a friday night no less. I really appreciate that! rinks: 0. Can you explain this in layman terms? Ultimately, I want to graph N.gamma as a function of r and rho, where I have 10 set values of r and 10 corresponding functions plotted with rho along xaxis. 1. OK. I was just curious if they were different functions, like max/Max for example, which I have found can make or break a book. 2. thanks again for the explanation. Cheers




Rank: Advanced Member Groups: Registered
Joined: 04/07/2015(UTC) Posts: 5,056 Was thanked: 809 time(s) in 642 post(s)

Originally Posted by: hanskl I need the output X to be a function of r and rho. Is that possible? Sorry, don't understand. If you want to see the evolution of the solution, use a while loop, SQRT(x) in Samples is a good example. Visit Golden ratio as well. problem_Kr Plot What.sm (18kb) downloaded 10 time(s).




Rank: Advanced Member Groups: Registered
Joined: 25/10/2017(UTC) Posts: 109 Was thanked: 7 time(s) in 5 post(s)

Originally Posted by: Jean Giraud Originally Posted by: hanskl I need the output X to be a function of r and rho. Is that possible? Sorry, don't understand. If you want to see the evolution of the solution, use a while loop, SQRT(x) in Samples is a good example. Visit Golden ratio as well. problem_Kr Plot What.sm (18kb) downloaded 10 time(s). Hi Jean, Thanks for the sample. I think I got my calculation working now with (a lot) of help from Martin. I cant plot N.gamma(r, rho) because it obviously takes an enormous amount of computational power. However, the formulas seem to return the correct value when you compare to the diagrams in this sheet (the benchmark for my calculation). PS. Is there some way to speed up the calculation of N.gamma? For some values of r and rho it can take Almost a minute to compute. Cheers problem_solved.sm (512kb) downloaded 10 time(s).Edited by user 03 October 2020 02:36:43(UTC)
 Reason: Not specified




Rank: Advanced Member Groups: Registered
Joined: 04/07/2015(UTC) Posts: 5,056 Was thanked: 809 time(s) in 642 post(s)

Originally Posted by: hanskl I cant plot N.gamma(r, rho) because it obviously takes an enormous amount of computational power. Get a quick insight ... that you can mesh finer. problem_solved More Insight.sm (518kb) downloaded 6 time(s).




Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,647 Was thanked: 915 time(s) in 575 post(s)

Originally Posted by: hanskl I think I got my calculation working now with (a lot) of help from Martin.
I cant plot N.gamma(r, rho) because it obviously takes an enormous amount of computational power.
However, the formulas seem to return the correct value when you compare to the diagrams in this sheet (the benchmark for my calculation).
PS. Is there some way to speed up the calculation of N.gamma? For some values of r and rho it can take Almost a minute to compute.
Cheers
There is no need to define the quantities as explicit functions if you just want to evaluate an expression. Removing the arguments together with setting the whole algorithm to numeric optimization in the context mnenu made a speedup from 30s to 3s. problem_solved_Kr.sm (510kb) downloaded 13 time(s). 




Rank: Advanced Member Groups: Registered
Joined: 25/10/2017(UTC) Posts: 109 Was thanked: 7 time(s) in 5 post(s)

Originally Posted by: mkraska Originally Posted by: hanskl I think I got my calculation working now with (a lot) of help from Martin.
I cant plot N.gamma(r, rho) because it obviously takes an enormous amount of computational power.
However, the formulas seem to return the correct value when you compare to the diagrams in this sheet (the benchmark for my calculation).
PS. Is there some way to speed up the calculation of N.gamma? For some values of r and rho it can take Almost a minute to compute.
Cheers
There is no need to define the quantities as explicit functions if you just want to evaluate an expression. Removing the arguments together with setting the whole algorithm to numeric optimization in the context mnenu made a speedup from 30s to 3s. Can you help me understand understand the difference between having the quantities as explicit functions like I did and your approach, and what will the advantange/disadvantages be for one or the other? I want to try to understand the best way to structure my calculations for the result I want. BTW, just installed the latest SMath version, what a tremendous speedup that gave on this sheet. My oldere version would go on for minutes without being able to plot N.gamma. Thanks




Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,647 Was thanked: 915 time(s) in 575 post(s)

Originally Posted by: hanskl Can you help me understand understand the difference between having the quantities as explicit functions like I did and your approach, and what will the advantange/disadvantages be for one or the other? I want to try to understand the best way to structure my calculations for the result I want.
BTW, just installed the latest SMath version, what a tremendous speedup that gave on this sheet. My oldere version would go on for minutes without being able to plot N.gamma.
Thanks
I am on very thin ice when making statements on that, because it is just based on observations, not on conceptual insight. If you define functions, these are stored as is and evaluated when the function is called. If define a variable, the expression is simplified using the current context (available definitions of variables) You can see that in the debugger window. You always see what you store in the dynamic callout when hovering over the definition. I know that in the end you want a quantity depending on r and rho, thus a function. Yet, inside the algorithm r and rho are fixed and not subject to change for a given call of that algorithm. So it is straightforward to use the global value of these variables throughout the algorithm, specified as arguments when calling the algorithm. 

1 user thanked mkraska for this useful post.



Rank: Advanced Member Groups: Registered
Joined: 04/07/2015(UTC) Posts: 5,056 Was thanked: 809 time(s) in 642 post(s)

Originally Posted by: hanskl My oldere version would go on for minutes without being able to plot N.gamma. My steam engine laptop SS 6179 runs Martin in 1.8 min. It runs the solver 168*10= 1680 for all 10 plots Set j=9 ... constrain gamma 0.9 reluctant Set the solver in eval(,) mode. Cheers ... Jean problem_solved_Kr Sanity.sm (511kb) downloaded 12 time(s).




Rank: Advanced Member Groups: Registered
Joined: 04/07/2015(UTC) Posts: 5,056 Was thanked: 809 time(s) in 642 post(s)





Rank: Advanced Member Groups: Registered
Joined: 04/07/2015(UTC) Posts: 5,056 Was thanked: 809 time(s) in 642 post(s)

N.γ(r,ρ) should plot in ~ 30 s on modern machines. Saved in my Project file ... any abstract about the project ? Project Hanski Unknown.sm (514kb) downloaded 14 time(s).




Rank: Administration Groups: Developers, Registered, Knovel Developers, Administrators, Advanced Member Joined: 11/07/2008(UTC) Posts: 1,475 Was thanked: 1742 time(s) in 590 post(s)

Originally Posted by: Jean Giraud N.γ(r,ρ) should plot in ~ 30 s on modern machines. Saved in my Project file ... any abstract about the project ? Project Hanski Unknown.sm (514kb) downloaded 14 time(s). Calculated in about 5 seconds on the latest stable SMath Studio. Originally Posted by: Jean Giraud Originally Posted by: hanskl My oldere version would go on for minutes without being able to plot N.gamma. My steam engine laptop SS 6179 runs Martin in 1.8 min. It runs the solver 168*10= 1680 for all 10 plots Set j=9 ... constrain gamma 0.9 reluctant Set the solver in eval(,) mode. Cheers ... Jean problem_solved_Kr Sanity.sm (511kb) downloaded 12 time(s). Calculated in about 2 seconds on the latest stable SMath Studio. Edited by user 04 October 2020 00:35:02(UTC)
 Reason: Not specified




Rank: Advanced Member Groups: Registered
Joined: 04/07/2015(UTC) Posts: 5,056 Was thanked: 809 time(s) in 642 post(s)

Originally Posted by: Andrey Ivashov Thanks Andrey ... very impressive ! What about Julia set ? Cheers ... Jean. Pattern Julia Set [PUBLISH].sm (173kb) downloaded 20 time(s).




Rank: Advanced Member Groups: Registered
Joined: 25/10/2017(UTC) Posts: 109 Was thanked: 7 time(s) in 5 post(s)

Originally Posted by: mkraska Originally Posted by: hanskl Can you help me understand understand the difference between having the quantities as explicit functions like I did and your approach, and what will the advantange/disadvantages be for one or the other? I want to try to understand the best way to structure my calculations for the result I want.
BTW, just installed the latest SMath version, what a tremendous speedup that gave on this sheet. My oldere version would go on for minutes without being able to plot N.gamma.
Thanks
I am on very thin ice when making statements on that, because it is just based on observations, not on conceptual insight. If you define functions, these are stored as is and evaluated when the function is called. If define a variable, the expression is simplified using the current context (available definitions of variables) You can see that in the debugger window. You always see what you store in the dynamic callout when hovering over the definition. I know that in the end you want a quantity depending on r and rho, thus a function. Yet, inside the algorithm r and rho are fixed and not subject to change for a given call of that algorithm. So it is straightforward to use the global value of these variables throughout the algorithm, specified as arguments when calling the algorithm. Thanks, that actually made a lot of sense! Will apply that to other sheets i have that struggle with Slow computation. Cheers!




Rank: Advanced Member Groups: Registered
Joined: 25/10/2017(UTC) Posts: 109 Was thanked: 7 time(s) in 5 post(s)

Originally Posted by: Jean Giraud N.γ(r,ρ) should plot in ~ 30 s on modern machines. Saved in my Project file ... any abstract about the project ? Project Hanski Unknown.sm (514kb) downloaded 14 time(s). The functions calculate soil bearing capacity factors for drained conditions. N.q is used internationally, but N.gamma is according to Norwegian recommendations. Cheers.




Rank: Administration Groups: Developers, Registered, Knovel Developers, Administrators, Advanced Member Joined: 11/07/2008(UTC) Posts: 1,475 Was thanked: 1742 time(s) in 590 post(s)

Originally Posted by: Jean Giraud Originally Posted by: Andrey Ivashov Thanks Andrey ... very impressive ! What about Julia set ? Cheers ... Jean. Pattern Julia Set [PUBLISH].sm (173kb) downloaded 20 time(s). 55 seconds. As I see it was about 14 minutes before...




Rank: Advanced Member Groups: Registered
Joined: 25/10/2017(UTC) Posts: 109 Was thanked: 7 time(s) in 5 post(s)

I have tinkered a bit because there seems to be a mismatch between my graph and the diagram for some values of r and rho. The plots of N.γ(0.9,atan(x)) and N.γ(0.95,atan(x)) has an odd shape when they approach x=1. It turns out there are some complex solutions in this area. Is there anything to be done about this, to make the graph closer resemble the diagram? Cheers problem_complex.sm (261kb) downloaded 12 time(s).




Forum Jump
You cannot post new topics in this forum.
You cannot reply to topics in this forum.
You cannot delete your posts in this forum.
You cannot edit your posts in this forum.
You cannot create polls in this forum.
You cannot vote in polls in this forum.