Message Boards Message Boards

0
|
16835 Views
|
15 Replies
|
5 Total Likes
View groups...
Share
Share this post:

Can somebody please tell me why the Wolfram example code will not compile?

Paul T
Posted 11 years ago
Hi,

I am attempting to compile the example code at Create Standalone Executables Using Compiled Functions, but I'm getting errors suggesting, among other things that:
  • the example code is not correct,
  • that there are unknown variable types,
  • and that 'something' is missing before the unknown variable type.
I am fairly conversant with Mathematica, but C is new to me. I've just been attempting to learn it for the past few weeks or so. I've tried various compilers and their consoles (or lack thereof), and IDEs (or, again, the lack thereof) such as MinGW, Pelles C, and Cygwin.

My main compiler, and the only one the CCompiler function seems to recognize, is the MinGW one. This is what the CCompiler[ ] function gives me:
In[10]:= CCompilers[]

Out[10]= {{"Name" -> "MinGW",
  "Compiler" -> CCompilerDriver`MinGWCompiler`MinGWCompiler,
  "CompilerInstallation" -> "C:\\MinGW\\bin\\gcc.exe",
  "CompilerName" -> Automatic}}
These are my specific errors and the command(s) that gave them:
 In[8]:= Needs["CCompilerDriver`"];
 
 lopassExe =
  CreateExecutable[{fnSource, lopassmainSrcFile}, "lowpass",
   "TargetDirectory" -> targetDir,
   "Libraries" -> "WolframRTL_Static_Minimal"]
 
 During evaluation of In[8]:= CreateExecutable::cmperr: Compile error: C:\Program Files\Wolfram Research\Mathematica\9.0\SystemFiles\IncludeFiles\C/WolframLibrary.h:31:15: error: expected '=', ',', ';', 'asm' or '__attribute__' before 'int' >>
 
During evaluation of In[8]:= CreateExecutable::cmperr: Compile error: C:\Program Files\Wolfram Research\Mathematica\9.0\SystemFiles\IncludeFiles\C/WolframLibrary.h:60:2: error: unknown type name 'mint' >>

During evaluation of In[8]:= CreateExecutable::cmperr: Compile error: C:\Program Files\Wolfram Research\Mathematica\9.0\SystemFiles\IncludeFiles\C/WolframLibrary.h:141:2: error: expected specifier-qualifier-list before 'mint' >>

During evaluation of In[8]:= General::stop: Further output of CreateExecutable::cmperr will be suppressed during this calculation. >>

Out[9]= $Failed
As you can see, this is simply part of the example code on the Wolfram page along with my errors.

Can anybody come up with any simple reasons (or at least explain the ones listed above in a bit more depth) the example code won't compile, and what I might do to resolve the problem?

Thank you.

P.S. Since this is my first post, I'm somewhat at sea as to where I should post. Also, it would seem to me that somebody else ought to have had the problem since it's just at the copy and paste level, but even though I did searches, I haven't seen anybody else with the problem, nor does this seem really relevant to the forum group classification system, so please forgive me if this is out of place.
POSTED BY: Paul T
15 Replies
Sort By:
Paul T
Paul T
Posted 11 years ago
Well now this seems to work:
 In[1]:= << CCodeGenerator`
 
 In[7]:= a = Compile[{{x, _Integer}}, x^2]
 
 Out[7]= \!\(
 TagBox[
 RowBox[{"CompiledFunction", "[",
 RowBox[{
 RowBox[{"{", "x", "}"}], ",",
SuperscriptBox["x", "2"], ",", "\<\"-CompiledCode-\"\>"}], "]"}],
False,
Editable->False]\)

In[9]:= b = SymbolicCGenerate[a, "fun"]

Out[9]= CProgram[{CInclude["math.h"]}, {CInclude[
   "WolframRTL.h"]}, {CDeclare[{"static",
    "WolframCompileLibrary_Functions"}, "funStructCompile"],
  CDeclare[{"static", "mbool"}, CAssign["initialize", 1]]},
CInclude["fun.h"],
CFunction[{"DLLEXPORT", "int"},
  "Initialize_fun", {{"WolframLibraryData", "libData"}},
  CBlock[{CIf["initialize",
     CBlock[{{CAssign["funStructCompile",
         CPointerMember["libData",
          "compileLibraryFunctions"]]}, {}, {}, {},
       CAssign["initialize", 0]}]], CReturn[0]}]],
CFunction[{"DLLEXPORT", "void"},
  "Uninitialize_fun", {{"WolframLibraryData", "libData"}},
  CBlock[CIf[COperator[Not, {"initialize"}],
    CBlock[{}, CAssign["initialize", 1]]]]], {},
CFunction[{"DLLEXPORT", "int"},
  "fun", {{"WolframLibraryData", "libData"}, {"mint", "A1"}, {"mint",
    CDereference["Res"]}},
  CBlock[{{CDeclare["mint", "I0_0"],
     CDeclare["mint", "I0_1"]}, {}, {CAssign["I0_0",
      "A1"]}, {{CAssign["I0_1",
       COperator[Times, {"I0_0", "I0_0"}]], {}}, {}}, {}, {{},
     CAssign[CDereference["Res"], "I0_1"], {}, {},
     CCall[CPointerMember["funStructCompile",
       "WolframLibraryData_cleanUp"], {"libData", 1}],
     CReturn[0]}}]], {}]

In[10]:= b // ToCCodeString

Out[10]= "#include \"math.h\"

#include \"WolframRTL.h\"

static WolframCompileLibrary_Functions funStructCompile;

static mbool initialize = 1;

#include \"fun.h\"

DLLEXPORT int Initialize_fun(WolframLibraryData libData)
{
if( initialize)
{
funStructCompile = libData->compileLibraryFunctions;
initialize = 0;
}
return 0;
}

DLLEXPORT void Uninitialize_fun(WolframLibraryData libData)
{
if( !initialize)
{
initialize = 1;
}
}

DLLEXPORT int fun(WolframLibraryData libData, mint A1, mint *Res)
{
mint I0_0;
mint I0_1;
I0_0 = A1;
I0_1 = I0_0 * I0_0;
*Res = I0_1;
funStructCompile->WolframLibraryData_cleanUp(libData, 1);
return 0;
}

"
I'm guessing it's not my Mathematica 9 then, but something about MinGW. What could it be?
POSTED BY: Paul T
The type "mint" is defined by Mathematica in a .h file -- I thought it was mathlink.h, but it might be WolframLibrary.h.  There are a few #ifdef's I don't want to track down.  Here's a guide to C development:

http://reference.wolfram.com/mathematica/tutorial/MathLinkDeveloperGuide-Windows.html

You should look at your WolframLibrary.h and see if it looks ok.  Maybe it got trashed.  Mine looks fine (and I get no errors).
POSTED BY: Michael Rogers
Paul T
Paul T
Posted 11 years ago
Thank you for that. I will study these things, although it will take me some time since I'm such a noob to making this part of it work.

I have one of the sales reps words for it to me that even if a month elapses from the date of my purchase (12/31/2013), if Mathematica 9 still isn't performing for me, then I can get consideration for extended time to make it work before crunch time comes and I have to either keep it or ask for a refund. I don't want to give up on Mathematica 9. I want to make it work the way it's supposed to.

Thank you again for your help.

Oh, wait, I have a couple more questions. One is about mint: why would Mathematica 9 fail to recognize its own data type? What's up with that? I didn't think of that just now either, it's been in the back of my mind all along.

The other question is about what I tried to do earlier in this thread when I attempted to force Mathematica to acknowledge the path to some other compiler than MinGW. It's as if I went fishing and not only baited the hook, but grabbed the fish out of the water, shoved the barbed hook through the side of its mouth, and tossed it back into the drink in the hope that I would get to see the rod, reel, and fishline haul it out of the water, only to find that it had taken my bait and then slipped the hook. How and why does Mathematica simply ignore all of my attempts to force it to acknowledge the presence of other compilers? While it agrees that there is a new default compiler if I say there is, it then simply ignores everything else. That just seems weird to me, but then I don't know what mistakes I might be making in my attempts to do that either . . .

I think Wolfram needs to be a whole lot more explicit by way of example in these matters. So far, the examples are proving to be much too scanty.
POSTED BY: Paul T
Paul T
Paul T
Posted 11 years ago
I've just barely gotten started poking around, but I've already found something a bit odd in one of the files.

I can't say that I know what looks good or what looks bad. To me, everything seemed OK, but what do I know?

I avoided poking around in most of the Wolfram files so as to avoid making associations about how the files ought to be read. I know you can remove unwanted associations, but since my system is Windows XP, if you want to completely remove an unwanted association rather than change them from one to another, you have to do a command prompt ASSOC filetype =  and leave nothing after the = sign. It isn't very hard and I've done it in the past, but the less fooling around the better, so I was cautious about attempting to read files.

Having said that, I did go along this path: C:\Program Files\Wolfram Research\Mathematica\9.0\SystemFiles\IncludeFiles\C and I looked at all of the stuff in the C folder. I felt OK about doing so because those can all be opened in notepad without making any undue associations.

There were indeed, as you said Michael, a bunch of #ifdefs, typedefs, etc., and quite a bit about mint which I didn't fully grasp. This is what struck me as a bit odd to be leaving in a distribution to the general public though:

          /*
          TODO these two definitions look strange.
          */

That was in the WolframCompileLibrary near the bottom of the #ifdef. That seems like an oversight on the programmer's part. I know it's just a comment, but it was the content of the comment that struck me. I would have thought they'd have gotten rid of stuff like that before they sold it to the general public. While I doubt it has anything to do with my particular problem, it does seem a bit odd to me.

    
POSTED BY: Paul T
Paul T
Paul T
Posted 11 years ago
I found the answer: I needed a shorter directory name. I suspected that earlier, but when I got rid of the CreateDirectory thing and forced a short name like "gloop" on the tag end of Temp, it still didn't run. I thought sure I was short enough with the "gloop" thing, but I guess not. Finally, when I used "C:/Temp" for targetDir though, it worked.

Ultimately, shortening file names has been the answer to some other things I wanted to do recently too, most notably when I was attempting to make some CUDA stuff work. I tried all of the stuff in this thread and none of it worked, up until I tried what was suggested in post #9.

I tend to prefer long descriptive inclusive directory and file names, but I've found in the past that my preference for long names gets me into trouble Then I end up having to do a lot of shortening and renaming in order to be able to open even simple text and picture files.

In conclusion, I'd say there is a subtle flaw in the example as given, in that the temporary file names Mathematica 9 likes to create are too long for some of the older systems.
POSTED BY: Paul T
Paul T
Paul T
Posted 11 years ago
Not to spam or belabor the issue, but my own much simpler example gives something similar:
 In[1]:= << CCodeGenerator`
 
 In[6]:= a = Compile[{{x, _Integer}}, x^2, CompilationTarget -> "C"]
 
 During evaluation of In[6]:= CreateLibrary::cmperr: Compile error: C:\Program Files\Wolfram Research\Mathematica\9.0\SystemFiles\IncludeFiles\C/WolframLibrary.h:31:15: error: expected '=', ',', ';', 'asm' or '__attribute__' before 'int' >>
 
 During evaluation of In[6]:= CreateLibrary::cmperr: Compile error: C:\Program Files\Wolfram Research\Mathematica\9.0\SystemFiles\IncludeFiles\C/WolframLibrary.h:60:2: error: unknown type name 'mint' >>
 
 During evaluation of In[6]:= CreateLibrary::cmperr: Compile error: C:\Program Files\Wolfram Research\Mathematica\9.0\SystemFiles\IncludeFiles\C/WolframLibrary.h:141:2: error: expected specifier-qualifier-list before 'mint' >>

During evaluation of In[6]:= General::stop: Further output of CreateLibrary::cmperr will be suppressed during this calculation. >>

During evaluation of In[6]:= Compile::nogen: A library could not be generated from the compiled function. >>

Out[6]= \!\(
TagBox[
RowBox[{"CompiledFunction", "[",
RowBox[{
RowBox[{"{", "x", "}"}], ",",
SuperscriptBox["x", "2"], ",", "\<\"-CompiledCode-\"\>"}], "]"}],
False,
Editable->False]\)
I wonder, is something wrong with my Mathematica 9? I just bought it not too long ago.

