Error

 9 Pages<12345>»
 Previous Topic Next Topic
 mkraska #41 Posted : 05 August 2013 21:19:02(UTC) Rank: Advanced MemberGroups: Registered Joined: 15/04/2012(UTC)Posts: 1,935Was thanked: 1094 time(s) in 698 post(s) I played around a bit with the examples by Radovan.My observations:- using stack() or not is irrelevant in the single first order equation case.- D and initial condition may be scalar or vector-valued. - if() and Max() and eval() work, if D is defined inline (in the solver argument list), they do not work, if D is defined outside.- the different behaviour of Rkadapt is due to poor convergence of the adaptivity algorithm.- if rkfixed (and possibly the other solvers as well) fails, then you can evaluate the expression using = and get an error- if rkfixed succeeds, there is no way to evaluate it using = (key is ignored), you have to use assignment insteadEdited by user 05 August 2013 21:29:55(UTC)  | Reason: Not specifiedFile Attachment(s): if-max-eval.sm (13kb) downloaded 64 time(s). Martin Kraska Pre-configured portable distribution of SMath Studio: https://smath.com/wiki/SMath_with_Plugins.ashx 1 user thanked mkraska for this useful post. on 05/08/2013(UTC)
 uni #42 Posted : 05 August 2013 21:37:42(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,469Was thanked: 1256 time(s) in 731 post(s) Updated. Source code refactoring.The problem here is not the solvers, the problem is in the values ​​of the function parameters. Russia ☭ forever Viacheslav N. Mezentsev 1 user thanked uni for this useful post. on 05/08/2013(UTC)
 uni #43 Posted : 05 August 2013 21:57:23(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,469Was thanked: 1256 time(s) in 731 post(s) Originally Posted by: mkraska I played around a bit with the examples by Radovan.My observations:- using stack() or not is irrelevant in the single first order equation case.- D and initial condition may be scalar or vector-valued. ...Do not recommend using scalars. The point is that the ode functions in SMath Studio are described as follows:Quote: new ArgumentInfo[] { new ArgumentInfo( ArgumentSections.ColumnVector ), new ArgumentInfo( ArgumentSections.RealNumber ), new ArgumentInfo( ArgumentSections.RealNumber ), new ArgumentInfo( ArgumentSections.RealNumber ), new ArgumentInfo( ArgumentSections.Function ) } );For now there is no strict monitoring parameter types, but it can occur in the future. I can make support for scalar parameters, but then for every function I need to write it analog with different set of types of parameters. Russia ☭ forever Viacheslav N. Mezentsev 1 user thanked uni for this useful post. on 05/08/2013(UTC)
 mkraska #44 Posted : 05 August 2013 22:41:23(UTC) Rank: Advanced MemberGroups: Registered Joined: 15/04/2012(UTC)Posts: 1,935Was thanked: 1094 time(s) in 698 post(s) Originally Posted by: uni Do not recommend using scalars. The point is that the ode functions in SMath Studio are described as follows:Quote: new ArgumentInfo[] { new ArgumentInfo( ArgumentSections.ColumnVector ), new ArgumentInfo( ArgumentSections.RealNumber ), new ArgumentInfo( ArgumentSections.RealNumber ), new ArgumentInfo( ArgumentSections.RealNumber ), new ArgumentInfo( ArgumentSections.Function ) } );For now there is no strict monitoring parameter types, but it can occur in the future. I can make support for scalar parameters, but then for every function I need to write it analog with different set of types of parameters.Ok, I see. How about a generic ODESolve(method, y0, t0, t1, n, D), where method is some optional string indicating the numerical procedure. Then, whatever convenience UI stuff you add goes to all solvers. Similar to the FindRoot function you might adopt the embedded assignment option, e.g. by specifying the initial conditions and iniial value of the independent variable in equation form like ODESolve(mat(v≡1,s≡0,2,3,1),t≡0,1,10,D(t,y))That would create vectors v, s and t, ready for plotting or further processing. In case of two requested vectors you might augment them into something like v_s which would immediately be good for plotting.Edited by user 05 August 2013 22:45:24(UTC)  | Reason: Not specified Martin Kraska Pre-configured portable distribution of SMath Studio: https://smath.com/wiki/SMath_with_Plugins.ashx
 ignasi #45 Posted : 06 August 2013 11:53:41(UTC) Rank: NewbieGroups: Registered Joined: 04/07/2013(UTC)Posts: 6Location: StockholmWas thanked: 7 time(s) in 3 post(s) Dear uni and omorrI can not express in words how thankful I am and my admiration for your work. SMath studio will be my preferred tool from now on.To help others I have the intention of making a “stirred tank reactor” example with microorganisms growing and consuming a substrate; I already have such an example in Mathcad from the time I was teaching. There is already a similar example posted in http://smath.info/wiki/Examples.ashx, but an example with the new plugin functions would be useful I believe.Thanks again! 2 users thanked ignasi for this useful post. on 06/08/2013(UTC),  on 07/08/2013(UTC)
 omorr #46 Posted : 06 August 2013 12:28:59(UTC) Rank: AdministrationGroups: Registered, Advanced MemberJoined: 23/06/2009(UTC)Posts: 1,740Was thanked: 318 time(s) in 268 post(s) You are welcome and thank you I am looking forward to your examples.Regards,Radovan When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
 mkraska #47 Posted : 07 August 2013 01:40:12(UTC) Rank: Advanced MemberGroups: Registered Joined: 15/04/2012(UTC)Posts: 1,935Was thanked: 1094 time(s) in 698 post(s) I tried to summarize the results of the recent discussions as to how to use the ODE solvers. Up to now I assume that everything might work with any of the solvers but that is not tested thoroughly. It would be desirable to have some sort of benchmark problems and to compare the performance of the individual solvers. That might result in sort of usage recommendations, if not in a generic solver just like FindRoot.Feel free to challenge and correct the statements in the attached handbook pages.Edit: Minor corrections in the introduction (example description had a sign error) and in the requirements for D (there was a spurious result display for an expression with embedded :=, which I could not reproduce)Edited by user 07 August 2013 12:52:16(UTC)  | Reason: Not specifiedFile Attachment(s): section math ode.zip (67kb) downloaded 77 time(s). Martin Kraska Pre-configured portable distribution of SMath Studio: https://smath.com/wiki/SMath_with_Plugins.ashx 3 users thanked mkraska for this useful post. on 07/08/2013(UTC),  on 07/08/2013(UTC),  on 07/08/2013(UTC)
 omorr #48 Posted : 07 August 2013 08:49:07(UTC) Rank: AdministrationGroups: Registered, Advanced MemberJoined: 23/06/2009(UTC)Posts: 1,740Was thanked: 318 time(s) in 268 post(s) Martin,Once more, all my respect to you for having a good will and patience to extract all these things and put them together. I have to admit that sometimes I envy you for that Regards,RadovanEdited by user 07 August 2013 14:42:25(UTC)  | Reason: Not specified When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" 1 user thanked omorr for this useful post. on 07/08/2013(UTC)
 uni #49 Posted : 14 August 2013 05:51:32(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,469Was thanked: 1256 time(s) in 731 post(s) Updated. Added functions from the Intel® Ordinary Differential Equations Solver Library: rkm9st(), mk52lfn(), mk52lfa(), rkm9mkn(), rkm9mka().mk52lfa() not finished yet.Edited by user 14 August 2013 05:57:35(UTC)  | Reason: Not specifiedFile Attachment(s): ODE_manual.pdf (240kb) downloaded 86 time(s). ODESolvers. Van der Pol oscillator.sm (11kb) downloaded 82 time(s). Russia ☭ forever Viacheslav N. Mezentsev 2 users thanked uni for this useful post. on 14/08/2013(UTC),  on 14/08/2013(UTC)
 omorr #50 Posted : 14 August 2013 11:16:00(UTC) Rank: AdministrationGroups: Registered, Advanced MemberJoined: 23/06/2009(UTC)Posts: 1,740Was thanked: 318 time(s) in 268 post(s) I have never imagined that so many ODE solvers would be available in SMath. I wish you to stay inspired because you are keeping our positive spirit high Many thanks for that.By the way, are you going to to introduce some additional parameters like RelTol,AbsTol in DotNumerics. As I could see from the Reference Manuel for these solvers, among many other parameters (I suppose you set them by default values) there are parameters like ep, tr for absolute and relative tolerance. Regards,RadovanEdited by user 14 August 2013 11:30:48(UTC)  | Reason: Not specified When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" 1 user thanked omorr for this useful post. on 14/08/2013(UTC)
 uni #51 Posted : 14 August 2013 11:52:10(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,469Was thanked: 1256 time(s) in 731 post(s) Originally Posted by: omorr I have never imagined that so many ODE solvers would be available in SMath. I wish you to stay inspired because you are keeping our positive spirit high Many thanks for that.By the way, are you going to to introduce some additional parameters like RelTol,AbsTol in DotNumerics. As I could see from the Reference Manuel for these solvers, among many other parameters (I suppose you set them by default values) there are parameters like ep, tr for absolute and relative tolerance. Regards,RadovanYou can use RelTol and AbsTol in the same way. I included them in to the plugin. This code is running in the unmanaged space. Therefore, it is platform-specific. During loading the plugin from its resources on disk is copied library for specific platform (iode.dll). I'm wondering to know does it work in 64-bit system?Now I can use any function that is written in C/C++ or other languages. Russia ☭ forever Viacheslav N. Mezentsev
 Davide Carpi #52 Posted : 14 August 2013 12:37:41(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 13/01/2012(UTC)Posts: 2,613Location: ItalyWas thanked: 1318 time(s) in 869 post(s) Works on x64 systems (at least on my notebook), except for the mk52lfa(Y,t0,tmax,n,D(t,y,4)) (Y not defined)best regards,Davide If you like my plugins consider to support SMath Studio buying a plan; to offer me a coffee: paypal.me/dcprojects 1 user thanked Davide Carpi for this useful post. on 14/08/2013(UTC)
 uni #53 Posted : 14 August 2013 12:46:46(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,469Was thanked: 1256 time(s) in 731 post(s) Originally Posted by: w3b5urf3r_reloaded Works on x64 systems (at least on my notebook), except for the mk52lfa(Y,t0,tmax,n,D(t,y,4)) (Y not defined)best regards,DavideThank you, Davide. Unfortunately, this function requires the Jacobian. It must be specified explicitly. Some other functions also require it, but it can be omitted. I'll have to change the number of parameters to a function call. Russia ☭ forever Viacheslav N. Mezentsev
 omorr #54 Posted : 14 August 2013 14:35:37(UTC) Rank: AdministrationGroups: Registered, Advanced MemberJoined: 23/06/2009(UTC)Posts: 1,740Was thanked: 318 time(s) in 268 post(s) Originally Posted by: uni You can use RelTol and AbsTol in the same way. I included them in to the plugin. This code is running in the unmanaged space. Therefore, it is platform-specific. During loading the plugin from its resources on disk is copied library for specific platform (iode.dll). I'm wondering to know does it work in 64-bit system?You are right. RelTol and AbsTol worked with these functions as well . I tried it and can confirm that mk52lfn() and rkm9mkn() are quite fast and worked as expected with Abstol=Reltol=10^-7 . There is no secret which one of the examples i tried . This worked on my 64-bit laptop. I also could see the iode.dll in the plugins folder. I used few SMath SMath versions (including few ones by Martin) and am not sure if 32-bit is enabled or not. Anyway, I think it is not so important at the moment Originally Posted by: uni Now I can use any function that is written in C/C++ or other languages.All my respect for that . I hope that you will find a nice "candy store" now and allow us to taste it as well Regards,Radovan When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
 mkraska #55 Posted : 14 August 2013 17:42:48(UTC) Rank: Advanced MemberGroups: Registered Joined: 15/04/2012(UTC)Posts: 1,935Was thanked: 1094 time(s) in 698 post(s) Here is a nice little toy for exploring precision control. Have fun!Edited by user 14 August 2013 17:52:58(UTC)  | Reason: Not specifiedFile Attachment(s): Section Math ODE precision.sm (18kb) downloaded 84 time(s).mkraska attached the following image(s): Martin Kraska Pre-configured portable distribution of SMath Studio: https://smath.com/wiki/SMath_with_Plugins.ashx 2 users thanked mkraska for this useful post. on 15/08/2013(UTC),  on 15/08/2013(UTC)
 mkraska #56 Posted : 15 August 2013 09:44:08(UTC) Rank: Advanced MemberGroups: Registered Joined: 15/04/2012(UTC)Posts: 1,935Was thanked: 1094 time(s) in 698 post(s) Originally Posted by: uni Unfortunately, this function requires the Jacobian. It must be specified explicitly. Some other functions also require it, but it can be omitted. I'll have to change the number of parameters to a function call.Does that mean that currently, the routines use an internal approximation of the Jacobi matrix, whenever they need it? Could the matrix be generated analytically or numerically by the functions from Davide's Nonlinear Solvers? Would that require user input? Martin Kraska Pre-configured portable distribution of SMath Studio: https://smath.com/wiki/SMath_with_Plugins.ashx
 ignasi #57 Posted : 21 August 2013 09:54:08(UTC) Rank: NewbieGroups: Registered Joined: 04/07/2013(UTC)Posts: 6Location: StockholmWas thanked: 7 time(s) in 3 post(s) Dear colleaguesI made an example on how SMath ODE solvers may be used.Comments and improvements are welcomeIf you thing it would be useful, feel free to post it or use it wherever it is best suited.I am not sure of what would be appropriate. File Attachment(s): Box_models.sm (109kb) downloaded 99 time(s). 4 users thanked ignasi for this useful post. on 21/08/2013(UTC),  on 21/08/2013(UTC),  on 21/08/2013(UTC),  on 21/08/2013(UTC)
 omorr #58 Posted : 21 August 2013 15:19:26(UTC) Rank: AdministrationGroups: Registered, Advanced MemberJoined: 23/06/2009(UTC)Posts: 1,740Was thanked: 318 time(s) in 268 post(s) Nice example As a chemical engineer I am quite "attached" to chemical, biochemical reactors . That is one of the reason I was nagging before about ode solvers (at the moment we have many of them - thanking to uni .Quite a chance that this model will fail with some nonstiff solvers. Fortunately, and thanking to uni, there are three plugins (at the moment) with ode solvers (OdeSolvers, DotNumerics and Mathlab C++ Math library). You can try all of them. Regards,Radovan When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
 madan26 #59 Posted : 15 May 2015 11:36:35(UTC) Rank: NewbieGroups: Registered Joined: 15/05/2015(UTC)Posts: 2 I have solved one equation in sage math but not able understand how to solve in smath anybody can guide me on this. below is my equationexercise1r = var('r'PI = pi.n()rc = 281rw = 262Bc = 36/180*PIF = (r^2-rw^2)*sqrt(rc^2-r^2*sin(Bc)^2)/sqrt((rc^2-rw^2)^2*rc^2*cos(Bc)^2-(r^2-rw^2)^2*(rc^2-r^2*sin(Bc)^2))F.show() z = function('z',r)d = desolve_rk4(diff(z,r)-F,z,ics=[rc,0],end_points=255,step=-0.01)list_plot(d,plotjoined=True, aspect_ratio=1) madan26 attached the following image(s):
 Jean Giraud #60 Posted : 08 October 2015 06:11:23(UTC) Rank: GuestGroups: Registered Joined: 04/07/2015(UTC)Posts: 6,792Was thanked: 974 time(s) in 802 post(s) From the start, you are all right and all wrong !1. You have to understand what you are doing and how it can be done numerically.2. No need for all these Moulton gadget, Fhelberg ... etc. For instance: Fhelbergextends the RK domain where it has no meaning [Trust my years old testing].3. All those librairies, proven ? only few of them [NAG, ACM... OK ? !]4. My Smath version came with Maxima plugins: some work some don't. On that, Smathis not an integrated version. Your clients will be soon fed up of shoping for zillionsof unproven plugins and eventual incompatibility between versions.5. On the other hand, as long as Smath is 32 bits it is not really a CAS.6. Too much automation in solvers may exclude finding a solution. What can be done with little manual pushing may take days to de-automate. In Mathcad 11, Levenberg-Marquardtworks generally well but many fits have to start manually, then Newton, then CG, refined LM. 6. All those 3D plots from "Uni" make me envious, but why should I scratch my right ear with my left hand or ask someone else to do it [unless paralysed]. The GNUplot works, but that'sthe only one. Martin couldn't put color in my famous "Breather". Very interesting.Cheers, Jean
 Users browsing this topic Guest
 9 Pages<12345>»
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.