/*!////////////////////////////////////////////////////////////////////////////////////////////////////////////////*/ //!////////////////////////////////////////////////////////////////////////////////////////////////////////////// /// /*!________________________________________________________________________________________________________________*/ /// mathSuite v7.10 --- by Marco Chiarelli aka DekraN aka Wesker013 /// /// LAST UPDATE: 16:00 of 16/11/2014 by Myself. /// /// This program is protected by Creative Commons CC BY-SA 2.0 License. For more informations contact me. /// /// You can contact me at: marco_chiarelli@yahoo.it or on the secundary mail: /// /// marcochiarelli.nextgenlab@gmail.com in order to report a bug or simply for sending me an advice that could be /// /// useful or could improve the speed or optimize my MSCenv System. /// //!-------------------------------------------------------------------------------------------------------------- /// //!______________________________________________________________________________________________________________ /// //!-------------------------------------------------------------------------------------------------------------- /// /// contact me at: marco_chiarelli@yahoo.it or marcochiarelli.nextgenlab@gmail.com /// /// I'll be glad to fix your scripts or simply to take away your doubts about the program /// //!-------------------------------------------------------------------------------------------------------------- /// //!-------------------------------------------------------------------------------------------------------------- /// /// Thanks to giggikr: http://forum.html.it/forum/showthread/t-1374455.html for his own function, cambiabase, /// /// which I renamed, modified and adapted to this program. Thanks to: /// /// http://elite.polito.it/files/courses/12BHD/progr/Esercizi-C-v2_01.pdf for some of their scripts. /// /// Thanks to Bibek Subedi, for his invertMatrix function, which I renamed, modified and adapted to this program. /// /// Link Source: http://programming-technique.blogspot.it/2011/09/numerical-methods-inverse-of-nxn-matrix.html /// /// Thanks to Paul Bourke:https://www.cs.rochester.edu/u/brown/Crypto/assts/projects/adj.html for his CoFactor fnc/// /// which I modified and adapted to this program. /// /// Thanks to my University friends: Dino Sbarro, Gabriele Accarino and Giampiero d'Autilia /// /// for inspiring me the implementation of some functions and for the everyday support. /// /// Thanks to W. Cochran wcochran@vancouver.wsu.edu for his Strassen Algorithm Implementation, which I renamed, /// /// adapted and modified to this program. Thanks also to: Computer Science Division | EECS at UC Berkeley for /// /// some notions about Matrix Multiplication Optimizations Techniques: www.cs.berkeley.edu/~knight/cs267/hw1.html /// /// Thanks to: http://cap-lore.com/MathPhys/eigen/j.c for the actual _matrixEigenValues function. /// /// Massive thanks to Brian Allen Vanderburg II for his fabulous C parser and inline functions solver, EXPREVAL, /// /// which elegantly gave in theory infinite functionalities and potential to my program. That's the project link /// /// with Online Documentation: http://expreval.sourceforge.net/ Thanks to: http://www.cprogramming.com/tips/ and /// /// http://stackoverflow.com/questions/599365/what-is-your-favorite-c-programming-trick /// /// http://stackoverflow.com/questions/132241/hidden-features-of-c /// /// that are some websites in which I found a lot of useful C tips and tricks, and they were an important /// /// checkpoint for resources retrieving in order to speed-up and optimize my program. Still greatly thanks to /// /// Bibek Subedi for his website: http://www.programming-techniques.com/ which put in front of my eyes a new /// /// world of C programming. I also recently renewed the program code by improving a lot of his C tricks and tips. /// /// For example, the upper-triangular Matrixes conversion, which was useful to enhance some functions like det(), /// /// sgeqsolver ExprEval inline command, etc. /// /// Greatly thanks to Daniel Veillard for his fabulous XML Parser, LIBXML2. /// /// Greatly thanks to vict85 of matematicamente.it Network, for informing me about the benefits of using /// /// generally a single reference for the Matrix Type, like LAPACK and the other Numeric Calculus Environments. /// /// Thanks to Francesco Palma for reporting me some bugs, and finally, massive thanks to my /// /// Informatic Fundaments Teacher, Mario Alessandro Bochicchio, which gave me a lot of C advices and some general /// /// tricks and tips, that enlarged my professional informatic horizonts. That's all... /// /*!________________________________________________________________________________________________________________*/ //!////////////////////////////////////////////////////////////////////////////////////////////////////////////////// /*!////////////////////////////////////////////////////////////////////////////////////////////////////////////////*/ #include "dutils.h" // DA RENDERE VISIBILE SIA AL COMPILATORE CHE AL LINKER #ifdef STACKALLOC struct program suite = { SYSTEM_LOG, INIT_EXPRTYPE, INIT_MATRIXOBJ, INIT_MATRIXOBJ, INIT_LOGOBJ, INIT_LOGOBJ, INIT_LAYOUTOBJ, { { NULL, 0 }, { NULL, 0 }, { NULL, 0 }, { NULL, 0 } }, { { { INIT_VARLIST, INIT_VARLIST }, MAIN_ENV, STARTING_ENVSNO }, { { INIT_MATLIST, INIT_MATLIST }, MAIN_MAT, STARTING_MATNO }, { { INIT_LOGSLIST, INIT_LOGSLIST }, MAIN_LOG, STARTING_LOGSNO }, { { INIT_LAYOUTLIST, INIT_LAYOUTLIST }, MAIN_LAYOUT, STARTING_LAYOUTNO } }, INIT_FUNCLIST, INIT_CONSTLIST, DEFAULT_COLORS_PATH, INITIALIZING_RANDOM_SEED, { INIT_COLOR, INIT_COLOR, INIT_COLOR, INIT_COLOR, INIT_COLOR, INIT_COLOR }, PROGRAM_BUSY, INVALID_EXITHANDLE, false, false }; #else struct program * suite; #endif const struct prog_constants suite_c = { { "Varlist", "Matrix", "Log", "Layout" }, { "Show Programs Descriptions", "Set Item as Current at Creation", "Items Selection By Path", "Lists Items Autosaving", #ifdef WINOS "System Log Security Check", #endif "System INPUT Parser", "Advanced Calculator Parser", "Matrices Items Inserting Parser", "Reset Consts List and Functions List", "Save last Results", "Show Varlist after Operations", "Show Average Time after Operations", "Show Average Time after Advanced Calculator Operations", "Show Average Time after Linear Algebra Operations", "Show Average Time after Programs Executions", "Show Date and Time in Textline", "Print Rows Labels", "Domains Checks", "Matrices Chars Insert", "Matrices Automatic Inserting", "Lazy Execution", "Inverse Operations", "Inverse Operations Automatic Deactivator", "Degrees Entering on Trigonometric Operations", "Request Programs Repetitions", "Strassen Optimization" }, { { BITMASK_SHOWDESCRIPTION, DEFAULT_SHOWDESCRIPTION }, { BITMASK_AUTOSETCURITEM, DEFAULT_AUTOSETCURITEM }, { BITMASK_ITEMSSELECTBYPATH, DEFAULT_ITEMSSELECTBYPATH }, { BITMASK_ITEMSAUTOSAVING, DEFAULT_ITEMSAUTOSAVING }, #ifdef WINOS { BITMASK_SYSLOGSECURITYCHECK, DEFAULT_SYSLOGSECURITYCHECK }, #endif // WINOS { BITMASK_SYSTEMPARSING, DEFAULT_SYSTEMPARSING }, { BITMASK_ADVCALCPARSING, DEFAULT_ADVCALCPARSING }, { BITMASK_MATRIXPARSING, DEFAULT_MATRIXPARSING }, { BITMASK_RESETLISTS, DEFAULT_RESETLISTS }, { BITMASK_SAVERESULTS, DEFAULT_SAVERESULTS }, { BITMASK_SHOWVARLIST, DEFAULT_SHOWVARLIST }, { BITMASK_SHOWDIFFTIME, DEFAULT_SHOWDIFFTIME }, { BITMASK_SHOWDIFFTIMEADVCALC, DEFAULT_SHOWDIFFTIMEADVCALC }, { BITMASK_SHOWDIFFTIMEALGOPS, DEFAULT_SHOWDIFFTIMEALGOPS }, { BITMASK_SHOWEXECTIME, DEFAULT_SHOWEXECTIME }, { BITMASK_PRINTTIME, DEFAULT_PRINTTIME }, { BITMASK_PRINTROWSLABELS, DEFAULT_PRINTROWSLABELS, }, { BITMASK_DOMAINCHECK, DEFAULT_DOMAINCHECK }, { BITMASK_CHARINSERT, DEFAULT_CHARINSERT }, { BITMASK_INSERTMODE, DEFAULT_INSERTMODE }, { BITMASK_LAZYEXECUTION, DEFAULT_LAZYEXECUTION }, { BITMASK_INVERSEOPERATIONS, DEFAULT_INVERSEOPERATIONS }, { BITMASK_AUTOTURNBACK, DEFAULT_AUTOTURNBACK }, { BITMASK_DEGREESENTERING, DEFAULT_DEGREESENTERING }, { BITMASK_PROGREPEATCHECK, DEFAULT_PROGREPEATCHECK }, { BITMASK_STRASSENOPTIMIZATION, DEFAULT_STRASSENOPTIMIZATION } }, { MAX_FIBONACCI_MEMOIZABLE_INDEX, MAX_FATTORIALE_MEMOIZABLE_INDEX, MAX_EVEN_SFATTORIALE_MEMOIZABLE_INDEX, MAX_ODD_SFATTORIALE_MEMOIZABLE_INDEX }, { "Fibonacci", "Factorial", "Even SemiFactorial", "Odd SemiFactorial", "Every" }, { "Real Numbers", "Complex Numbers", "Quaternions", "Octonions", "Sedenions" }, { { REAL_UNIT_NAME }, { REAL_UNIT_NAME, "i" }, { REAL_UNIT_NAME, "i", "j", "k" }, { REAL_UNIT_NAME, "e1", "e2", "e3", "e4", "e5", "e6", "e7" }, { REAL_UNIT_NAME, "e1", "e2", "e3", "e4", "e5", "e6", "e7", "e8", "e9", "e10", "e11", "e12", "e13", "e14", "e15" } }, { "Default COLOR", "Errors COLOR", "Credits and Directives COLOR", "User Executions COLOR", "System COLOR", "Signatures/Decorations COLOR" }, { "BLACK", "BLUE", "GREEN", "AZURE", "DARK RED", "PURPLE", "YELLOW", "WHITE", "GREY", "LUM BLUE", "LUM GREEN", "CRYSTAL", "LUM RED", "LUM PURPLE", "LUM YELLOW", "WHITEPLUS" #ifdef WINOS , UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, "LUM WHITE", UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, "LUM AZURE", UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING, UNKNOWN_STRING #endif } }; sprog main_menu[MAX_PROGRAMMI] = { [MAIN_BCALCULATOR] = { "Basic Calculator", CMD_BCALC, USAGE_BCALC, basicCalculator, BY_USER, CHILD }, [MAIN_ADVANCEDCALCULATOR] = { "Advanced Calculator", CMD_ADVCALC, USAGE_ADVCALC, calcolatoreAvanzato, AUTOMATIC, FATHER }, [MAIN_ALGEBRAOPERATIONS] = { "Linear Algebra Operations", CMD_LINEARALGEBRAOPERATIONS, USAGE_LINEARALGEBRAOPERATIONS, algebraOperations, AUTOMATIC, FATHER }, #ifdef ALLOW_MSSMANAGER [MAIN_MSSMANAGER] = { "PROGRAM Scripts Manager", CMD_MSSMANAGER, USAGE_MSSMANAGER, mssManager, AUTOMATIC, FATHER }, #endif [MAIN_CHANGESETTINGS] = { "Edit PROGRAM Settings", CMD_CHANGESETTINGS, USAGE_CHANGESETTINGS, changeProgramSettings, AUTOMATIC, FATHER } }; int main(int argc, char **argv) { #ifndef STACKALLOC suite = malloc(sizeof(struct program)); errMem(suite, HEAPALLOC_ERROR); // INITIALIZING suite vars strcpy(access(sysLogPath), SYSTEM_LOG); access(exprVars) = INIT_EXPRTYPE; access(curMatrix) = INIT_MATRIXOBJ; access(lmpMatrix) = INIT_MATRIXOBJ; access(curLog) = INIT_LOGOBJ; access(sysLog) = INIT_LOGOBJ; access(curLayout) = INIT_LAYOUTOBJ; #pragma omp parallel for num_threads(MAX_MEMOIZABLE_FUNCTIONS) for(dim_typ i=0; i 1 && ((extension = strrchr(argv[1], '.')+1) == NULL || (!strcmp(extension, DEFAULT_PATHLIST_FILE_EXTENSION)))) strcpy(path, argv[1]); else sprintf(path, "%s.%s", DEFAULT_MAIN_PATH, DEFAULT_PATHLIST_FILE_EXTENSION); dim_typ i; const sel_typ mode = extension ? checkItemTypeByExtension(extension) : MAX_LISTS; _lfLoader(path); if((mode != LAYOUTS && !getItemsListNo(LAYOUTS)) || (mode == LAYOUTS && !file_exists(argv[1]))) { const char bools_identifiers[MAX_DIMENSIONS][SIGN_STRING] = { "false", "true" }; xmlInitParser(); LIBXML_TEST_VERSION; doc = xmlNewDoc(BAD_CAST "1.0"); root_node = xmlNewNode(NULL, BAD_CAST "settings"); xmlDocSetRootElement(doc, root_node); node = xmlNewNode(NULL, BAD_CAST "generalSettings"); xmlNewChild(node, NULL, BAD_CAST "exitChar", BAD_CAST EXIT_CHAR); char tmp_string[INFO_STRING] = NULL_CHAR; sprintf(tmp_string, "%hu", DEFAULT_PRECISION); xmlNewChild(node, NULL, BAD_CAST "programPrecision", BAD_CAST tmp_string); sprintf(tmp_string, "%hu", DEFAULT_STABILIZER_FACTOR); xmlNewChild(node, NULL, BAD_CAST "stabilizerFactor", BAD_CAST tmp_string); sprintf(tmp_string, "%hu", MIN_STIRLING_NUMBER); xmlNewChild(node, NULL, BAD_CAST "minStirlingNumber", BAD_CAST tmp_string); sprintf(tmp_string, "%hu", DEFAULT_ALGEBRA); xmlNewChild(node, NULL, BAD_CAST "algebra", BAD_CAST tmp_string); sprintf(tmp_string, "%.*f", DEFAULT_PRECISION, OUTLIER_CONSTANT); xmlNewChild(node, NULL, BAD_CAST "outlierConstant", BAD_CAST tmp_string); xmlAddChild(root_node, node); xmlFreeNode(node); node = xmlNewNode(NULL, BAD_CAST "matricesOptions"); sprintf(tmp_string, "%hu", MAX_RIGHE); xmlNewChild(node, NULL, BAD_CAST "maxRows", BAD_CAST tmp_string); sprintf(tmp_string,"%hu", MAX_COLONNE); xmlNewChild(node, NULL, BAD_CAST "maxColumns", BAD_CAST tmp_string); sprintf(tmp_string, "%hu", _BLOCK_SIZE); xmlNewChild(node, NULL, BAD_CAST "blockSize", BAD_CAST tmp_string); sprintf(tmp_string,"%hu", STARTING_MINOSMMDIM); xmlNewChild(node, NULL, BAD_CAST "minOSMMDim", BAD_CAST tmp_string); sprintf(tmp_string,"%hu", STARTING_MINSTRASSENDIM); xmlNewChild(node, NULL, BAD_CAST "minStrassenDim", BAD_CAST tmp_string); sprintf(tmp_string, "%hu", MAX_DSVD_ITERATIONS); xmlNewChild(node, NULL, BAD_CAST "maxDSVDIterations", BAD_CAST tmp_string); sprintf(tmp_string, "%hu", MAX_SIMPLEXMETHOD_ITERATIONS); xmlNewChild(node, NULL, BAD_CAST "maxSimplexIterations", BAD_CAST tmp_string); xmlAddChild(root_node, node); xmlFreeNode(node); node = xmlNewNode(NULL, BAD_CAST "memoizerOptions"); #ifdef WINOS #pragma omp parallel for num_threads(MAX_MEMOIZABLE_FUNCTIONS) #endif for(i=0; i 2) _handleCmdLine(argc-2, &argv[2]); progInfo(WITHOUT_DESCRIPTION); updInfo(); operationsGroupMenu(MAX_PROGRAMMI, main_menu, NULL_CHAR, BY_CHARS); prepareToExit(); #ifdef WINOS (void) system("PAUSE"); // if working on Windows Environment... #else printf("\e[0m"); #endif return NOERROR_EXIT; }