Or perhaps something is wrong with MinGW? MinGW was mentioned as a decent free compiler to use. Perhaps it's not.
POSTED BY: Paul T
Paul T
Paul T
Posted 11 years ago
Just for yuks, here's another fruitless attempt to change things a bit:
 In[1]:= << CCompilerDriver`
 
 In[2]:= CCompilers[]
 
 Out[2]= {{"Name" -> "MinGW",
   "Compiler" -> CCompilerDriver`MinGWCompiler`MinGWCompiler,
   "CompilerInstallation" -> "C:\\MinGW\\bin\\gcc.exe",
   "CompilerName" -> Automatic}}
 
In[3]:= $CCompiler = {"Compiler" ->
   CCompilerDriver`CygwinGCC`CygwinGCC,
  "CompileInstallation" -> "C:\\cygwin\\bin\\gcc"}

Out[3]= {"Compiler" -> CCompilerDriver`CygwinGCC`CygwinGCC,
"CompileInstallation" -> "C:\\cygwin\\bin\\gcc"}

In[4]:= DefaultCCompiler[]

Out[4]= CCompilerDriver`CygwinGCC`CygwinGCC

In[5]:= CCompilers[Full]

Out[5]= {{"Name" -> "Visual Studio",
  "Compiler" ->
   CCompilerDriver`VisualStudioCompiler`VisualStudioCompiler,
  "CompilerInstallation" -> None,
  "CompilerName" -> Automatic}, {"Name" -> "MinGW",
  "Compiler" -> CCompilerDriver`MinGWCompiler`MinGWCompiler,
  "CompilerInstallation" -> "C:\\MinGW\\bin\\gcc.exe",
  "CompilerName" -> Automatic}, {"Name" -> "Cygwin GCC",
  "Compiler" -> CCompilerDriver`CygwinGCC`CygwinGCC,
  "CompilerInstallation" -> None,
  "CompilerName" -> Automatic}, {"Name" -> "Intel Compiler",
  "Compiler" -> CCompilerDriver`IntelCompiler`IntelCompiler,
  "CompilerInstallation" -> None,
  "CompilerName" -> Automatic}, {"Name" -> "Generic C Compiler",
  "Compiler" -> CCompilerDriver`GenericCCompiler`GenericCCompiler,
  "CompilerInstallation" -> None, "CompilerName" -> Automatic}}
Notice how I've attempted to change the default compiler from MinGW to CygwinGCC. Of course it didn't work at all. I suppose that what I'm just doing is just utter nonsense, although I have absolutely no idea of why it should be utter nonsense.

It seems to me like there ought to be some way to tell Mathematica to look elsewhere for other compilers, but there doesn't seem to be.

I've changed the environment variable path too, not that it did any good.
POSTED BY: Paul T
Your link is to a V8 reference.  The corresponding V9 page has another option set, "SystemLibraries" -> {}:
lopassExe =
CreateExecutable[{fnSource, lopassmainSrcFile}, "lowpass",
  "SystemLibraries" -> {}, "TargetDirectory" -> targetDir,
  "Libraries" -> "WolframRTL_Static_Minimal"]
I got different error messages (on a Mac) with your code than you did, but using the V9 example got rid of them.  Maybe it will work on your system.
POSTED BY: Michael Rogers
Paul T
Paul T
Posted 11 years ago
Thank you Michael, I will try that a little later on today.
POSTED BY: Paul T
Paul T
Paul T
Posted 11 years ago
Well, nice thought, nice try, but still no cigar:
  In[1]:= targetDir = CreateDirectory[]
  
  Out[1]= "C:\\Documents and Settings\\HP_Owner\\Local \
  Settings\\Temp\\m-8973967a-d7cd-4736-b0c9-a0939b2474e5"
  
  In[2]:= lopass =
    Compile[{{x, _Real, 1}, dt, RC},
     Module[{a = dt/(RC + dt), yprev = First[x], yi},
      Table[yi = a*x[[i]] + (1 - a)*yprev;
      yprev = yi;
      yi, {i, 1, Length[x]}]]];
 
 In[3]:= (*build the example in a specified directory*)
 fnSource = FileNameJoin[{targetDir, "lopass.c"}];
 Export[fnSource, lopass];
 FileNames["*", targetDir]
 
 Out[5]= {"C:\\Documents and Settings\\HP_Owner\\Local \
 Settings\\Temp\\m-8973967a-d7cd-4736-b0c9-a0939b2474e5\\lopass.c", \
 "C:\\Documents and Settings\\HP_Owner\\Local \
 Settings\\Temp\\m-8973967a-d7cd-4736-b0c9-a0939b2474e5\\lopass.h"}
 
 In[6]:= (*main function C source*)lopassmainSrc = "
     #include \"stdio.h\"
     #include \"stdlib.h\"
     #include \"lopass.h\"
     #include \"WolframRTL.h\"
    
     static WolframLibraryData libData = 0;
    
     int main()
     {
       int err = 0;
       mint i, type, rank, nelems, *dims;
       double *data;
       MTensor x, y;
       double dt;
       double RC;
       libData = WolframLibraryData_new(WolframLibraryVersion);
    
       /* read x */
       type = MType_Real;
       rank = 1;
       dims = (mint*)malloc(rank * sizeof(mint));
       scanf(\" %d\", &nelems);
       dims[0] = nelems;
    
       err = (*(libData->MTensor_new))(type, rank, dims, &x);
       if (err) return 1;
    
       free(dims);
    
       data = (*(libData->MTensor_getRealData))(x);
       for(i = 0; i < nelems; i++) {
         scanf(\" %lf\", &(data[i]));
       }
    
       /* read dt */
       scanf(\" %lf\", &dt);
       /* read RC */
       scanf(\" %lf\", &RC);
    
       err = Initialize_lopass(libData);
    
       y = 0;
       err = lopass(libData, x, dt, RC, &y);
       printf(\"%d\\n\", err);
       if(0 == err){
         dims = (mint*)libData->MTensor_getDimensions(y);
         nelems = dims[0];
         data = (*(libData->MTensor_getRealData))(y);
         printf(\"%d\\n\", nelems);
         for(i = 0; i < nelems; i++)
           printf(\"%f\\n\", data[i]);
       }
    
       Uninitialize_lopass(libData);
       return 0;
     }
     ";
 
 In[7]:= (*create the main source file*)lopassmainSrcFile =
  FileNameJoin[{targetDir, "lopassMain.c"}]
 Export[lopassmainSrcFile, lopassmainSrc, "Text"];
 
 Out[7]= "C:\\Documents and Settings\\HP_Owner\\Local \
 Settings\\Temp\\m-8973967a-d7cd-4736-b0c9-a0939b2474e5\\lopassMain.c"
 
 In[9]:= Needs["CCompilerDriver`"];
 lopassExe =
  CreateExecutable[{fnSource, lopassmainSrcFile}, "lowpass",
   "SystemLibraries" -> {}, "TargetDirectory" -> targetDir,
   "Libraries" -> "WolframRTL_Static_Minimal"]
 
 During evaluation of In[9]:= CreateExecutable::cmperr: Compile error: C:\Program Files\Wolfram Research\Mathematica\9.0\SystemFiles\IncludeFiles\C/WolframLibrary.h:31:15: error: expected '=', ',', ';', 'asm' or '__attribute__' before 'int' >>
 
 During evaluation of In[9]:= CreateExecutable::cmperr: Compile error: C:\Program Files\Wolfram Research\Mathematica\9.0\SystemFiles\IncludeFiles\C/WolframLibrary.h:60:2: error: unknown type name 'mint' >>
 
 During evaluation of In[9]:= CreateExecutable::cmperr: Compile error: C:\Program Files\Wolfram Research\Mathematica\9.0\SystemFiles\IncludeFiles\C/WolframLibrary.h:141:2: error: expected specifier-qualifier-list before 'mint' >>

During evaluation of In[9]:= General::stop: Further output of CreateExecutable::cmperr will be suppressed during this calculation. >>

Out[10]= $Failed
What strikes me as most odd, not that this is the only thing wrong with it, but what strikes me as most odd is that an unknown data type named mint keeps cropping up. I'm still a novice C programmer, but I've now been exposed to most of the data types a number of times, and I don't ever recall a type known as mint. I've seen int, long, double, float, long int, short, char, etc., but never a mint. What in the heck is a mint, I mean other than a candy or tea flavor or a big pile of money???

Here's Wikipedia on the subject of C data types: C data types. Where is there a mint in that list?
POSTED BY: Paul T
Paul T
Paul T
Posted 11 years ago
OK, I figured out what the problem was. I'm not exactly sure what to call it, but Micheal Rogers was right when he gave me the link to the Mathlink tutorial. I don't know exactly what to call it, configuring, set-up, or whatever, but I needed to link Mathematica 9 with one of my C compilers. Today I got the freely downloadable Microsoft C++ 2008 Version 9.0 of  Visual Studio. I spent the afternoon familiarizing myself with the IDE (I'm still a little unsure of how to make it display what I want, but that's sort of beside the point here).

Anyway, in the Mathlink tutorial it said to get the bin, lib and include files from the folder at the end of a certain Mathematica 9 path: C:\Program Files\Wolfram Research\Mathematica\9.0\SystemFiles\Links\MathLink\DeveloperKit\Windows\CompilerAdditions\MLDev32, and then hold down the Control button while doing a drag and drop of those files into the same files of Microsoft Visual Studio 2008. I did that and after a while I came back to this and ran the example again, and this time it worked!

Here's proof:
 In[1]:= targetDir = CreateDirectory[]
 
 Out[1]= "C:\\Documents and Settings\\HP_Owner\\Local \
 Settings\\Temp\\m-6535021e-b3c9-4783-b3a4-6b6d2714640f"
 
 In[2]:= lopass =
   Compile[{{x, _Real, 1}, dt, RC},
    Module[{a = dt/(RC + dt), yprev = First[x], yi},
     Table[yi = a*x[[i]] + (1 - a)*yprev;
     yprev = yi;
     yi, {i, 1, Length[x]}]]];

In[3]:= (*build the example in a specified directory*)
fnSource = FileNameJoin[{targetDir, "lopass.c"}];
Export[fnSource, lopass];
FileNames["*", targetDir]

Out[5]= {"C:\\Documents and Settings\\HP_Owner\\Local \
Settings\\Temp\\m-6535021e-b3c9-4783-b3a4-6b6d2714640f\\lopass.c", \
"C:\\Documents and Settings\\HP_Owner\\Local \
Settings\\Temp\\m-6535021e-b3c9-4783-b3a4-6b6d2714640f\\lopass.h"}

