Skip to main content

mathlink or wstp - How to implement LinkLaunch-like functionality with timeout?


tl;dr: What is the best way (and correct sequence of operations) to connect two programs using MathLink while anticipating failures / timeouts?



Launch program and connect, with timeout


I have a program P that I want to connect to a Mathematica kernel K through MathLink/WSTP. P may take a long time to launch, and it may also fail completely (depending on the computer where it is being run).



The expected behaviour is this:



  • K starts up P, and waits for a connection.

  • If P does not connect back within time t, K should fail gracefully and report an error.

  • If P took longer to start than K expected, and at that point it can no longer connect, it should also fail gracefully.

  • They should not leave any broken links behind if the connection process fails. Preferably, no shared memory link names should become unusable in case of failure.


K will use Wolfram Language functions. P will use the C MathLink API.


I believe I should rely on LinkReadyQ for the timeout (see its second argument), but it is not clear what is the correct sequence of operations to perform if robustness is critical.


It is also not clear to me which of K and P should use the Connect and Listen link modes (I think K should use Listen and P should use Connect).




Connect two already running programs


I would also like to implement a second mode of operation for K and P:


Suppose that both K and P are already running. I would like to implement a user-friendly method to connect them together. Ideally, the user would issue a command both in K and P:


connectToP[] (* Issued in K. *)

connectToK[] (* Issued in P.  P is not Mathematica, but for this example, it's okay to pretend it is. *)

Preferably, the order in which these commands are issued should not matter, and the connection would be set up automatically. This is for the sake of user-friendliness. Again, it is preferable to make the commands abortable without bad side effects (or to allow them to time out).


Use case: Connect Mathematica with another scientific computing system. Suppose that you have already done some work in the other system and have a workspace full of variables. You realize that you want to transfer some of these to Mathematica without having to restart the system.



The key question here is if it is possible to implement a robust connection method which does not care about the order in which the commands are issued. Alternatively, there could be a user-friendly error message (and no bad side effects, such as broken links) when the commands are issued in the wrong order.



This is related to Why does LinkRead block when LinkReadyQ returns True?.




Comments

Post a Comment

Popular posts from this blog

plotting - Plot 4D data with color as 4th dimension

I have a list of 4D data (x position, y position, amplitude, wavelength). I want to plot x, y, and amplitude on a 3D plot and have the color of the points correspond to the wavelength. I have seen many examples using functions to define color but my wavelength cannot be expressed by an analytic function. Is there a simple way to do this? Answer Here a another possible way to visualize 4D data: data = Flatten[Table[{x, y, x^2 + y^2, Sin[x - y]}, {x, -Pi, Pi,Pi/10}, {y,-Pi,Pi, Pi/10}], 1]; You can use the function Point along with VertexColors . Now the points are places using the first three elements and the color is determined by the fourth. In this case I used Hue, but you can use whatever you prefer. Graphics3D[ Point[data[[All, 1 ;; 3]], VertexColors -> Hue /@ data[[All, 4]]], Axes -> True, BoxRatios -> {1, 1, 1/GoldenRatio}]
Post a Comment
Read more

plotting - Filling between two spheres in SphericalPlot3D

Manipulate[ SphericalPlot3D[{1, 2 - n}, {θ, 0, Pi}, {ϕ, 0, 1.5 Pi}, Mesh -> None, PlotPoints -> 15, PlotRange -> {-2.2, 2.2}], {n, 0, 1}] I cant' seem to be able to make a filling between two spheres. I've already tried the obvious Filling -> {1 -> {2}} but Mathematica doesn't seem to like that option. Is there any easy way around this or ... Answer There is no built-in filling in SphericalPlot3D . One option is to use ParametricPlot3D to draw the surfaces between the two shells: Manipulate[ Show[SphericalPlot3D[{1, 2 - n}, {θ, 0, Pi}, {ϕ, 0, 1.5 Pi}, PlotPoints -> 15, PlotRange -> {-2.2, 2.2}], ParametricPlot3D[{ r {Sin[t] Cos[1.5 Pi], Sin[t] Sin[1.5 Pi], Cos[t]}, r {Sin[t] Cos[0 Pi], Sin[t] Sin[0 Pi], Cos[t]}}, {r, 1, 2 - n}, {t, 0, Pi}, PlotStyle -> Yellow, Mesh -> {2, 15}]], {n, 0, 1}]
Post a Comment
Read more

plotting - Mathematica: 3D plot based on combined 2D graphs

I have several sigmoidal fits to 3 different datasets, with mean fit predictions plus the 95% confidence limits (not symmetrical around the mean) and the actual data. I would now like to show these different 2D plots projected in 3D as in but then using proper perspective. In the link here they give some solutions to combine the plots using isometric perspective, but I would like to use proper 3 point perspective. Any thoughts? Also any way to show the mean points per time point for each series plus or minus the standard error on the mean would be cool too, either using points+vertical bars, or using spheres plus tubes. Below are some test data and the fit function I am using. Note that I am working on a logit(proportion) scale and that the final vertical scale is Log10(percentage). (* some test data *) data = Table[Null, {i, 4}]; data[[1]] = {{1, -5.8}, {2, -5.4}, {3, -0.8}, {4, -0.2}, {5, 4.6}, {1, -6.4}, {2, -5.6}, {3, -0.7}, {4, 0.04}, {5, 1.0}, {1, -6.8}, {2, -4.7}, {3, -1....
Post a Comment
Read more