I am solving pde. For the post processing bvp, there is solution given in output which is the solution from mathematica version 10.0.1.
But when the similar code executed in mathematica 12.1, gives error "General:: Exp[-717.401] is too small to represent as a normalized machine number; precision may be lost." and does not give any results. As I checked it, i tried to solve the problem by changing the exponential fuction in "SMTAddInitialBoundary". So some changes are made with the fuction by removing "-Exp" as shown ,
By doing some changes in the function, I could able to understand that the problem is with "SMTAddInitialBoundary" syntax.
But at the end question is why first code does not give the solution in mathematica 12.1? I tried to find the problem but could not able to reach.
Is there a problem with mathematica version or with "SMTAddInitialBoundary"?