nm

11353 Reputation

20 Badges

13 years, 12 days

MaplePrimes Activity


These are replies submitted by nm

@sursumCorda 

Thanks. Unfortunately this does not work all the time. Here is an example

s:="A& {test} & \\textrm{The fundamental matrix has } &C & {test} ":

StringTools:-RegSubs("(.*)\\\\textrm\\{(.*)\\}(.*)" = "\\1\\begin{minipage}{\\linewidth}\\textrm{\\2}\\end{minipage}\\3", s);

It gives

It did not match the first closing }.   

May be this is what you meant by the follow up you had where you said there is a problem.  

I need to do this inside Maple. So using Python afterwords will not really work for my workflow.

 

@ecterrab 

Ok, that is fair enough for this one example 

ode:=diff(y(x),x)=(y(x)-3)^2;  

But then why did Maple then say y=3 is singular solution when using Lie method?

Also for the second example I shown

ode:=diff(y(x),x)=2*x*sqrt(1-y(x)^2);

with y=1, y=-1. These can not be obtained from the general solution for any value of c. 

I have dozens more examples I can show if you like where 

'singsol'=all

does not produce singular solutions. So I am not sure now if one needs to use Lie method always to see singular solution or not. I also have examples where Maple gives singular solution without the option 'singsol' at all. Here is one

dsolve(ode,y(x))
dsolve(ode,y(x),'singsol'=all)
dsolve(ode,'Lie','singsol'=all)

All the above produce singular solutions y=r, y=-r. So the first command should not have produced these, since the option singsol was not used.

It is not possible to obtain y=r,y=-r from general solution for any value of c.

So for me, it is just confusing how and when singsol is supposed to work in Maple.

@vv 

I do not understand. if do as you suggested, I get wrong solution like I did before

ode:=x*y(x)*diff(y(x), x) = (x + 1)*(y(x) + 1);
ic:=y(1) = 1;
sol_no_IC:=dsolve(ode,'implicit');
eq:= eval(sol_no_IC,[y(x)=1,x=1]);
PDEtools:-Solve(eq,c__1);
eval(sol_no_IC,%);
PDEtools:-Solve(%,y(x));
odetest(%,[ode,ic])

 

I am looking for an algorithmic approach, not trial and error. I know Maple used -1 branch when plugging in the constant of integration. I just do not know how to decide to use -1 or 1 or 0 for LambertW after finding the constant of integration and before plugging the constant back into the generate solution. i.e. Why Maple decided to use -1 branch when doing so?

well, Maple is made of the front end process and the engine process. These are completely separate processes. They communicate using interprocess communication. The engine is what server.exe is. The memory shown on the botton of the worksheet is for the frontend (how much ram used for content of the worksheet itself and for the front end process itself).  

So these ofcourse will be different memory as they are different processes. In the windows task manger, the front end process is the one that has the little Java icon next to it. The engine is called server.exe. If you have many worksheets open and have set the option to use separate engine for each, then you will see many server.exe's running but one front end process.

@ecterrab 

Thanks for finding where the bug is. Hopefully this will be fixed in some future version of Maple. 

I never said this was solution generated by Maple. I use Maple's odetest to test different types of solutions not all generated by dsolve. Some of these solutions can ofcourse be wrong or not valid in some way or another.

@Preben Alsholm 

I also use windows 10, pro version with 128 GB. This is strange that C_R does not see the crash and I see it and we are both on windows 10. Something else is different.

But good you are able to reproduce it, even if on windows 11. so it is not just me. I do not think the windows version makes a difference. There is something else going on.

I noticed using Maple 2024, I get many more server crashes on same PC than I had on Maple 2023.2.1 and I have not made any changes really to my PC or the software I use. I do not have any special anti-virus software, other than what windows 10 itself has.  (This is different from the random hangs which happens much more frequent).

I wish Maple would generate a log file when the server crashes. May be a core dump somewhere to use for analysis to see the cause of the crash?

But nothing is generated as far as I know. Even the pop-up message does not give any indication of the cause of the server crash. This is not good.  

@C_R 

Thanks for trying, This is very strange, Why some get this crash and some do not.  I have no idea what is going on. May be someone will figure it out.

I'll try to turn on tracing and see if I can find something new.

@C_R 

I noticed you do not have Physics updates installed on your system. This could be the difference., I do not know. Try with Physics updates installed with latest version and see.

@ecterrab 

FYI, I've added two movies in my post.  I see you are using mac OS. I am on windows. This could be the difference.

my ini file has only one command which just prints the process id. Nothing else.

Please see updated movies. Another person above was also able to reproduce the crash on windows.

@acer 

" Does it crash if you force only method=contour? "

Yes, It still crashes.

u:=Int(2/3/alpha^(3/2)*sin(1/2*3^(1/2)*ln(alpha))*sin(alpha)*3^(1/2),alpha=0..x,method=contour);
value(u);

#server crash

@acer 

How does one check "without Physics update"? Other than removing it from libname, how to check? I tried below removing Physics update lib from libname, but crash still happens.

 


 

26028