In[6]:= (*main function C source*)lopassmainSrc = "
    #include \"stdio.h\"
    #include \"stdlib.h\"
    #include \"lopass.h\"
    #include \"WolframRTL.h\"
   
    static WolframLibraryData libData = 0;
   
    int main()
    {
      int err = 0;
      mint i, type, rank, nelems, *dims;
      double *data;
      MTensor x, y;
      double dt;
      double RC;
      libData = WolframLibraryData_new(WolframLibraryVersion);
   
      /* read x */
      type = MType_Real;
      rank = 1;
      dims = (mint*)malloc(rank * sizeof(mint));
      scanf(\" %d\", &nelems);
      dims[0] = nelems;
   
      err = (*(libData->MTensor_new))(type, rank, dims, &x);
      if (err) return 1;
   
      free(dims);
   
      data = (*(libData->MTensor_getRealData))(x);
      for(i = 0; i < nelems; i++) {
        scanf(\" %lf\", &(data[i]));
      }
   
      /* read dt */
      scanf(\" %lf\", &dt);
      /* read RC */
      scanf(\" %lf\", &RC);
   
      err = Initialize_lopass(libData);
   
      y = 0;
      err = lopass(libData, x, dt, RC, &y);
      printf(\"%d\\n\", err);
      if(0 == err){
        dims = (mint*)libData->MTensor_getDimensions(y);
        nelems = dims[0];
        data = (*(libData->MTensor_getRealData))(y);
        printf(\"%d\\n\", nelems);
        for(i = 0; i < nelems; i++)
          printf(\"%f\\n\", data[i]);
      }
   
      Uninitialize_lopass(libData);
      return 0;
    }
    ";

In[7]:= (*create the main source file*)lopassmainSrcFile =
FileNameJoin[{targetDir, "lopassMain.c"}]
Export[lopassmainSrcFile, lopassmainSrc, "Text"];

Out[7]= "C:\\Documents and Settings\\HP_Owner\\Local \
Settings\\Temp\\m-6535021e-b3c9-4783-b3a4-6b6d2714640f\\lopassMain.c"

In[9]:= Needs["CCompilerDriver`"];
lopassExe =
CreateExecutable[{fnSource, lopassmainSrcFile}, "lowpass",
  "SystemLibraries" -> {}, "TargetDirectory" -> targetDir,
  "Libraries" -> "WolframRTL_Static_Minimal"]


Out[10]= "C:\\Documents and Settings\\HP_Owner\\Local \
Settings\\Temp\\m-6535021e-b3c9-4783-b3a4-6b6d2714640f\\lowpass.exe"
See, no failure, no unknown variable type named mint, just an .exe file!

Wow! After all that messing around and a couple of weeks of effort, it all came down to three drag-and-drops.

The instructions for the 2010 version they spoke of in the Mathlink tutorial were suitable for the 2008 version. You just find the relevant bin, lib and include files in the folder at the end of this path: C:\Program Files\Microsoft Visual Studio 9.0\VC -- the VC folder and then do the Control drag-and-drop.

Also, while I haven't attempted this yet, I noticed there was a special cygwin folder in the same folder as the mldev32 folder. I'm going to attempt the same sort of thing on cygwin, MinGW, and Pelles C compilers and see if I can get it to work with those. If I can, then I'll know that my problem was simply that I had to do this manually, and that the most inexpensive version of Mathematica 9, the Home Edition, may need a little help from the user before it will work properly.

Oh, and a few days ago, I realized there was another language I might potentially try these things with, the R language. I know nothing about it, and I'd probably have to compose my own examples, but I wouldn't be a bit surprised if something similar to this same procedure is probably the way to link up Mathematica with all of those other languages in Mathlink.

OK then, for now, solved, thanks to Michael Rogers!

Thank you Michael !!!  :^)

Edit: well, I have been looking at the RLink stuff and I can see I was off base there. On the otherhand, I don't think there is any reason I shouldn't try to get the other C compilers going.
POSTED BY: Paul T
Paul T
Paul T
Posted 11 years ago
Well, I guess I should count my blessings; one out of four of my compilers functions with Mathematica 9.

Why won't MinGW work? I don't know. Attempting to make a connection by putting the Mathematica bin, lib and include in folders of the same name is a non-starter. On the other hand, I've seen that the same three folders exist inside other folders, so there may be more things I could try, but so far I've only put the Mathematica folders into the folders 'closest to the surface' since that was how it worked with Microsoft Visual Studio.

As for the Cygwin and Pelles C compilers, it still doesn't even know where they are, and I still haven't figured out how to tell it where to look.

One out of four so far for establishing a Mathlink connection . . . not very good odds. And here I was so certain it would work . . .

Oh well, I suppose I could try the Generic compiler option (bottom of the page) next, but the more I look at it, the more I suspect that that too will also be unsuitable.

Why are the odds of success so low? What is it that I'm missing here? I got one to work, why not more than one? I do not know, I just do not know. I am more mystified than ever . . .
POSTED BY: Paul T
Paul T
Paul T
Posted 11 years ago
I thought I might have to do this, but I don't because the Mathematica 9 installation already did it:
All devices and libraries in the SystemAdditions folder should be copied to the Windows system directory. On some 32-bit versions of Microsoft Windows, the Windows system directory is C:\WINDOWS\SYSTEM, whereas on other versions it is C:\WINNT\SYSTEM32 or C:\WINDOWS\SYSTEM32.
I got this from MathLink Development in C (Windows).
POSTED BY: Paul T
Paul T
Paul T
Posted 11 years ago
OK, well, ultimately, it still doesn't work. I added a line at the end so that you can sort of see where it fails, but I don't know why it does. There is a certain point in the Import function where it says to Stringjoin the expression with this " < ", and that leads to the empty set. I've tried putting in " << " and "<" all of which led to the empty set. I also tried putting in " ", but that led to an extremely large output and a graph with millions of zeros, something like 4,666,464 if memory serves.

It's sad that this example has to have graphs in the output, because that leads to such ridiculously large nonsensical looking outputs that I fear people will be put off from looking at it. I will repeat the code in another post but leave out the ugly "graphical" outputs.

Here is my code:
  In[1]:= targetDir = CreateDirectory[]
  
  Out[1]= "C:\\Documents and Settings\\HP_Owner\\Local \
  Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8"
  
  In[2]:= lopass = Compile[{{x, _Real, 1}, dt, RC},
     Module[{a = dt/(RC + dt), yprev = First[x], yi},
      Table[yi = a*x[[i]] + (1 - a)*yprev;
       yprev = yi;
      yi, {i, 1, Length[x]}
      ]
     ]
    ];
 
 In[3]:= (* build the example in a specified directory *)
 fnSource = FileNameJoin[{targetDir, "lopass.c"}];
 Export[fnSource, lopass];
 FileNames["*", targetDir]
 
 Out[5]= {"C:\\Documents and Settings\\HP_Owner\\Local \
 Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8\\lopass.c", \
 "C:\\Documents and Settings\\HP_Owner\\Local \
 Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8\\lopass.h"}
 
 In[6]:= (* main function C source *)
 lopassmainSrc = "
   #include \"stdio.h\"
   #include \"stdlib.h\"
   #include \"lopass.h\"
   #include \"WolframRTL.h\"
  
   static WolframLibraryData libData = 0;
  
   int main()
   {
     int err = 0;
     mint i, type, rank, nelems, *dims;
     double *data;
     MTensor x, y;
     double dt;
     double RC;
     libData = WolframLibraryData_new(WolframLibraryVersion);
  
     /* read x */
     type = MType_Real;
     rank = 1;
     dims = (mint*)malloc(rank * sizeof(mint));
     scanf(\" %d\", &nelems);
     dims[0] = nelems;
  
     err = (*(libData->MTensor_new))(type, rank, dims, &x);
     if (err) return 1;
  
     free(dims);
  
     data = (*(libData->MTensor_getRealData))(x);
     for(i = 0; i < nelems; i++) {
       scanf(\" %lf\", &(data[i]));
     }
  
     /* read dt */
     scanf(\" %lf\", &dt);
     /* read RC */
     scanf(\" %lf\", &RC);
  
     err = Initialize_lopass(libData);
  
     y = 0;
     err = lopass(libData, x, dt, RC, &y);
     printf(\"%d\\n\", err);
     if(0 == err){
       dims = (mint*)libData->MTensor_getDimensions(y);
       nelems = dims[0];
       data = (*(libData->MTensor_getRealData))(y);
       printf(\"%d\\n\", nelems);
       for(i = 0; i < nelems; i++)
         printf(\"%f\\n\", data[i]);
     }
  
     Uninitialize_lopass(libData);
     return 0;
   }
   ";
 
 In[7]:= (* create the main source file *)
 lopassmainSrcFile = FileNameJoin[{targetDir, "lopassMain.c"}]
 Export[lopassmainSrcFile, lopassmainSrc, "Text"];
 
 Out[7]= "C:\\Documents and Settings\\HP_Owner\\Local \
 Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8\\lopassMain.c"
 
 In[9]:= Needs["CCompilerDriver`"];
 lopassExe = CreateExecutable[{fnSource, lopassmainSrcFile}, "lowpass",
   "SystemLibraries" -> {}, "TargetDirectory" -> targetDir,
   "Libraries" -> "WolframRTL_Static_Minimal"]
 
 Out[10]= "C:\\Documents and Settings\\HP_Owner\\Local \
 Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8\\lowpass.exe"
 
In[11]:= dt = 0.01;
RC = 0.3;
input = Table[Sin[x] + Sin[x*10] + Sin[x*50], {x, 0, 10, dt}];
ListPlot[input, ImageSize -> Medium, Joined -> True]

Out[14]= \!\(\*
GraphicsBox[{{},
{RGBColor[0.24720000000000017`, 0.24, 0.6], LineBox[CompressedData["
1:eJxVW2lYzN/7TrZkC0X2sZa1yBKiGyFbEloI06ZVpn2vmZqZppoWJBUxfS2V
JdmTaAoJrSpKxUhU1koqS/z6Xz2fF39vXI8z53Oe7X62c0yxPGRkIy8nJ9fc
V07u//7u/dOiK/f//+DvjhW6i2vq6d8VoPli6uadHczvlNC66tmaoXN+EK2M
uIG1isaODK2KjuSRCQWD24iegAT/T5UV3CaiWYhI29xwZPVLoqdiQjjbe03m
WaKnY3PT2139b2bl9NIz0bIn71Kr4hGi1RFppRt5tOk5/X42JiWnRf28xXx/
LmJ3z3lf9bWd6PlY4uQh+xP8h2hN5K5Tt229LIdeegHKn2bdSdLoQ/RCrD9a
2OBs9o9+rwWb9OhP2oo/iV4E1tFXRU9/MPpYjDPVl++9HMzoawnenB8fGHek
gOilMM2ZKvqPfYJobWSsVzI5M+MO0cvwdBl/+NhxVUQvx2hjUdr4NkaeFRiW
PfFbkMp3onVg/V+L0KKzk+iVuNVfIjqezPC3Ch1uemujzJl1XVz6OOTSMOtW
ooEAh+ruIdff9tJcYOqePdr2OhLS72oMWWird8b9bS/NXQ2LLx4qOus+0/oa
GOgrt3/2/kbra2A8RO2g5aePtL4Wvk522y6YyGh9LQ6s//Yyc9ttWtfD6ksb
s0tGlNP5emgqGy3W8SB+5NYh8YeR1wRv0id3HQ6OsTibVMv4y3qU7k4Qv9t+
kr6/HuefXf2aPY74lduA1JuF10t8Wmh9A2YmaCfV7eqgdX1cnGl5vTizi9b1
sfx1X98xip20vhHVM8xFZ8+30vpGWN4K9t4Q+47WN0FhTq123sQbtL4J6WFx
zZOUGPtthvGuFQ831TcQ/5sRXOjs1LbjA61vwQBWTsNKoYzWtyDVak19Yt0j
Wt8KI+Q/XjasmL6/FY2tV/qPPVxP5xtAVemiZ7/1TbRuAD1JyVx2wAda3wbr
sJfWG3OraX0bxtw+csVcMZO+b4j54xee/DqK+IEhLoxYtGn6KsIL1xBr9EY1
G+j97qWlhiiNq+oe9Owv7d8OV7MXdn13E43taI7qUhY7/6L925ERPya1rzl9
T7od5fecBj779on2G0H/88Whw7aQfXukVdRXHhosqKP9RnC9/K15xj1alxph
oK6QM+n+F9q/A6VjGlbU5HTQ/h1wOr6vxBSEV+4OHBnRFl23tC969+9ATvDn
Exo7+vfScjsxP06sn+ZJNHbCqHx5tkkS/Z67E4rt4on9llN8kO7Ec/us0I1L
uuj8Xcg8d/V3KYfwj11QKy89VVX1ns7fBY2+m3Rzb5M/S3chfmWbUehUWpcz
hk5D4uHqiq+03xjejTN01v+k+Mk1xnNhleqkfRSvpMZ4YJw6//M+Jl6ZoHhA
QHdoDdEwweaOWzZed8geXBNUHCh/PWI28Ss1QXm/ZeFP4ug8OVOIPo3qfHqP
8ARTJDaFH/GNyyF/MUWe50V9ZeuaXlpqimPNW+RKTN6Qf5nhSIsodfDkil4a
ZuDWuFvs/pFM55thwDvfEWIlsqfUDOvvfnt/cvtHOn83RmckRf0XxuhvNxad
rnasMSaauxs5xW+0TUsp/kl3Q2HTvne5E17Q/j0onidUCRycT+fvwYRZ8fuc
h7wn/vfgdP+ZAy+8p/gk3YOopmrzI9mEZzlzCHM+rTD48JX2m6O/kVGZqTbh
m2uO5w8K0jeo3aH95kjQ3Ln//fIaOn8vlmq1DObqkj9jL0Z+5lbfsaL4yt2L
8IVOagY2REv34s6vH0vmPmL8fx8Kz6/9tOnHG9q/Dy3dvzO0Ft+m/ftQUaUz
38SP8C/dhwHzb39/n1JF/O+H7fatKdE3nxP/+yE5oO/8opvyJ3c/fvQ/Pqam
jPG//ag9kH2ky5XJx2zILY0+4bCD8M1iw2wFJ8VcxvgTG0EpUYV5evK9NJuN
qJZt9UuuUH7kspGrX54T00D+JmHj6RDbWB0x429s3MrPG920mM6TsVH+O9lG
6MXIbwEN7X4O93+SfVkWmBU675nqUAYPFjgfMiVneivFD7YFplhfGtXQRXjg
WkAx6ukH8/PEn8QCnSM5SQ3/CM9SC2T5yqc1Sgf20jILfMjrXODLJVrOEovc
qyMjDOj3LEtsUMlT7JNB8sESLj79jdz/UDxjWyLzRuHyBivKv1xLJJ2/0PI3
4TPJb4noVB33LYebSX5L7LPH6ln/KF7JLJF80Fy+U52pR6wwQKV9/JJs+j7L
Cve7tOY32DH6t4LeTZ3xaVsZ/fes7/k7rdKBaK4V5sRPDXAQ0+8lVnCdeXY1
u43yv9QKAiXdsWkPCE8yKzRqGEbMukp4lLMGe37B+eHzBL3+w7IG3LMldcMZ
f7KGID/EtetmeS/NtkbD7DUWW7lUj3Gt4TPL67SOpJTkt4ZqP5O33UKKb1Jr
JL6784e/ijnfGiNrfNf3+cX4nw0eVN/KKBtENMsGiwdHrxa/pvgAGwgbuLf9
bCifsm2w1r4tvSotj863wbSW8E/JVYRXiQ1GatyJOF78hfBig8EJIyYWjqf8
L7PBz0cVRzUETH1yAHs7ugyGTK8l+Q9Ab0aeXE3cVTr/AFSHeiUkBFD+YR8A
+8gNj/aR38j+B3Crdl1zVjzxL+n5vem7wHOvGbwfwK6ii7qZfPJ32QFcWNdq
X6T6muS3RYz58F1+fa+Q/LY4Wey8qlCulPRvi9kRak73bF6S/m0xteksv92A
qQdscfDl56RJWlI63xY7wvaeTBzC6N8WgxpdvxZ2k7/KbGGw+D/FY95M/WuH
xr6hyWWTyd9Zdrh8uezvuynkX7CDw6BJ1pFMPcy2g3fR+kCTvd0kv11PfnAx
HPeE8q/EDuK0sBKHcLK31A75nL/9mmsJ3zI7BNl9y37g2Ejn2+Pyjf5/zSMI
Hyx73N6+XblyKOU/2GP92nPHh86l89g96wPF94cmE39ce0zZ31aw7ibhV2KP
G9s3pp4fQ/iW2mOYTZtttOuAXlpmj8QHzcU3PlJ+l3MAp8vQUX004YflgCKf
lTN2XaF6GQ7wPbh0D9+A5GE74Jduif5FY+Kf6wDHTv3TN2xI3xIHeHmbX9xX
QP4rdcDJP7nv1pUx8c8BB/Wb9paaMfW5Iz77h79pGEH1CssRoeJhSJjIxB9H
GKtVzomsJv7YjvioMcPv5nJG/46ovp/gPi2U9CVxhPKWVRmr5zL50hGZ8Ve1
povz6XxHZPW/U/H2VBn5vxMOXJffFutI9TnLCQrvOgbzxhENJ3wb/vN8O5t+
z3aC4YuOgd/XZ9H5Tsgb86u/oRbVrxIn1F4+3MwBI78TUhPD2+5ZUDyXOeFY
8NLOY4OY+H8Qy0PU2iojaT/rIM7mxxi9W0X1Aw4ib07ZTatmqjfYB2FlneSj
lkb5mnsQTZ5GmkPLvxP+D6JI4eGoqoPthP+DULN6/0x/IOV72UF8Mc7Od/zQ
QPI744ml6uWgsRRPWM74M8jgidCZ6gs4I4ztvr3tKtXHbGdEPh00n2dD/HOd
cahkldkPIelb4oyNKzcMMnAgfEud0W6meNstnOKJzBlx6UF5fU9V0vmH0HS0
7+WXg4kf1iFUlKStCU6h+h2H0Gxm971xD/Uz7EPYLvJPjdV8QPIfglTWd9Lo
tdV0/iGs/5dkZmBK8Ul6CNdMBt1PukP+LDuE6EHKkbkvGfxz4HpUv9hjPOVv
JQ7CF/fvEImYeoCDiliVK7XDCd+aPbTBgzOjU5j8zMGk55sX2gdRfDTkYObL
bdtleqQ/NgdNGx0FqyWVvTSHg3/D6gKM5Jh+joN475VKK/WI3xgOfm3dcTm5
hPiVcNBeYbk5VUL4yOBAOuvQmk1J8uiVj4OsrYqZKv8Rnks56Awv7vp1ktZl
HNz8My/HKor2t3AwpFXB4Xkk0w+7wLxvjfqZcSSPkgs+KJb/mSijeozlAklA
se6MrGKS3wVypbPXjjj2jOR3wfj21JORga9Ifhd8N+C9t0wgf2C7YGBX0rL7
Dwn/HBdMVFKN+zuO6iOuC6w9AteZnaX8H+OC0D3Jp4vYFB8kLvCc+exUnxVU
L2S4YMe1yH7PJjD4dsGdc1o2244+6KVLXaA3N6plUfIr8ncX1JX/alN5S/mu
xQVh68Jf+JyjfCjnir7To9cHFFB9rOQK0Up/vaDARvJHV7gb2A46ZfOil9Z0
xZ+dNuHmWudJfle0JF6xtFUk+xq6wuvmupaP10kf7J7vhc/PSltTRvK74vmj
VUU6b6n/57pCZ6Wnup8K5e8YV1gpebSfHdxGeHaF3WC7wtMjfvbSGa5of3Kk
6LdcN+HbFYq/tHPGdfzppUtdEeuyxFOlm/p5mSteKO+OyR9G32txxemdMzI+
pzD9uxv4K8q1LT7nkvxu8N3SnSg8/5js7wZlyVTDZ6E0X9J0w4zl++XeOD4k
+d0wZKGD3oLoh737Dd3QKjOYdC6I+nO2G561bfo7XJP0y3HD4CEll28+onjF
dYPTuvfdI70oXsW4Ydh8Db+vmhSvJG442Bbm5lND8SHDDUGyL5z/FDJJfjes
OrG4YFc7xZtSN5i6uvUblM3EWzecWv9ZYFFN+b/FDbltd50e2jP1pzscRatO
ntWhekXJHbFl2jqnCinesdyhsy50Z6gHI7879qqcmX979HWKT+44ubBs/z1X
yg+G7jDILwyyW1ZC8rtDKfG6XYxFDNnfHXz2mqrsf7Xk/+74Ne0nyrqI3xh3
yLsOGR6uwODfHdFdrurvxlC8ynCH+pLA4AAO5T+pO46M3VeecJvwU+qOeXbG
/16Hknyynu/vjjIZe5zw0uKOZWkhnrGTS0h+D6StrfdbYkH9o5IH7pdrF9qZ
0byE5YHdNq8sr9pSv6XpAeOK6v0quyj+wgOjnIvf355aSP7vgdbuQW8Tq96R
/3vg6s01d2blU33D8UC9zuHCynKmv/XAuclDy/oz6zEeUJ/rObJrI8VHiQcU
I3Ln+umeIfk9sPHn6+UHJpN/ST1gvdxn9/gG8pfSnu935AYPm8b4vwe+fFy4
9NOlX+T/HggKW+A+Rp/wJOeJFe1lszXSyf+UPLF/wkXHfXdpvsfyRIywse6L
J9UDmp5wfRBZq3yQqU89cWBPRtxlq6dkf0+sC9mQ9HYC6Y/tiWz/F7OVJzaT
/3vCvlugtLGYzuN64hQvwVdxPuE5xhOxXQdPL6rvI+2V3xOzU1yWbArt20tn
eGLc2PFKTmpESz3RoquiElgn10uXemK1aMSuSYdIXpknLt2p0oU6g39PzMw9
5znyCs0H5LyQalLodfAM2VfJCx+qd1RV2NM8g+UFzpLUK9f+vCb5vbAm+UJj
njvFU3ghZsa2C2c6CM+GXrjuOCXu4iA6n+2FZ2ryfdUuk3wcL4zf43s2P51o
rhfcFpxbNi6b7BHjhVetL7nnvZj45wWtCvspI8+SvTN6+Ik9pp++S0T294JC
TkUAj0V4KvWC6sD2m36VVI/KvBD2fWtD+lWiW7zwsZAlKT1Lv5fzhvz4+DmT
xlJ9ouSN4z/NTi+Jp3jP8sa8mDrRbTWqRzS9Ucwb4Fu18RPJ7w33el/vSX3I
vobe2PduSYl8F+mL7Q1ZfqHjtbLjJL83rrVtUJG/RfmV643MeRVzlJ0Z/Hvj
5tqdGdvzKF5JvPH6YERljSeT/7xhdLsoWG8a0195Y/SNHWZnR5J8pd74/e3S
9Vz3XJLfGyNnqM77pUj2bfHGwlHzPM1taV4q54Px2yr3Js+kelLJB/POwr9k
O+P/Ppi5ZY1x4VvKF5o+KJAa7n7/tYjk94Hmm66/96ooHxj6YHT0i6aEl6Rf
tg++pow4c3sezV84PfsDW289mEL9LNcHkk9CR68rEWR/HzzcHdiRrkj2lvhA
1lp53+Q0+UOGD0Y93auQOOI32d8HFwNSlg41ZfzfB1JFtaQluoQfmQ/anF6k
bF1C6y0+4Nc/cmftpv1yvkjRfqpq9JrBvy+SMkUOJx5T/8zyxYYTXpdZ+xn7
++JfnN7klHrCD3xRO/r51aGtpC9DX/wsLRKsY/0g+/tCq2142UIvxv998dHh
De9eoXwvP1xfrGwyf+Bd3r+XjvGFhsrP65oTBvbSEl/EKPz4eXHPgF46wxf7
jWQnPHIZ/PviYr7/Eb+qfxT/evhz6ZwwRUjze5kvRLMmnstj5v8tvsiNSFr2
I4axvx9eJI/4IjlO9Y6SH84O4r1rqmbk98Ps93drGi+QfjT9YDyl22IYoz/4
oWP6j7XPhKRfQz9Eq0tM3BVIPrYfqjtPPbiSRvbg+OH6a23Vpot/Cf9+OK1+
N77rHd1HxPjB5aLWQLNhhC+JH6xv+Ux3HUrzmAw/mG15q5Dim0f+74eZuxLH
+PKp3i71Q7Fk/Obt3+g+ReYHORvb2sl1Fwn/fnjfNGpBejDTf/jD+MLVI4Jq
Rn5/WMaavGp7S/pi+SN/1qDaC/nkf5r+yGN1PdrZROvwRwInNGjPcNpv6I9l
7rHKOZtpPsH2x6DJQ9Lb1lD9wPGHbuDLh4K5TP/sj6/hdXJ6q5n63x8B+weP
SfEgWuKPx6vr92zLot9n+GNg0dXpCdkkn9Qfi+a4rtAXPSP7+yPZv3TiBaZ/
kvnjv3UjBxm2Mvb3x2urgxNcVCleywXg46rI9q0+lD+VAtA+eETJGg7FM1YA
DLMmTJCNoHymGQA5a+uOtuKlvecjAFE7zAQNf2j+ZRiAPxppSwU7Kgj/AdDu
HGazQi+H5A+A8TLjPk+/kD25AfiSuSl3tBPZOyYAuz1Nrt9SJ/+VBCDOrFaw
fxr5d0YALB4YPzvG+Je0h3+1xk07IokuDUB05fsdT+Yy+A/A0htVbRllVB+3
9PB7uqZGX5O5/wpE/MMlU5ZZkn6UAsHvr7Zo0GLKb6xAHC0Y+d+CYqI1A+EZ
vb1eyZn0hUBUDZ+5vciK6g3DQEgeS0+r7SM8sAOhz75r8XZdv16aE4hEOZ2V
m/8N6KW5gSgUVyrWvSe8xwTiRtPKNkEerUsCYTx5ucnUaib/B+JwfUPi0yuk
D2kgONX+juFmpK/SQJg6/ci9EsfM+wJRO+bfX9tHZM+WQCwVujx98JZouSCk
+25Nb+Uy8gdhzfovrXsOkX5YQZif+PyCZTfpTzMI3UPPVh+bTPpFEI5879vs
u4n0bxiELQM7LtaKaJ0dhJ0ai2+PuMXEvyD85fzc2ceF4iM3CN8UQhKzCgk/
MUGwnbYx4N6GArJ/EFxSquP2tD8l/w+ClUTRWjuY8oc0CIeMn/dxzSE8lAYB
it2rOptPEv6DcKlOo0wrjOaHLUFQe15RP9GSiX9cCENsOhROk/wKXBy7xHK/
WczggYsqC8n5HHvKj6pcTP+r13HBk/ITi4vNEdopdt7Xeml1Lhbkm24tf075
T7Nnv3i33Ylqmpdoc+GUnWZQ5kL5G1yk/t4Z2PSd6mF9LrS8t7/cNJvm54Zc
rPjTvPYen/pLUy4K5+8duDX5McUXLswWaPWpaqP+yI6L2QLPpk+txC+Hi9J+
Ziul70g+by7k2ate8WZSfOdysefe8pJEpr4ScRHomtkWbPyI7MGFSLHlT8Ur
6vfjuZCcVFvivInm5xIuTlbWVBfYUP+VysXiYyclWw2pv8jg4oLzyLPH9j/p
/V4mFxPU1tnPiyD9S7l4cmnWj5Rg8ocCLmqtb2/U+cX0s1xMTdirmh1PeKri
Ys2Ylrvpd4mW9Zyfsklj4XDKJ01ceL0//x9rMuGxhYtJYVkZO2zJnl1cGF0/
K6g4T/4mx8N6j7ONafIkvwIPD7IGHY65RfWOUs/6vmFHQph+QJWHLTasy+XH
iV8WDxrvBg04/578W50HxZKf1QcvEh40eVg+5sms96Mov2vzcOrOo4mfjxMN
Hr6q/H7QrUzxQZ+H+glXDIrqmXzKg7Jep/UlP6qPTXlot0/deXwhycPmQe+o
TJgVRfWIHQ+Lol9oXram+pPDg7Nw84skW+pfvHnov60kpY85c1/Ow60Zlg0P
11J+E/FgdCMu2u486S+Gh/3RfxaJnEm+eB7sMoa3zlem+CPhgZPEzT1XQvZK
5WHFDdFQnWcUPzJ44H+OlEvXIv1l8rBguaXazufUL0t50K5q7NchpfxRwMOb
lrSkPSvpvUApDw3GA48WH6H+tYqHhyU1rw6MZep7HuQervVcVE35v4kH7uWO
zSrZ58j+PGDCt1cDX9M8qIsHseSekd0ipv8JRmvKuV/f51G+VAhGzurX4j5z
ad6qFAw5BXel40J636DaQ3d2GVVpEb5ZwRAI37yKuET1u3owRDe/akel0zxa
MxhNG+b5LkymeZd2MCZPXWIoL2Dm68HQGbctIvYczcf0g3HXSss2wYLmD4bB
yPNu/Ciyo3mMaTDE53d7dr6g/pMdjE8Z4r335tN81C4YM+v0NlsnkD9zgrHo
uwMPVpS/vYORqtfw5qQn6YcbjPNZf54etqZ4KgrGZrPEDeajSN8xwfDkx/xS
LCP54oPBz3iV/t8soiXBeGvxweTqGbJHajAKTQYsXX2LidfBWHZw/zWZO8Xf
zGBcy8z6kZvD3E8FY2lbxo03zxn8B2P9uRCvbhPy99Jg2Hx+w/r0nvyxKhj/
GsaqtDVRPSrr4edw3prR2+h7TcG4Y/h2xbZ8sl9LML5+1D1qcYjiWVcwOCNC
V86UxJP9Q1D7/N2Ff6vvkf1DsOFvc2WGMd33K4VA+2bM9IXnGPyH4HP/ypOz
jjD5MQT8ibGxHf8x+A9B+OHk3MUfCb+aIai0G7lllRLlQ+0QiD3/jr/hSPEK
IQhd46raHEH5Wz8Ep8quy9rVKH4bhuBKS81S+0fEj2kI4ubNzhFz0sn+Icj7
KT7UvJ3uI+xCcImX/OHL1qRemhMCh9s6SeciaR7lHQLlIuV7CnsoX3BDcLEq
RfUVU/+IQvAjY1eO7jVmHheCyG+mQpN0sk98CDJKB3Bn32L68xAsLzr559Ny
iiepIbAp6759t4Dyd0YIwvKPL+bNp3loZghGBLw0r86meY80BIP8nHZP1KX+
tiAEU78ujAt9SPgpDYG13MGVl9ypH64Kwd+Hd3CRT/fNsh59qwVWjVpL/XRT
CEZ/faN8w4LuA1tCsOWH/8MbHYTfrhCcnq5nNFyX8rUcH/K8Ci35tXQ/pMDH
6UktFXO9aX6lxMeOkuGXvZcSPlX5eHd8v8cwI+a+mo9d+hMWt6rS/ZQ6H9rs
o5dSdOn+SpMPlZgPdwPNaV2bj81x9hO5zHsD8OH+z9jjxUmST5+PjBMxVkcq
Gfzz4ZBUN0Ioof7BlI/SE9oPbBOpn2Hzod59ZK1CcRrZn48BYRz16auZeTMf
hzYGDXyXRPHUm4/D7TvLjJWZ9xl8PN++tSrzOt1viPjYMm+S+ugQmr/H8FHw
yLl4K5/iVzwfx+xGCpeeZeYTfGwt0Hnqy6X+JJWP9i03oi6LaB6fwUeom+0t
ewuKP5l82N0PmGl1nOpPKR/ZMbPqTu2g+XABH1abkv8Omk54Lu2Rz9tlXeoC
8tcqPu7UXRvqrULzBxkf3Aui8yuSKT429cj7QuXTiEXkDy183J+LwMifNE/s
4mN2+tYxt48x9/8COIkT4qxe3CP7C9ARoadgHEvxSkmAJTd1wn7NZ/AvAMdt
3n+yp4QPlgDOJcf9iqwpPqkLELPjZ/SpvkRrCqDuIxScuUn9lrYApxrypq79
RHiBAKOSz4fqOhE+9QWYdED5wYJIqp8MBagz9XNxfk78mwpgkXkzJuo6cx8n
gGu7h9aWkeSvdgJ8k25rmill+k0BVpwzm5uic5fwL8C5scWCex1UH3AFGPvk
3ChFLaqHRAJIgt/eWPWVqccFyO4XKHqbRfVDvACf47aO0j5MeJIIUJwXdmL2
I9J3qgD9Fm4W7fMi/8oQoHOA5XF7c/KnTAF+mT+Wa1Kl+zapAPKRK/p5NxNd
IED863cp0m7CR6kAt9apBL12JP+sEqCq5bHCqxfUH8sE2JPv4TH1Ofl3kwCp
ExfwjRdRPdoi6MEr+8OZc3S/2iWARmtjQlcs1ddyQtSOyy8/pUzxR0GIQY1o
3RZP83YlIX4ob+jjriWPXvsLcU087afNq369NEuIeNkbg7qzdP+uLkSqS8e+
6AO0rinE7qtnh9QcpHigLcTT33Yf8ZbiB4Son3GyJDGHyf9CqCwVzpS/RvnX
UIgv8e5VGYsY+wtx+wbfU6THzNeFMFssn3BuMfUXdkL0042JmlxC+OQIEXLj
eEGSF53nLYTTKvE28wvM+0AhRq1e3fK2meKTSIgPJ9ePnjyU6BghHuamuOyN
pf3xQmSIxipDj3nvIURS4su8pGs030sVQr8m45RTw2GyvxC84BOGkx7TPDVT
iPTCrOfXyyk/S4VY9iC0onYWxYcCIbZUr1J5LrtF+Bfi9bjjGg95NE+oEmKq
hYtWQCf5v0yINRpti08kkz80CdEd4FpTV0b9VEuPPCGe/3JGUL7oEmLVgGHj
6srJX+VCwX2VOH3wFHfK/6HY5/XV+7Qr+btSKIo2PL6lMY3ikWoowqWYFXWC
6m9WKOqvVP8LP0H5Uz0UXksOWY08SvjRDEXBgxGja4dSvasdirH3aixu1FO9
i1DEDts3VeJG8VM/FOlfGmI0jhH/hqHIu8EycykmPJmG4vTq+VH1aVR/sUMx
YIve5bBh5P92oRj5be+FpzLq1zihmPTjjlOMEt3XeIdi3GEjwaSlFN+4ocj9
06lSaEb2EYViy+k4NeV7pP+YnvOG/Tbsv4Tqw/hQsANfTyjNZ+q/UKScMdrR
coXyQ2oo+JMO/b0wkPwnIxT2j/WDcxTle+nMHv2eCrs6ehTR0lCsfGRUo76V
fl8Qirv633L3pjH3WaHo34fvMXMC4bOqR76fQ4JcM5n3ZKHwz7/2wruD8NLU
w3/b2donpwkPLaGYOzat3wMp+WtXKAL3Lrn3ZSN9X06EgvQ791X2ET4VROjq
3zl4xCjCs5IIlU3v1zQeoPc0qiIk2J+pWVtMNEuEisrCmMgF9L5GXYS0GbdZ
F9tIPk0RHsuH/J11ge7rtEVgn3Q5Mfkk875HhDNaT/YOmEb+qy/CVYePT7Wu
M/2/CGEr8o7+riR5TEWQGWQpjt9A8rBFePK2a39BAunfToQPDkOv6QwlfHNE
eO8+tbXpNcnnLcKVUo3jSfHEH1eEOtXkuL9mtC4SwdYtr2yNkOJfjAgP719d
NKqI4km8CC1dubGDjpO/SXroN2VNI3fSfC9VhAG1NRMrjQt7/SdDhMXm1lkL
m+k+KVOETdn2Rsd2U30oFeFr7V3RsP50X1HQI49wzJa/FhRPSkUYvjeku8GD
7F0lwvSkayp7fZn3RCIole2IRS7RTSIUffIpvcon/2zp+b1OsqrORvpelwj8
jYEatnfiyP5h8OrW9P1UTflQIQz7/0RoVu2neaRSGPzvPFmTO4fwrRqGk9c5
jevmMPPxMGg9kpWdSaJ+Sz0Mpgvr2SMO0HtZzTCIb9WsnTCK3stqh+HDgKdX
FS4w72XDYBazvZQ/gIn/YRhnGZSdzSPaMAxbH4wJMc8heUzDsCjv9CHn9Uz8
D4N748wrRgVU/9iFod8qn3efo6nf4YTB87JVl04a4ds7DDPU9z1x5NN7Rm4Y
WEaV8f2koYT/MGxr7O4e40/5MyYMJnEj//C6KR7Fh2HevD7Zud3M+4wwNLU+
ydCwJvymhmHLqpeTTwwk/8oIg/zhwqbJzHvZzDAsydvaOGwp+ac0DH891Fo0
W8h/C8KQsGfJcnc2yV8aBrellUEplfSeqioMlxWfbP48jOKlLAwD+szxflJE
601hyNy76oVWPu1v6eFHc0PfV4GE964wnPgxwHiOLfP/W8KhY7BHY3g95WuF
cKxICenjWkd4VgrHj9rb9flviFYNxyLftGjDD0z+DwevOvpDoQ19Tz0ckxxi
7r7cQ+dphkNTXJ094SDxox2O7PSKx41iwjfCUV6klr5tPvXP+uF4LeIFfjFj
8n84jjuPGxP7lvRvGg7Lzn1Kf9fS/IAdjjmvd20UulF8sQtHRO1X1pRnZA9O
OML3W45oncngv0e+9zOOGg5j8n84nnk2yjWA9C8Kh2rAiiPyJsx9ZDjmDa7d
NzyU5ivxPfu1uy9vtqZ8IglHos+dYpsrdb10ag9/2ukeTkNpHpHR8/2runuN
rJn+PxwszTfP0qcnEv7Dwf5hEVvpQvmlIBxPOD1N+RGyb2k4NhUVcD6+YPAf
jpThnZJ7X5n+r0fe08Ynmg4z9V84lrwsC8p1fkL2D8diLyUb/kK67+kKh9V+
jv4p5v5aLgJz35SXWBRQP6sQgbWjYoM4V5j7wAhkzG1vDgqhddUITFV+ub/4
JHM/EgH7D+xGr0zCk3oEJKnN360H0PxDMwJ+guHj1l4hvGpHIPdE3/pLmcz7
3gj8u2myib2CaP0IDAneMXxzNvOeLAK7uzP1JqnS90wj8Pfbsq4tErrvYUfg
g/+9usRbzPwnAhMKY/vYp9I8ixOB+XHjD+s30jzTOwKbkhb8u9zK3L9EwPzq
fa0p/aheEEXg2jE9tw3+5J8xEdA/6t3qZU/5Kj4CXVwvbn4i1eeSCNwcMkV7
8XSq31MjsHhA7NecLOa9SgRC9OKiLwwk/8rs4W9d4eq5gcz9dQT4LV36B/6R
/xdE4Dx/jFl9LL3nK42AIa9gZ9Vspv6LgLdzh5bKe+rfZBFotdCQawxj8B+B
GtWNJSkzid+WCFwfdMMufzPhoysC90VpUwc1MPgXo8Tm1NoOIb2fUxCjtHNM
mlSNaCUxVBVHf4moIjypitHS0t9iEZj+X4w0VrHs9lWSR12MmX6C2Tpfyd6a
Ynw7ve5J2yi6n9AWI9t0mJuyHvknxCiefeTodFWa7+uLwa82ueSdTP5sKMav
Baw+lvcY/IshW8YaW3aG+iG2GP1Kh4UtW0r82Imhdudo3wKG5ohxPIL3s3Q1
/d5bjFSf+gadYsITVwzj29PfX9UkfIvEcB9penO3C/WjMWKsrVVf+mUV1bPx
Yqyc1vBj0zDKhxIxvsQs+HonguhUMeRX7brhlEX4yhCjquLD/fQLlF8zxTjC
l1+Ra0Tvp6RijPAY7bBzPzP/FWNARMM7J+Y9YakY7Ss3ThE3kn6qxBD9uOcZ
LD5H9hdDL/HsHmVTmo81ibHTvPRbcinNz1rE2Ku6b+6fsTRf6+rR3/b9PH11
mlfLRWLt89Ef2E40v1OIhOZsD479GXq/ohQJ8dnVQ14GUjxQjcS1D55zJrZR
/86KxIdnK1alvqP6Xj0SySYp9dwdNI/TjMRXq5y89mHUj2tHYsMtmz93HtN7
KUTC/0THKPYPmsfqR6LdOdh35Bp6j28YiePLuTqr3akeMY1E2+RDic6PKB6x
I6Forikrukj82UUiIcfpqyiF1jmRqHxU/1T3KcUr70h8OylXpP2e+OVGwvbQ
7PcqoeSfokjALt7u4iHCU0wksg+f5vWroPwVH4ldl6If1l6gelwSicZaO626
RFpPjeyx3/KzIY20PyMSz1v11oweTPVPZiS0Jy+Ztn8q3c9JIyGffSPpfhbN
cwoicT819M7C18Rfac/5QVOUAzqoH6mKhLp2RfIhLeZ+v4f/VSHTHl6heNXU
Y5936c46z4m/lkh0FmzTdn/J5P9ImFg5LzrgTvFALgq74j7kdfaheKUQhXu/
7Du/FRFelKLg1/3ZoHk45UPVKHw9r/t4cx3FJ1YUMmYYN88W0TxWPQqmujYn
Pzkz70Wj8O35SX6aCd0/aEfBaPSt3McgfCEKW4ZYterdo//vpB+Fu15q8z8W
Un9qGIVA/VVza6cx9X8UttWlbuCeYN7jR+HAeQXuo7fUP9hFYffNsyETrlB/
zInC9EWuK36tI3/z7jn/z8O/5Z+Y/i8K3ycHG1w8THgRRUEgu2S1OIR5jxGF
qEucYM5WwkN8FLZeTXi+i0W0JArSFWpbjRzo96lRKPpZ7s9zo3iQEQX1d8tK
OlVpfpcZhVHbm2rGHaJ+QBqF9J0NjokadD9YEAX21tZNs25TPi2NQtYBjQE/
Qpn5XxRa5w1dP2Uk+bcsCpZzvzf4niJ+mqIQc3MMv+I3894nCgYRBpdGL6X7
864ojHMuC/wSSLRcNE6sXqp3O5/m8wrRuNM293JWO8UHpWhsVwoo/PCYwX80
vJ/Usnc0M/k/GpO+urRnzaT6Rz0a898acTJGk3ya0bjw4beTthpj/2hUzHy3
8yOfed8ejTHlh6s7ntE8QT8adY5JEvamrpz/ARysT5U=
"]]}, {}},
AspectRatio->0.6180339887498948,
Axes->True,
AxesLabel->{None, None},
AxesOrigin->{0, 0.},
ImageSize->Medium,
Method->{},
PlotRange->{{0, 1001.}, {-2.98723801101878, 2.9757153662914786`}},
PlotRangeClipping->True,
PlotRangePadding->{{20.02, 20.02}, {0.11925906754620517`,
    0.11925906754620517`}}]\)

In[15]:= inputText =
  StringJoin[ToString[Length[input]], "\n", Riffle[ToString /@ input, "\n"],
   ToString[dt], ToString[RC]];
inputFile = FileNameJoin[{targetDir, "input.txt"}]
Export[inputFile, inputText, "Text"];

Out[16]= "C:\\Documents and Settings\\HP_Owner\\Local \
Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8\\input.txt"

In[18]:= SetDirectory[targetDir]; outputLines =
Import["!" <> FileNameTake[lopassExe, -1] <> " < " <> inputFile, "Lines"];
ResetDirectory[];

In[20]:= output = ToExpression /@ outputLines;
output =
  If[MatchQ[output, {0, x_ /; x + 2 === Length[output], ___}],
   Drop[output, 2],
   $Failed
   ];

In[22]:= ListPlot[{input, output}, ImageSize -> Medium, Joined -> True]

Out[22]= \!\(\*
GraphicsBox[{{}, {{}, {},
{RGBColor[0.24720000000000017`, 0.24, 0.6], LineBox[CompressedData["
1:eJxVW2lYzN/7TrZkC0X2sZa1yBKiGyFbEloI06ZVpn2vmZqZppoWJBUxfS2V
JdmTaAoJrSpKxUhU1koqS/z6Xz2fF39vXI8z53Oe7X62c0yxPGRkIy8nJ9fc
V07u//7u/dOiK/f//+DvjhW6i2vq6d8VoPli6uadHczvlNC66tmaoXN+EK2M
uIG1isaODK2KjuSRCQWD24iegAT/T5UV3CaiWYhI29xwZPVLoqdiQjjbe03m
WaKnY3PT2139b2bl9NIz0bIn71Kr4hGi1RFppRt5tOk5/X42JiWnRf28xXx/
LmJ3z3lf9bWd6PlY4uQh+xP8h2hN5K5Tt229LIdeegHKn2bdSdLoQ/RCrD9a
2OBs9o9+rwWb9OhP2oo/iV4E1tFXRU9/MPpYjDPVl++9HMzoawnenB8fGHek
gOilMM2ZKvqPfYJobWSsVzI5M+MO0cvwdBl/+NhxVUQvx2hjUdr4NkaeFRiW
PfFbkMp3onVg/V+L0KKzk+iVuNVfIjqezPC3Ch1uemujzJl1XVz6OOTSMOtW
ooEAh+ruIdff9tJcYOqePdr2OhLS72oMWWird8b9bS/NXQ2LLx4qOus+0/oa
GOgrt3/2/kbra2A8RO2g5aePtL4Wvk522y6YyGh9LQ6s//Yyc9ttWtfD6ksb
s0tGlNP5emgqGy3W8SB+5NYh8YeR1wRv0id3HQ6OsTibVMv4y3qU7k4Qv9t+
kr6/HuefXf2aPY74lduA1JuF10t8Wmh9A2YmaCfV7eqgdX1cnGl5vTizi9b1
sfx1X98xip20vhHVM8xFZ8+30vpGWN4K9t4Q+47WN0FhTq123sQbtL4J6WFx
zZOUGPtthvGuFQ831TcQ/5sRXOjs1LbjA61vwQBWTsNKoYzWtyDVak19Yt0j
Wt8KI+Q/XjasmL6/FY2tV/qPPVxP5xtAVemiZ7/1TbRuAD1JyVx2wAda3wbr
sJfWG3OraX0bxtw+csVcMZO+b4j54xee/DqK+IEhLoxYtGn6KsIL1xBr9EY1
G+j97qWlhiiNq+oe9Owv7d8OV7MXdn13E43taI7qUhY7/6L925ERPya1rzl9
T7od5fecBj779on2G0H/88Whw7aQfXukVdRXHhosqKP9RnC9/K15xj1alxph
oK6QM+n+F9q/A6VjGlbU5HTQ/h1wOr6vxBSEV+4OHBnRFl23tC969+9ATvDn
Exo7+vfScjsxP06sn+ZJNHbCqHx5tkkS/Z67E4rt4on9llN8kO7Ec/us0I1L
uuj8Xcg8d/V3KYfwj11QKy89VVX1ns7fBY2+m3Rzb5M/S3chfmWbUehUWpcz
hk5D4uHqiq+03xjejTN01v+k+Mk1xnNhleqkfRSvpMZ4YJw6//M+Jl6ZoHhA
QHdoDdEwweaOWzZed8geXBNUHCh/PWI28Ss1QXm/ZeFP4ug8OVOIPo3qfHqP
8ARTJDaFH/GNyyF/MUWe50V9ZeuaXlpqimPNW+RKTN6Qf5nhSIsodfDkil4a
ZuDWuFvs/pFM55thwDvfEWIlsqfUDOvvfnt/cvtHOn83RmckRf0XxuhvNxad
rnasMSaauxs5xW+0TUsp/kl3Q2HTvne5E17Q/j0onidUCRycT+fvwYRZ8fuc
h7wn/vfgdP+ZAy+8p/gk3YOopmrzI9mEZzlzCHM+rTD48JX2m6O/kVGZqTbh
m2uO5w8K0jeo3aH95kjQ3Ln//fIaOn8vlmq1DObqkj9jL0Z+5lbfsaL4yt2L
8IVOagY2REv34s6vH0vmPmL8fx8Kz6/9tOnHG9q/Dy3dvzO0Ft+m/ftQUaUz
38SP8C/dhwHzb39/n1JF/O+H7fatKdE3nxP/+yE5oO/8opvyJ3c/fvQ/Pqam
jPG//ag9kH2ky5XJx2zILY0+4bCD8M1iw2wFJ8VcxvgTG0EpUYV5evK9NJuN
qJZt9UuuUH7kspGrX54T00D+JmHj6RDbWB0x429s3MrPG920mM6TsVH+O9lG
6MXIbwEN7X4O93+SfVkWmBU675nqUAYPFjgfMiVneivFD7YFplhfGtXQRXjg
WkAx6ukH8/PEn8QCnSM5SQ3/CM9SC2T5yqc1Sgf20jILfMjrXODLJVrOEovc
qyMjDOj3LEtsUMlT7JNB8sESLj79jdz/UDxjWyLzRuHyBivKv1xLJJ2/0PI3
4TPJb4noVB33LYebSX5L7LPH6ln/KF7JLJF80Fy+U52pR6wwQKV9/JJs+j7L
Cve7tOY32DH6t4LeTZ3xaVsZ/fes7/k7rdKBaK4V5sRPDXAQ0+8lVnCdeXY1
u43yv9QKAiXdsWkPCE8yKzRqGEbMukp4lLMGe37B+eHzBL3+w7IG3LMldcMZ
f7KGID/EtetmeS/NtkbD7DUWW7lUj3Gt4TPL67SOpJTkt4ZqP5O33UKKb1Jr
JL6784e/ijnfGiNrfNf3+cX4nw0eVN/KKBtENMsGiwdHrxa/pvgAGwgbuLf9
bCifsm2w1r4tvSotj863wbSW8E/JVYRXiQ1GatyJOF78hfBig8EJIyYWjqf8
L7PBz0cVRzUETH1yAHs7ugyGTK8l+Q9Ab0aeXE3cVTr/AFSHeiUkBFD+YR8A
+8gNj/aR38j+B3Crdl1zVjzxL+n5vem7wHOvGbwfwK6ii7qZfPJ32QFcWNdq
X6T6muS3RYz58F1+fa+Q/LY4Wey8qlCulPRvi9kRak73bF6S/m0xteksv92A
qQdscfDl56RJWlI63xY7wvaeTBzC6N8WgxpdvxZ2k7/KbGGw+D/FY95M/WuH
xr6hyWWTyd9Zdrh8uezvuynkX7CDw6BJ1pFMPcy2g3fR+kCTvd0kv11PfnAx
HPeE8q/EDuK0sBKHcLK31A75nL/9mmsJ3zI7BNl9y37g2Ejn2+Pyjf5/zSMI
Hyx73N6+XblyKOU/2GP92nPHh86l89g96wPF94cmE39ce0zZ31aw7ibhV2KP
G9s3pp4fQ/iW2mOYTZtttOuAXlpmj8QHzcU3PlJ+l3MAp8vQUX004YflgCKf
lTN2XaF6GQ7wPbh0D9+A5GE74Jduif5FY+Kf6wDHTv3TN2xI3xIHeHmbX9xX
QP4rdcDJP7nv1pUx8c8BB/Wb9paaMfW5Iz77h79pGEH1CssRoeJhSJjIxB9H
GKtVzomsJv7YjvioMcPv5nJG/46ovp/gPi2U9CVxhPKWVRmr5zL50hGZ8Ve1
povz6XxHZPW/U/H2VBn5vxMOXJffFutI9TnLCQrvOgbzxhENJ3wb/vN8O5t+
z3aC4YuOgd/XZ9H5Tsgb86u/oRbVrxIn1F4+3MwBI78TUhPD2+5ZUDyXOeFY
8NLOY4OY+H8Qy0PU2iojaT/rIM7mxxi9W0X1Aw4ib07ZTatmqjfYB2FlneSj
lkb5mnsQTZ5GmkPLvxP+D6JI4eGoqoPthP+DULN6/0x/IOV72UF8Mc7Od/zQ
QPI744ml6uWgsRRPWM74M8jgidCZ6gs4I4ztvr3tKtXHbGdEPh00n2dD/HOd
cahkldkPIelb4oyNKzcMMnAgfEud0W6meNstnOKJzBlx6UF5fU9V0vmH0HS0
7+WXg4kf1iFUlKStCU6h+h2H0Gxm971xD/Uz7EPYLvJPjdV8QPIfglTWd9Lo
tdV0/iGs/5dkZmBK8Ul6CNdMBt1PukP+LDuE6EHKkbkvGfxz4HpUv9hjPOVv
JQ7CF/fvEImYeoCDiliVK7XDCd+aPbTBgzOjU5j8zMGk55sX2gdRfDTkYObL
bdtleqQ/NgdNGx0FqyWVvTSHg3/D6gKM5Jh+joN475VKK/WI3xgOfm3dcTm5
hPiVcNBeYbk5VUL4yOBAOuvQmk1J8uiVj4OsrYqZKv8Rnks56Awv7vp1ktZl
HNz8My/HKor2t3AwpFXB4Xkk0w+7wLxvjfqZcSSPkgs+KJb/mSijeozlAklA
se6MrGKS3wVypbPXjjj2jOR3wfj21JORga9Ifhd8N+C9t0wgf2C7YGBX0rL7
Dwn/HBdMVFKN+zuO6iOuC6w9AteZnaX8H+OC0D3Jp4vYFB8kLvCc+exUnxVU
L2S4YMe1yH7PJjD4dsGdc1o2244+6KVLXaA3N6plUfIr8ncX1JX/alN5S/mu
xQVh68Jf+JyjfCjnir7To9cHFFB9rOQK0Up/vaDARvJHV7gb2A46ZfOil9Z0
xZ+dNuHmWudJfle0JF6xtFUk+xq6wuvmupaP10kf7J7vhc/PSltTRvK74vmj
VUU6b6n/57pCZ6Wnup8K5e8YV1gpebSfHdxGeHaF3WC7wtMjfvbSGa5of3Kk
6LdcN+HbFYq/tHPGdfzppUtdEeuyxFOlm/p5mSteKO+OyR9G32txxemdMzI+
pzD9uxv4K8q1LT7nkvxu8N3SnSg8/5js7wZlyVTDZ6E0X9J0w4zl++XeOD4k
+d0wZKGD3oLoh737Dd3QKjOYdC6I+nO2G561bfo7XJP0y3HD4CEll28+onjF
dYPTuvfdI70oXsW4Ydh8Db+vmhSvJG442Bbm5lND8SHDDUGyL5z/FDJJfjes
OrG4YFc7xZtSN5i6uvUblM3EWzecWv9ZYFFN+b/FDbltd50e2jP1pzscRatO
ntWhekXJHbFl2jqnCinesdyhsy50Z6gHI7879qqcmX979HWKT+44ubBs/z1X
yg+G7jDILwyyW1ZC8rtDKfG6XYxFDNnfHXz2mqrsf7Xk/+74Ne0nyrqI3xh3
yLsOGR6uwODfHdFdrurvxlC8ynCH+pLA4AAO5T+pO46M3VeecJvwU+qOeXbG
/16Hknyynu/vjjIZe5zw0uKOZWkhnrGTS0h+D6StrfdbYkH9o5IH7pdrF9qZ
0byE5YHdNq8sr9pSv6XpAeOK6v0quyj+wgOjnIvf355aSP7vgdbuQW8Tq96R
/3vg6s01d2blU33D8UC9zuHCynKmv/XAuclDy/oz6zEeUJ/rObJrI8VHiQcU
I3Ln+umeIfk9sPHn6+UHJpN/ST1gvdxn9/gG8pfSnu935AYPm8b4vwe+fFy4
9NOlX+T/HggKW+A+Rp/wJOeJFe1lszXSyf+UPLF/wkXHfXdpvsfyRIywse6L
J9UDmp5wfRBZq3yQqU89cWBPRtxlq6dkf0+sC9mQ9HYC6Y/tiWz/F7OVJzaT
/3vCvlugtLGYzuN64hQvwVdxPuE5xhOxXQdPL6rvI+2V3xOzU1yWbArt20tn
eGLc2PFKTmpESz3RoquiElgn10uXemK1aMSuSYdIXpknLt2p0oU6g39PzMw9
5znyCs0H5LyQalLodfAM2VfJCx+qd1RV2NM8g+UFzpLUK9f+vCb5vbAm+UJj
njvFU3ghZsa2C2c6CM+GXrjuOCXu4iA6n+2FZ2ryfdUuk3wcL4zf43s2P51o
rhfcFpxbNi6b7BHjhVetL7nnvZj45wWtCvspI8+SvTN6+Ik9pp++S0T294JC
TkUAj0V4KvWC6sD2m36VVI/KvBD2fWtD+lWiW7zwsZAlKT1Lv5fzhvz4+DmT
xlJ9ouSN4z/NTi+Jp3jP8sa8mDrRbTWqRzS9Ucwb4Fu18RPJ7w33el/vSX3I
vobe2PduSYl8F+mL7Q1ZfqHjtbLjJL83rrVtUJG/RfmV643MeRVzlJ0Z/Hvj
5tqdGdvzKF5JvPH6YERljSeT/7xhdLsoWG8a0195Y/SNHWZnR5J8pd74/e3S
9Vz3XJLfGyNnqM77pUj2bfHGwlHzPM1taV4q54Px2yr3Js+kelLJB/POwr9k
O+P/Ppi5ZY1x4VvKF5o+KJAa7n7/tYjk94Hmm66/96ooHxj6YHT0i6aEl6Rf
tg++pow4c3sezV84PfsDW289mEL9LNcHkk9CR68rEWR/HzzcHdiRrkj2lvhA
1lp53+Q0+UOGD0Y93auQOOI32d8HFwNSlg41ZfzfB1JFtaQluoQfmQ/anF6k
bF1C6y0+4Nc/cmftpv1yvkjRfqpq9JrBvy+SMkUOJx5T/8zyxYYTXpdZ+xn7
++JfnN7klHrCD3xRO/r51aGtpC9DX/wsLRKsY/0g+/tCq2142UIvxv998dHh
De9eoXwvP1xfrGwyf+Bd3r+XjvGFhsrP65oTBvbSEl/EKPz4eXHPgF46wxf7
jWQnPHIZ/PviYr7/Eb+qfxT/evhz6ZwwRUjze5kvRLMmnstj5v8tvsiNSFr2
I4axvx9eJI/4IjlO9Y6SH84O4r1rqmbk98Ps93drGi+QfjT9YDyl22IYoz/4
oWP6j7XPhKRfQz9Eq0tM3BVIPrYfqjtPPbiSRvbg+OH6a23Vpot/Cf9+OK1+
N77rHd1HxPjB5aLWQLNhhC+JH6xv+Ux3HUrzmAw/mG15q5Dim0f+74eZuxLH
+PKp3i71Q7Fk/Obt3+g+ReYHORvb2sl1Fwn/fnjfNGpBejDTf/jD+MLVI4Jq
Rn5/WMaavGp7S/pi+SN/1qDaC/nkf5r+yGN1PdrZROvwRwInNGjPcNpv6I9l
7rHKOZtpPsH2x6DJQ9Lb1lD9wPGHbuDLh4K5TP/sj6/hdXJ6q5n63x8B+weP
SfEgWuKPx6vr92zLot9n+GNg0dXpCdkkn9Qfi+a4rtAXPSP7+yPZv3TiBaZ/
kvnjv3UjBxm2Mvb3x2urgxNcVCleywXg46rI9q0+lD+VAtA+eETJGg7FM1YA
DLMmTJCNoHymGQA5a+uOtuKlvecjAFE7zAQNf2j+ZRiAPxppSwU7Kgj/AdDu
HGazQi+H5A+A8TLjPk+/kD25AfiSuSl3tBPZOyYAuz1Nrt9SJ/+VBCDOrFaw
fxr5d0YALB4YPzvG+Je0h3+1xk07IokuDUB05fsdT+Yy+A/A0htVbRllVB+3
9PB7uqZGX5O5/wpE/MMlU5ZZkn6UAsHvr7Zo0GLKb6xAHC0Y+d+CYqI1A+EZ
vb1eyZn0hUBUDZ+5vciK6g3DQEgeS0+r7SM8sAOhz75r8XZdv16aE4hEOZ2V
m/8N6KW5gSgUVyrWvSe8xwTiRtPKNkEerUsCYTx5ucnUaib/B+JwfUPi0yuk
D2kgONX+juFmpK/SQJg6/ci9EsfM+wJRO+bfX9tHZM+WQCwVujx98JZouSCk
+25Nb+Uy8gdhzfovrXsOkX5YQZif+PyCZTfpTzMI3UPPVh+bTPpFEI5879vs
u4n0bxiELQM7LtaKaJ0dhJ0ai2+PuMXEvyD85fzc2ceF4iM3CN8UQhKzCgk/
MUGwnbYx4N6GArJ/EFxSquP2tD8l/w+ClUTRWjuY8oc0CIeMn/dxzSE8lAYB
it2rOptPEv6DcKlOo0wrjOaHLUFQe15RP9GSiX9cCENsOhROk/wKXBy7xHK/
WczggYsqC8n5HHvKj6pcTP+r13HBk/ITi4vNEdopdt7Xeml1Lhbkm24tf075
T7Nnv3i33Ylqmpdoc+GUnWZQ5kL5G1yk/t4Z2PSd6mF9LrS8t7/cNJvm54Zc
rPjTvPYen/pLUy4K5+8duDX5McUXLswWaPWpaqP+yI6L2QLPpk+txC+Hi9J+
Ziul70g+by7k2ate8WZSfOdysefe8pJEpr4ScRHomtkWbPyI7MGFSLHlT8Ur
6vfjuZCcVFvivInm5xIuTlbWVBfYUP+VysXiYyclWw2pv8jg4oLzyLPH9j/p
/V4mFxPU1tnPiyD9S7l4cmnWj5Rg8ocCLmqtb2/U+cX0s1xMTdirmh1PeKri
Ys2Ylrvpd4mW9Zyfsklj4XDKJ01ceL0//x9rMuGxhYtJYVkZO2zJnl1cGF0/
K6g4T/4mx8N6j7ONafIkvwIPD7IGHY65RfWOUs/6vmFHQph+QJWHLTasy+XH
iV8WDxrvBg04/578W50HxZKf1QcvEh40eVg+5sms96Mov2vzcOrOo4mfjxMN
Hr6q/H7QrUzxQZ+H+glXDIrqmXzKg7Jep/UlP6qPTXlot0/deXwhycPmQe+o
TJgVRfWIHQ+Lol9oXram+pPDg7Nw84skW+pfvHnov60kpY85c1/Ow60Zlg0P
11J+E/FgdCMu2u486S+Gh/3RfxaJnEm+eB7sMoa3zlem+CPhgZPEzT1XQvZK
5WHFDdFQnWcUPzJ44H+OlEvXIv1l8rBguaXazufUL0t50K5q7NchpfxRwMOb
lrSkPSvpvUApDw3GA48WH6H+tYqHhyU1rw6MZep7HuQervVcVE35v4kH7uWO
zSrZ58j+PGDCt1cDX9M8qIsHseSekd0ipv8JRmvKuV/f51G+VAhGzurX4j5z
ad6qFAw5BXel40J636DaQ3d2GVVpEb5ZwRAI37yKuET1u3owRDe/akel0zxa
MxhNG+b5LkymeZd2MCZPXWIoL2Dm68HQGbctIvYczcf0g3HXSss2wYLmD4bB
yPNu/Ciyo3mMaTDE53d7dr6g/pMdjE8Z4r335tN81C4YM+v0NlsnkD9zgrHo
uwMPVpS/vYORqtfw5qQn6YcbjPNZf54etqZ4KgrGZrPEDeajSN8xwfDkx/xS
LCP54oPBz3iV/t8soiXBeGvxweTqGbJHajAKTQYsXX2LidfBWHZw/zWZO8Xf
zGBcy8z6kZvD3E8FY2lbxo03zxn8B2P9uRCvbhPy99Jg2Hx+w/r0nvyxKhj/
GsaqtDVRPSrr4edw3prR2+h7TcG4Y/h2xbZ8sl9LML5+1D1qcYjiWVcwOCNC
V86UxJP9Q1D7/N2Ff6vvkf1DsOFvc2WGMd33K4VA+2bM9IXnGPyH4HP/ypOz
jjD5MQT8ibGxHf8x+A9B+OHk3MUfCb+aIai0G7lllRLlQ+0QiD3/jr/hSPEK
IQhd46raHEH5Wz8Ep8quy9rVKH4bhuBKS81S+0fEj2kI4ubNzhFz0sn+Icj7
KT7UvJ3uI+xCcImX/OHL1qRemhMCh9s6SeciaR7lHQLlIuV7CnsoX3BDcLEq
RfUVU/+IQvAjY1eO7jVmHheCyG+mQpN0sk98CDJKB3Bn32L68xAsLzr559Ny
iiepIbAp6759t4Dyd0YIwvKPL+bNp3loZghGBLw0r86meY80BIP8nHZP1KX+
tiAEU78ujAt9SPgpDYG13MGVl9ypH64Kwd+Hd3CRT/fNsh59qwVWjVpL/XRT
CEZ/faN8w4LuA1tCsOWH/8MbHYTfrhCcnq5nNFyX8rUcH/K8Ci35tXQ/pMDH
6UktFXO9aX6lxMeOkuGXvZcSPlX5eHd8v8cwI+a+mo9d+hMWt6rS/ZQ6H9rs
o5dSdOn+SpMPlZgPdwPNaV2bj81x9hO5zHsD8OH+z9jjxUmST5+PjBMxVkcq
Gfzz4ZBUN0Ioof7BlI/SE9oPbBOpn2Hzod59ZK1CcRrZn48BYRz16auZeTMf
hzYGDXyXRPHUm4/D7TvLjJWZ9xl8PN++tSrzOt1viPjYMm+S+ugQmr/H8FHw
yLl4K5/iVzwfx+xGCpeeZeYTfGwt0Hnqy6X+JJWP9i03oi6LaB6fwUeom+0t
ewuKP5l82N0PmGl1nOpPKR/ZMbPqTu2g+XABH1abkv8Omk54Lu2Rz9tlXeoC
8tcqPu7UXRvqrULzBxkf3Aui8yuSKT429cj7QuXTiEXkDy183J+LwMifNE/s
4mN2+tYxt48x9/8COIkT4qxe3CP7C9ARoadgHEvxSkmAJTd1wn7NZ/AvAMdt
3n+yp4QPlgDOJcf9iqwpPqkLELPjZ/SpvkRrCqDuIxScuUn9lrYApxrypq79
RHiBAKOSz4fqOhE+9QWYdED5wYJIqp8MBagz9XNxfk78mwpgkXkzJuo6cx8n
gGu7h9aWkeSvdgJ8k25rmill+k0BVpwzm5uic5fwL8C5scWCex1UH3AFGPvk
3ChFLaqHRAJIgt/eWPWVqccFyO4XKHqbRfVDvACf47aO0j5MeJIIUJwXdmL2
I9J3qgD9Fm4W7fMi/8oQoHOA5XF7c/KnTAF+mT+Wa1Kl+zapAPKRK/p5NxNd
IED863cp0m7CR6kAt9apBL12JP+sEqCq5bHCqxfUH8sE2JPv4TH1Ofl3kwCp
ExfwjRdRPdoi6MEr+8OZc3S/2iWARmtjQlcs1ddyQtSOyy8/pUzxR0GIQY1o
3RZP83YlIX4ob+jjriWPXvsLcU087afNq369NEuIeNkbg7qzdP+uLkSqS8e+
6AO0rinE7qtnh9QcpHigLcTT33Yf8ZbiB4Son3GyJDGHyf9CqCwVzpS/RvnX
UIgv8e5VGYsY+wtx+wbfU6THzNeFMFssn3BuMfUXdkL0042JmlxC+OQIEXLj
eEGSF53nLYTTKvE28wvM+0AhRq1e3fK2meKTSIgPJ9ePnjyU6BghHuamuOyN
pf3xQmSIxipDj3nvIURS4su8pGs030sVQr8m45RTw2GyvxC84BOGkx7TPDVT
iPTCrOfXyyk/S4VY9iC0onYWxYcCIbZUr1J5LrtF+Bfi9bjjGg95NE+oEmKq
hYtWQCf5v0yINRpti08kkz80CdEd4FpTV0b9VEuPPCGe/3JGUL7oEmLVgGHj
6srJX+VCwX2VOH3wFHfK/6HY5/XV+7Qr+btSKIo2PL6lMY3ikWoowqWYFXWC
6m9WKOqvVP8LP0H5Uz0UXksOWY08SvjRDEXBgxGja4dSvasdirH3aixu1FO9
i1DEDts3VeJG8VM/FOlfGmI0jhH/hqHIu8EycykmPJmG4vTq+VH1aVR/sUMx
YIve5bBh5P92oRj5be+FpzLq1zihmPTjjlOMEt3XeIdi3GEjwaSlFN+4ocj9
06lSaEb2EYViy+k4NeV7pP+YnvOG/Tbsv4Tqw/hQsANfTyjNZ+q/UKScMdrR
coXyQ2oo+JMO/b0wkPwnIxT2j/WDcxTle+nMHv2eCrs6ehTR0lCsfGRUo76V
fl8Qirv633L3pjH3WaHo34fvMXMC4bOqR76fQ4JcM5n3ZKHwz7/2wruD8NLU
w3/b2donpwkPLaGYOzat3wMp+WtXKAL3Lrn3ZSN9X06EgvQ791X2ET4VROjq
3zl4xCjCs5IIlU3v1zQeoPc0qiIk2J+pWVtMNEuEisrCmMgF9L5GXYS0GbdZ
F9tIPk0RHsuH/J11ge7rtEVgn3Q5Mfkk875HhDNaT/YOmEb+qy/CVYePT7Wu
M/2/CGEr8o7+riR5TEWQGWQpjt9A8rBFePK2a39BAunfToQPDkOv6QwlfHNE
eO8+tbXpNcnnLcKVUo3jSfHEH1eEOtXkuL9mtC4SwdYtr2yNkOJfjAgP719d
NKqI4km8CC1dubGDjpO/SXroN2VNI3fSfC9VhAG1NRMrjQt7/SdDhMXm1lkL
m+k+KVOETdn2Rsd2U30oFeFr7V3RsP50X1HQI49wzJa/FhRPSkUYvjeku8GD
7F0lwvSkayp7fZn3RCIole2IRS7RTSIUffIpvcon/2zp+b1OsqrORvpelwj8
jYEatnfiyP5h8OrW9P1UTflQIQz7/0RoVu2neaRSGPzvPFmTO4fwrRqGk9c5
jevmMPPxMGg9kpWdSaJ+Sz0Mpgvr2SMO0HtZzTCIb9WsnTCK3stqh+HDgKdX
FS4w72XDYBazvZQ/gIn/YRhnGZSdzSPaMAxbH4wJMc8heUzDsCjv9CHn9Uz8
D4N748wrRgVU/9iFod8qn3efo6nf4YTB87JVl04a4ds7DDPU9z1x5NN7Rm4Y
WEaV8f2koYT/MGxr7O4e40/5MyYMJnEj//C6KR7Fh2HevD7Zud3M+4wwNLU+
ydCwJvymhmHLqpeTTwwk/8oIg/zhwqbJzHvZzDAsydvaOGwp+ac0DH891Fo0
W8h/C8KQsGfJcnc2yV8aBrellUEplfSeqioMlxWfbP48jOKlLAwD+szxflJE
601hyNy76oVWPu1v6eFHc0PfV4GE964wnPgxwHiOLfP/W8KhY7BHY3g95WuF
cKxICenjWkd4VgrHj9rb9flviFYNxyLftGjDD0z+DwevOvpDoQ19Tz0ckxxi
7r7cQ+dphkNTXJ094SDxox2O7PSKx41iwjfCUV6klr5tPvXP+uF4LeIFfjFj
8n84jjuPGxP7lvRvGg7Lzn1Kf9fS/IAdjjmvd20UulF8sQtHRO1X1pRnZA9O
OML3W45oncngv0e+9zOOGg5j8n84nnk2yjWA9C8Kh2rAiiPyJsx9ZDjmDa7d
NzyU5ivxPfu1uy9vtqZ8IglHos+dYpsrdb10ag9/2ukeTkNpHpHR8/2runuN
rJn+PxwszTfP0qcnEv7Dwf5hEVvpQvmlIBxPOD1N+RGyb2k4NhUVcD6+YPAf
jpThnZJ7X5n+r0fe08Ynmg4z9V84lrwsC8p1fkL2D8diLyUb/kK67+kKh9V+
jv4p5v5aLgJz35SXWBRQP6sQgbWjYoM4V5j7wAhkzG1vDgqhddUITFV+ub/4
JHM/EgH7D+xGr0zCk3oEJKnN360H0PxDMwJ+guHj1l4hvGpHIPdE3/pLmcz7
3gj8u2myib2CaP0IDAneMXxzNvOeLAK7uzP1JqnS90wj8Pfbsq4tErrvYUfg
g/+9usRbzPwnAhMKY/vYp9I8ixOB+XHjD+s30jzTOwKbkhb8u9zK3L9EwPzq
fa0p/aheEEXg2jE9tw3+5J8xEdA/6t3qZU/5Kj4CXVwvbn4i1eeSCNwcMkV7
8XSq31MjsHhA7NecLOa9SgRC9OKiLwwk/8rs4W9d4eq5gcz9dQT4LV36B/6R
/xdE4Dx/jFl9LL3nK42AIa9gZ9Vspv6LgLdzh5bKe+rfZBFotdCQawxj8B+B
GtWNJSkzid+WCFwfdMMufzPhoysC90VpUwc1MPgXo8Tm1NoOIb2fUxCjtHNM
mlSNaCUxVBVHf4moIjypitHS0t9iEZj+X4w0VrHs9lWSR12MmX6C2Tpfyd6a
Ynw7ve5J2yi6n9AWI9t0mJuyHvknxCiefeTodFWa7+uLwa82ueSdTP5sKMav
Baw+lvcY/IshW8YaW3aG+iG2GP1Kh4UtW0r82Imhdudo3wKG5ohxPIL3s3Q1
/d5bjFSf+gadYsITVwzj29PfX9UkfIvEcB9penO3C/WjMWKsrVVf+mUV1bPx
Yqyc1vBj0zDKhxIxvsQs+HonguhUMeRX7brhlEX4yhCjquLD/fQLlF8zxTjC
l1+Ra0Tvp6RijPAY7bBzPzP/FWNARMM7J+Y9YakY7Ss3ThE3kn6qxBD9uOcZ
LD5H9hdDL/HsHmVTmo81ibHTvPRbcinNz1rE2Ku6b+6fsTRf6+rR3/b9PH11
mlfLRWLt89Ef2E40v1OIhOZsD479GXq/ohQJ8dnVQ14GUjxQjcS1D55zJrZR
/86KxIdnK1alvqP6Xj0SySYp9dwdNI/TjMRXq5y89mHUj2tHYsMtmz93HtN7
KUTC/0THKPYPmsfqR6LdOdh35Bp6j28YiePLuTqr3akeMY1E2+RDic6PKB6x
I6Forikrukj82UUiIcfpqyiF1jmRqHxU/1T3KcUr70h8OylXpP2e+OVGwvbQ
7PcqoeSfokjALt7u4iHCU0wksg+f5vWroPwVH4ldl6If1l6gelwSicZaO626
RFpPjeyx3/KzIY20PyMSz1v11oweTPVPZiS0Jy+Ztn8q3c9JIyGffSPpfhbN
cwoicT819M7C18Rfac/5QVOUAzqoH6mKhLp2RfIhLeZ+v4f/VSHTHl6heNXU
Y5936c46z4m/lkh0FmzTdn/J5P9ImFg5LzrgTvFALgq74j7kdfaheKUQhXu/
7Du/FRFelKLg1/3ZoHk45UPVKHw9r/t4cx3FJ1YUMmYYN88W0TxWPQqmujYn
Pzkz70Wj8O35SX6aCd0/aEfBaPSt3McgfCEKW4ZYterdo//vpB+Fu15q8z8W
Un9qGIVA/VVza6cx9X8UttWlbuCeYN7jR+HAeQXuo7fUP9hFYffNsyETrlB/
zInC9EWuK36tI3/z7jn/z8O/5Z+Y/i8K3ycHG1w8THgRRUEgu2S1OIR5jxGF
qEucYM5WwkN8FLZeTXi+i0W0JArSFWpbjRzo96lRKPpZ7s9zo3iQEQX1d8tK
OlVpfpcZhVHbm2rGHaJ+QBqF9J0NjokadD9YEAX21tZNs25TPi2NQtYBjQE/
Qpn5XxRa5w1dP2Uk+bcsCpZzvzf4niJ+mqIQc3MMv+I3894nCgYRBpdGL6X7
864ojHMuC/wSSLRcNE6sXqp3O5/m8wrRuNM293JWO8UHpWhsVwoo/PCYwX80
vJ/Usnc0M/k/GpO+urRnzaT6Rz0a898acTJGk3ya0bjw4beTthpj/2hUzHy3
8yOfed8ejTHlh6s7ntE8QT8adY5JEvamrpz/ARysT5U=
"]]}}, {}},
AspectRatio->0.6180339887498948,
Axes->True,
AxesLabel->{None, None},
AxesOrigin->{0, 0.},
ImageSize->Medium,
Method->{},
PlotRange->{{0, 1001.}, {-2.98723801101878, 2.9757153662914786`}},
PlotRangeClipping->True,
PlotRangePadding->{{20.02, 20.02}, {0.11925906754620517`,
    0.11925906754620517`}}]\)

In[23]:= DeleteDirectory[targetDir, DeleteContents -> True]

In[24]:= If[output === $Failed \[And] outputLines === {}, {output, outputLines,
  targetDir}, Null]

Out[24]= {$Failed, {}, "C:\\Documents and Settings\\HP_Owner\\Local \
Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8"}
POSTED BY: Paul T
Paul T
Paul T
Posted 11 years ago
Here is my code minus the graphical outputs:
  In[1]:= targetDir = CreateDirectory[]
  
  Out[1]= "C:\\Documents and Settings\\HP_Owner\\Local \
  Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8"
  
  In[2]:= lopass = Compile[{{x, _Real, 1}, dt, RC},
     Module[{a = dt/(RC + dt), yprev = First[x], yi},
      Table[yi = a*x[[i]] + (1 - a)*yprev;
       yprev = yi;
      yi, {i, 1, Length[x]}
      ]
     ]
    ];
 
 In[3]:= (* build the example in a specified directory *)
 fnSource = FileNameJoin[{targetDir, "lopass.c"}];
 Export[fnSource, lopass];
 FileNames["*", targetDir]
 
 Out[5]= {"C:\\Documents and Settings\\HP_Owner\\Local \
 Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8\\lopass.c", \
 "C:\\Documents and Settings\\HP_Owner\\Local \
 Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8\\lopass.h"}
 
 In[6]:= (* main function C source *)
 lopassmainSrc = "
   #include \"stdio.h\"
   #include \"stdlib.h\"
   #include \"lopass.h\"
   #include \"WolframRTL.h\"
  
   static WolframLibraryData libData = 0;
  
   int main()
   {
     int err = 0;
     mint i, type, rank, nelems, *dims;
     double *data;
     MTensor x, y;
     double dt;
     double RC;
     libData = WolframLibraryData_new(WolframLibraryVersion);
  
     /* read x */
     type = MType_Real;
     rank = 1;
     dims = (mint*)malloc(rank * sizeof(mint));
     scanf(\" %d\", &nelems);
     dims[0] = nelems;
  
     err = (*(libData->MTensor_new))(type, rank, dims, &x);
     if (err) return 1;
  
     free(dims);
  
     data = (*(libData->MTensor_getRealData))(x);
     for(i = 0; i < nelems; i++) {
       scanf(\" %lf\", &(data[i]));
     }
  
     /* read dt */
     scanf(\" %lf\", &dt);
     /* read RC */
     scanf(\" %lf\", &RC);
  
     err = Initialize_lopass(libData);
  
     y = 0;
     err = lopass(libData, x, dt, RC, &y);
     printf(\"%d\\n\", err);
     if(0 == err){
       dims = (mint*)libData->MTensor_getDimensions(y);
       nelems = dims[0];
       data = (*(libData->MTensor_getRealData))(y);
       printf(\"%d\\n\", nelems);
       for(i = 0; i < nelems; i++)
         printf(\"%f\\n\", data[i]);
     }
  
     Uninitialize_lopass(libData);
     return 0;
   }
   ";
 
 In[7]:= (* create the main source file *)
 lopassmainSrcFile = FileNameJoin[{targetDir, "lopassMain.c"}]
 Export[lopassmainSrcFile, lopassmainSrc, "Text"];
 
 Out[7]= "C:\\Documents and Settings\\HP_Owner\\Local \
 Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8\\lopassMain.c"
 
 In[9]:= Needs["CCompilerDriver`"];
 lopassExe = CreateExecutable[{fnSource, lopassmainSrcFile}, "lowpass",
   "SystemLibraries" -> {}, "TargetDirectory" -> targetDir,
   "Libraries" -> "WolframRTL_Static_Minimal"]
 
 Out[10]= "C:\\Documents and Settings\\HP_Owner\\Local \
 Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8\\lowpass.exe"
 
In[11]:= dt = 0.01;
RC = 0.3;
input = Table[Sin[x] + Sin[x*10] + Sin[x*50], {x, 0, 10, dt}];
ListPlot[input, ImageSize -> Medium, Joined -> True]

In[15]:= inputText =
  StringJoin[ToString[Length[input]], "\n", Riffle[ToString /@ input, "\n"],
   ToString[dt], ToString[RC]];
inputFile = FileNameJoin[{targetDir, "input.txt"}]
Export[inputFile, inputText, "Text"];

Out[16]= "C:\\Documents and Settings\\HP_Owner\\Local \
Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8\\input.txt"

In[18]:= SetDirectory[targetDir]; outputLines =
Import["!" <> FileNameTake[lopassExe, -1] <> " < " <> inputFile, "Lines"];
ResetDirectory[];

In[20]:= output = ToExpression /@ outputLines;
output =
  If[MatchQ[output, {0, x_ /; x + 2 === Length[output], ___}],
   Drop[output, 2],
   $Failed
   ];

In[22]:= ListPlot[{input, output}, ImageSize -> Medium, Joined -> True]

In[23]:= DeleteDirectory[targetDir, DeleteContents -> True]

In[24]:= If[output === $Failed \[And] outputLines === {}, {output, outputLines,
  targetDir}, Null]

Out[24]= {$Failed, {}, "C:\\Documents and Settings\\HP_Owner\\Local \
Settings\\Temp\\m-7dc3e81d-6370-48ba-93f4-809adf3127a8"}
One last thing, I do get the .exe, both of the .c's, the input file, and the .h in the Temporary directory. I think I'm getting all of the stuff I'm supposed to get, but for some reason the Import can't pick it up, and I don't know why.
POSTED BY: Paul T
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Tag limit exceeded
Note: Only the first five people you tag will receive an email notification; the other tagged names will appear as links to their profiles.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract