Skip to main content
41 events
when toggle format what by license comment
Dec 2, 2019 at 4:50 comment added xzczd @HabibUllah As indicated by the warning, your PC doesn't have a C compiler installed or the C compiler isn't properly configured. You can simply remove the CompilationTarget -> C option, but I do recommend you to install a C compiler (or configure it properly, if you've already installed one), that'll make the code faster.
Dec 1, 2019 at 9:25 comment added Mathematicain Sir when is i run the FDM-based solution it gives the error (CCompilerDriver`CreateLibrary::nocomp: A C compiler cannot be found on your system. Please consult the documentation to learn how to set up suitable compilers. >> Compile::nogen: A library could not be generated from the compiled function. >>)
Nov 2, 2018 at 7:27 comment added xzczd @zhk It's possible to modify the program of course, but given the current version is already fast, I suggest you to simply modify the end time and use the obtained data to build a new InterpolatingFunction.
Nov 2, 2018 at 7:19 comment added zhk I meant to say that if the end time is 5 then I’m unable to see what happened at t=0.5. What I’m looking for is to see variations in the dependent variables wrt t, x, y. Let say, I want to visualise u vs t, x , y using parametricplot3d.
Nov 2, 2018 at 7:13 comment added xzczd @zhk What do you mean by "vary time"?
Nov 2, 2018 at 5:55 comment added zhk It is possible to vary time as we are varying x and y in the FDM-based solution?
Aug 30, 2018 at 14:36 history edited xzczd CC BY-SA 4.0
added 8 characters in body
Aug 30, 2018 at 12:20 history edited xzczd CC BY-SA 4.0
added 1 character in body
Oct 8, 2017 at 12:17 history edited xzczd CC BY-SA 3.0
Typo fixed
Apr 13, 2017 at 12:55 history edited CommunityBot
replaced http://mathematica.stackexchange.com/ with https://mathematica.stackexchange.com/
Mar 24, 2017 at 7:22 comment added zhk Let us continue this discussion in chat.
Mar 24, 2017 at 7:15 comment added xzczd @Maple solver[5, 1/150, #, 6] & /@ {2, 4, 6} triggers CompiledFunction::cfne but solver[5, 1/150, 1, 6] works well, so this isn't a coding issue, I'm afraid.
Mar 24, 2017 at 7:07 comment added zhk @xzczd K0 is a parameter, I kept it zero, {gm = 2, N1 = 3, S1 = 0, Sc = 6, K0 = 0} the new dependent variable is c[t,x,y]. With the conditions of the old problem it is working but if i change it then its not working.
Mar 24, 2017 at 7:05 comment added xzczd @maple If K0 is another independent variable, add it to the first argument of Compile; c isn't initialized, it need to be first defined in the same way as u, v, T, etc.; your PC doesn't have a C compiler installed, or your C compiler isn't successfully configured, if speed isn't that important, simply take away the CompilationTarget option, or search in this site to find the proper way to configure a C compiler.
Mar 24, 2017 at 6:50 comment added zhk @xzczd Please have a look, I don't know where I am making a mistake. dropbox.com/s/otzgxjczzl6untx/FDM.nb?dl=0
Mar 24, 2017 at 5:35 comment added xzczd @MapleSE-Area51Proposal Just modify {u = table, v = table, T = table} to {u = Table[1., {index@x}, {index@y}], v = Table[0., {index@x}, {index@y}], T = Table[1., {index@x}, {index@y}]}. (The code can be modified in a conciser way, but I think this modification is easier to understand. )
Mar 24, 2017 at 4:36 comment added zhk @xzczd How to change the boundary conditions? From u=T=0 at y=0 to u=T=1 at y=0 and v=0 remain the same.
Jan 3, 2017 at 13:11 history edited xzczd CC BY-SA 3.0
remove redundant HoldPattern; add another useful link.
Dec 24, 2016 at 13:51 history edited xzczd CC BY-SA 3.0
modify the naming a little
Dec 23, 2016 at 13:48 history edited xzczd CC BY-SA 3.0
deleted 4 characters in body
Dec 23, 2016 at 13:41 history edited xzczd CC BY-SA 3.0
FDM solution added
Dec 22, 2016 at 6:06 comment added zhk Let us continue this discussion in chat.
Dec 21, 2016 at 13:34 comment added xzczd @mmm Of course, because ContourPlot doesn't accept this syntax. Just compare the result of ContourPlot[Cos[x] Cos[y], {x, -3, 3}, {y, -3, 3}] and ContourPlot[Sin[x] Sin[y], {x, -3, 3}, {y, -3, 3}] and ContourPlot[{Sin[x] Sin[y], Cos[x] Cos[y]}, {x, -3, 3}, {y, -3, 3}]. For more details, check the document of ContourPlot carefully. As to the ListLinePlot issue, the step size is too large, try something like {n, 1, 5, 1/5}, you may also need to adjust the InterpolationOrder option.
Dec 21, 2016 at 13:30 comment added zhk I was looking for some sort of smooth curve but with this ListLinePlot[ Table[Abs@D[solfunc[5, 4, 1, 1, n][[1]][1, 1, y], y]^n /. y -> 0, {n, 1, 5, 1}]] I am not getting it.
Dec 21, 2016 at 13:22 comment added zhk I am using this ContourPlot[ sollst[#][[1]][1, x, y] & /@ {1, 2, 4} // Evaluate, {x, 0, 10}, {y, 0, 10}, PlotPoints -> 180, Frame -> True but the contours does not make any sense.
Dec 21, 2016 at 13:13 comment added xzczd @mmm Well, I think it's quite straightforward. What have you tried?
Dec 21, 2016 at 12:57 comment added zhk I still didn't get the contour plot. Any suggestions?
Dec 21, 2016 at 4:07 comment added xzczd @mmm Try a larger eps e.g. eps =10^2
Dec 21, 2016 at 3:56 comment added zhk Contour of u[5,x,y] vs (x,y) for different values of n.
Dec 20, 2016 at 18:23 history edited xzczd CC BY-SA 3.0
added 493 characters in body
Dec 20, 2016 at 14:35 comment added xzczd @J.M. You're right, though this circumvention triggers some warning, it's feasible. Check my edit for more details.
Dec 20, 2016 at 14:33 history edited xzczd CC BY-SA 3.0
added 1 character in body
Dec 20, 2016 at 14:27 history edited xzczd CC BY-SA 3.0
added 1 character in body
Dec 20, 2016 at 14:24 comment added xzczd @mmm 1) Check my edit. 2) Notice the output of solfunc is a list of solution, the 1st one is the solution for u, the 2nd one for v, the 3rd one for T, so e.g. solfunc[5, 4, 1, 1, n][[1]] is a solution for u. To obtain a list of (Abs[D[u[1, 1, y],y]])^(n) at y=0 vs n, you can e.g.: Table[Abs@D[solfunc[5, 4, 1, 1, n][[1]][1, 1, y], y]^n /. y -> 0, {n, 1, 2}], then you can use ListLinePlot to get the picture. 3) Contour of which variable?
Dec 20, 2016 at 14:11 history edited xzczd CC BY-SA 3.0
added 1114 characters in body
Dec 20, 2016 at 12:54 vote accept zhk
Dec 20, 2016 at 12:38 comment added zhk First of all, I really appreciate your efforts. I have couple questions. 1) When I try some other values of n=0.5, 1.5, 3 I am getting an error. 2) How to plot (Abs[D[u[1, 1, y],y]])^(n) at y=0 vs n? 3) How to plot contours of `u[t,x,y]' for various values of n?
Dec 20, 2016 at 7:19 history edited xzczd CC BY-SA 3.0
add another remark
Dec 18, 2016 at 8:57 history edited xzczd CC BY-SA 3.0
added 522 characters in body
Dec 18, 2016 at 8:42 comment added J. M.'s missing motivation Actually, why not replace Abs[f] with Sqrt[f^2] and sidestep the piecewise conversion altogether?
Dec 18, 2016 at 8:26 history answered xzczd CC BY-SA 3.0