** This is part of the CTSim program
** Copyright (C) 1983-2000 Kevin Rosenberg
**
-** $Id: ctsimtext.cpp,v 1.26 2003/04/01 18:56:59 kevin Exp $
+** $Id$
**
** This program is free software; you can redistribute it and/or modify
** it under the terms of the GNU General Public License (version 2) as
// If linked to ctsimtext, but executed as another name, eg pjrec, then program will use that
// linked name as name of function.
-static const char* const g_szIdStr = "$Id: ctsimtext.cpp,v 1.26 2003/04/01 18:56:59 kevin Exp $";
+static const char* const g_szIdStr = "$Id$";
static const char* const s_szProgramName = "ctsimtext";
static const char* const s_szProgramName2 = "ctsimtext.exe";
static const char* const s_szProgramName3 = "ctsimtext-lam";
static void convertStringToArgcv (char* szLine, int* piArgc, char*** pppArgv);
-void
+void
ctsimtext_usage (const char *program)
{
std::cout << "usage: " << fileBasename(program) << " ctsim-function-name ctstim-function-parameters...\n";
std::cout << " phm2if Convert a geometric phantom into an imagefile\n";
std::cout << " phm2pj Take projections of a phantom object\n";
std::cout << " phm2helix Take projections of a phantom object\n";
- std::cout << " pjHinterp Interpolate helical projections of a phantom object\n";
+ std::cout << " pjHinterp Interpolate helical projections of a phantom object\n";
std::cout << " linogram Print linogram sampling\n";
}
-void
+void
interactive_usage ()
{
std::cout << "usage: function-name parameters...\n";
static bool s_bInteractive = false;
-int
+int
ctsimtext_main (int argc, char * argv[])
{
- int iReturn = 0;
-
+ int iReturn = 0;
+
if (argc > 1 && (strcmp(s_szProgramName, fileBasename (argv[0])) == 0 || strcmp(s_szProgramName2, fileBasename (argv[0])) == 0 || strcmp(s_szProgramName3, fileBasename (argv[0])) == 0)) {
argv++;
argc--;
std::cout << ", Version " << VERSION;
#endif
std::cout << " (Type \"quit\" to end)\n\n";
-
+
while (1) {
#ifdef HAVE_READLINE
char* pszInputLine = readline (szPrompt);
break;
if (*pszInputLine != EOS)
add_history (pszInputLine);
-
+
#else // DONT_HAVE_READLINE
-
+
static const int s_MaxLineLength = 1024;
char* pszInputLine = new char [s_MaxLineLength+1];
std::cout << szPrompt;
std::cin.getline (pszInputLine, s_MaxLineLength);
-
+
#ifdef DEBUG
std::cout << "#" << pszInputLine << "#\n";
#endif
-
+
std::cout << std::flush;
std::cout << "\n";
#endif // DONT_HAVE_READLINE
-
- if (strncasecmp (pszInputLine, "quit", 4) == 0)
+
+ if (strncasecmp (pszInputLine, "quit", 4) == 0)
break;
-
+
convertStringToArgcv (pszInputLine, &argc, &argv);
#ifdef DEBUG
for (int i = 0; i < argc; i++)
std::cout << "Token " << i << ": " << argv[i] << "\n";
#endif
iReturn = processCommand (argc, argv);
-
+
delete pszInputLine;
}
}
-
+
return iReturn;
}
char* pCurrentPos = pszLine;
int nTokens = 0;
std::vector<char*> vecpszToken;
-
+
// Process line
bool bInDoubleQuote = false;
bool bInSingleQuote = false;
bool bInToken = false;
-
+
while (*pCurrentPos) {
if (isspace (*pCurrentPos)) {
if (! bInToken)
nTokens++;
vecpszToken.push_back (pCurrentPos);
}
-
+
pCurrentPos++;
}
-
+
*piArgc = nTokens;
if (nTokens > 0) {
*pppArgv = new char* [nTokens];
*pppArgv = NULL;
}
-static int
+static int
processCommand (int argc, char* const argv[])
{
if (argc < 1)
return 1;
-
+
const char* const pszFunction = fileBasename (argv[0]);
- try {
+ try {
if (strcasecmp (pszFunction, "if1") == 0)
return if1_main (argc, argv);
else if (strcasecmp (pszFunction, "if2") == 0)
std::cerr << "Exception: " << e.what() << std::endl;
} catch (...) {
std::cerr << "Unknown exception caught\n";
- }
+ }
return 1;
}
-int
+int
main (int argc, char* argv[])
{
- int retval = 1;
-
- retval = ctsimtext_main(argc, argv);
-
+#if defined(HAVE_FFTW) && defined(HAVE_GETENV)
+ const char* const pszHome = getenv("HOME");
+ char* pszWisdom = NULL;
+
+ if (pszHome) {
+ const char szFileBase[] = ".fftw3-wisdom";
+ int nHome = strlen(pszHome);
+ int nBase = strlen(szFileBase);
+ int len = nHome + nBase + 1;
+ pszWisdom = new char [ len + 1 ];
+ strcpy (pszWisdom, pszHome);
+ pszWisdom[nHome] = '/';
+ strcpy(pszWisdom+nHome+1,szFileBase);
+ pszWisdom[nHome+nBase+2] = 0;
+
+ FILE *wisdom = fopen(pszWisdom,"r");
+ if (wisdom) {
+ fftw_import_wisdom_from_file(wisdom);
+ fclose(wisdom);
+ }
+ }
+#endif
+
+ int retval = ctsimtext_main(argc, argv);
+
+#if defined(HAVE_FFTW) && defined(HAVE_GETENV)
+ if (pszWisdom) {
+ FILE* wisdom = fopen(pszWisdom,"w+");
+ if (wisdom) {
+ fftw_export_wisdom_to_file(wisdom);
+ fclose(wisdom);
+ delete [] pszWisdom;
+ }
+ }
+#endif
+
return retval;
}