Error

 3 Pages<123>
 Previous Topic Next Topic
 Jean Giraud #21 Posted : 05 October 2020 03:30:34(UTC) Rank: Advanced MemberGroups: Registered Joined: 04/07/2015(UTC)Posts: 5,322Was thanked: 842 time(s) in 672 post(s) Originally Posted by: hanskl 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?This is a good observation, I did catch it but didn't mention.Yes the system generates some complex values, can't track it backin the generating algo. Will have another look.The atan(,) from numerical libraries does not generate complex.That's why you had no plot initially, complex blinds the Quick plot.Cheers ... Jean
 hanskl #22 Posted : 05 October 2020 09:09:48(UTC) Rank: Advanced MemberGroups: Registered Joined: 25/10/2017(UTC)Posts: 143Was thanked: 7 time(s) in 5 post(s) Originally Posted by: Jean Giraud This is a good observation, I did catch it but didn't mention.Yes the system generates some complex values, can't track it backin the generating algo. Will have another look.The atan(,) from numerical libraries does not generate complex.That's why you had no plot initially, complex blinds the Quick plot.Cheers ... Jean Another interesting observation is that j.lim in the even Numbers generates values in the areas mentiod that are far too big, while odd j.lims generate conservatively small values. I really appreciate all the help I can get with this. Cheers.
 hanskl #23 Posted : 05 October 2020 10:16:25(UTC) Rank: Advanced MemberGroups: Registered Joined: 25/10/2017(UTC)Posts: 143Was thanked: 7 time(s) in 5 post(s) Originally Posted by: Jean Giraud Is there anything to be done about this, to make the graph closer resemble the diagram?This is a good observation, I did catch it but didn't mention.Yes the system generates some complex values, can't track it backin the generating algo. Will have another look.The atan(,) from numerical libraries does not generate complex.That's why you had no plot initially, complex blinds the Quick plot.Cheers ... Jean Take a look at this sheet, where every factor of the equation has been ploted separately.It seems the complex numbers starts with X.c before the iterations start, due to c going negative.EDIT: It seems in order to produce a plot similar to that of the diagram, dc/dx can not be allowed to be negative. How would I be able to manipulate this?Cheers. problem_complex_separated_5.sm (320kb) downloaded 9 time(s).Edited by user 05 October 2020 15:05:10(UTC)  | Reason: New attachment
 Jean Giraud #24 Posted : 05 October 2020 18:55:36(UTC) Rank: Advanced MemberGroups: Registered Joined: 04/07/2015(UTC)Posts: 5,322Was thanked: 842 time(s) in 672 post(s) From your last report, you are right and wrong.In many more projects, residuals complex have been observedwhere the maths do not generate complex. The numerical Smathengine is in cause. In other applications Smath 32 bits is too granular.Your project is a two time constants XFR [Transfer Function],associated wrt the composite process gain Kp. That is much differentthan XFR in Process Control Techniques. Maybe the model is ill posed.In the mean time, by pushing j 1..99 we get a much better recovery.Cheers ... Jean Project Hanski Soil Drainage.sm (517kb) downloaded 12 time(s).
 Jean Giraud #25 Posted : 05 October 2020 23:12:00(UTC) Rank: Advanced MemberGroups: Registered Joined: 04/07/2015(UTC)Posts: 5,322Was thanked: 842 time(s) in 672 post(s) ... this version is suggestive of completing the project. Project Hanski Soil Drainage.sm (518kb) downloaded 10 time(s).
 hanskl #26 Posted : 06 October 2020 02:01:07(UTC) Rank: Advanced MemberGroups: Registered Joined: 25/10/2017(UTC)Posts: 143Was thanked: 7 time(s) in 5 post(s) Originally Posted by: Jean Giraud ... this version is suggestive of completing the project. Project Hanski Soil Drainage.sm (518kb) downloaded 10 time(s).JMG - great tip to put eval() on all the lines inside the function, sped up the computation considerably, far beyond just right clicking and choosing numeric evaluation.I think unless someone on the forum can find a way to differentiate c, I will stop pursuing this further.I have attached a new sheet with my final findings of this function. In my furter application of N.gamma i will have the function output only real values of the solution, and use deltaX = 0,0001 and j.lim = 99.Thanks to everyone for the help.Cheers. problem_complex_separated_6.sm (281kb) downloaded 13 time(s).
 Jean Giraud #27 Posted : 06 October 2020 20:56:19(UTC) Rank: Advanced MemberGroups: Registered Joined: 04/07/2015(UTC)Posts: 5,322Was thanked: 842 time(s) in 672 post(s) This project is nicely concluded.j.lim 39 is as good as j.lim 99 ... much faster.Should be same in your machine and SS version.Cheers, thanks ... Jean. Project Hanski Soil Drainage Ng.sm (296kb) downloaded 15 time(s).
 hanskl #28 Posted : 25 October 2020 01:58:07(UTC) Rank: Advanced MemberGroups: Registered Joined: 25/10/2017(UTC)Posts: 143Was thanked: 7 time(s) in 5 post(s) Hi guys,Back at it again, I am now trying to plot these functions in Zedgraph with a log scale.I've found that "if" within the function is not supported by Zedgraph ("x is not defined" and "cases" breaks the iteration process.Is there any way to make zedgraph plot my function?
 Jean Giraud #29 Posted : 25 October 2020 21:03:09(UTC) Rank: Advanced MemberGroups: Registered Joined: 04/07/2015(UTC)Posts: 5,322Was thanked: 842 time(s) in 672 post(s) Originally Posted by: hanskl Is there any way to make zedgraph plot my function?log scale 'y' deteriorates the reality. The project is in two different plots.Easy to add labels for the part of large deviations, like this demo.IMHO, much better for publishing the reality. You can label part two as well.Cheers ... Jean Project Hanski Soil Drainage Ng.sm (314kb) downloaded 9 time(s). 1 user thanked Jean Giraud for this useful post. on 25/10/2020(UTC)
 hanskl #30 Posted : 25 October 2020 21:58:57(UTC) Rank: Advanced MemberGroups: Registered Joined: 25/10/2017(UTC)Posts: 143Was thanked: 7 time(s) in 5 post(s) Originally Posted by: Jean Giraud Originally Posted by: hanskl Is there any way to make zedgraph plot my function?log scale 'y' deteriorates the reality. The project is in two different plots.Easy to add labels for the part of large deviations, like this demo.IMHO, much better for publishing the reality. You can label part two as well.Cheers ... JeanThanks Jean for your reply, but it doesnt quite solve my problem.Attached is a sheet showing what I want to achieve, and explaining the troubles I have trying to do that.Best regards, problem_log.sm (692kb) downloaded 15 time(s).
 Jean Giraud #31 Posted : 26 October 2020 03:40:07(UTC) Rank: Advanced MemberGroups: Registered Joined: 04/07/2015(UTC)Posts: 5,322Was thanked: 842 time(s) in 672 post(s) X_Y plots reject conditional. Try to remove them in the program.You can log scale the calculated ordinates, example attached.Cheers ... Jean Project Hanski Soil Drainage Ng.sm (322kb) downloaded 12 time(s).
 hanskl #32 Posted : 26 October 2020 10:35:08(UTC) Rank: Advanced MemberGroups: Registered Joined: 25/10/2017(UTC)Posts: 143Was thanked: 7 time(s) in 5 post(s) Originally Posted by: Jean Giraud X_Y plots reject conditional. Try to remove them in the program.You can log scale the calculated ordinates, example attached.Cheers ... JeanBut the program needs the conditional due to the iteration process I think.I will have a look into Maxima to see if that gives me what I want.On the other hand, I have found another method to calculate Ngamma. I have tried to write a similar program to this method, but when I plot this method it differs vastly from the first method/original diagram. Can someone please have a look and check if the setup for method 2 is correct?Thanks. Method 1 vs Method 2.sm (460kb) downloaded 14 time(s).
 Jean Giraud #33 Posted : 26 October 2020 21:00:10(UTC) Rank: Advanced MemberGroups: Registered Joined: 04/07/2015(UTC)Posts: 5,322Was thanked: 842 time(s) in 672 post(s) Please, visit what I just posted under "Uni" ZedGraph. Project Hanski Soil Drainage Ng User.sm (313kb) downloaded 13 time(s).
 hanskl #34 Posted : 28 October 2020 00:01:17(UTC) Rank: Advanced MemberGroups: Registered Joined: 25/10/2017(UTC)Posts: 143Was thanked: 7 time(s) in 5 post(s) Originally Posted by: Jean Giraud Please, visit what I just posted under "Uni" ZedGraph.Thanks Jean, and everyone else who has spent their time trying to figure this out with me.I think I have come as far as I can on this project, and a lot of new cool smath stuff has been learned.What an amazing program and community!Cheers.
 Jean Giraud #35 Posted : 28 October 2020 03:03:56(UTC) Rank: Advanced MemberGroups: Registered Joined: 04/07/2015(UTC)Posts: 5,322Was thanked: 842 time(s) in 672 post(s) Originally Posted by: hanskl What an amazing program and community!Please visit the just posted latest brew in ZedGraph thread.For our education: what are r, ρ, Nγ ???'r' percolation factor ?'ρ' the unary radius ?'Nγ' some spread ?No such thing as soil drainage function in the web !Thanks, cheers ... Jean 1 user thanked Jean Giraud for this useful post. on 28/10/2020(UTC)
 hanskl #36 Posted : 28 October 2020 11:48:13(UTC) Rank: Advanced MemberGroups: Registered Joined: 25/10/2017(UTC)Posts: 143Was thanked: 7 time(s) in 5 post(s) Originally Posted by: Jean Giraud Please visit the just posted latest brew in ZedGraph thread.For our education: what are r, ρ, Nγ ???'r' percolation factor ?'ρ' the unary radius ?'Nγ' some spread ?No such thing as soil drainage function in the web !Thanks, cheers ... JeanNγ is a factor used in the following equation to calculate soil bearing capacity.σ.v + a = N.q*(p' + a) + 0,5*N.γ*γ'*B.0r = is roughness of the foundation surface, and ρ = internal friction angle of the soil.https://www.geostru.eu/b...for-shallow-foundations/Cheers.
 Jean Giraud #37 Posted : 28 October 2020 16:40:57(UTC) Rank: Advanced MemberGroups: Registered Joined: 04/07/2015(UTC)Posts: 5,322Was thanked: 842 time(s) in 672 post(s) Thanks for the details. As is, the project is complete.You can create an interpolated blue track trace, read Nγ on the plot.I'm not completely happy of the Y log, however immaterial for application.Thanks Hanski for enriching the Smath CommunityJean ... last brew attached. Project Hanski Soil Drainage Ng ZedGraph.sm (400kb) downloaded 8 time(s).
 hanskl #38 Posted : 28 October 2020 17:41:18(UTC) Rank: Advanced MemberGroups: Registered Joined: 25/10/2017(UTC)Posts: 143Was thanked: 7 time(s) in 5 post(s) Originally Posted by: Jean Giraud Thanks for the details. As is, the project is complete.You can create an interpolated blue track trace, read Nγ on the plot.I'm not completely happy of the Y log, however immaterial for application.Thanks Hanski for enriching the Smath CommunityJean ... last brew attached.The whole business of ploting the function is just for documentation purposes, to make sure the function is correct (compared to the referenced diagrams).In practical application, only the function will be used.I think the whole problems is ill posed, as you suggested earlier, as evidenced by the complex results. I have made a modification to the function to combat this in my final version of this function.Attached is the document I will be storing on my server for future reference/use on projects.Best regards, project_complete.sm (944kb) downloaded 16 time(s).
 Jean Giraud #39 Posted : 28 October 2020 22:03:36(UTC) Rank: Advanced MemberGroups: Registered Joined: 04/07/2015(UTC)Posts: 5,322Was thanked: 842 time(s) in 672 post(s) Originally Posted by: hanskl I have made a modification to the function to combat this in my final version of this function.I will visit later tomorrow.In the mean time, if you like this ZedGraph plot ... attached.Cheers ... Jean Project Hanski Soil Drainage Ng ZedGraph.sm (520kb) downloaded 9 time(s).
 Jean Giraud #40 Posted : 29 October 2020 02:37:46(UTC) Rank: Advanced MemberGroups: Registered Joined: 04/07/2015(UTC)Posts: 5,322Was thanked: 842 time(s) in 672 post(s) Originally Posted by: hanskl Attached is the document I will be storing on my server for future reference/use on projects.Your latest "complete" does not calculate same values ?Much shorter in real values. For the remaining of the project, all red SS 6179.Can't even copy/paste algo in same document, neither export.
 Users browsing this topic
 3 Pages<123>
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.