interface(version);

`Standard Worksheet Interface, Maple 2024.0, Windows 10, March 01 2024 Build ID 1794891`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1715 and is the same as the version installed in this computer, created 2024, April 3, 20:27 hours Pacific Time.`

restart;

24896

libname

"C:\Users\Owner\maple\toolbox\2024\Physics Updates\lib", "C:\Program Files\Maple 2024\lib"

libname:="C:\\Program Files\\Maple 2024\\lib"

"C:\Program Files\Maple 2024\lib"

u:=Int(2/3/alpha^(3/2)*sin(1/2*3^(1/2)*ln(alpha))*sin(alpha)*3^(1/2),alpha = 0 .. x)

Int((2/3)*sin((1/2)*3^(1/2)*ln(alpha))*sin(alpha)*3^(1/2)/alpha^(3/2), alpha = 0 .. x)

value(u);


 

Download another_server_crash_on_int_maple_2024_april_4_2024.mw

@acer 

That is correct. verified on Maple 2024

 

integrand:=(2*x^2022+1)/(x^2023+x);
convert(integrand,parfrac);  #kernel crash

windows 10

@ecterrab 

Ah!  I found the problem. few days ago I asked about expansion of exp here Why-Maple-Likes-To-Write-Exp2a-As-expa2-

So I used the suggestion in one of the answers

"There's also a way to suppress the expansion of exp (or any other function) that maintains the suppression until you turn it off; so it's suitable for use in an initialization file:

expand(expandoff()): expandoff(exp):"

And edited my Maple ini file and added the above. Now that you mentioned ini file, I removed the above and the problem is now gone. 

So that was the cause. Here is worksheet that shows it is fixed now after removing the above from my ini file

restart;

142864

interface(version);

`Standard Worksheet Interface, Maple 2024.0, Windows 10, March 01 2024 Build ID 1794891`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1708 and is the same as the version installed in this computer, created 2024, March 27, 16:20 hours Pacific Time.`

expand(expandoff()): expandoff(exp):

integrand:=(d*x)^m/(a+b*arctanh(c*x^n))^2;

(d*x)^m/(a+b*arctanh(c*x^n))^2

int(integrand,x);

Error, (in int/gparse/gmon) too many levels of recursion

restart;

142864

integrand:=(d*x)^m/(a+b*arctanh(c*x^n))^2;

(d*x)^m/(a+b*arctanh(c*x^n))^2

int(integrand,x);

2*x*(c*exp(n*ln(x))-1)*(c*exp(n*ln(x))+1)*exp(m*(ln(d)+ln(x)-((1/2)*I)*Pi*csgn(I*d*x)*(-csgn(I*d*x)+csgn(I*d))*(-csgn(I*d*x)+csgn(I*x))))/(b*c*n*exp(n*ln(x))*(-b*ln(1-c*exp(n*ln(x)))+b*ln(c*exp(n*ln(x))+1)+2*a))+int(-2*exp(m*(ln(d)+ln(x)-((1/2)*I)*Pi*csgn(I*d*x)*(-csgn(I*d*x)+csgn(I*d))*(-csgn(I*d*x)+csgn(I*x))))*(c^2*m*(exp(n*ln(x)))^2+c^2*n*(exp(n*ln(x)))^2+(exp(n*ln(x)))^2*c^2-m+n-1)/(b*c*n*exp(n*ln(x))*(-b*ln(1-c*exp(n*ln(x)))+b*ln(c*exp(n*ln(x))+1)+2*a)), x)

 


Download int_gparse_gmon_error_FIXED_maple_2024_march_27_2024.mw

I had no idea that disabling expansion of exp will cause such internal issue. Thanks for everyone's comments and answers.

@Carl Love 

opps. fixed it now. Was sleepy and did not notice.

@Ronan 

Thanks. But I tried RealDomain and it did not work for me. The code is still on my screen. Here it is 

35788

restart;

35788

n:=3;m:=2;
eqx:=x^(n/m)=a;
use RealDomain in PDEtools:-Solve(eqx,x); end use;
F:=map(X->eval(eqx,X),[%]);
map(X->evalb(X),F);

3

2

x^(3/2) = a

x = a^(2/3), x = (1/4)*a^(2/3)*(1+I*3^(1/2))^2, x = (1/4)*a^(2/3)*(I*3^(1/2)-1)^2

[a = a, (1/16)*4^(1/2)*(a^(2/3)*(1+I*3^(1/2))^2)^(3/2) = a, (1/16)*4^(1/2)*(a^(2/3)*(I*3^(1/2)-1)^2)^(3/2) = a]

[true, false, false]

 

Download why_real_domain_does_not_work_for_me.mw

What is the difference between what I did above and what you did for it to have worked with you?

The strange thing, is that in Mathematica by default it works in complex domain. Maple should be the same. But Mathematica gave x=a^(2/3) as only solution without saying real domain.

It looks like when I used solve() with real domain, it works. But when using PDEtools:-Solve with real domain, it does not work!  is this a bug?

First 19 20 21 22 23 24 25 Last Page 21 of 91