// Run as mos-script using OpenModelica 1.9.1+dev (r18507) (RML version)

"Modelica.Blocks.Examples.InverseModel";
echo(false);
loadString("
function cell
  input Real time;
  input Boolean greenOnOk := false;
  input String greenColor := \"#00FF00\";
  output String str;
algorithm
  str := if time==-1.0 then \"<td bgcolor=\\\"#FF0000\\\">&nbsp;</td>\" else ((if greenOnOk then \"<td bgcolor=\\\"\" + greenColor + \"\\\">\" else \"<td>\") + 
    OpenModelica.Scripting.Internal.Time.readableTime(if time < 1e-5 then 0 else time)+\"</td>\");
end cell;
function length input String str; output Integer len; external \"C\" len=ModelicaStrings_length(str); end length;
function substring input String str; input Integer i1,i2; output String out; external \"C\" out=ModelicaStrings_substring(str,i1,i2); end substring;
function firstPart input String str; input Integer maxLength := 4096; output String out; algorithm out := if length(str) < maxLength then str else (substring(str,1,maxLength) + \"
... (output truncated, was \" + String(length(str)) + \" bytes)\"); end firstPart;
function csvFileToVariable
  input String str;
  output String out;
protected
  String matches[2];
algorithm
  (,matches) := OpenModelica.Scripting.regex(str,\"^.*[.]diff[.](.*)$\",2);
  out := matches[2];
end csvFileToVariable;
function testcase
  input Real time;
  input Boolean success;
  input String errXMLLong;
  output String str;
protected
  String errXML;
algorithm
  errXML := firstPart(errXMLLong,16*1024); // 16 kB log per test should be enough; increase if there is a good reason to do so
  str:=\"<testcase name=\\\"Modelica.Blocks.Examples.InverseModel\\\" time=\\\"\"+String(time)+\"\\\">\"+(
    if not success then \"<error type=\\\"Failure\\\" />\"
    else \"\"
  )+(
    if errXML <> \"\" then (\"<system-out>\" + errXML + \"</system-out>\")
    else \"\"
  )+
  \"</testcase>
\";
end testcase;
");getErrorString();
statFile := "Modelica.Blocks.Examples.InverseModel.stat";
writeFile("BuildModelRecursive.html","Modelica.Blocks.Examples.InverseModel
",append=true);
writeFile(statFile,"<tr><td bgcolor=\"#ff0000\">Modelica.Blocks.Examples.InverseModel</td></tr>");getErrorString();
writeFile(statFile + ".flat", testcase(0,false,"killed"));getErrorString();
writeFile(statFile + ".compile", testcase(0,false,"killed"));getErrorString();
writeFile(statFile + ".sim", testcase(0,false,"killed"));getErrorString();
writeFile(statFile + ".verify", testcase(0,false,"killed"));getErrorString();
outputFormat:="csv";

runScript("CustomCommands.mos");
loadModel(Modelica,{"default"});

reference_reltol:=0.003;
reference_reltolDiffMinMax:=0.003;
reference_rangeDelta:=0.001;
(startTime,stopTime,tolerance):=getSimulationOptions(Modelica.Blocks.Examples.InverseModel,defaultTolerance=1e-06);

referenceOK := false;
referenceFiles := "/home/martin/trunk/Examples/tmp2/refs";
referenceCell := if referenceFiles == "" then "" else "<td>&nbsp;</td>";
reference := "/home/martin/trunk/Examples/tmp2/refs/Modelica.Blocks.Examples.InverseModel.csv";
referenceExists := referenceFiles <> "" and regularFileExists(reference);
variableFilter := if referenceExists
  then sum((stringReplace(stringReplace(s,"[","."),"]",".") + "|") for s in readSimulationResultVars(reference))
  else ".*";

res:=buildModel(Modelica.Blocks.Examples.InverseModel,tolerance=tolerance,outputFormat=outputFormat,variableFilter=variableFilter);
greenColor := "#00FF00";
errFile:="files/Modelica.Blocks.Examples.InverseModel.err";
simFile:="files/Modelica.Blocks.Examples.InverseModel.sim";
err:=getErrorString();
system("rm -f " + errFile);
writeFile(simFile,"");
if err <> "" then
  writeFile(errFile,err);
end if;

frontend :=OpenModelica.Scripting.Internal.Time.timerTock(OpenModelica.Scripting.Internal.Time.RT_CLOCK_FRONTEND);
backend  :=OpenModelica.Scripting.Internal.Time.timerTock(OpenModelica.Scripting.Internal.Time.RT_CLOCK_BACKEND);
simcode  :=OpenModelica.Scripting.Internal.Time.timerTock(OpenModelica.Scripting.Internal.Time.RT_CLOCK_SIMCODE);
templates:=OpenModelica.Scripting.Internal.Time.timerTock(OpenModelica.Scripting.Internal.Time.RT_CLOCK_TEMPLATES);
total    :=OpenModelica.Scripting.Internal.Time.timerTock(OpenModelica.Scripting.Internal.Time.RT_CLOCK_SIMULATE_TOTAL);
build    :=OpenModelica.Scripting.Internal.Time.timerTock(OpenModelica.Scripting.Internal.Time.RT_CLOCK_BUILD_MODEL);
frontend :=if backend <> -1.0 then frontend-backend else frontend;
backend  :=if simcode <> -1.0 then backend-simcode else backend;
simcode  :=if templates <> -1.0 then simcode-templates else simcode;
templates:=if build <> -1.0 then templates-build else templates;
timeDiff := -1.0;

OpenModelica.Scripting.Internal.Time.timerTick(OpenModelica.Scripting.Internal.Time.RT_CLOCK_USER_RESERVED);
buildRes := res[1] <> "";
simRes  := if not buildRes then false else 0 == system("ulimit -t 480 ; ./Modelica.Blocks.Examples.InverseModel > "+simFile+" 2>&1");

system("sed -i '300,$ d' '" + simFile + "'"); // Only keep the top 300 lines
timeSim := OpenModelica.Scripting.Internal.Time.timerTock(OpenModelica.Scripting.Internal.Time.RT_CLOCK_USER_RESERVED);
if simRes then
  system("touch Modelica.Blocks.Examples.InverseModel.simsuccess");
  prefix := "files/Modelica.Blocks.Examples.InverseModel.diff";
  if referenceExists then
    OpenModelica.Scripting.Internal.Time.timerTick(OpenModelica.Scripting.Internal.Time.RT_CLOCK_USER_RESERVED);
    (referenceOK,diffVars) := diffSimulationResults("Modelica.Blocks.Examples.InverseModel_res."+outputFormat,reference,prefix,relTol=reference_reltol,relTolDiffMinMax=reference_reltolDiffMinMax,rangeDelta=reference_rangeDelta);
    if referenceOK then
      system("touch Modelica.Blocks.Examples.InverseModel.verifysuccess");
    end if;
    timeDiff := OpenModelica.Scripting.Internal.Time.timerTock(OpenModelica.Scripting.Internal.Time.RT_CLOCK_USER_RESERVED);
    diffFiles := {prefix + "." + var for var in diffVars};
    // Create a file containing only the calibrated variables, for easy display
    if not referenceOK then
      timeDiff := OpenModelica.Scripting.Internal.Time.timerTock(OpenModelica.Scripting.Internal.Time.RT_CLOCK_USER_RESERVED);
      referenceCell := "<td bgcolor=\"#FF0000\">"+OpenModelica.Scripting.Internal.Time.readableTime(timeDiff)+", <a href=\"files/Modelica.Blocks.Examples.InverseModel.diff.html\">"+String(size(diffFiles,1))+" failed</a></td>";
      writeFile("files/Modelica.Blocks.Examples.InverseModel.diff.html","<html><body><h1>Modelica.Blocks.Examples.InverseModel differences from the reference file</h1><p>startTime: "+String(startTime)+"</p><p>stopTime: "+String(stopTime)+"</p><p>Simulated using tolerance: "+String(tolerance)+"</p><ul>" + sum("<li>"+csvFileToVariable(file)+" <a href=\""+OpenModelica.Scripting.basename(file)+".html\">(javascript)</a> <a href=\""+OpenModelica.Scripting.basename(file)+".csv\">(csv)</a></li>" for file in diffFiles) + "</ul></body></html>");
      {writeFile(prefix + "." + var + ".html","<html>
<head>
<script type=\"text/javascript\" src=\"dygraph-combined.js\"></script>
    <style type=\"text/css\">
    #graphdiv {
      position: absolute;
      left: 10px;
      right: 10px;
      top: 40px;
      bottom: 10px;
    }
    </style>
</head>
<body>
<div id=\"graphdiv\"></div>
<p><input type=checkbox id=\"0\" checked onClick=\"change(this)\">
<label for=\"0\">reference</label>
<input type=checkbox id=\"1\" checked onClick=\"change(this)\">
<label for=\"1\">actual</label>
<input type=checkbox id=\"2\" checked onClick=\"change(this)\">
<label for=\"2\">high</label>
<input type=checkbox id=\"3\" checked onClick=\"change(this)\">
<label for=\"3\">low</label>
<input type=checkbox id=\"4\" checked onClick=\"change(this)\">
<label for=\"4\">error</label>
<input type=checkbox id=\"5\" onClick=\"change(this)\">
<label for=\"5\">actual (original)</label>
Parameters used for the comparison: Relative tolerance "+String(reference_reltol)+" (local), "+String(reference_reltolDiffMinMax)+" (relative to max-min). Range delta "+String(reference_rangeDelta)+".</p>
<script type=\"text/javascript\">
g = new Dygraph(document.getElementById(\"graphdiv\"),
                 \""+OpenModelica.Scripting.basename(prefix + "." + var+".csv")+"\",{title: '"+var+"',
  legend: 'always',
  connectSeparatedPoints: true,
  xlabel: ['time'],
  y2label: ['error'],
  series : { 'error': { axis: 'y2' } },
  colors: ['blue','red','teal','lightblue','orange','black'],
  visibility: [true,true,true,true,true,false]
});
function change(el) {
  g.setVisibility(parseInt(el.id), el.checked);
}
</script>
</body>
</html>") for var in diffVars};

    else
      referenceCell := "<td bgcolor=\"#00FF00\">"+OpenModelica.Scripting.Internal.Time.readableTime(timeDiff)+"</td>";
    end if;
  end if;
end if;

sim:=readFile(simFile);
str:="<tr><td>" + (if err <> "" then "<a href=\""+errFile+"\">Modelica.Blocks.Examples.InverseModel</a>" else "Modelica.Blocks.Examples.InverseModel") +
  (if sim <> "" then " (<a href=\""+simFile+"\">sim</a>)" else "") +
  "</td>" + referenceCell +
  cell(if simRes then timeSim else -1.0,greenOnOk=true,greenColor=greenColor) +
  "<td bgcolor=\""+(if res[1]<>"" then greenColor else "#FF0000")+"\">"+OpenModelica.Scripting.Internal.Time.readableTime(total)+"</td>" +
  sum(cell(d,greenColor=greenColor) for d in {frontend,backend,simcode,templates,build}) + "</tr>\n";getErrorString();
writeFile(statFile,str);getErrorString();
errXML:=escapeXML(err);
skip:="<testcase name=\"Modelica.Blocks.Examples.InverseModel\"><skipped /></testcase>";
str:=testcase(frontend,backend <> -1.0,errXML);
writeFile(statFile + ".flat", str);getErrorString();
str:=if backend <> -1.0 then testcase(total,buildRes,errXML) else skip;
writeFile(statFile + ".compile", str);getErrorString();
str:=if buildRes then testcase(timeSim,simRes,escapeXML(sim)) else skip;
writeFile(statFile + ".sim", str);getErrorString();
str:=if simRes then testcase(timeDiff,referenceOK,if referenceOK then "" else "See the generated report for details") else skip;
writeFile(statFile + ".verify", str);getErrorString();