/*
* receives a SPICE string from the <inputfile>.* with #<text># and replaces
* it with the correct parameters. Depending on the optimize variable value:
* optimize=0 : we are initializing, just replace fixed variables
* optimize=1 : " " optimizing, replace remaining optimizing variables
*/
void ReplaceSymbol(char *ret, int optimize)
{
int i, k, ii;
char laux[LONGSTRINGSIZE], inlinecomment[LONGSTRINGSIZE] /*, lxp[LONGSTRINGSIZE]*/;
k=inlinestrpos(ret);
inlinecomment[0]='\0';
if (k) {
strsub(inlinecomment, ret, k, (int)strlen(ret)); /*copies the in-line comment*/
ret[k-1]='\0';
}
ii=1;
ReadSubKey(laux, ret, &ii, '#', '#', 5);
while (ii<=(int)strlen(ret)) {
/*Str2Lower(laux);*/
i=0;
while (strcmp(parameters[i].symbol, laux)) {
if (parameters[i].symbol[0] == '\0') { /*if symbol is not found*/
printf("initialize.c - ReplaceSymbol -- Symbol in <inputfile>.* not found in <inputfile>.cfg: %s\n", laux);
exit(EXIT_FAILURE);
}
i++;
}
// strsub(laux, ret, ii+1, (int)strlen(ret)-ii); /*copies the last part of the string to laux*/
// /*ret[strpos2(ret, "#", 1)-1]='\0';*/ /*properly finishes the string*/
// ret[ii-(int)strlen(parameters[i].symbol)-2]='\0'; /*properly finishes the string*/
if (optimize==0) { /*optimize=0 : we are initializing*/
if (parameters[i].optimize==0)
sprintf(laux, "%E", parameters[i].value); /*writes the value*/
else
sprintf(laux, "#%s#",parameters[i].symbol); /*writes the #<symbol># back again*/
} else { /*optimize=1 : we are optimizing*/
sprintf(laux, "%E", parameters[i].value); /*writes the value*/
}
InsertString(ret, laux, (ii-(int)strlen(parameters[i].symbol)-1), ii+1);
// strcat(ret, lxp);
// strcat(ret, laux);
ii++; /*because 'ReadSubKey' starts at [ii-1] instead of [ii]*/
ReadSubKey(laux, ret, &ii, '#', '#', 0);
}
strcat(ret, inlinecomment); /*concatenates the in-line comment*/
}
void SimpleParametersCategory(int num_measures, char *llog, statistics stats, FILE **fOut)
{
int i, j;
char lkk1[LONGSTRINGSIZE], lkk2[LONGSTRINGSIZE], lkk3[LONGSTRINGSIZE];
llog[0]=58; /* ':' */
j=1;
for (i = 1; i <= (num_measures); i=i+2) {
ReadSubKey(lkk1, llog, &j, ':', ':', 0);
StripSpaces(lkk1);
ReadSubKey(lkk2, llog, &j, ':', ':', 0);
StripSpaces(lkk2);
sprintf(lkk3, "%s:#%s#:%s:%E:%E:LIN_DOUBLE:OPT", lkk1, lkk1, lkk2, stats.min[i], stats.max[i]);
fprintf(*fOut, "%s\n", lkk3);
}
}
/*
* receives a string from the extract/ directory with #SYMBOL# and/or #NODE# text
* requiring a measurement extraction and replaces it with the appropriate value
*/
void DecodeSymbolNode(char *ret, int i)
{
int ii;
char laux[LONGSTRINGSIZE];
ii=1;
ReadSubKey(laux, ret, &ii, '#', '#', 0);
while (ii<=(int)strlen(ret)) {
Str2Lower(laux);
if (!strcmp(laux, "symbol"))
ii=ii+1000; /*it is a "symbol" and encode this information by adding '1000'*/
else {
if (strcmp(laux, "node")) { /*if it is not "node" then exit*/
printf("initialize.c - DecodeSymbolNode -- Unrecognized option: %s\n", laux);
exit(EXIT_FAILURE);
}
}
if (ii>1000)
strsub(laux, ret, ii-1000+1, (int)strlen(ret)-(ii-1000)); /*copies the last part of the string to laux*/
else
strsub(laux, ret, ii+1, (int)strlen(ret)-ii); /*copies the last part of the string to laux*/
ret[strpos2(ret, "#", 1)-1]='\0'; /*properly finishes string */
if (ii>1000) {
ii=ii-1000;
strcat(ret, UNIQUECHAR); /* unique sequence added to the symbol*/
strcat(ret, measurements[i].meas_symbol); /* adds the symbol or... */
/* strcat(ret, " "); not necessary anymore (*) */
} else
strcat(ret, measurements[i].node); /* ...adds the node */
strcat(ret, laux);
ii=1;
ReadSubKey(laux, ret, &ii, '#', '#', 0);
}
}
void ComplexParametersCategory(char *llog, statistics stats, FILE **fOut, FILE **fcfg)
{
int i, j;
char lkk1[LONGSTRINGSIZE], lkk2[LONGSTRINGSIZE];
fgets2(lkk, LONGSTRINGSIZE, *fcfg);
while ((lkk[0] != '#') && (lkk[0] != '\0') && (!feof(*fcfg))) {
i=strpos2(lkk, ":OPT", 1);
if (lkk[0]=='*') /*if the line is comment out, simply ignore it*/
i=0;
if (i) { /* line with OPT */
j=1;
ReadSubKey(lkk1, lkk, &j, '#', '#', 5); /* lkk1=symbol */
lkk2[0]='\0';
i=0;
j=1;
while ( (strcmp (lkk1, lkk2)) && (j<strlen(llog) ) ) {
i++;
ReadSubKey(lkk2, llog, &j, ':', ':', 0);
StripSpaces(lkk2);
}
if (j==strlen(llog)) {
printf("auxfunc_log.c - ComplexParametersCategory -- Wrong config opened.\n");
exit(EXIT_FAILURE);
}
ReadSubKey(lkk2, llog, &j, ':', ':', 0); /*the measured value for the previously text*/
i++; /*the correct position of the value*/
j=1;
ReadSubKey(lkk1, lkk, &j, ':', ':', 0);
strsub(lkk1, lkk, 1, j);
sprintf(lkk2, "%s:%E:%E", lkk2, stats.min[i], stats.max[i]);
strcat(lkk1, lkk2);
ReadSubKey(lkk2, lkk, &j, ':', ':', 0);
ReadSubKey(lkk2, lkk, &j, ':', ':', 0);
ReadSubKey(lkk2, lkk, &j, ':', ':', 0);
strsub(lkk2, lkk, j, (int)strlen(lkk));
strcat(lkk1, lkk2);
fprintf(*fOut, "%s\n", lkk1);
} else { /* line with --- */
fprintf(*fOut, "%s\n", lkk);
}
fgets2(lkk, LONGSTRINGSIZE, *fcfg);
}
}
/*
* create the file 'alter.inc' to be included in <inputfile>.* file. Depending
* the value of the 'append' variable, the data is append or a new file is
* created.
*/
void CreateALTERinc(char *ConfigFile, char *OutputFile, int append)
{
int i, j, k, kk, add, jj, l;
char data[LONGSTRINGSIZE];
FILE *fsweepINI, *falterINC;
alter_line alter[ALTERLINES]; /*ALTERLINES in 'sweep.ini' file*/
int sweep[SWEEPLINES]; /*SWEEPLINES variables: [2123] */
int order[SWEEPLINES]; /*SWEEPLINES variables: [1111] [1112] [1113] [1121] ...*/
int ptr, num_lines, index, alter_times;
/**/
/*Step0: zero variables*/
for (i = 0; i <= ALTERLINES-1; i++) {
*alter[i].text = '\0';
for (j = 0; j <= SWEEPVARS-1; j++)
*alter[i].data[j] = '\0';
}
for (i = 0; i <= SWEEPLINES-1; i++) {
sweep[i] = 0;
order[i] = 0;
}
/**/
/*Step1: open input and output files*/
if ((fsweepINI=fopen(ConfigFile,"rt")) == 0) {
printf("auxfunc_alter.c - CreateALTERinc -- Cannot open config file: %s\n", ConfigFile);
exit(EXIT_FAILURE);
}
if (append) {
if ((falterINC=fopen(OutputFile,"r+t")) == 0) { /*append to a file*/
printf("auxfunc_alter.c - CreateALTERinc -- Cannot write to output file: %s\n", OutputFile);
exit(EXIT_FAILURE);
}
} else {
if ((falterINC=fopen(OutputFile,"wt")) == 0) { /*create an empty file for writing*/
printf("auxfunc_alter.c - CreateALTERinc -- Cannot write to output file: %s\n", OutputFile);
exit(EXIT_FAILURE);
}
}
/**/
/**/
ReadKey(lkk, "# ALTER #", fsweepINI);
if (!lkk[0]) {
printf("auxfunc_alter.c - CreateALTERinc -- # ALTER # key not found\n");
exit(EXIT_FAILURE);
} else {
/**/
/*Step2: read the options in the configuration file*/
fgets2(lkk, LONGSTRINGSIZE, fsweepINI);
switch(spice) { /*Spectre special case*/
case 1: /*Eldo*/
case 2: /*HSPICE*/
break;
case 3: /*LTspice*/
printf("auxfunc_alter.c - CreateALTERinc -- LTspice not supported\n");
exit(EXIT_FAILURE);
case 4: /*Spectre*/
if (!strpos2(lkk, " altergroup", 1)) {
printf("auxfunc_alter.c - CreateALTERinc -- 'altergroup' must exist in Spectre mode\n");
exit(EXIT_FAILURE);
}
break;
case 50: /*Qucs*/
printf("auxfunc_alter.c - CreateALTERinc -- Qucs not supported\n");
exit(EXIT_FAILURE);
case 100: /*general*/
printf("auxfunc_alter.c - CreateALTERinc -- GENERAL not supported\n");
exit(EXIT_FAILURE);
default:
printf("auxfunc_alter.c - CreateALTERinc -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
i = 1;
ptr = 0;
while (lkk[0] != '#') {
if (lkk[0] != '*') {
if (i > (ALTERLINES)) {
printf("auxfunc_alter.c - CreateALTERinc -- Maximum number of %d lines reached. Increase ALTERLINES in auxfunc_alter.h\n", ALTERLINES);
exit(EXIT_FAILURE);
}
j = strpos2(lkk, "$", 1); /*This will skip the characters after '$', the inline comment used by the sweep tools*/
if (j != 0)
lkk[j-1]='\0'; /* sprintf(lkk, "%.*s", (int)(j - 1), strcpy(data, lkk)); */
StripSpaces(lkk);
/*****/
/*The next block reads the data to sweep*/
j = strpos2(lkk, "[", 1);
if (j == 0 || lkk[0] == '*') {
sprintf(alter[i - 1].text, "%.*s", (int)strlen(lkk), lkk);
i++;
}
else {
if (strpos2(lkk, " ", j)) {
printf("auxfunc_alter.c - CreateALTERinc -- More than one space exist in: %s\n", strsub(data, lkk, j, LONGSTRINGSIZE));
exit(EXIT_FAILURE);
}
if (strpos2(lkk, "[ ", j)) {
printf("auxfunc_alter.c - CreateALTERinc -- Space exist after '[' in: %s\n", strsub(data, lkk, j, LONGSTRINGSIZE));
exit(EXIT_FAILURE);
}
if (strpos2(lkk, " ]", j)) {
printf("auxfunc_alter.c - CreateALTERinc -- Space exist before ']' in: %s\n", strsub(data, lkk, j, LONGSTRINGSIZE));
exit(EXIT_FAILURE);
}
sprintf(alter[i - 1].text, "%.*s", (int)(j - 1), lkk);
/*begin*/
k = 0;
ReadSubKey(data, lkk, &j, '[', ' ', 4);
if (*data != '\0') {
k++;
read_sweep_vars(alter[i - 1].data[0], data);
}
l=1;
while (l<SWEEPVARS) {
ReadSubKey(data, lkk, &j, ' ', ' ', 0);
if (*data != '\0') {
k++;
read_sweep_vars(alter[i - 1].data[l], data);
} else
break;
l++;
}
sweep[ptr] = k;
order[ptr] = 1;
ptr++;
/*end;*/
i++;
}
/*****/
}
fgets2(lkk, LONGSTRINGSIZE, fsweepINI);
}
num_lines = i - 1;
/**/
/*Step3: */
alter_times = 1; /*alter_times: how many simulations?*/
kk = 0;
for (i = 0; i < num_lines; i++) {
if (*alter[i].data[0] != '\0') {
kk++;
alter_times *= sweep[kk - 1];
}
}
if (alter_times > 65536L || kk > SWEEPLINES) {
if (kk > SWEEPLINES) {
printf("auxfunc_alter.c - CreateALTERinc -- Maximum number of %d variables to sweep reached. Increase SWEEPLINES in auxfunc_alter.h\n", SWEEPLINES);
exit(EXIT_FAILURE);
} else {
printf("auxfunc_alter.c - CreateALTERinc -- More than 65536 simulations reached\n");
exit(EXIT_FAILURE);
}
} else {
if (append) { /*Assumes that is comming from errfunc.c*/
switch(spice) {
case 1: /*Eldo*/
case 2: /*HSPICE*/
fseek(falterINC, -5, SEEK_END); /*properly position the pointer*/
break;
case 3: /*LTspice*/
printf("auxfunc_alter.c - CreateALTERinc -- LTspice not supported\n");
exit(EXIT_FAILURE);
case 4: /*Spectre*/
break;
case 50: /*Qucs*/
printf("auxfunc_alter.c - CreateALTERinc -- Qucs not supported\n");
exit(EXIT_FAILURE);
case 100: /*general*/
printf("auxfunc_alter.c - CreateALTERinc -- GENERAL not supported\n");
exit(EXIT_FAILURE);
default:
printf("auxfunc_alter.c - CreateALTERinc -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
#ifdef __MINGW32__
fseek(falterINC, -1, SEEK_CUR); /*properly position the pointer*/
#endif
}
switch(spice) { /*Spectre special case*/
case 1: /*Eldo*/
case 2: /*HSPICE*/
break;
case 3: /*LTspice*/
break;
printf("auxfunc_alter.c - CreateALTERinc -- LTspice not supported\n");
case 4: /*Spectre*/
fprintf(falterINC, "simulator lang=spectre\n\n");
break;
case 50: /*Qucs*/
printf("auxfunc_alter.c - CreateALTERinc -- Qucs not supported\n");
break;
case 100: /*general*/
printf("auxfunc_alter.c - CreateALTERinc -- GENERAL not supported\n");
break;
default:
printf("auxfunc_alter.c - CreateALTERinc -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
/**/
/*Step4*/
for (index = 1; index <= alter_times; index++) {
add = 1; /*build ( [1111] [1112] [1113] [1121] ...*/
for (jj = kk - 1; jj >= 0; jj--) {
if (sweep[jj] != 0) {
if (order[jj] != sweep[jj]) {
if (index != 1)
order[jj] += add;
add = 0;
} else {
if (add == 1)
order[jj] = 1;
}
}
}
/**/
/*Step5: print the header with the information about the sweeped vars*/
j = 1;
switch(spice) {
case 1: /*Eldo*/
case 2: /*HSPICE*/
fprintf(falterINC, ".ALTER @%d ->", index);
break;
case 3: /*LTspice*/
printf("auxfunc_alter.c - CreateALTERinc -- LTspice not supported\n");
exit(EXIT_FAILURE);
case 4: /*Spectre*/
fprintf(falterINC, "//.ALTER @%d ->", index);
break;
case 50: /*Qucs*/
printf("auxfunc_alter.c - CreateALTERinc -- Qucs not supported\n");
exit(EXIT_FAILURE);
case 100: /*general*/
printf("auxfunc_alter.c - CreateALTERinc -- GENERAL not supported\n");
exit(EXIT_FAILURE);
default:
printf("auxfunc_alter.c - CreateALTERinc -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
for (k = 0; k < num_lines; k++) {
if (*alter[k].data[0] != '\0') {
strcpy(data, alter[k].text + 0); /*previously it had '1', but with '0' nothing is removed*/
sprintf(data, " %s", alter[k].text + 0);
fputs(data, falterINC);
l = order[j - 1];
j++;
fprintf(falterINC, "%s;", alter[k].data[l - 1]);
}
}
putc('\n', falterINC);
/**/
/*Step6: print each one of the alter lines*/
i = 0; /*For Spectre only: in the following block, 'i' is for the indication of EndOfAlterGroup*/
j = 1;
for (k = 0; k < num_lines; k++) {
if (*alter[k].data[0] != '\0') {
fputs(alter[k].text, falterINC);
l = order[j - 1];
j++;
fprintf(falterINC, "%s\n", alter[k].data[l - 1]);
} else {
if (*alter[k].text != '\0') {
switch(spice) { /*Spectre altergroup handling*/
case 1: /*Eldo*/
case 2: /*HSPICE*/
fprintf(falterINC, "%s\n", alter[k].text);
break;
case 3: /*LTspice*/
printf("auxfunc_alter.c - CreateALTERinc -- LTspice not supported\n");
exit(EXIT_FAILURE);
case 4: /*Spectre*/
if ( (k==0) || (i==1) ){ /*add index to altergroup line*/
strcpy(lkk, alter[k].text);
l = strpos2(lkk, " ", 1);
sprintf(data, "%d", index);
InsertString(lkk, data, l, l);
fprintf(falterINC, "%s\n", lkk);
} else
fprintf(falterINC, "%s\n", alter[k].text);
if (!strcmp(alter[k].text, "}")) /*Detects the end of altergroup as soon as it finds '}'*/
i=1; /* and stores this information in variable 'i' */
l = strpos2(alter[k].text, "}", 1); /*Verification that detection went well*/
if (l && (i==0) ) {
printf("auxfunc_alter.c - CreateALTERinc -- Failed the detection of the end of altergroup\n");
exit(EXIT_FAILURE);
}
break;
case 50: /*Qucs*/
printf("auxfunc_alter.c - CreateALTERinc -- Qucs not supported\n");
exit(EXIT_FAILURE);
case 100: /*general*/
printf("auxfunc_alter.c - CreateALTERinc -- GENERAL not supported\n");
exit(EXIT_FAILURE);
default:
printf("auxfunc_alter.c - CreateALTERinc -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
}
}
}
putc('\n', falterINC);
}
}
}
if (fsweepINI != NULL)
fclose(fsweepINI);
if (falterINC != NULL)
fclose(falterINC);
} /*CreateALTERinc*/
/*
* initialize variables and necessary files to optimize the file in the variable "filename"
* 1: Optimization Flow: should Alter and/or MonteCarlo be done?
* 2: read "# Parameters #"
* 3: read "# Measurements #"
* 4: create <hostname>.tmp file; add measurements and replace symbols
* 5: erase from memory parameters that will not be optimized
* 6: use initial values stored in the <inputfile>.cfg file
*/
int initialize(char *filename) /* , double *x) */
{
int i, ii;
int ccode;
char laux[LONGSTRINGSIZE], laux2[SHORTSTRINGSIZE], hostname[SHORTSTRINGSIZE] = "0";
/* <inputfile>.* <inputfile>.cfg <hostname>.tmp /extract/<file> */
FILE *fspice_source, *fspice_cfg, *fspice_tmp, *fextract;
/**/
/*Step1: Optimization Flow: should Alter and/or MonteCarlo be done?*/
sprintf(laux, "%s%s", filename, ".cfg");
if ((fspice_cfg =fopen(laux ,"rt")) == 0) {
printf("initialize.c - Step1 -- Cannot open config file: %s\n", laux);
exit(EXIT_FAILURE);
}
ReadKey(lkk, "#Optimization Flow#", fspice_cfg);
if (strcmp(lkk, "#Optimization Flow#")) {
printf("INFO: initialize.c - Step1 -- #Optimization Flow# key not found\n");
fflush(stdout);
} else {
fgets2(lkk, LONGSTRINGSIZE, fspice_cfg); /*should Alter be done?*/
Str2Lower(lkk);
ii=1;
ReadSubKey(laux, lkk, &ii, ':', ' ', 4);
if (!strcmp(laux, "yes")) /* Alter==yes */
AlterMC+=2;
else {
if (strcmp(laux, "no")) { /* Alter!=no */
printf("initialize.c - Step1 -- Incorrect line format: %s\n", lkk);
exit(EXIT_FAILURE);
}
}
fgets2(lkk, LONGSTRINGSIZE, fspice_cfg); /*should MonteCarlo be done?*/
Str2Lower(lkk);
ii=1;
ReadSubKey(laux, lkk, &ii, ':', ' ', 4);
if (!strcmp(laux, "yes")) /*MonteCarlo==yes*/
AlterMC+=1;
else {
if (strcmp(laux, "no")) { /*MonteCarlo!=no*/
printf("initialize.c - Step1 -- Incorrect line format: %s\n", lkk);
exit(EXIT_FAILURE);
}
}
fgets2(lkk, LONGSTRINGSIZE, fspice_cfg);
ii=1;
ReadSubKey(laux, lkk, &ii, ':', ' ', 4);
AlterMCcost=asc2real(laux, 1, (int)strlen(laux));
if (AlterMCcost < 0) {
printf("initialize.c - Step1 -- Minumum cost=%f, should be >= 0\n", AlterMCcost);
exit(EXIT_FAILURE);
}
#ifdef DEBUG
AlterMCcost=1.7976931348623157e+308; /*DBL_MAX from <float.h>*/
#endif
fgets2(lkk, LONGSTRINGSIZE, fspice_cfg); /*should RF parasitics be added?*/
Str2Lower(lkk);
ii=1;
ReadSubKey(laux, lkk, &ii, ':', ' ', 4);
if (!strcmp(laux, "yes")) /*RF==yes*/
ExecuteRF=1;
else {
if (strcmp(laux, "no")) { /*RF!=no*/
printf("initialize.c - Step1 -- Incorrect line format: %s\n", lkk);
exit(EXIT_FAILURE);
}
}
}
/**/
/*Step2: read "# Parameters #"*/
fseek(fspice_cfg, 0, SEEK_SET);
ReadKey(lkk, "# Parameters #", fspice_cfg); /*configuration parameters*/
if (strcmp(lkk, "# Parameters #")) {
printf("INFO: initialize.c - Step2 -- No parameters in config file\n");
fflush(stdout);
} else {
i=0;
fgets2(lkk, LONGSTRINGSIZE, fspice_cfg);
while ((lkk[0] != '#') && (lkk[0] != '\0') && (!feof(fspice_cfg))) {
if (lkk[0] != '*') {
ii=strpos2(lkk, ":", 1);
if (lkk[ii-1] != ':') {
printf("initialize.c - Step2 -- Incorrect line format: %s\n", lkk);
exit(EXIT_FAILURE);
}
strsub(parameters[i].name, lkk, 1, ii-1); /*name */
if (!(int)strlen(parameters[i].name)) {
printf("initialize.c - Step2 -- Name size equal to zero not allowed: %s\n", lkk);
exit(EXIT_FAILURE);
}
ReadSubKey(parameters[i].symbol, lkk, &ii, '#', '#', 5); /*symbol */
if (!(int)strlen(parameters[i].symbol)) {
printf("initialize.c - Step2 -- Symbol size equal to zero not allowed: %s\n", lkk);
exit(EXIT_FAILURE);
}
ReadSubKey(laux, lkk, &ii, ':', ':', 5);
parameters[i].value=asc2real(laux, 1, (int)strlen(laux)); /*value */
ReadSubKey(laux, lkk, &ii, ':', ':', 5);
parameters[i].minimum=asc2real(laux, 1, (int)strlen(laux)); /*minimum */
ReadSubKey(laux, lkk, &ii, ':', ':', 5);
parameters[i].maximum=asc2real(laux, 1, (int)strlen(laux)); /*maximum */
if (parameters[i].minimum > parameters[i].maximum) { /*just to help*/
printf("initialize.c - Step2 -- Minimum is larger than Maximum in line: %s\n", lkk);
exit(EXIT_FAILURE);
}
ReadSubKey(laux, lkk, &ii, ':', ':', 5); /*format */
if (laux[0] == 'E') {
char IEC60063[] = " 3 6 12 24 48 96 192 ";
strsub(laux2, laux, 2, 4); /*To keep with existing format of xxx_ while adding support for IEC 60063*/
sprintf(laux, "EEE_%s", laux2); /*preferred number series for resistors, capacitors and inductors */
if (!strpos2(IEC60063, laux2, 1)) { /*Only E3, E6, E12, E24, E48, E96 and E192 are allowed.*/
printf("initialize.c - Step2 -- Unrecognized option: %s\n", laux);
exit(EXIT_FAILURE);
}
}
if (laux[3] != 95) { /* 95="_" */
printf("initialize.c - Step2 -- Unrecognized option: %s\n", laux);
exit(EXIT_FAILURE);
}
parameters[i].format=-1; /*number '-1' is used for posterior value validation*/
strsub(laux2, laux, 1, 3);
if (!strcmp(laux2, "LIN"))
parameters[i].format=0;
if (!strcmp(laux2, "LOG"))
parameters[i].format=2;
if (!strcmp(laux2, "EEE"))
parameters[i].format=8;
if (parameters[i].format==-1) { /*validation*/
printf("initialize.c - Step2 -- Unrecognized option: %s\n", laux);
exit(EXIT_FAILURE);
}
strsub(laux2, laux, 5, (int)strlen(laux));
if (!strcmp(laux2, "DOUBLE"))
parameters[i].format=parameters[i].format+1;
if (!strcmp(laux2, "INT")) {
parameters[i].format=parameters[i].format+2;
/*Enforce INT format*/
if ((parameters[i].value - (int)parameters[i].value) > 0) {
printf("initialize.c - Step2 -- Number defined as INT is not integer in line: '%s'\n", lkk);
exit(EXIT_FAILURE);
}
if ((parameters[i].minimum - (int)parameters[i].minimum) > 0) {
printf("initialize.c - Step2 -- Number defined as INT is not integer in line: '%s'\n", lkk);
exit(EXIT_FAILURE);
}
if ((parameters[i].maximum - (int)parameters[i].maximum) > 0) {
printf("initialize.c - Step2 -- Number defined as INT is not integer in line: '%s'\n", lkk);
exit(EXIT_FAILURE);
}
/* if ((parameters[i].maximum-parameters[i].minimum) <= 0) {
printf("INFO: initialize.c - Step2 -- Minimum and Maximum are equal in line: '%s'\n", lkk);
} */
}
if (parameters[i].format==8)
parameters[i].format=parameters[i].format+(((int)asc2real(laux2, 1, (int)strlen(laux2)))<<8);
if (parameters[i].format==0) { /*validation*/
printf("initialize.c - Step2 -- Unrecognized option: %s\n", laux);
exit(EXIT_FAILURE);
}
ReadSubKey(laux, lkk, &ii, ':', ':', 4); /*optimize */
if (!strcmp(laux, "OPT"))
parameters[i].optimize=1; /*is it "OPT"? */
else {
if (!strcmp(laux, "---")) { /*if it is "---", then */
parameters[i].optimize=0; /*do not optimize */
parameters[i].minimum=parameters[i].value; /*furthermore, if it is just to define a quantity */
parameters[i].maximum=parameters[i].value; /*then make minimum=maximum=value */
} else {
printf("initialize.c - Step2 -- Unrecognized option: %s\n", laux);
exit(EXIT_FAILURE);
}
}
i++;
if (i > MAXPARAMETERS) {
printf("initialize.c - Step2 -- Maximum number of parameter reached (>%d). Increase MAXPARAMETERS in initialize.h\n",MAXPARAMETERS);
exit(EXIT_FAILURE);
}
}
fgets2(lkk, LONGSTRINGSIZE, fspice_cfg);
}
}
/**/
/*Step3: read "# Measurements #"*/
fseek(fspice_cfg, 0, SEEK_SET);
ReadKey(lkk, "# Measurements #", fspice_cfg); /*configuration measurements*/
if (strcmp(lkk, "# Measurements #")) {
printf("INFO: initialize.c - Step3 -- No measurements in config file\n");
fflush(stdout);
} else {
i=0;
fgets2(lkk, LONGSTRINGSIZE, fspice_cfg);
while ((lkk[0] != '#') && (lkk[0] != '\0') && (!feof(fspice_cfg))) {
if (lkk[0] != '*') {
ii=strpos2(lkk, ":", 1);
if (lkk[ii-1] != ':') {
printf("initialize.c - Step3 -- Incorrect line format: %s\n", lkk);
exit(EXIT_FAILURE);
}
strsub(measurements[i].meas_symbol, lkk, 1, ii-1); /*meas_symbol*/
sprintf(laux, "%i", i);
strcat(measurements[i].meas_symbol, laux); /* add number to symbol (*) */
ReadSubKey(measurements[i].node, lkk, &ii, ':', ':', 5); /*node*/
ReadSubKey(laux, lkk, &ii, ':', ':', 5); /*objective_constraint*/
measurements[i].objective_constraint=0; /*number '0' is used for posterior value validation*/
if (!strcmp(laux, "MIN"))
measurements[i].objective_constraint=1;
if (!strcmp(laux, "MAX"))
measurements[i].objective_constraint=2;
if (!strcmp(laux, "MON"))
measurements[i].objective_constraint=3;
if (!strcmp(laux, "LE"))
measurements[i].objective_constraint=4;
if (!strcmp(laux, "GE"))
measurements[i].objective_constraint=5;
if (!strcmp(laux, "EQ"))
measurements[i].objective_constraint=6;
if (measurements[i].objective_constraint==0) { /*validation*/
printf("initialize.c - Step3 -- Unrecognized option: %s\n", laux);
exit(EXIT_FAILURE);
}
ReadSubKey(laux, lkk, &ii, ':', ':', 4);
measurements[i].constraint_value=asc2real(laux, 1, (int)strlen(laux)); /*constraint_value*/
i++;
if (i > MAXMEASUREMENTS-1) {
printf("initialize.c - Step3 -- Maximum number of measurements reached (>%d). Increase MAXMEASUREMENTS in initialize.h\n", MAXMEASUREMENTS);
exit(EXIT_FAILURE);
}
}
fgets2(lkk, LONGSTRINGSIZE, fspice_cfg);
}
measurements[i].meas_symbol[0]='\0'; /*just in case... zero the first remaining entry*/
}
fclose(fspice_cfg);
/**/
/*Step4: create <hostname>.tmp file; add measurements and replace symbols*/
switch(spice) {
case 1: /*Eldo*/
sprintf(laux, "%s%s", filename, ".cir");
break;
case 2: /*HSPICE*/
sprintf(laux, "%s%s", filename, ".sp");
break;
case 3: /*LTspice*/
sprintf(laux, "%s%s", filename, ".net");
break;
case 4: /*Spectre*/
sprintf(laux, "%s%s", filename, ".scs");
break;
case 50: /*Qucs*/
sprintf(laux, "%s%s", filename, ".txt");
break;
case 51: /*ngspice*/
sprintf(laux, "%s%s", filename, ".sp");
break;
case 100: /*general*/
sprintf(laux, "%s%s", filename, ".txt");
break;
default:
printf("initialize.c - Step4 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
if ((fspice_source=fopen(laux, "rt")) == 0) { /*source netlist*/
printf("initialize.c - Step4 -- Cannot open input file: %s\n", laux);
exit(EXIT_FAILURE);
}
/**/
if ((ccode = gethostname(hostname, sizeof(hostname))) != 0) { /* !=0 can most likelly be deleted from all lines */
printf("initialize.c - Step4 -- gethostname failed, ccode = %d\n", ccode);
exit(EXIT_FAILURE);
}
/* printf("host name: %s\n", hostname); */
ii=strpos2(hostname, ".", 1);
if (ii) /* hostname is "longmorn.xx.xx.xx" */
hostname[ii-1]='\0';
sprintf(lkk, "%s%s", hostname, ".tmp"); /* hostname is "longmorn" */
if ((fspice_tmp =fopen(lkk ,"wt")) == 0) { /* netlist to simulate given by "hostname" */
printf("initialize.c - Step4 -- Cannot write to tmp file: %s\n", lkk);
exit(EXIT_FAILURE);
}
/*Step4.1: ".end" not yet found*/
fgets2(lkk, LONGSTRINGSIZE, fspice_source); /*read and */
fprintf(fspice_tmp, "%s\n", lkk); /*write the first line */
while (!P_eof(fspice_source)) {
fgets2(lkk, LONGSTRINGSIZE, fspice_source);
strcpy(laux, lkk); /*detect ".end", ".END", ".End", ... */
Str2Lower(laux);
StripSpaces(laux); /* avoid spaces after the command ".end" */
if (!strcmp(laux, ".end"))
break;
/***** -------------- *********** -------------- *****/
/***** -------------- ** BEGIN ** -------------- *****/
if (lkk[0]!='*') {
i=inlinestrpos(lkk);
ii=1;
ReadSubKey(laux, lkk, &ii, '#', '#', 0);
if ( (laux[0]=='\0') || (ii>(int)strlen(lkk)) || ((i<ii) && (i!=0)) ) { /*does it contains #<text>#? */
if ((int)strlen(lkk) && (!RFModule(lkk, 0, fspice_tmp)) )
fprintf(fspice_tmp, "%s\n", lkk); /* no: write line to <hostname>.tmp */
} else { /* yes: replace #<text># in this line */
if (!RFModule(lkk, 0, fspice_tmp)) { /* -Symbol replaced in the RFmodule*/
ReplaceSymbol(lkk, 0); /* -Symbol yet to be replaced */
fprintf(fspice_tmp, "%s\n", lkk); /* write line to <hostname>.tmp */
}
}
}
/***** -------------- ** END ** -------------- *****/
/***** -------------- *********** -------------- *****/
}
switch(spice) {
case 1: /*Eldo*/
if (strcmp(laux, ".end")) { /*Exit if ".end" is not found*/
printf("initialize.c - Step4.1 -- End not found in %s.cir\n", filename);
exit(EXIT_FAILURE);
}
break;
case 2: /*HSPICE*/
if (strcmp(laux, ".end")) { /*Exit if ".end" is not found*/
printf("initialize.c - Step4.1 -- End not found in %s.sp\n", filename);
exit(EXIT_FAILURE);
}
break;
case 3: /*LTspice*/
if (strcmp(laux, ".end")) { /*Exit if ".end" is not found*/
printf("initialize.c - Step4.1 -- End not found in %s.net\n", filename);
exit(EXIT_FAILURE);
}
break;
case 4: /*Spectre*/ /* ".end" does not exist in Spectre syntax */
break;
case 50: /*Qucs*/ /* ".end" does not exist in Qucs syntax */
break;
case 51: /*ngspice*/
if (strcmp(laux, ".end")) { /*Exit if ".end" is not found*/
printf("initialize.c - Step4.1 -- End not found in %s.sp\n", filename);
exit(EXIT_FAILURE);
}
break;
case 100: /*general*/
break;
default:
printf("initialize.c - Step4.1 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
/*Special case to deal with Spectre MDL*/
switch(spice) {
case 1: /*Eldo*/
break;
case 2: /*HSPICE*/
break;
case 3: /*LTspice*/
break;
case 4: /*Spectre*/
fclose(fspice_tmp);
sprintf(lkk, "%s%s", hostname, ".mdl"); /* hostname is "longmorn" */
if ((fspice_tmp =fopen(lkk ,"wt")) == 0) { /* netlist to simulate given by "hostname" */
printf("initialize.c - Step4.1 -- Cannot write to tmp file: %s\n", lkk);
exit(EXIT_FAILURE);
}
fseek(fspice_source, 0, SEEK_SET);
while (!P_eof(fspice_source)) {
fgets2(lkk, LONGSTRINGSIZE, fspice_source);
if ( (lkk[0] != '*') && (lkk[0] != '\0') && (!strpos2(lkk, "//", 1)) ) {
strcpy(laux, lkk); /*detect "TRAN", "tran", "TRan", ... */
Str2Lower(laux);
StripSpaces(laux); /*avoid spaces after the command*/
if (strpos2(laux, " dc ", 1)) {
fprintf(fspice_tmp, "alias measurement dc_run {\n");
i=strpos2(laux, " dc ", 1);
strsub(laux, lkk, 1, i);
StripSpaces(laux);
fprintf(fspice_tmp, "run %s\n", laux);
}
if (strpos2(laux, " ac ", 1)) {
fprintf(fspice_tmp, "alias measurement ac_run {\n");
i=strpos2(laux, " ac ", 1);
strsub(laux, lkk, 1, i);
StripSpaces(laux);
fprintf(fspice_tmp, "run %s\n", laux);
}
if (strpos2(laux, " tran ", 1)) {
fprintf(fspice_tmp, "alias measurement tran_run {\n");
i=strpos2(laux, " tran ", 1);
strsub(laux, lkk, 1, i);
StripSpaces(laux);
fprintf(fspice_tmp, "run %s\n", laux);
}
}
}
break;
case 50: /*Qucs*/
break;
case 51: /*ngspice*/
break;
case 100: /*general*/
break;
default:
printf("initialize.c - Step4.1 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
/*Special case to deal with Spectre MDL*/
/*Step4.2: End of file is found and now add measurements*/
i=0;
fprintf(fspice_tmp, "\n");
while (strcmp(measurements[i].meas_symbol,"\0") ) { /*until the end of all symbols*/
sprintf(lkk, "%i", i); /* to remove integer added in (*) */
ii=(int)strlen(lkk); /* to remove integer added in (*) */
sprintf(lkk, "%s%s", "extract/", measurements[i].meas_symbol); /* to remove integer added in (*) */
lkk[(int)strlen(lkk)-ii]='\0'; /* to remove integer added in (*) */
if ((fextract =fopen(lkk ,"rt")) == 0) {
printf("initialize.c - Step4.2 -- Cannot find measurement file: %s\n", lkk);
exit(EXIT_FAILURE);
}
strcpy(laux,lkk);
ReadKey(lkk, "# Commands #", fextract);
if (strcmp(lkk, "# Commands #")) { /*finds "# Commands #" and writes to <hostname>.tmp until the end of file*/
printf("initialize.c - Step4.2 -- Wrong format in file: %s\n", laux);
exit(EXIT_FAILURE);
}
/*Step4.2.1: Add the measurement to <hostname>.tmp*/
/*---------------------------------------------------------------*/
sprintf(lkk, "%i", i); /* to remove integer added in (*) */
ii=(int)strlen(lkk); /* to remove integer added in (*) */
strcpy(lkk, measurements[i].meas_symbol); /* to remove integer added in (*) */
lkk[(int)strlen(lkk)-ii]='\0'; /* to remove integer added in (*) */
switch(spice) {
case 1: /*Eldo*/
fprintf(fspice_tmp, "* %i) Extract \'%s\'\n", i+1, lkk);
break;
case 2: /*HSPICE*/
fprintf(fspice_tmp, "* %i) Extract \'%s\'\n", i+1, lkk);
break;
case 3: /*LTspice*/
fprintf(fspice_tmp, "* %i) Extract \'%s\'\n", i+1, lkk);
break;
case 4: /*Spectre*/
fprintf(fspice_tmp, "// %i) Extract \'%s\'\n", i+1, lkk);
break;
case 50: /*Qucs*/
fprintf(fspice_tmp, "# %i) Extract \'%s\'\n", i+1, lkk);
break;
case 51: /*ngspice*/
fprintf(fspice_tmp, "* %i) Extract \'%s\'\n", i+1, lkk);
break;
case 100: /*general*/
fprintf(fspice_tmp, "* %i) Extract \'%s\'\n", i+1, lkk);
break;
default:
printf("initialize.c - Step4.2.1 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
fgets2(lkk, LONGSTRINGSIZE, fextract); /* reads from directory "extract/" */
while ((lkk[0] != '#') && (lkk[0] != '\0') && (!feof(fextract))) {
switch(spice) {
case 1: /*Eldo*/
strcpy(laux,lkk);
StripSpaces(laux);
Str2Lower(laux);
if (strpos2(laux, ".meas ", 1)) {
sprintf(lkk, "%i", i); /* to remove integer added in (*) */
ii=(int)strlen(lkk); /* to remove integer added in (*) */
sprintf(lkk, "%s%s", "extract/", measurements[i].meas_symbol); /* to remove integer added in (*) */
lkk[(int)strlen(lkk)-ii]='\0'; /* to remove integer added in (*) */
printf("initialize.c - Step4.2.1 -- .MEAS not supported in file: %s. Use .EXTRACT instead.\n", lkk);
exit(EXIT_FAILURE);
}
break;
case 2: /*HSPICE*/
break;
case 3: /*LTspice*/
break;
case 4: /*Spectre*/
break;
case 50: /*Qucs*/
break;
case 51: /*ngspice*/
break;
case 100: /*general*/
break;
default:
printf("initialize.c - Step4.2.1 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
ii=1;
ReadSubKey(laux, lkk, &ii, '#', '#', 0);
if (laux[0]!='\0') {
DecodeSymbolNode(lkk, i); /*has to replace #<text># in this line*/
}
if (lkk[0]!='#') { /*if end of block has not been reached*/
fprintf(fspice_tmp, "%s\n", lkk); /* writes to <hostname>.tmp */
}
fgets2(lkk, LONGSTRINGSIZE, fextract); /*reads from directory "extract/"*/
}
fprintf(fspice_tmp, "\n");
/*Step4.2.2: Add entry to variable 'measure[i].var_name' and 'measure[i].search'*/
/*---------------------------------------------------------------*/
fseek(fextract, 0, SEEK_SET);
ReadKey(lkk, "MEASURE_VAR", fextract);
{
measure[0].search[0]=10;/*"ProcessOutputFile", line 1231: makes data to read not found on measure[0]*/
}
ii=1;
while ((*measure[ii].var_name) != '\0') /* finds the proper entry */
ii++; /* place. Store in 'ii' */
if (ii > MAXMEAS-1) {
printf("initialize.c - Step4.2.2 -- Maximum number of measurements reached (>%d). Increase MAXMEAS in auxfunc_measurefromlis.h\n", MAXMEAS);
exit(EXIT_FAILURE);
}
if (strcmp((sprintf(laux, "%.11s", lkk), laux), "MEASURE_VAR")) { /*general case, if it's not a "MEASURE_VAR"; data from "/extract/<file>"*/
sprintf(lkk, "%s%s", UNIQUECHAR,measurements[i].meas_symbol);
strcpy(measure[ii].var_name, lkk); /*measure[ii].var_name*/
switch(spice) { /* the format in which the variables are written in the output file: fast read */
case 1: /*Eldo*/
sprintf(lkk, " %s%s =", UNIQUECHAR, measurements[i].meas_symbol);
Str2Upper(lkk);
break;
case 2: /*HSPICE*/
sprintf(lkk, " %s%s=", UNIQUECHAR, measurements[i].meas_symbol);
Str2Lower(lkk);
break;
case 3: /*LTspice*/
sprintf(lkk, "%s%s:", UNIQUECHAR, measurements[i].meas_symbol);
Str2Lower(lkk);
break;
case 4: /*Spectre*/
sprintf(lkk, "%s%s", UNIQUECHAR, measurements[i].meas_symbol);
ccode=(int)strlen(lkk);
while (ccode<18) {
strcat(lkk, " ");
ccode++;
}
strcat(lkk, " =");
break;
case 50: /*Qucs*/
break;
case 51: /*ngspice*/
sprintf(lkk, "%s%s", UNIQUECHAR, measurements[i].meas_symbol);
Str2Lower(lkk);
ccode=(int)strlen(lkk);
while (ccode<20) {
strcat(lkk, " ");
ccode++;
}
strcat(lkk, "=");
break;
case 100: /*general*/
break;
default:
printf("initialize.c - Step4.2.2 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
strcpy(measure[ii].search, lkk); /*measure[ii].var_name*/
} else { /*a line with "MEASURE_VAR" exist in "/extract/<file>"*/
while (!strcmp((sprintf(laux, "%.11s", lkk), laux), "MEASURE_VAR")) {
DecodeSymbolNode(lkk, i);
#ifndef __MINGW32__
ii=ProcessMeasureVar(lkk, ii, "/dev/null");
#else
ii=ProcessMeasureVar(lkk, ii, "NUL");
#endif
switch(spice) {
case 1: /*Eldo*/
ccode=strpos2(measure[ii].search, UNIQUECHAR, 1);
//if (ccode) {
// Str2Upper(measure[ii].search);
//}
break;
case 2: /*HSPICE*/
ccode=strpos2(measure[ii].search, UNIQUECHAR, 1);
//if (ccode) {
// Str2Lower(measure[ii].search);
//}
break;
case 3: /*LTspice*/
ccode=strpos2(measure[ii].search, UNIQUECHAR, 1);
//if (ccode) {
// Str2Lower(measure[ii].search);
//}
break;
case 4: /*Spectre*/
break;
case 50: /*Qucs*/
break;
case 51: /*ngspice*/
ccode=strpos2(measure[ii].search, "MATH", 1); /*DO NOT change MATH lines character case*/
if (!ccode){
ccode=strpos2(measure[ii].search, UNIQUECHAR, 1);
if (ccode) {
Str2Lower(measure[ii].search);
}
}
break;
case 100: /*general*/
break;
default:
printf("initialize.c - Step4.2.2 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
ReadKey(lkk, "MEASURE_VAR", fextract);
if ((int)strlen(lkk)) {
ii++;
if (ii > MAXMEAS-1) {
printf("initialize.c - Step4.2.2 -- Maximum number of measurements reached (>%d). Increase MAXMEAS in auxfunc_measurefromlis.h\n", MAXMEAS);
exit(EXIT_FAILURE);
}
}
}
}
/*---------------------------------------------------------------*/
/*Step4.2.3: Continues ... */
fclose(fextract);
i++;
}
/*Step4.3: Add ".end" where required*/
switch(spice) {
case 1: /*Eldo*/
fprintf(fspice_tmp, "%s\n", ".end");
break;
case 2: /*HSPICE*/
fprintf(fspice_tmp, "%s\n", ".end");
break;
case 3: /*LTspice*/
fprintf(fspice_tmp, "%s\n", ".end");
break;
case 4: /*Spectre*/
/*Special case to deal with Spectre MDL*/
ccode=0; /* at this momment, only one measurement can exist in the <inputfile>.scs */
fseek(fspice_source, 0, SEEK_SET);
fprintf(fspice_tmp, "}\n\n");
while (!P_eof(fspice_source)) {
fgets2(lkk, LONGSTRINGSIZE, fspice_source);
if ( (lkk[0] != '*') && (lkk[0] != '\0') && (!strpos2(lkk, "//", 1)) ) {
strcpy(laux, lkk); /*detect "TRAN", "tran", "TRan", ... */
Str2Lower(laux);
StripSpaces(laux); /*avoid spaces after the command*/
if (strpos2(laux, " dc ", 1)) {
if (ccode!=0) {
printf("initialize.c - Step4.3 -- Only one type of simulation is implememted at this time!\n");
exit(EXIT_FAILURE);
}
ccode++;
i=strpos2(laux, " dc ", 1);
strsub(laux, lkk, 1, i);
StripSpaces(laux);
fprintf(fspice_tmp, "run dc_run as dc1\n");
}
if (strpos2(laux, " ac ", 1)) {
if (ccode!=0) {
printf("initialize.c - Step4.3 -- Only one type of simulation is implememted at this time!\n");
exit(EXIT_FAILURE);
}
ccode++;
i=strpos2(laux, " ac ", 1);
strsub(laux, lkk, 1, i);
StripSpaces(laux);
fprintf(fspice_tmp, "run ac_run as ac1\n");
}
if (strpos2(laux, " tran ", 1)) {
if (ccode!=0) {
printf("initialize.c - Step4.3 -- Only one type of simulation is implememted at this time!\n");
exit(EXIT_FAILURE);
}
ccode++;
i=strpos2(laux, " tran ", 1);
strsub(laux, lkk, 1, i);
StripSpaces(laux);
fprintf(fspice_tmp, "run tran_run as tran1\n");
}
}
}
/*Special case to deal with Spectre MDL*/
break;
case 50: /*Qucs*/
break;
case 51: /*ngspice*/
fprintf(fspice_tmp, "%s\n", ".end");
break;
case 100: /*general*/
break;
default:
printf("initialize.c - Step4.3 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
fclose(fspice_source);
fclose(fspice_tmp);
/**/
/*Step5: erase from memory parameters that will not be optimized*/
ii=0;
for (i = 0; i < MAXPARAMETERS; i++) {
if ((parameters[i].optimize == 1) && (i>ii)) { /*if ... then move*/
strcpy (parameters[ii].name, parameters[i].name);
strcpy (parameters[ii].symbol,parameters[i].symbol);
parameters[ii].value = parameters[i].value;
parameters[ii].minimum = parameters[i].minimum;
parameters[ii].maximum = parameters[i].maximum;
parameters[ii].format = parameters[i].format;
parameters[ii].optimize = parameters[i].optimize;
ii++;
} else
if (parameters[i].optimize)
ii++;
}
for (i = ii; i < MAXPARAMETERS; i++) { /*just in case... zero all remaining entries*/
parameters[i].name[0] ='\0';
parameters[i].symbol[0]='\0';
parameters[i].value = 0;
parameters[i].minimum = 0;
parameters[i].maximum = 0;
parameters[i].format = 0;
parameters[i].optimize = 0;
}
/**/
/*Initialization6: use initial values stored in the <inputfile>.cfg file */
/*for (ii = 0; ii < MAXPARAMETERS; ii++) { */
/* x[ii] = scaleto(parameters[ii].value, parameters[ii].minimum, parameters[ii].maximum, -10, +10); */
/*} */
/**/
/*Initialization7*/
return EXIT_SUCCESS;
}
/*
* This is the main function
* 1: scale input parameters from [-10, +10] to fit within desired range values
* 2: read from <hostname>.tmp and write to <hostname>.*
* 3: execute 'SPICE <hostname>.*'
* 4: read output <hostname>.out
* 5: find cost
* 6: save log information to file <hostname>.log
* 7: ALTER and MC simulation
* 8: return cost
*/
double errfunc(char *filename, double *x)
{
/*double currentcost;*/ /*total cost*/
int i, ii;
int ccode;
char laux[LONGSTRINGSIZE], hostname[SHORTSTRINGSIZE] = "0", filename_x[SHORTSTRINGSIZE] = "0";
/* <hostname>.* <hostname>.out <hostname>.tmp <hostname>.log*/
FILE *fspice_input, *fspice_output, *fspice_tmp, *fspice_log;
/**/
/*Step1: scale input parameters from [-10, +10] to fit within desired range values*/
#ifdef DEBUG
printf("DEBUG: errfunc.c - Step1\n");
#endif
for (i = 0; i < MAXPARAMETERS; i++) {
if (parameters[i].format) {
parameters[i].value = scaleto(x[i], -10, +10, parameters[i].minimum, parameters[i].maximum, parameters[i].format);
if (parameters[i].value<0) {
ii=1;
/*printf("INFO: errfunc.c - Step1 -- Negative values\n");*/
}
}
}
/**/
/*Step2: read from <hostname>.tmp and write to <hostname>.**/
/* */
/* It is not necessary to check every time for the possibility to read/write. Move to initicialization if possible!! */
/* */
#ifdef DEBUG
printf("DEBUG: errfunc.c - Step2\n");
#endif
if ((ccode = gethostname(hostname, sizeof(hostname))) != 0) {
printf("errfunc.c - Step2 -- gethostname failed, ccode = %d\n", ccode);
exit(EXIT_FAILURE);
}
/* printf("host name: %s\n", hostname); */
ii=strpos2(hostname, ".", 1);
if (ii) /* hostname is "longmorn.xx.xx.xx" */
hostname[ii-1]='\0';
sprintf(lkk, "%s%s", hostname, ".tmp"); /* hostname is "longmorn" */
if ((fspice_tmp =fopen(lkk ,"rt")) == 0) { /* netlist to simulate given by "hostname" */
printf("errfunc.c - Step2 -- Cannot read from tmp file: %s\n", lkk);
exit(EXIT_FAILURE);
}
switch(spice) {
case 1: /*Eldo*/
sprintf(lkk, "%s%s", hostname, ".cir");
break;
case 2: /*HSPICE*/
sprintf(lkk, "%s%s", hostname, ".sp");
break;
case 3: /*LTspice*/
sprintf(lkk, "%s%s", hostname, ".net");
break;
case 4: /*Spectre*/
sprintf(lkk, "%s%s", hostname, ".scs");
break;
case 50: /*Qucs*/
sprintf(lkk, "%s%s", hostname, ".txt");
break;
case 51: /*ngspice*/
sprintf(lkk, "%s%s", hostname, ".sp");
break;
case 100: /*general*/
sprintf(lkk, "%s%s", hostname, ".txt");
break;
default:
printf("errfunc.c - Step2 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
if ((fspice_input =fopen(lkk ,"wt")) == 0) { /*netlist to simulate given by 'hostname'*/
printf("errfunc.c - Step2 -- Cannot write to output file: %s\n", lkk);
exit(EXIT_FAILURE);
}
/*Step2.1: ".end" not yet found*/
fgets2(lkk, LONGSTRINGSIZE, fspice_tmp); /*read and*/
fprintf(fspice_input, "%s\n", lkk); /*write the first line*/
while (!P_eof(fspice_tmp)) {
fgets2(lkk, LONGSTRINGSIZE, fspice_tmp);
strcpy(laux, lkk); /*detect ".end", ".END", ".End", ... */
Str2Lower(laux);
StripSpaces(laux); /* avoid spaces after the command ".end" */
if (!strcmp(laux, ".end"))
break;
/***** -------------- *********** -------------- *****/
/***** -------------- ** BEGIN ** -------------- *****/
if (lkk[0]!='*') {
i=inlinestrpos(lkk);
ii=1;
ReadSubKey(laux, lkk, &ii, '#', '#', 0);
if ( (laux[0]=='\0') || (ii>(int)strlen(lkk)) || ((i<ii) && (i!=0)) ) { /* does it contains #<text>#? */
if ((int)strlen(lkk) && (!RFModule(lkk, 1, fspice_input)) )
fprintf(fspice_input, "%s\n", lkk); /*no, write line to <hostname>.* */
} else { /*yes, replace #<text># in this line */
if (!RFModule(lkk, 1, fspice_input)) {
ReplaceSymbol(lkk, 1);
fprintf(fspice_input, "%s\n", lkk); /* write line to <hostname>.* */
}
}
}
/***** -------------- ** END ** -------------- *****/
/***** -------------- *********** -------------- *****/
}
switch(spice) {
case 1: /*Eldo*/
if (strcmp(laux, ".end")) { /*Exit if ".end" is not found*/
printf("errfunc.c - Step2.1 -- End not found in %s.cir\n", filename);
exit(EXIT_FAILURE);
}
break;
case 2: /*HSPICE*/
if (strcmp(laux, ".end")) { /*Exit if ".end" is not found*/
printf("errfunc.c - Step2.1 -- End not found in %s.sp\n", filename);
exit(EXIT_FAILURE);
}
break;
case 3: /*LTspice*/
if (strcmp(laux, ".end")) { /*Exit if ".end" is not found*/
printf("errfunc.c - Step2.1 -- End not found in %s.net\n", filename);
exit(EXIT_FAILURE);
}
break;
case 4: /*Spectre*/ /* ".end" does not exist in Spectre syntax */
break;
case 50: /*Qucs*/ /* ".end" does not exist in Qucs syntax */
break;
case 51: /*ngspice*/
if (strcmp(laux, ".end")) { /*Exit if ".end" is not found*/
printf("errfunc.c - Step2.1 -- End not found in %s.sp\n", filename);
exit(EXIT_FAILURE);
}
break;
case 100: /*general*/
break;
default:
printf("errfunc.c - Step2.1 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
/*Step2.3: Add ".end" where required*/
switch(spice) {
case 1: /*Eldo*/
fprintf(fspice_input, "%s\n", ".end");
break;
case 2: /*HSPICE*/
fprintf(fspice_input, "%s\n", ".end");
break;
case 3: /*LTspice*/
fprintf(fspice_input, "%s\n", ".end");
break;
case 4: /*Spectre*/
break;
case 50: /*Qucs*/
break;
case 51: /*ngspice*/
fprintf(fspice_input, "%s\n", ".end");
break;
case 100: /*general*/
break;
default:
printf("errfunc.c - Step2.3 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
fclose(fspice_input);
fclose(fspice_tmp);
/**/
/*Step3: execute 'SPICE <hostname>.*'*/
#ifdef DEBUG
printf("DEBUG: errfunc.c - Step3\n");
#endif
switch(spice) {
case 1: /*Eldo*/
sprintf(lkk, "nice -n 19 eldo -noconf -i %s.cir > %s.out", hostname, hostname);
break;
case 2: /*HSPICE*/
#ifndef __MINGW32__
sprintf(lkk, "nice -n 19 hspice -i %s.sp -o %s.lis > /dev/null", hostname, hostname);
#else
sprintf(lkk, "hspice -i %s.sp -o %s.lis > NUL", hostname, hostname);
#endif
break;
case 3: /*LTspice*/
#ifndef __MINGW32__
sprintf(lkk, "nice -n 19 ltspice -b %s.net > /dev/null", hostname);
#else
sprintf(lkk, "ltspice -b %s.net > NUL", hostname);
#endif
break;
case 4: /*Spectre*/
sprintf(lkk, "nice -n 19 spectremdl -batch %s.mdl -design %s.scs > /dev/null", hostname, hostname);
break;
case 50: /*Qucs*/
#ifndef __MINGW32__
sprintf(lkk, "nice -n 19 qucsator -i %s.txt -o %s.dat > /dev/null", hostname, hostname);
#else
sprintf(lkk, "qucsator -i %s.txt -o %s.dat > NUL", hostname, hostname);
#endif
break;
case 51: /*ngspice*/
#ifndef __MINGW32__
sprintf(lkk, "nice -n 19 ngspice -b -o %s.out %s.sp > /dev/null 2>&1", hostname, hostname);
#else
sprintf(lkk, "ngspice -o %s.out %s.sp > NUL", hostname, hostname);
#endif
break;
case 100: /*general*/
#ifndef __MINGW32__
sprintf(lkk, "nice -n 19 ./general.sh %s %s", hostname, hostname);
#else
sprintf(lkk, "./general.sh %s %s", hostname, hostname);
#endif
break;
default:
printf("errfunc.c - Step3 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
ii=system(lkk);
#ifndef __MINGW32__
if (WIFSIGNALED(ii) && (WTERMSIG(ii) == SIGINT || WTERMSIG(ii) == SIGQUIT)) {
printf("errfunc.c - Step3 -- Ctrl-C key pressed. Exiting optimization loop.\n");
fflush(stdout);
#ifdef MPI
/* exit(EXIT_FAILURE); */
#else
return(0); /*returned simulation cost is zero; Ctrl-C detection*/
#endif
}
#endif
/**/
/*Step4: read output <hostname>.out, find maxcost and save log information to file <hostname>.log*/
#ifdef DEBUG
printf("DEBUG: errfunc.c - Step4\n");
#endif
switch(spice) {
case 1: /*Eldo*/
sprintf(lkk, "%s.out", hostname);
break;
case 2: /*HSPICE*/
sprintf(lkk, "%s.lis", hostname);
break;
case 3: /*LTspice*/
sprintf(lkk, "%s.log", hostname);
break;
case 4: /*Spectre*/
sprintf(lkk, "%s.measure", hostname);
break;
case 50: /*Qucs*/
sprintf(lkk, "%s.dat", hostname);
break;
case 51: /*ngspice*/
sprintf(lkk, "%s.out", hostname);
break;
case 100: /*general*/
sprintf(lkk, "%s.out", hostname);
break;
default:
printf("errfunc.c - Step4 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
maxcost=0;
#ifdef DEBUG
ProcessOutputFile(lkk, 3); /* =>3:mem+file */
#else
ProcessOutputFile(lkk, 1); /* =>1:mem */
#endif
/**/
/*Step5: find cost*/
/* #ifdef DEBUG */
/* printf("DEBUG: errfunc.c - Step5\n"); */
/* #endif */
/* currentcost=CostFunction(); */
/* if (maxcost<currentcost) */
/* maxcost=currentcost; */
/**/
/*Step6: save log information to file <hostname>.log*/
/* #ifdef DEBUG */
/* printf("DEBUG: errfunc.c - Step6\n"); */
/* #endif */
/* if (LOG) { */
/* LogtoFile(currentcost); */
/* } */
/**/
/*Step7: ALTER and MC simulation*/
#ifdef DEBUG
printf("DEBUG: errfunc.c - Step7\n");
#endif
if ( (((AllConstraintsMet()) && (AlterMCcost > maxcost)) || (AlterMCcost > maxcost)) && (AlterMC)) {
switch (AlterMC) {
case 1: /*Monte Carlo simulation*/
/* bla bla bla*/
printf("INFO: errfunc.c - Step7 -- altermc=%d\n", AlterMC);
fflush(stdout);
if (LOG) {
sprintf(laux, "%s.log", hostname);
if ((fspice_log=fopen(laux,"at")) == 0) {
printf("errfunc.c - Step7 -- Cannot open log file: %s\n", laux);
exit(EXIT_FAILURE);
}
fprintf(fspice_log, "Alter: 1 - altermc=%d\n", AlterMC);
fclose(fspice_log);
}
strcpy (filename_x, filename);
strcat(filename_x, ".cfg");
sprintf(lkk, "%s%s", hostname, ".tmp"); /* hostname is "longmorn" */
sprintf(laux, "%s%s", hostname, ".mc"); /* hostname is "longmorn" */
MonteCarlo(filename_x, lkk, laux); /* execute the 'monte' tool */
sprintf(lkk, "%s%s", hostname, ".mc"); /* hostname is "longmorn" */
if ((fspice_input =fopen(lkk ,"rt")) == 0) { /*netlist to simulate given by "hostname" */
printf("errfunc.c - Step7 -- Cannot read from mc file: %s\n", lkk);
exit(EXIT_FAILURE);
}
sprintf(lkk, "%s%s", hostname, ".tmp");
if ((fspice_output =fopen(lkk ,"wt")) == 0) { /*netlist to simulate given by 'hostname' */
printf("errfunc.c - Step7 -- Cannot write to tmp file: %s\n", lkk);
exit(EXIT_FAILURE);
}
while (!P_eof(fspice_input)) { /* reads from <hostname>.mc and writes to <hostname>.tmp */
fgets2(lkk, LONGSTRINGSIZE, fspice_input); /* read line from <hostname>.mc */
strcpy(laux, lkk);
Str2Lower(laux);
fprintf(fspice_output, "%s\n", lkk); /* write line to <hostname>.tmp */
}
fclose(fspice_input);
fclose(fspice_output);
AlterMC=AlterMC-1; /*decrease by one; nothing more will be executed*/
break;
case 2: /*Alter simulation - variable 'AlterMC'=2*/
case 3: /*Alter simulation - variable 'AlterMC'=3 and will execute a MonteCarlo simulation afterwards*/
/* bla bla bla*/
printf("INFO: errfunc.c - Step7 -- altermc=%d\n", AlterMC);
fflush(stdout);
if (LOG) {
sprintf(laux, "%s.log", hostname);
if ((fspice_log=fopen(laux,"at")) == 0) {
printf("errfunc.c - Step7 -- Cannot open log file: %s\n", laux);
exit(EXIT_FAILURE);
}
fprintf(fspice_log, "Alter: 2 - altermc=%d\n", AlterMC);
fclose(fspice_log);
}
sprintf(lkk, "%s%s", hostname, ".tmp"); /*hostname is "longmorn"*/
if ((fspice_tmp =fopen(lkk ,"r+t")) == 0) { /*netlist to simulate given by "hostname"*/
printf("errfunc.c - Step7 -- Cannot read from tmp file: %s\n", lkk);
exit(EXIT_FAILURE);
}
if ((fspice_log =fopen("alter.inc" ,"rt")) == 0) { /*If file 'alter.inc' does not exist */
strcpy (filename_x, filename);
strcat(filename_x, ".cfg");
CreateALTERinc(filename_x, lkk, 1); /*execute the 'alter' tool*/
fseek(fspice_tmp, 0, SEEK_END); /*properly position the pointer*/
} else { /* file 'alter.inc', so use it instead of the 'alter' tool */
#ifndef __MINGW32__
fseek(fspice_tmp, -5, SEEK_END); /*properly position the pointer*/
#else
fseek(fspice_tmp, -6, SEEK_END); /*properly position the pointer*/
#endif
fprintf(fspice_tmp, ".INCLUDE alter.inc\n"); /*write*/
fclose(fspice_log);
}
/* fseek(fspice_tmp, 0, SEEK_END); */ /*properly position the pointer*/
fprintf(fspice_tmp, ".end\n"); /*add ".end"*/
fclose(fspice_tmp);
AlterMC=AlterMC-2; /*decrease by two; it will not execute a Monte Carlo simulation next time if equal to zero*/
break;
default:
printf("errfunc.c - Step7 -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
}
/**/
/*Step8: return maximum cost (ALTER on MONTECARLO)*/
#ifdef DEBUG
printf("DEBUG: errfunc.c - Step8\n");
#endif
return (maxcost);
}
void WriteToMem(int num_measures)
{
int i, ii, j;
char laux[LONGSTRINGSIZE];
double currentcost;
i=1; /*It must start with 1*/
while (measure[i].search[0] != '\0') { /*LTspice special case*/
switch(spice) {
case 1: /*Eldo*/
break;
case 2: /*HSPICE*/
break;
case 3: /*LTspice*/
/*Str2Lower(laux);*/
ii=strpos2(measure[i].data, "=", 1); /* LTspice measurement format is different than Eldo and HSPICE */
if (ii) { /* so specific code is added in here to cope with the difference */
strsub(laux, measure[i].data, ii+1, LONGSTRINGSIZE); /* */
strcpy(measure[i].data, laux); /* */
if (strpos2(laux, "dB,", 1)) { /*Is the measurement a complex AC signal?*/
ii=1;
if (strpos2(measure[i].search, "phase", 1))
ReadSubKey(laux, measure[i].data, &ii, ',', 0xB0, 0); /*measure the AC phase */
else
ReadSubKey(laux, measure[i].data, &ii, '(', 'd', 0); /*measure the AC magnitude*/
strcpy(measure[i].data, laux);
}
}
break;
case 4: /*Spectre*/
break;
case 50: /*Qucs*/
break;
case 51: /*ngspice*/
break;
case 100: /*general*/
break;
default:
printf("errfunc.c - WriteToMem -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
i++;
}
DoMath(num_measures); /*'MATH=&...'; information is in the variable 'measured_data'*/
i=1; /*It must start with 1*/
while (measure[i].search[0] != '\0') {
strcpy(laux, UNIQUECHAR);
ii=0;
ii= strpos2(measure[i].var_name, laux, 1); /*if UNIQUECHAR is inexisting in measure[i].var_name and in Str2Lower(UNIQUECHAR) */
Str2Lower(laux); /*as well, this means we are processing a "MEASURE_VAR" line with intermediate data. */
ii=ii + strpos2(measure[i].var_name, laux, 1); /*ALL OTHERS MUST enter the following 'if' to have its measure[i].data processed */
if (ii >1) {
ii=0; /*UNIQUECHAR must be the first character. If not, this is an intermediate "MEASURE_VAR"*/
}
if (ii) {
strcpy(laux, UNIQUECHAR);
ii=(int)strlen(laux);
strsub(laux, measure[i].var_name, ii+1, LONGSTRINGSIZE); /*1- find output measure variable */
/* by copying the last part of the string starting from UNIQUECHAR*/
j=0;
while (strcasecmp(laux, measurements[j].meas_symbol)) { /*2- get its order in measurements[i] */
if (measurements[j].meas_symbol[0] == '\0') { /*if symbol is not found*/
printf("errfunc.c - WriteToMem -- Symbol from SPICE simulation not found in memory: %s\n", laux);
exit(EXIT_FAILURE);
}
j++;
}
if ((int)strlen(measure[i].data)) {
strcpy(laux, measure[i].data); /*3- read measurement */
StripSpaces(laux);
} else {
#ifdef DEBUG
i--; /* to keep the same code from initialize.c */
sprintf(lkk, "%i", i); /* to remove integer added in (*) */
ii=(int)strlen(lkk); /* to remove integer added in (*) */
strcpy(lkk, measurements[i].meas_symbol); /* to remove integer added in (*) */
lkk[(int)strlen(lkk)-ii]='\0'; /* to remove integer added in (*) */
printf("errfunc.c - WriteToMem -- Data not read for measurement in *.cfg: %s\n", lkk);
fflush(stdout); /*Sometimes the simulation does not converge, and so do not exit*/
i++;
#endif
}
measurements[j].measured_value=asc2real(laux, 1, (int)strlen(laux)); /*4- convert it to double */
if (!fcmp(measurements[j].measured_value, 0)) { /*5- check NaN and other text strings */
if ((laux[0] < 43) || (laux[0] > 57) ) /*if its text*/
if (measurements[j].objective_constraint == 1) /*1=MIN*/
measurements[j].measured_value=+1e+30; /*so that a large cost is later assigned*/
if (measurements[j].objective_constraint == 2) /*2=MAX*/
measurements[j].measured_value=-1e+30; /*so that a large cost is later assigned*/
if (measurements[j].objective_constraint == 4) /*4=LE*/
measurements[j].measured_value=+1e+30; /*so that a large cost is later assigned*/
if (measurements[j].objective_constraint > 4) /*5=GE, 6=EQ*/
measurements[j].measured_value=-1e+30; /*so that a large cost is later assigned*/
}
}
i++;
}
/*----------------------------------------------------------------------------*/
/**/
/*Step5: find cost*/
currentcost=CostFunction();
if (maxcost<currentcost)
maxcost=currentcost;
/**/
/*Step6: save log information to file <hostname>.log*/
if (LOG) {
LogtoFile(currentcost);
}
/*----------------------------------------------------------------------------*/
}
/*
* update the *.lis file with the state of the transistor
*/
void UpdateLIS(char *ConfigFile, char *InputFile)
{
int i, j, k;
char ltitle[LONGSTRINGSIZE]; /*title: found .alter @*/
char lelement[LONGSTRINGSIZE];
/*lmodel,*/ /*transistor/model name*/
char lVGS[LONGSTRINGSIZE], lVDS[LONGSTRINGSIZE], lVth[LONGSTRINGSIZE], lVDSAT[LONGSTRINGSIZE], fileLJR[LONGSTRINGSIZE];
FILE *fLJR; /**.ljr*/
FILE *fLIS, *fNoSat;
FILE *fsweepINI;
/*skip: array [1..100] of string [20];*/ /*100 transistor to skip*/
char laux[LONGSTRINGSIZE];
double Vovd, Voff, Vdst; /*overdrive, off and 'Vds-Vdsat' voltage read from sweep.ini file*/
ThreeLines stats;
fNoSat = NULL;
fLIS = NULL;
fLJR = NULL;
if ((fLIS=fopen(InputFile,"rt")) == 0) {
printf("auxfunc_updatelis.c - Cannot open input file: %s\n", InputFile);
exit(EXIT_FAILURE);
}
sprintf(fileLJR, "%.*sjr", (int)(strlen(InputFile) - 2), InputFile);
if ((fLJR=fopen(fileLJR,"wt")) == 0) {
printf("auxfunc_updatelis.c - Cannot open input file: %s\n", fileLJR);
exit(EXIT_FAILURE);
}
if ((fNoSat=fopen(NoSat,"wt")) == 0) {
printf("auxfunc_updatelis.c - Cannot open input file: %s\n", NoSat);
exit(EXIT_FAILURE);
}
for (i = 0; i < SKIPTRAN; i++) /*INITIALIZE: read transistors to skip*/
*skip[i] = '\0';
if ((fsweepINI=fopen(ConfigFile,"rt")) == 0) {
printf("auxfunc_updatelis.c - Cannot open input file: %s\n", ConfigFile);
exit(EXIT_FAILURE);
}
ReadKey(lkk, "SKIP_NOSAT", fsweepINI);
j = strpos2(lkk, "$", 1);
if (j != 0) /*This will skip the characters after '$', the inline comment used by the sweep tools*/
lkk[j-1]='\0';
StripSpaces(lkk);
k = 1;
if (!lkk[0])
printf("INFO: auxfunc_updatelis.c - 'SKIP_NOSAT' not found\n");
else {
while (!strcmp((sprintf(laux, "%.10s", lkk), laux), "SKIP_NOSAT")) { /*read SKIP_NOSAT from file*/
i = 1;
lkk[10] = ';';
while (i < (int)strlen(lkk)) {
ReadSubKey(laux, lkk, &i, ';', ';', 0);
StripSpaces(laux);
strcpy(skip[k - 1], laux);
k++;
if (k == SKIPTRAN) {
printf("auxfunc_updatelis.c - Maximum number of %d transistors to skip reached. Increase SKIPTRAN in auxfunc_updatelis\n", SKIPTRAN);
exit(EXIT_FAILURE);
}
}
ReadKey(lkk, "SKIP_NOSAT", fsweepINI);
j = strpos2(lkk, "$", 1);
if (j != 0) /*This will skip the characters after '$', the inline comment used by the sweep tools*/
lkk[j-1]='\0';
StripSpaces(lkk);
}
}
Vovd = 0.05; /*default OVERDRIVE VOLTAGE value is 50mV*/
fseek(fsweepINI, 0, SEEK_SET);
ReadKey(lkk, "VOVD", fsweepINI);
if (strcmp((sprintf(laux, "%.4s", lkk), laux), "VOVD"))
printf("INFO: auxfunc_updatelis.c - No Vovd, default value (50mV) will be used instead\n");
else {
i = 1;
i = (sscanf(ReadSubKey(laux, lkk, &i, ':', 'm', 5), "%lg", &Vovd) == 0);
if (i) { /*if (i=0) then laux contains a number and not text*/
printf("auxfunc_updatelis.c - Incorrect format: number must exist in '%s'\n", lkk);
exit(EXIT_FAILURE);
}
Vovd /= 1000;
/*the following line will print information about the symbols in 'fNoSat' file*/
fprintf(fNoSat, "` `:weak inversion if Vgs<Vth+%smV\n", laux);
}
Voff = 0.1; /*default OFF VOLTAGE value is 100mV*/
fseek(fsweepINI, 0, SEEK_SET);
ReadKey(lkk, "VOFF", fsweepINI);
if (strcmp((sprintf(laux, "%.4s", lkk), laux), "VOFF"))
printf("INFO: auxfunc_updatelis.c - No Voff: default value (100mV) will be used instead\n");
else {
i = 1;
i = (sscanf(ReadSubKey(laux, lkk, &i, ':', 'm', 5), "%lg", &Voff) == 0);
if (i) { /*if (i=0) then laux contains a number and not text*/
printf("auxfunc_updatelis.c - Incorrect format: number must exist in '%s'\n", lkk);
exit(EXIT_FAILURE);
}
Voff /= 1000;
}
Vdst = 0.0; /*default DST VOLTAGE value is 000mV*/
fseek(fsweepINI, 0, SEEK_SET);
ReadKey(lkk, "VDST", fsweepINI);
if (strcmp((sprintf(laux, "%.4s", lkk), laux), "VDST"))
printf("INFO: auxfunc_updatelis.c - No Vdst: default value (0mV) will be used instead\n");
else {
i = 1;
i = (sscanf(ReadSubKey(laux, lkk, &i, ':', 'm', 5), "%lg", &Vdst) == 0); /*read 'Vds-Vdsat' from file*/
if (i) {/*if (i=0) then laux contains a number and not text*/
printf("auxfunc_updatelis.c - Incorrect format: number must exist in '%s'\n", lkk);
exit(EXIT_FAILURE);
}
Vdst /= 1000;
/*the following lines will print information about the symbols in 'fNoSat' file*/
if (Vdst != 0)
fprintf(fNoSat, "`#`:Vds<Vdsat+%smV\n", laux);
fprintf(fNoSat, "`*`:Vds<Vdsat\n\n");
}
if (fsweepINI != NULL)
fclose(fsweepINI);
fsweepINI = NULL;
/* */
/* */
*ltitle = '\0';
while (!P_eof(fLIS)) {
switch(spice) {
case 1: /* Eldo */
while ((strcmp((sprintf(laux, "%.37s", lkk), laux), "0**** OPERATING POINT") != 0) & (!P_eof(fLIS))) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
StripSpaces(lkk); /*required due to Solaris OS*/
//if (lkk[0] == '@' && strcmp(lkk, ltitle)) {
if (strpos2(lkk, " .ALTER @", 1) !=0) {
fprintf(fNoSat, "%s | ", lkk);
strcpy(ltitle, lkk);
}
}
sprintf(laux, "%.37s", lkk);
if (!strcmp(laux, "0**** OPERATING POINT")) {
for (i = 1; i <= 6; i++) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
}
while ( (strpos2(lkk, " M", 1) != 0) || (strpos2(lkk, " X", 1) != 0)) { /*find operating region for all transistors*/
strcpy(lelement, lkk);
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
//if (strpos2(lkk, "region", 1) !=0) { /* Due to HSPICE 2001.2 line */
// fgets2(lkk, LONGSTRINGSIZE, fLIS);
// fprintf(fLJR, "%s\n", lkk); /* with the operation region */
//}
/* fgets2(lkk, LONGSTRINGSIZE, fLIS); */
/* fprintf(fLJR, "%s\n", lkk); */
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
fgets2(lVGS, LONGSTRINGSIZE, fLIS); /*Vgs*/
fprintf(fLJR, "%s\n", lVGS);
fgets2(lVDS, LONGSTRINGSIZE, fLIS); /*Vds*/
fprintf(fLJR, "%s\n", lVDS);
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
fgets2(lVth, LONGSTRINGSIZE, fLIS); /*Vth*/
fprintf(fLJR, "%s\n", lVth);
fgets2(lVDSAT, LONGSTRINGSIZE, fLIS); /*Vdsat*/
fprintf(fLJR, "%s\n", lVDSAT);
for (i = 1; i <= 24; i++) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
}
DoCalculations(lelement, lVGS, lVDS, lVth, lVDSAT, Vovd, Voff, Vdst, stats, &fNoSat); /*gets three lines with operating region*/
fprintf(fLJR, "%s\n", stats[0]);
fprintf(fLJR, "%s\n", stats[1]);
fprintf(fLJR, "%s\n\n\n", stats[2]);
for (i = 1; i <= 2; i++) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
if (i > 1)
fprintf(fLJR, "%s\n", lkk);
}
}
putc('\n', fNoSat);
}
break;
case 2: /* HSPICE */
i=0;
while ((strcmp((sprintf(laux, "%.12s", lkk), laux), "**** mosfets") != 0) & (!P_eof(fLIS))) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
StripSpaces(lkk); /*required due to Solaris OS*/
if (lkk[0] == '@' && strcmp(lkk, ltitle)) {
if (i>0) {
fprintf(fNoSat, "\n");
}
i++;
fprintf(fNoSat, "%s | ", lkk);
strcpy(ltitle, lkk);
}
}
sprintf(laux, "%.12s", lkk);
if (!strcmp(laux, "**** mosfets")) {
for (i = 1; i <= 4; i++) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
}
while (strpos2(lkk, "element ", 1) != 0) { /*find operating region for all transistors*/
strcpy(lelement, lkk);
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
if (strpos2(lkk, "region", 1) !=0) { /* Due to HSPICE 2001.2 line */
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk); /* with the operation region */
}
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
fgets2(lVGS, LONGSTRINGSIZE, fLIS); /*Vgs*/
fprintf(fLJR, "%s\n", lVGS);
fgets2(lVDS, LONGSTRINGSIZE, fLIS); /*Vds*/
fprintf(fLJR, "%s\n", lVDS);
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
fgets2(lVth, LONGSTRINGSIZE, fLIS); /*Vth*/
fprintf(fLJR, "%s\n", lVth);
fgets2(lVDSAT, LONGSTRINGSIZE, fLIS); /*Vdsat*/
fprintf(fLJR, "%s\n", lVDSAT);
for (i = 1; i <= 12; i++) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
}
DoCalculations(lelement, lVGS, lVDS, lVth, lVDSAT, Vovd, Voff, Vdst, stats, &fNoSat); /*gets three lines with operating region*/
fprintf(fLJR, "%s\n", stats[0]);
fprintf(fLJR, "%s\n", stats[1]);
fprintf(fLJR, "%s\n\n\n", stats[2]);
for (i = 1; i <= 4; i++) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
if (i > 2)
fprintf(fLJR, "%s\n", lkk);
}
}
putc('\n', fNoSat);
}
break;
case 3: /* LTspice */
while ((strcmp((sprintf(laux, "%.21s", lkk), laux), "--- BSIM3 MOSFETS ---") != 0) & (!P_eof(fLIS))) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
StripSpaces(lkk); /*required due to Solaris OS*/
if (lkk[0] == '@' && strcmp(lkk, ltitle)) {
fprintf(fNoSat, "%s | ", lkk);
strcpy(ltitle, lkk);
}
}
sprintf(laux, "%.21s", lkk);
if (!strcmp(laux, "--- BSIM3 MOSFETS ---")) {
for (i = 1; i <= 1; i++) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
}
while (strpos2(lkk, "Name: ", 1) != 0) { /*find operating region for all transistors*/
strcpy(lelement, lkk);
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
//fgets2(lkk, LONGSTRINGSIZE, fLIS);
//fprintf(fLJR, "%s\n", lkk);
//if (strpos2(lkk, "region", 1) !=0) { /* Due to HSPICE 2001.2 line */
// fgets2(lkk, LONGSTRINGSIZE, fLIS);
// fprintf(fLJR, "%s\n", lkk); /* with the operation region */
//}
//fgets2(lkk, LONGSTRINGSIZE, fLIS);
//fprintf(fLJR, "%s\n", lkk);
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
fgets2(lVGS, LONGSTRINGSIZE, fLIS); /*Vgs*/
fprintf(fLJR, "%s\n", lVGS);
fgets2(lVDS, LONGSTRINGSIZE, fLIS); /*Vds*/
fprintf(fLJR, "%s\n", lVDS);
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
fgets2(lVth, LONGSTRINGSIZE, fLIS); /*Vth*/
fprintf(fLJR, "%s\n", lVth);
fgets2(lVDSAT, LONGSTRINGSIZE, fLIS); /*Vdsat*/
fprintf(fLJR, "%s\n", lVDSAT);
for (i = 1; i <= 17; i++) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
}
DoCalculations(lelement, lVGS, lVDS, lVth, lVDSAT, Vovd, Voff, Vdst, stats, &fNoSat); /*gets three lines with operating region*/
fprintf(fLJR, "%s\n", stats[0]);
fprintf(fLJR, "%s\n", stats[1]);
fprintf(fLJR, "%s\n\n", stats[2]);
for (i = 1; i <= 2; i++) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
if (i > 1)
fprintf(fLJR, "%s\n", lkk);
}
}
putc('\n', fNoSat);
}
break;
case 4: /* Spectre */
// in input.scs put "dcOp dc oppoint=logfile"
// run spectre with "spectre input.scs =log input.lis"
//
// Code specific for "Primitive: bsim3v3" and is quite possible
// not to work on future Spectre versions and other transistor models
//
while ((strcmp((sprintf(laux, "%.31s", lkk), laux), "operating-point information `dc") != 0) & (!P_eof(fLIS))) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
StripSpaces(lkk); /*required due to Solaris OS*/
if (lkk[0] == '@' && strcmp(lkk, ltitle)) {
fprintf(fNoSat, "%s | ", lkk);
strcpy(ltitle, lkk);
}
Str2Lower(lkk); //required because versions until 2007 have "... Information" and after have "... information"
}
if (P_eof(fLIS))
break;
char l1[LONGSTRINGSIZE], l2[LONGSTRINGSIZE]; /*l1 and l2 are hold places for previous lines*/
fgets2(l1, LONGSTRINGSIZE, fLIS); fprintf(fLJR, "%s\n", l1);
fgets2(l2, LONGSTRINGSIZE, fLIS); fprintf(fLJR, "%s\n", l2);
/*bsim3v3*/ while ((strcmp((sprintf(laux, "%.18s", lkk), laux), "Primitive: bsim3v3") != 0) & (!P_eof(fLIS))) {
/*bsim4*/ // while ((strcmp((sprintf(laux, "%.16s", lkk), laux), "Primitive: bsim4") != 0) & (!P_eof(fLIS))) {
strcpy(l1, l2);
strcpy(l2, lkk);
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
StripSpaces(lkk);
if (lkk[0] == '@' && strcmp(lkk, ltitle)) {
fprintf(fNoSat, "%s | ", lkk);
strcpy(ltitle, lkk);
}
}
/*bsim3v3*/ sprintf(laux, "%.18s", lkk);
/*bsim3v3*/ if (!strcmp(laux, "Primitive: bsim3v3")) {
/*bsim4*/ // sprintf(laux, "%.16s", lkk);
/*bsim4*/ // if (!strcmp(laux, "Primitive: bsim4")) {
for (i = 1; i <= 0; i++) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
}
/*bsim3v3*/ while (strpos2(lkk, "Primitive: bsim3v3", 1) != 0) { /*find operating region for all transistors*/
/*bsim4*/ // while (strpos2(lkk, "Primitive: bsim4", 1) != 0) { /*find operating region for all transistors*/
strcpy(lelement, l1);
if (strpos2(lelement, " of ", 1)) { /*to remove the library name if and when it appears*/
lelement[strpos2(lelement, " of ", 1)]='\0';
}
/*bsim3v3*/ for (i = 1; i <= 9; i++) {
/*bsim4*/ // for (i = 1; i <= 7; i++) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
}
fgets2(lVGS, LONGSTRINGSIZE, fLIS); /*Vgs*/
fprintf(fLJR, "%s\n", lVGS);
{
i=1; ReadSubKey(laux, lVGS, &i, '=', 'V', 5);
i=strpos2(laux, " ", 2);
if ((laux[i]==109)) { /* 109=m; mV */
laux[i-1]='\0';
strcat(laux, "e-3");
}
if ((laux[i]==117)) { /* 117=u; uV */
laux[i-1]='\0';
strcat(laux, "e-6");
}
if ((laux[i]==110)) { /* 110=n; nV */
laux[i-1]='\0';
strcat(laux, "e-9");
}
sprintf(lVGS, " %s", laux); /* Puts the value on the expected column for DoCalculations */
}
fgets2(lVDS, LONGSTRINGSIZE, fLIS); /*Vds*/
fprintf(fLJR, "%s\n", lVDS);
{
i=1; ReadSubKey(laux, lVDS, &i, '=', 'V', 5);
i=strpos2(laux, " ", 2);
if ((laux[i]==109)) { /* 109=m; mV */
laux[i-1]='\0';
strcat(laux, "e-3");
}
if ((laux[i]==117)) { /* 117=u; uV */
laux[i-1]='\0';
strcat(laux, "e-6");
}
if ((laux[i]==110)) { /* 110=n; nV */
laux[i-1]='\0';
strcat(laux, "e-9");
}
sprintf(lVDS, " %s", laux); /* Puts the value on the expected column for DoCalculations */
}
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
fgets2(lVth, LONGSTRINGSIZE, fLIS); /*Vth*/
fprintf(fLJR, "%s\n", lVth);
{ //begin -- required for versions after the 2nd half of 2006 who have 3 more lines in between "vds" and "vth"
if (!strpos2(lVth, "vth =", 1)) {
fgets2(lVth, LONGSTRINGSIZE, fLIS); /*Vth*/
fprintf(fLJR, "%s\n", lVth);
fgets2(lVth, LONGSTRINGSIZE, fLIS); /*Vth*/
fprintf(fLJR, "%s\n", lVth);
fgets2(lVth, LONGSTRINGSIZE, fLIS); /*Vth*/
fprintf(fLJR, "%s\n", lVth);
}
if (!strpos2(lVth, "vth =", 1)) { /*Double-check just to be sure that vth was found*/
printf("auxfunc_updatelis.c - Spectre netlist format not supported\n");
exit(EXIT_FAILURE);
}
} //end -- required for versions after the 2nd half of 2006 who have 3 more lines in between "vds" and "vth"
{
i=1; ReadSubKey(laux, lVth, &i, '=', 'V', 5);
i=strpos2(laux, " ", 2);
if ((laux[i]==109)) { /* 109=m; mV */
laux[i-1]='\0';
strcat(laux, "e-3");
}
if ((laux[i]==117)) { /* 117=u; uV */
laux[i-1]='\0';
strcat(laux, "e-6");
}
if ((laux[i]==110)) { /* 110=n; nV */
laux[i-1]='\0';
strcat(laux, "e-9");
}
sprintf(lVth, " %s", laux); /* Puts the value on the expected column for DoCalculations */
}
fgets2(lVDSAT, LONGSTRINGSIZE, fLIS); /*Vdsat*/
fprintf(fLJR, "%s\n", lVDSAT);
{
i=1; ReadSubKey(laux, lVDSAT, &i, '=', 'V', 5);
i=strpos2(laux, " ", 2);
if ((laux[i]==109)) { /* 109=m; mV */
laux[i-1]='\0';
strcat(laux, "e-3");
}
if ((laux[i]==117)) { /* 117=u; uV */
laux[i-1]='\0';
strcat(laux, "e-6");
}
if ((laux[i]==110)) { /* 110=n; nV */
laux[i-1]='\0';
strcat(laux, "e-9");
}
sprintf(lVDSAT, " %s", laux); /* Puts the value on the expected column for DoCalculations */
}
for (i = 1; i <= 45; i++) {
fgets2(lkk, LONGSTRINGSIZE, fLIS);
fprintf(fLJR, "%s\n", lkk);
}
{ //begin -- required for versions after the 2nd half of 2006 who have 3 more lines in between "vds" and "vth"
fgets2(lkk, LONGSTRINGSIZE, fLIS);
//fprintf(fLJR, "%s\n", lkk);
//if ((int)(strlen(InputFile))) { /*empty line has been found, goes ahead*/
// i=0;
//}
while (strpos2(lkk, " = ", 1)) { /*empty line has been found, goes ahead*/
fprintf(fLJR, "%s\n", lkk);
fgets2(lkk, LONGSTRINGSIZE, fLIS);
}
} //end -- required for versions after the 2nd half of 2006 who have 3 more lines in between "vds" and "vth"
DoCalculations(lelement, lVGS, lVDS, lVth, lVDSAT, Vovd, Voff, Vdst, stats, &fNoSat); /*gets three lines with operating region*/
fprintf(fLJR, "%s\n", stats[0]);
fprintf(fLJR, "%s\n", stats[1]);
fprintf(fLJR, "%s\n\n", stats[2]);
for (i = 1; i <= 3; i++) {
strcpy(l1, l2);
strcpy(l2, lkk);
fgets2(lkk, LONGSTRINGSIZE, fLIS);
if (i > 0)
fprintf(fLJR, "%s\n", lkk);
}
}
putc('\n', fNoSat);
}
break;
case 50: /*Qucs*/
printf("auxfunc_updatelis.c - UpdateLIS -- Updatelis not implemente for Qucs\n");
exit(EXIT_FAILURE);
case 100: /* rosen */
printf("auxfunc_updatelis.c - UpdateLIS -- Updatelis not implemente for rosen\n");
exit(EXIT_FAILURE);
break;
default:
printf("auxfunc_updatelis.c - UpdateLIS -- Something unexpected has happened!\n");
exit(EXIT_FAILURE);
}
}
if (fLIS != NULL)
fclose(fLIS);
fLIS = NULL;
if (fLJR != NULL)
fclose(fLJR);
fLJR = NULL;
if (fNoSat != NULL)
fclose(fNoSat);
fNoSat = NULL;
if (fLJR != NULL)
fclose(fLJR);
if (fLIS != NULL)
fclose(fLIS);
if (fNoSat != NULL)
fclose(fNoSat);
} /*UpdateLIS*/
void CreateStatistics(char *InputFile, char *OutputFile)
{
int i, j, k;
char lkk1[LONGSTRINGSIZE], llog[LONGSTRINGSIZE];
statistics stats;
FILE *fIn, *fOut, *fcfg;
double aux;
int num_measures;
/* ------------------------------------------------------------------
------------------------------------------------------------------ */
/*create summary: table version*/
for (i = 0; i <= (MAXMEAS-1); i++) {
stats.avg[i] = 0.0; /*initialization of statistics variables*/
stats.sig[i] = 0.0;
stats.max[i] = -1.7e38;
stats.min[i] = 1.7e38;
}
if ((fIn=fopen(InputFile,"rt")) == 0) {
printf("auxfunc_log.c - Cannot open input file: %s\n", InputFile);
exit(EXIT_FAILURE);
}
if ((fOut=fopen(OutputFile,"wt")) == 0) {
printf("auxfunc_log.c - Cannot open output file: %s\n", OutputFile);
exit(EXIT_FAILURE);
}
fgets2(lkk, LONGSTRINGSIZE, fIn);
if (P_eof(fIn)) /*if there is only 1 line*/
fseek(fIn, 0, SEEK_SET);
strsub(lkk1, lkk, 1, strpos2(lkk, ":", 1)-1);
i=0; j=1; k=0;
while (!P_eof(fIn)) {
aux = asc2real(lkk1, 1, (int)strlen(lkk1));
stats.avg[i] += aux;
stats.sig[i] += aux * aux;
/*variancia=sum(x^2)/n - (avg(x))^2 # sigma=sqrt(variancia)*/
if (aux > stats.max[i])
stats.max[i] = aux;
if (aux < stats.min[i])
stats.min[i] = aux;
if (j==(int)strlen(lkk)) {
num_measures=i;
i=0;j=1;k++;
fgets2(lkk, LONGSTRINGSIZE, fIn);
strsub(lkk1, lkk, 1, strpos2(lkk, ":", 1)-1);
} else {
ReadSubKey(lkk1, lkk, &j, ':', ':', 1);
i++; j--;
}
}
WriteStats(num_measures-1, k, stats, &fOut);
/* -- */
i=0;
dd_ffblk fb;
char *mask="*.cfg";
if (!dd_findfirst( mask, &fb, DD_DIREC )) { /*Find first file having extension .cfg */
if (!strcmp(fb.dd_name, "rfmodule.cfg")) /*and make sure that it is not rfmodule.cfg*/
i=dd_findnext(&fb); /*otherwise find next .cfg file */
if (i==0)
printf("Opening config file: %s\n", fb.dd_name);
}
/* -- */
fprintf(fOut, "\n\n# Parameters #\n");
/*--------*/
strcpy(llog,lkk); /*last line from log file*/
fcfg=NULL;
if ((fcfg=fopen(fb.dd_name,"rt")) == 0) {
if (fb.dd_name) {
printf("auxfunc_log.c - Cannot open config file: %s\n", fb.dd_name);
/*exit(EXIT_FAILURE);*/
}
SimpleParametersCategory(num_measures-1, llog, stats, &fOut);
} else {
ReadKey(lkk, "# Parameters #", fcfg); /*configuration parameters*/
if (strcmp(lkk, "# Parameters #")) {
printf("INFO: auxfunc_log.c - Step2 -- No parameters in config file\n");
SimpleParametersCategory(num_measures-1, llog, stats, &fOut);
} else {
ComplexParametersCategory(llog, stats, &fOut, &fcfg);
}
}
/*--------*/
fprintf(fOut, "#\n");
if (fOut != NULL)
fclose(fOut);
fOut = NULL;
if (fIn != NULL)
fclose(fIn);
fIn = NULL;
/* ------------------------------------------------------------------
------------------------------------------------------------------ */
}
