version 1.3, 1999/12/16 22:18:35
|
version 1.7, 2000/06/30 21:36:16
|
Line 1
|
Line 1
|
|
/* definition of all capa functions |
|
Copyright (C) 1992-2000 Michigan State University |
|
|
|
The CAPA system is free software; you can redistribute it and/or |
|
modify it under the terms of the GNU Library General Public License as |
|
published by the Free Software Foundation; either version 2 of the |
|
License, or (at your option) any later version. |
|
|
|
The CAPA system is distributed in the hope that it will be useful, |
|
but WITHOUT ANY WARRANTY; without even the implied warranty of |
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
|
Library General Public License for more details. |
|
|
|
You should have received a copy of the GNU Library General Public |
|
License along with the CAPA system; see the file COPYING. If not, |
|
write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, |
|
Boston, MA 02111-1307, USA. */ |
|
|
/* =||>>================================================================<<||= */ |
/* =||>>================================================================<<||= */ |
/* 45678901234567890123456789012345678901234567890123456789012345678901234567 */ |
/* 45678901234567890123456789012345678901234567890123456789012345678901234567 */ |
Line 14
|
Line 31
|
#include "capaCommon.h" |
#include "capaCommon.h" |
#include "ranlib.h" |
#include "ranlib.h" |
|
|
|
|
char Parse_class[QUARTER_K]; |
char Parse_class[QUARTER_K]; |
int Parse_set; |
int Parse_set; |
int Parse_section; |
int Parse_section; |
Line 33 extern int Current_line[MAX_OPEN
|
Line 49 extern int Current_line[MAX_OPEN
|
extern int Func_idx; |
extern int Func_idx; |
extern Symbol FuncStack[MAX_FUNC_NEST]; |
extern Symbol FuncStack[MAX_FUNC_NEST]; |
|
|
|
#ifdef TTH |
|
extern int textohtmldyn(char*,char**,char**,int); |
|
char *tth_err; |
|
#endif |
|
|
/* --------------------------------------------------------------------------- */ |
/* --------------------------------------------------------------------------- */ |
int |
int |
match_function(func, argc) char *func; int argc; |
match_function(func, argc) char *func; int argc; |
{ |
{ |
if( !strcmp(func,"random") ) return (((argc==2 || argc==3)? RANDOM_F : MIS_ARG_COUNT)); |
if( !strcmp(func,"random") ) return (((argc==2 || argc==3)? RANDOM_F : MIS_ARG_COUNT)); |
|
if( !strcmp(func,"random_normal") ) return ((argc==5)? RANDOM_NORMAL_F : MIS_ARG_COUNT); |
|
if( !strcmp(func,"random_beta") ) return ((argc==5)? RANDOM_BETA_F : MIS_ARG_COUNT); |
|
if( !strcmp(func,"random_gamma") ) return ((argc==5)? RANDOM_GAMMA_F : MIS_ARG_COUNT); |
|
if( !strcmp(func,"random_poisson") ) return ((argc==4)? RANDOM_POISSON_F : MIS_ARG_COUNT); |
|
if( !strcmp(func,"random_exponential") ) return ((argc==4)? RANDOM_EXPONENTIAL_F : MIS_ARG_COUNT); |
|
if( !strcmp(func,"random_chi") ) return ((argc==4)? RANDOM_CHI_F : MIS_ARG_COUNT); |
|
if( !strcmp(func,"random_noncentral_chi") ) return ((argc==5)? RANDOM_NONCENTRAL_CHI_F : MIS_ARG_COUNT); |
if( !strcmp(func,"choose") ) return (CHOOSE_F); |
if( !strcmp(func,"choose") ) return (CHOOSE_F); |
if( !strcmp(func,"tex") ) return (((argc==2)? TEX_F: MIS_ARG_COUNT)); |
if( !strcmp(func,"tex") ) return (((argc==2)? TEX_F: MIS_ARG_COUNT)); |
if( !strcmp(func,"var_in_tex") ) return (VAR_IN_TEX_F); |
if( !strcmp(func,"var_in_tex") ) return (VAR_IN_TEX_F); |
Line 60 match_function(func, argc) char *func; i
|
Line 88 match_function(func, argc) char *func; i
|
if( !strcmp(func,"strlen") ) return (((argc==1)? STRLEN_F: MIS_ARG_COUNT)); |
if( !strcmp(func,"strlen") ) return (((argc==1)? STRLEN_F: MIS_ARG_COUNT)); |
if( !strcmp(func,"get_seed") ) return (((argc==0)? GET_SEED_F: MIS_ARG_COUNT)); |
if( !strcmp(func,"get_seed") ) return (((argc==0)? GET_SEED_F: MIS_ARG_COUNT)); |
if( !strcmp(func,"set_seed") ) return (((argc==1)? SET_SEED_F: MIS_ARG_COUNT)); |
if( !strcmp(func,"set_seed") ) return (((argc==1)? SET_SEED_F: MIS_ARG_COUNT)); |
|
if( !strcmp(func,"init_array") ) return (((argc==1)? INIT_ARRAY_F: MIS_ARG_COUNT)); |
if( !strcmp(func,"array_index") ) return (((argc==1)? ARRAY_INDEX_F: MIS_ARG_COUNT)); |
if( !strcmp(func,"array_index") ) return (((argc==1)? ARRAY_INDEX_F: MIS_ARG_COUNT)); |
if( !strcmp(func,"array_sorted_index") ) return (((argc==2)? ARRAY_SORTED_INDEX_F: MIS_ARG_COUNT)); |
if( !strcmp(func,"array_sorted_index") ) return (((argc==2)? ARRAY_SORTED_INDEX_F: MIS_ARG_COUNT)); |
if( !strcmp(func,"array_max") ) return (((argc==1)? ARRAY_MAX_F: MIS_ARG_COUNT)); |
if( !strcmp(func,"array_max") ) return (((argc==1)? ARRAY_MAX_F: MIS_ARG_COUNT)); |
Line 275 ArgNode_t *argp;
|
Line 304 ArgNode_t *argp;
|
capa_msg(MESSAGE_ERROR,tmpS); |
capa_msg(MESSAGE_ERROR,tmpS); |
} |
} |
} break; |
} break; |
|
|
case CHOOSE_F: { int ii, pick=1; |
case CHOOSE_F: { int ii, pick=1; |
ArgNode_t *tmpArgp; |
ArgNode_t *tmpArgp; |
|
|
Line 746 ArgNode_t *argp;
|
Line 776 ArgNode_t *argp;
|
resultp->s_type = I_CONSTANT; |
resultp->s_type = I_CONSTANT; |
resultp->s_int = 0; |
resultp->s_int = 0; |
} break; |
} break; |
case ARRAY_MOMENTS_F: /* it */ |
/* generate random numbers according to a pre-defined distributions and a seed */ |
|
case RANDOM_NORMAL_F: /* random_normal(return_array,item_cnt,seed,av,std_dev) */ |
|
case RANDOM_BETA_F: /* random_beta(return_array,item_cnt,seed,aa,bb) */ |
|
case RANDOM_GAMMA_F: /* random_gamma(return_array,item_cnt,seed,a,r) */ |
|
case RANDOM_POISSON_F: /* random_poisson(return_array,item_cnt,seed,mu) */ |
|
case RANDOM_EXPONENTIAL_F: |
|
/* random_exponential(return_array,item_cnt,seed,av) */ |
|
case RANDOM_CHI_F: /* random_chi(return_array,item_cnt,seed,df) */ |
|
case RANDOM_NONCENTRAL_CHI_F: |
|
/* random_noncentral_chi(return_array,item_cnt,seed,df,xnonc) */ |
|
/* gen_random_by_selector(output_p,sel,seed,item_cnt,p1,p2) */ |
|
{ int sel, item_cnt, tmp_int; |
|
float para1, para2; |
|
char *tmp_str; |
|
long tmp_long; |
|
Symbol *r_p; |
|
|
|
switch(func) { /* assigns the function selector */ |
|
case RANDOM_NORMAL_F: sel = NORMAL_DIS; break; |
|
case RANDOM_BETA_F: sel = BETA_DIS; break; |
|
case RANDOM_GAMMA_F: sel = GAMMA_DIS; break; |
|
case RANDOM_POISSON_F: sel = POISSON_DIS; break; |
|
case RANDOM_EXPONENTIAL_F: sel = EXPONENTIAL_DIS; break; |
|
case RANDOM_CHI_F: sel = CHI_DIS; break; |
|
case RANDOM_NONCENTRAL_CHI_F: sel = NONCENTRAL_CHI_DIS; break; |
|
} |
|
switch(func) { |
|
case RANDOM_NORMAL_F: |
|
case RANDOM_BETA_F: |
|
case RANDOM_GAMMA_F: /* two-parameter functions */ |
|
case RANDOM_NONCENTRAL_CHI_F: |
|
{ errCode = 0; |
|
switch( FIRST_ARGTYPE(argp) ) { /* parameter two */ |
|
case I_VAR: case I_CONSTANT: |
|
para2 = (float)FIRST_ARGINT(argp); |
|
break; |
|
case R_VAR: case R_CONSTANT: |
|
para2 = (float)FIRST_ARGREAL(argp); |
|
break; |
|
case S_VAR: case S_CONSTANT: |
|
case IDENTIFIER: |
|
resultp->s_type = S_CONSTANT; |
|
resultp->s_str = strsave("<<LAST ARG. OF THIS FUNCTION MUST BE A NUMBER>>"); |
|
sprintf(tmpS,"%s()'s last arg. must be a number.\n",FuncStack[Func_idx].s_name); |
|
capa_msg(MESSAGE_ERROR,tmpS); |
|
errCode = 1; |
|
break; |
|
} |
|
if(errCode == 0 ) { |
|
switch( SECOND_ARGTYPE(argp) ) { /* parameter one */ |
|
case I_VAR: case I_CONSTANT: |
|
para1 = (float)SECOND_ARGINT(argp); |
|
break; |
|
case R_VAR: case R_CONSTANT: |
|
para1 = (float)SECOND_ARGREAL(argp); |
|
break; |
|
case S_VAR: case S_CONSTANT: |
|
case IDENTIFIER: |
|
resultp->s_type = S_CONSTANT; |
|
resultp->s_str = strsave("<<FOURTH ARG. OF THIS FUNCTION MUST BE A NUMBER>>"); |
|
sprintf(tmpS,"%s()'s fourth arg. must be a number.\n",FuncStack[Func_idx].s_name); |
|
capa_msg(MESSAGE_ERROR,tmpS); |
|
errCode = 1; |
|
break; |
|
} |
|
if(errCode == 0 ) { |
|
switch( THIRD_ARGTYPE(argp) ) { /* seed */ |
|
case I_VAR: case I_CONSTANT: |
|
tmp_str = (char *)capa_malloc(32,1); |
|
sprintf(tmp_str,"%ld",THIRD_ARGINT(argp) ); |
|
break; |
|
case R_VAR: case R_CONSTANT: |
|
tmp_long = (long)THIRD_ARGREAL(argp); |
|
tmp_str = (char *)capa_malloc(32,1); |
|
sprintf(tmp_str,"%ld",tmp_long); |
|
break; |
|
case S_VAR: case S_CONSTANT: |
|
tmp_str = strsave(THIRD_ARGSTR(argp)); |
|
break; |
|
case IDENTIFIER: |
|
resultp->s_type = S_CONSTANT; |
|
resultp->s_str = strsave("<<THIRD ARG. OF THIS FUNCTION MUST BE A NUMBER>>"); |
|
sprintf(tmpS,"%s()'s third arg. must be a number.\n",FuncStack[Func_idx].s_name); |
|
capa_msg(MESSAGE_ERROR,tmpS); |
|
errCode = 1; |
|
break; |
|
} |
|
if(errCode == 0 ) { |
|
switch( FOURTH_ARGTYPE(argp) ) { /* item_cnt */ |
|
case I_VAR: case I_CONSTANT: |
|
item_cnt = FOURTH_ARGINT(argp); |
|
break; |
|
case R_VAR: case R_CONSTANT: |
|
item_cnt = (int)FOURTH_ARGREAL(argp); |
|
break; |
|
case S_VAR: case S_CONSTANT: |
|
case IDENTIFIER: |
|
resultp->s_type = S_CONSTANT; |
|
resultp->s_str = strsave("<<SECOND ARG. OF THIS FUNCTION MUST BE A NUMBER>>"); |
|
sprintf(tmpS,"%s()'s second arg. must be a number.\n",FuncStack[Func_idx].s_name); |
|
capa_msg(MESSAGE_ERROR,tmpS); |
|
errCode = 1; |
|
break; |
|
} |
|
if(errCode == 0 ) { |
|
switch( FIFTH_ARGTYPE(argp) ) { /* array_name, clear the content of this array first */ |
|
case I_VAR: case I_CONSTANT: |
|
case R_VAR: case R_CONSTANT: |
|
resultp->s_type = S_CONSTANT; |
|
resultp->s_str = strsave("<<FIRST ARG. OF THIS FUNCTION MUST BE AN ARRAY NAME>>"); |
|
sprintf(tmpS,"%s()'s first arg. must be a name of an array.\n",FuncStack[Func_idx].s_name); |
|
capa_msg(MESSAGE_ERROR,tmpS); |
|
errCode = 1; |
|
break; |
|
case S_VAR: case S_CONSTANT: |
|
tmp_int = free_array(FIFTH_ARGSTR(argp)); |
|
r_p = gen_random_by_selector(FIFTH_ARGSTR(argp),sel,tmp_str,item_cnt,para1,para2); |
|
capa_mfree((char *)resultp); |
|
resultp = r_p; |
|
break; |
|
case IDENTIFIER: |
|
tmp_int = free_array(FIFTH_ARGNAME(argp)); |
|
r_p = gen_random_by_selector(FIFTH_ARGNAME(argp),sel,tmp_str,item_cnt,para1,para2); |
|
capa_mfree((char *)resultp); |
|
resultp = r_p; |
|
break; |
|
} |
|
} /* the fourth argument of this function (item_cnt) */ |
|
} /* the third argument of this function (seed) */ |
|
} /* the second argument of this function (paramenter one) */ |
|
} /* the first argument of this function (parameter two) */ |
|
|
|
} |
|
break; |
|
case RANDOM_POISSON_F: |
|
case RANDOM_EXPONENTIAL_F: |
|
case RANDOM_CHI_F: /* one parameter functions */ |
|
{ errCode = 0; |
|
switch( FIRST_ARGTYPE(argp) ) { /* parameter one */ |
|
case I_VAR: case I_CONSTANT: |
|
para1 = (float)FIRST_ARGINT(argp); |
|
break; |
|
case R_VAR: case R_CONSTANT: |
|
para1 = (float)FIRST_ARGREAL(argp); |
|
break; |
|
case S_VAR: case S_CONSTANT: |
|
case IDENTIFIER: |
|
resultp->s_type = S_CONSTANT; |
|
resultp->s_str = strsave("<<LAST ARG. OF THIS FUNCTION MUST BE A NUMBER>>"); |
|
sprintf(tmpS,"%s()'s last arg. must be a number.\n",FuncStack[Func_idx].s_name); |
|
capa_msg(MESSAGE_ERROR,tmpS); |
|
errCode = 1; |
|
break; |
|
} |
|
if(errCode == 0 ) { |
|
switch( SECOND_ARGTYPE(argp) ) { /* seed */ |
|
case I_VAR: case I_CONSTANT: |
|
tmp_str = (char *)capa_malloc(32,1); |
|
sprintf(tmp_str,"%ld",SECOND_ARGINT(argp) ); |
|
break; |
|
case R_VAR: case R_CONSTANT: |
|
tmp_long = (long)SECOND_ARGREAL(argp); |
|
tmp_str = (char *)capa_malloc(32,1); |
|
sprintf(tmp_str,"%ld",tmp_long); |
|
break; |
|
case S_VAR: case S_CONSTANT: |
|
tmp_str = strsave(SECOND_ARGSTR(argp)); |
|
break; |
|
case IDENTIFIER: |
|
resultp->s_type = S_CONSTANT; |
|
resultp->s_str = strsave("<<THIRD ARG. OF THIS FUNCTION MUST BE A NUMBER>>"); |
|
sprintf(tmpS,"%s()'s third arg. must be a number.\n",FuncStack[Func_idx].s_name); |
|
capa_msg(MESSAGE_ERROR,tmpS); |
|
errCode = 1; |
|
break; |
|
} |
|
if(errCode == 0 ) { |
|
switch( THIRD_ARGTYPE(argp) ) { /* item_cnt */ |
|
case I_VAR: case I_CONSTANT: |
|
item_cnt = THIRD_ARGINT(argp); |
|
break; |
|
case R_VAR: case R_CONSTANT: |
|
item_cnt = (int)THIRD_ARGREAL(argp); |
|
break; |
|
case S_VAR: case S_CONSTANT: |
|
case IDENTIFIER: |
|
resultp->s_type = S_CONSTANT; |
|
resultp->s_str = strsave("<<SECOND ARG. OF THIS FUNCTION MUST BE A NUMBER>>"); |
|
sprintf(tmpS,"%s()'s second arg. must be a number.\n",FuncStack[Func_idx].s_name); |
|
capa_msg(MESSAGE_ERROR,tmpS); |
|
errCode = 1; |
|
break; |
|
} |
|
if(errCode == 0 ) { |
|
switch( FOURTH_ARGTYPE(argp) ) { /* array_name, clear the content of this array first */ |
|
case I_VAR: case I_CONSTANT: |
|
case R_VAR: case R_CONSTANT: |
|
resultp->s_type = S_CONSTANT; |
|
resultp->s_str = strsave("<<FIRST ARG. OF THIS FUNCTION MUST BE AN ARRAY NAME>>"); |
|
sprintf(tmpS,"%s()'s first arg. must be a name of an array.\n",FuncStack[Func_idx].s_name); |
|
capa_msg(MESSAGE_ERROR,tmpS); |
|
errCode = 1; |
|
break; |
|
case S_VAR: case S_CONSTANT: |
|
tmp_int = free_array(FOURTH_ARGSTR(argp)); |
|
r_p = gen_random_by_selector(FOURTH_ARGSTR(argp),sel,tmp_str,item_cnt,para1,para2); |
|
capa_mfree((char *)resultp); |
|
resultp = r_p; |
|
break; |
|
case IDENTIFIER: |
|
tmp_int = free_array(FOURTH_ARGNAME(argp)); |
|
r_p = gen_random_by_selector(FOURTH_ARGNAME(argp),sel,tmp_str,item_cnt,para1,para2); |
|
capa_mfree((char *)resultp); |
|
resultp = r_p; |
|
break; |
|
} |
|
|
|
} /* the third argument of this function (seed) */ |
|
} /* the second argument of this function (paramenter one) */ |
|
} /* the first argument of this function (parameter two) */ |
|
} |
|
break; |
|
} /* end second switch */ |
|
} break; |
|
case ARRAY_MOMENTS_F: /* */ |
{ |
{ |
char *tmp_input; |
char *tmp_input; |
Symbol *r_p; |
Symbol *r_p; |
Line 819 ArgNode_t *argp;
|
Line 1073 ArgNode_t *argp;
|
resultp->s_str = strsave("NOT YET"); |
resultp->s_str = strsave("NOT YET"); |
} break; |
} break; |
|
|
|
case INIT_ARRAY_F: |
|
{ int rr; |
|
|
|
switch( FIRST_ARGTYPE(argp) ) { |
|
case I_VAR: case I_CONSTANT: |
|
case R_VAR: case R_CONSTANT: |
|
resultp->s_type = S_CONSTANT; |
|
resultp->s_str = strsave("<<ARG. OF THIS FUNCTION MUST BE AN ARRAY NAME>>"); |
|
sprintf(tmpS,"init_array()'s arg. must be a name of an array.\n"); |
|
capa_msg(MESSAGE_ERROR,tmpS); |
|
break; |
|
case S_VAR: case S_CONSTANT: /* allows the use of init_array(array[1]) which array[1]="another" */ |
|
rr = free_array(FIRST_ARGSTR(argp)); |
|
resultp->s_type = I_CONSTANT; |
|
resultp->s_int = rr; |
|
break; |
|
case IDENTIFIER: |
|
rr = free_array(FIRST_ARGNAME(argp)); |
|
resultp->s_type = I_CONSTANT; |
|
resultp->s_int = rr; |
|
break; |
|
} |
|
} break; |
case ARRAY_MAX_F: |
case ARRAY_MAX_F: |
case ARRAY_MIN_F: |
case ARRAY_MIN_F: |
{ int min; |
{ int min; |
Line 1014 ArgNode_t *argp;
|
Line 1291 ArgNode_t *argp;
|
} |
} |
break; |
break; |
case TEX_F: { if (Parsemode_f != TeX_MODE) { |
case TEX_F: { if (Parsemode_f != TeX_MODE) { |
|
#ifdef TTH |
|
#define CHARLEN 1024*1024 |
|
{ |
|
char *html; |
|
if ( (Parsemode_f==HTML_MODE) && |
|
((SECOND_ARGTYPE(argp) == S_VAR) || |
|
(SECOND_ARGTYPE(argp) == S_CONSTANT)) |
|
) { |
|
printf("Hi There %s\n",SECOND_ARGSTR(argp)); |
|
resultp->s_type = SECOND_ARGTYPE(argp); |
|
if(tth_err) { free(tth_err); tth_err=NULL; } |
|
textohtmldyn(SECOND_ARGSTR(argp),&html,&tth_err,CHARLEN); |
|
if(html) { |
|
resultp->s_str=strsave(html); |
|
capa_mfree(html); |
|
} else { |
|
resultp->s_str=strsave(""); |
|
} |
|
break; |
|
} |
|
} |
|
#undef CHARLEN |
|
#endif |
resultp->s_type = FIRST_ARGTYPE(argp); |
resultp->s_type = FIRST_ARGTYPE(argp); |
switch(FIRST_ARGTYPE(argp)) { |
switch(FIRST_ARGTYPE(argp)) { |
case I_VAR: |
case I_VAR: |