Message Boards

WOLFRAM COMMUNITY

10309 Views

6 Replies

2 Total Likes

View groups...

Follow this post

Share this post:

GROUPS:

Engineering Mathematics Mechanical Engineering Mathematica External Programs and Systems Wolfram Language Optimization

interface to external program / command line and optimisation [beginner]

geotrupes a

Posted 9 years ago

Trying to: -1- define a function which takes a number [step 1], write to a file [step 2], runs an external bash script [step 3], read a text file and return the result as a number [step 4] -2- trying to optimise this function with NMinimize Here is the code of the function: f[x_] := CompoundExpression[SetDirectory["~/Desktop/test"]; Print[x]; (this is for debugging) Export["in.data", N[x], "Text"]; (this saves the value x in a text file) Run["~/Desktop/test/run.sh"]; (this runs a shell script which \ takes the text file, read it and makes the square of the number in in, a place holder for more to come) N[ToExpression[ Import["out.data"]]](this reads the number from the script \ output and convert to number) ] If I use: f[10], all if fine. the result is 100. However if I use: NMinimize[{f[y], 0 < y < 10}, {y}, Method -> "SimulatedAnnealing"] The results is nonsense as the function writes "y" in the first text file as symbolic form. Any idea on how to make this work? Thanks G.

Trying to: -1- define a function which takes a number [step 1], write to a file [step 2], runs an external bash script [step 3], read a text file and return the result as a number [step 4] -2- trying to optimise this function with NMinimize

Here is the code of the function:

f[x_] := CompoundExpression[SetDirectory["~/Desktop/test"];

  Print[x]; (*this is for debugging*)

  Export["in.data", N[x], 
   "Text"]; (*this saves the value x in a text file*)
  Run["~/Desktop/test/run.sh"]; (*this runs a shell script which \
takes the text file, 
  read it and makes the square of the number in in, 
  a place holder for more to come*)

  N[ToExpression[
    Import["out.data"]]](*this reads the number from the script \
output and convert to number*) ]

If I use: f[10], all if fine. the result is 100.

However if I use:

NMinimize[{f[y], 0 < y < 10}, {y}, Method -> "SimulatedAnnealing"]

The results is nonsense as the function writes "y" in the first text file as symbolic form.

Any idea on how to make this work? Thanks G.

POSTED BY: geotrupes a

6 Replies

Sort By:

Bill Simpson

Posted 9 years ago

Perhaps this

POSTED BY: Bill Simpson

Bill Simpson

Posted 9 years ago

You think you are seeing this behaving as you expect In[1]:= f[x_?NumericQ] := RunProcess[{"echo", ToString[InputForm[x^2, NumberMarks->False]]},"StandardOutput"]; f[3] but if you then do this Print[FullForm[%]] I believe you will understand that your function is returning a string, not a numeric value, and trying to "minimize a string" isn't going to work. If you are running Windows then f[x_?NumericQ] := ToExpression[StringDrop[RunProcess[{"echo", ToString[InputForm[x^2, NumberMarks->False]]}, "StandardOutput"], -2]]; Print[FullForm[f[3]]] might be closer to what you want. And if you aren't running Windows then the StringDrop may need to be adjusted to match your local conventions. I think there should be a much better way of doing what you are trying to do.

POSTED BY: Bill Simpson

geotrupes a

Posted 9 years ago

Thank you Bill, tried and it works f[x_?NumericQ] := CompoundExpression[ SetDirectory["~/Desktop/test"]; Print[N[x]]; Export["in.data", N[x], "Text"]; Run["~/Desktop/test/run.sh"]; result = N[ToExpression[Import["out.data"]]] ; Run["rm *.data"]; Print[result] There are few issues to resolve but this is great. THANK YOU! G. PS do you happen to know if this is running in a multicore fashion?

Thank you Bill, tried and it works

f[x_?NumericQ] := CompoundExpression[
  SetDirectory["~/Desktop/test"];
  Print[N[x]]; 
  Export["in.data", N[x], "Text"];
  Run["~/Desktop/test/run.sh"]; 
  result = N[ToExpression[Import["out.data"]]] ;
  Run["rm *.data"];
  Print[result]

There are few issues to resolve but this is great.

THANK YOU! G.

PS do you happen to know if this is running in a multicore fashion?

POSTED BY: geotrupes a

geotrupes a

Posted 9 years ago

POSTED BY: geotrupes a

geotrupes a

Posted 9 years ago

Ok, I tried a simpler function. Using echo to print a square, made in Mathematica. See this: f[x_?NumericQ] := RunProcess[{"echo", ToString[InputForm[x^2, NumberMarks -> False]]}, "StandardOutput"]; It works for a single value of x, but then NMinimize complains that: NMinimize[{f[y], 1 < y < 10}, {y}, Method -> "RandomSearch"] NMinimize::nnum: The function value 1.9925755057065486` is not a number at {y} = {1.41159}. >> What is the issue? I must confess to be rather lost at this stage :-( Thanks G>

Ok, I tried a simpler function. Using echo to print a square, made in Mathematica.

See this:

f[x_?NumericQ] := 
  RunProcess[{"echo", ToString[InputForm[x^2, NumberMarks -> False]]}, 
   "StandardOutput"];

It works for a single value of x, but then NMinimize complains that:

NMinimize[{f[y], 1 < y < 10}, {y}, Method -> "RandomSearch"]


NMinimize::nnum: The function value 1.9925755057065486`
 is not a number at {y} = {1.41159}. >>

What is the issue? I must confess to be rather lost at this stage :-( Thanks G>

POSTED BY: geotrupes a

geotrupes a

Posted 9 years ago

The ToExpression solves the issue and all works!! I am very new to Mathematica an trying to run without knowing how to walk :-) If there is a better way and you have some idea please shoot! Thank you, G>

POSTED BY: geotrupes a

Reply to this discussion

Reply Preview

Attachments

Remove Add a file to this post

Follow this discussion

or Discard

Group Abstract

Feedback