Doxygen Source Code Documentation
Main Page Alphabetical List Data Structures File List Data Fields Globals Search
mri_symbolize.c
Go to the documentation of this file.00001 #include "mrilib.h"
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015 floatvecvec * SYM_expand_ranges( int nlast, int nrang, SYM_irange *rang, char *str )
00016 {
00017 floatvec *fv ;
00018 floatvecvec *fvv=NULL ;
00019 int rr , ii , ss , gg, *qlist , nvec=0 , iv ;
00020 NI_str_array *sar ;
00021 char qname[64] , *qstr , *qpt , *qls ;
00022 float fac ;
00023
00024 ENTRY("SYM_expand_ranges") ;
00025
00026 if( nlast < 0 ) RETURN(NULL) ;
00027
00028
00029
00030 if( nrang < 1 || rang == NULL || str == NULL || *str == '\0' ) RETURN(NULL) ;
00031
00032
00033
00034 for( ii=0 ; str[ii] != '\0' && isspace(str[ii]) ; ii++ ) ;
00035
00036 if( str[ii] == '\0' ||
00037 str[ii] == '#' ||
00038 (str[ii] == '/' && str[ii+1] == '/')
00039 ) RETURN(NULL) ;
00040
00041 fv = (floatvec *)malloc(sizeof(floatvec)) ;
00042 fv->nar = nlast+1 ;
00043 fv->ar = (float *)calloc(sizeof(float),nlast+1) ;
00044
00045
00046
00047 sar = NI_decode_string_list( str , "~" ) ;
00048 if( sar == NULL ){
00049 fvv = (floatvecvec *)malloc(sizeof(floatvecvec)) ;
00050 fvv->nvec = 1 ;
00051 fvv->fvar = fv ;
00052 RETURN(fvv) ;
00053 }
00054
00055
00056
00057 for( ss=0 ; ss < sar->num ; ss++ ){
00058 qstr = sar->str[ss] ;
00059 if( qstr == NULL || *qstr == '\0' ) continue ;
00060 if( *qstr == '#' ||
00061 (*qstr == '/' && *(qstr+1) == '/') ) break ;
00062
00063 qstr = strdup(sar->str[ss]) ;
00064 qls = strchr(qstr,'[') ;
00065 qlist = NULL ;
00066 if( qls != NULL ){
00067 *qls = '\0' ;
00068 qls++ ;
00069 }
00070
00071 qpt = strchr(qstr,'*') ;
00072 if( qpt != NULL ){
00073 fac = (float)strtod(qstr,NULL) ;
00074 if( fac == 0.0 && *qstr != '0' ) fac = 1.0 ;
00075 qpt++ ;
00076 } else if( *qstr == '+' ){
00077 qpt = qstr+1 ; fac = 1.0 ;
00078 } else if( *qstr == '-' ){
00079 qpt = qstr+1 ; fac = -1.0 ;
00080 } else {
00081 qpt = qstr ; fac = 1.0 ;
00082 }
00083
00084 for( rr=0 ; rr < nrang ; rr++ )
00085 if( strcmp(qpt,rang[rr].name) == 0 ) break ;
00086 if( rr == nrang ){
00087 fprintf(stderr,"** ERROR: can't match symbolic name '%s'\n",qpt) ;
00088 free((void *)qstr) ; continue ;
00089 }
00090
00091 if( qls != NULL ){
00092 MCW_intlist_allow_negative( (rang[rr].nbot < 0) ) ;
00093 qlist = MCW_get_intlist( rang[rr].ntop+1 , qls ) ;
00094
00095 if( qlist != NULL && *qls == '[' ){
00096 if( nvec == 0 ){
00097 nvec = qlist[0] ;
00098 fvv = (floatvecvec *)malloc(sizeof(floatvecvec)) ;
00099 fvv->nvec = nvec ;
00100 fvv->fvar = (floatvec *)calloc(sizeof(floatvec),nvec) ;
00101 for( iv=0 ; iv < nvec ; iv++ ){
00102 fvv->fvar[iv].nar = nlast+1 ;
00103 fvv->fvar[iv].ar = (float *)calloc(sizeof(float),nlast+1) ;
00104 }
00105 } else if( qlist[0] != nvec ){
00106 fprintf(stderr,"** ERROR: mismatch in use of -gltsym [[...]]: '%s'\n",
00107 sar->str[ss] ) ;
00108 free((void *)qlist) ; free((void *)qstr) ;
00109 continue ;
00110 }
00111 for( iv=0 ; iv < nvec ; iv++ ){
00112 gg = qlist[iv+1] - rang[rr].nbot + rang[rr].gbot ;
00113 if( gg >= 0 && gg <= nlast ) fvv->fvar[iv].ar[gg] = fac ;
00114 }
00115 free((void *)qlist) ; free((void *)qstr) ;
00116 continue ;
00117 }
00118 }
00119
00120 if( qlist == NULL ){
00121 qlist = (int *)malloc(sizeof(int)*(rang[rr].ntop-rang[rr].nbot+2)) ;
00122 qlist[0] = rang[rr].ntop-rang[rr].nbot+1 ;
00123 for( ii=0 ; ii < qlist[0] ; ii++ ) qlist[ii+1] = rang[rr].nbot+ii ;
00124 }
00125
00126
00127 for( ii=0 ; ii < qlist[0] ; ii++ ){
00128 if( qlist[ii+1] < rang[rr].nbot || qlist[ii+1] > rang[rr].ntop ){
00129 fprintf(stderr,"** ERROR: subscript %s[%d] out of range %d..%d\n",
00130 rang[rr].name , qlist[ii+1] , rang[rr].nbot,rang[rr].ntop ) ;
00131 continue ;
00132 }
00133 gg = qlist[ii+1] - rang[rr].nbot + rang[rr].gbot ;
00134 if( gg >= 0 && gg <= nlast ) fv->ar[gg] = fac ;
00135 }
00136
00137 free((void *)qlist) ; free((void *)qstr) ;
00138 }
00139 MCW_intlist_allow_negative(0) ;
00140
00141 NI_delete_str_array(sar);
00142
00143
00144
00145 if( nvec == 0 ){
00146 fvv = (floatvecvec *)malloc(sizeof(floatvecvec)) ;
00147 fvv->nvec = 1 ;
00148 fvv->fvar = fv ;
00149 } else {
00150 for( iv=0 ; iv < nvec ; iv++ ){
00151 for( gg=0 ; gg <= nlast ; gg++ ){
00152 if( fvv->fvar[iv].ar[gg] == 0.0f ) fvv->fvar[iv].ar[gg] = fv->ar[gg] ;
00153 }
00154 }
00155 KILL_floatvec(fv) ;
00156 }
00157
00158 RETURN(fvv) ;
00159 }