Guillermo P. Ortiz

Hello Robert, I appreciate your effort very much, but unfortunately this result has unexpected results also. The imaginary part should be positive. Let me ask for some guide to understand maxima syntax used in your work 1) What is the meaning of the first two lines? linel:100 and display2d:false 2) what is %emode and why is located there first as false and after as true 3)mathdeclare I suspect that is for aa. but you are using a lot of other function like lambda, constantp and apply1 I guess that...

Thank you very much for your patch, Robert. Just I test your version for mysqrt() and it seems working fine. But I am not sure i rectform are using the default sqrt instead of mysqrt, it could be? The question arise from the fact that I am not specting a negative sign of the imaginary part of the output as one obtained here https://sourceforge.net/p/maxima/bugs/_discuss/thread/e6e3db30f2/a456/attachment/dodier-epsMwithoutInt.dat Regards

Sorry, here the bug report related to https://sourceforge.net/p/maxima/bugs/3933/#15dd Please report bugs to: https://sourceforge.net/p/maxima/bugs To report a bug, you must have a Sourceforge account. Please include the following information with your bug report: Maxima version: "5.44.0" Maxima build date: "2021-04-24 14:52:58" Host type: "x86_64-pc-linux-gnu" Lisp implementation type: "GNU Common Lisp (GCL)" Lisp implementation version: "GCL 2.6.12" User dir: "/home/gortiz/.maxima" Temp dir: "/tmp"...

Hi Robert, I understand you. Here, for simplicity, consider the error output +++++++ Maxima encountered a Lisp error: Condition in MACSYMA-TOP-LEVEL [or a callee]: INTERNAL-SIMPLE-ERROR: The storage for CONS is exhausted. 1614811 pages allocated. Use ALLOCATE to expand the space. Automatically continuing. ++++++ obtained for the simpler case code that i post here https://sourceforge.net/p/maxima/bugs/3933/#3743 without intepolated data file issue. As you can see, it is not a misleading on how to...

Thanks to @dauti, I am noticing that environment variable can manage the memory used. By default the maxima setup this variable more aggressive in order to allows multiple task. May be I can probe to increase this values, but not sure. I have by default export GCL_GC_PAGE_THRESH=0.2 export GCL_GC_ALLOC_MIN=0.01 export GCL_GC_PAGE_MAX=0.05 and export GCL_MEM_MULTIPLE=0.2 what values for each would you suggest me? Regards

OK, let me share a code version similar to the above but without cspline. I am totally agree with you about numerical versus symbolic, the problem for me is I am not how to manage the sqrt redefinition according to the neeede branch cut that I said before. dynamalloc:true; eq1:ED1*exp(%i*k1*d1)-ED2*exp(%i*k2*d1)+EI1*exp(-%i*k1*d1)-EI2*exp(-%i*k2*d1)=DE0*exp(%i*k*d1); eq2:k1*ED1*exp(%i*k1*d1)-k2*ED2*exp(%i*k2*d1)-k1*EI1*exp(-%i*k1*d1)+k2*EI2*exp(-%i*k2*d1)=k*DE0*exp(%i*k*d1); eq3:ED1-ED2*exp(%i*(k2-k)*d)+EI1-EI2*exp(-%i*(k2+k)*d)=DE0;...

Hello Stavros, I think that main problem is with mysqrt definition. I have not found any other form to manage the selection of other branch cut for sqrt function defined by Maxima. I am trying a code version without cspline and seems after 3 hs with core i7 7ma generation processor and 16 GB RAM , that Maxima using 8GB is working to hard to manage it I will kill it and I could try with optimize before, thank you in advance Regards

Hello Robert, I have set dynamalloc to true and a lot enough of ram memory, the maxima process in question only take a half of it, but thinking about your comment I change to a non interpolating version of my code. One issue that I am not sure is how to manage the sqrt redefinition to mysqrt. It is needed because the default definition of the branch cut selected by maxima but not adjustable by user, is not it? As a mathematical task, briefly I want to solve a 4x4 linear eq. system from there I obtain...