PRO dustem_run_plugins, p_dim ,$ param_descs,$ param_values,$ param_func,$ scopes,$ st=st,$ avoid=avoid,$ dustem_run=dustem_run,$ help=help ;+ ; NAME: ; dustem_run_plugins ; ; PURPOSE: ; runs the plugins according to their scopes ; a N_ELEMENTS(PLUGINS)X N_ELEMENTS(SCOPES) matrix is created to achieve this ; ; CATEGORY: ; DustEMWrap, Distributed, Mid-Level, Plugins ; ; CALLING SEQUENCE: ; dustem_run_plugins,p_dim,param_descs,param_values,scopes,[,st][,avoid][,dustem_run][,/help] ; ; INPUTS: ; p_dim = parameter values ; param_descs = parameter description vector ; param_values = current parameter values ; param-func = plugin indices array ; scopes = scopes of plugins ; ; OPTIONAL INPUT PARAMETERS: ; avoid = scopes should be avoided ; dustem_run = Dustem should be run ; ; OUTPUTS: ; None ; ; OPTIONAL OUTPUT PARAMETERS: ; st = Dustem output structure ; ; ACCEPTED KEY-WORDS: ; help = if set, print this help ; ; COMMON BLOCKS: ; None ; ; SIDE EFFECTS: ; None ; ; RESTRICTIONS: ; The DustEM fortran code must be installed ; The DustEMWrap IDL code must be installed ; ; PROCEDURES AND SUBROUTINES USED: ; ; EXAMPLES ; dustem_activate_plugins, p_min, pram_descs, param_values, param_func, scopes, st=st, avoid=avoid, dustem_run=dustem_run ; ; MODIFICATION HISTORY: ; Written by IC ; Evolution details on the DustEMWrap gitlab. ; See http://dustemwrap.irap.omp.eu/ for FAQ and help. ;- ;Because I don't want to run this block each time the procedure is run I will test on the presence of all the scope tags in !dustem_plugin ;if all of them are present then this block has been run. ;stop dew_prefix='dustem_plugin_' ; we assume that this prefix is at the start of all plugin routines IF (tag_names(*!dustem_plugin))(0) EQ 'NONE' THEN goto, the_end ;Just to be sure if not keyword_set(avoid) then avoid=0 if not keyword_set(dustem_run) then dustem_run=0 ;stop test_all_scopes = 1. FOR i=0L,n_tags(*!dustem_plugin)-1 DO BEGIN IF isa(((*!dustem_plugin).(i).scope)) then test=1 else test=0 test_all_scopes*=test ENDFOR ;stop IF ~test_all_scopes THEN BEGIN ;FIRST LOOP TO RETRIEVE THE SCOPES OF ALL THE PLUGINS f=1 ; Initializing the index that is associated to each plugin FOR i=0L,n_elements(param_descs)-1 DO BEGIN parameter_type=dustem_parameter_description2type(param_descs[i],string_name=string_name) ; Looping over the parameter description vector (here saved in a system variable) to sort out its different elements with repect to their types IF parameter_type EQ 'PLUGIN' THEN BEGIN ; Selecting the plugins ftn = strmid(param_descs[i],0) ; String containing the name of the plugin and the keyword used (ie: dustem_create_continuum_2) ii = strsplit(ftn,'_',count=countx) & ii = ii[countx-1]-1 ; Locating the last underscore to automate the extraction of the plugin's keyword ftn = strmid(ftn,0,ii) ; String containing the name of the plugin without the associated keyword this_plugin = strupcase(strmid(ftn,strlen(dew_prefix))) ; end_char=strlen(this_plugin) ; k=where(tag_names(*!dustem_plugin) eq strmid(this_plugin,0,end_char),counte) ; Selecting a plugin through matching the string name of the plugin from the scope system variable with the one read from the parameter description vector k=where(tag_names(*!dustem_plugin) eq this_plugin,counte) ; Selecting a plugin through matching the string name of the plugin from the scope system variable with the one read from the parameter description vector tmp = where(param_func eq f, count) ; Array of indices helping with the sorting out of the different plugins (as a plugin can be called several times with different keywords so as to fit several quantities) PDO_tmp = param_descs(tmp) ; Array containing the different calls of each plugin at a time p_dim_tmp = param_values(tmp) ; Array containing the values of the fitted parameters vari = strmid(PDO_tmp(0),ii+1) ; Test variable to be used to alter the initialization of the index and value arrays below as a keyword can be a scalar or a string ; Initialization of the index and value arrays. These arrays will be used (below) when calling the plugins (with the default key=key and val=val keywords) ;JPB: I have seen cases where count=0 there, which makes the code fail. So including that test below IF count EQ 0 THEN BEGIN message, 'something went wrong, signé JPB',/continue stop ENDIF if strlen(strtrim(string(vari),2)) gt 1 then begin index = strarr(count) & value = fltarr(count) endif else begin index = fltarr(count) & value = fltarr(count) endelse ; Filling the index and value arrays for one plugin at a time FOR l=0, count-1 DO BEGIN index[l] = strmid(PDO_tmp[l],ii+1) value[l] = p_dim_tmp[l] ENDFOR ;==============Dry run of the plugins to obtain their scopes and run them accordingly============== str='toto='+ftn+'(scope=scope)' & str=str(0) toto=execute(str) ((*!dustem_plugin).(k).scope)=ptr_new(scope) str='toto='+ftn+'(paramtag=paramtag)' & str=str(0) toto=execute(str) ((*!dustem_plugin).(k).paramtag)=ptr_new(paramtag) f=f+1 & i=i+count-1 ; Incrementing the parameter and same-type plugin indices ENDIF ENDFOR ENDIF ;stop ;LOCATING THE PLUGINS THAT SHOULD BE RUN indices = fltarr(n_tags(*!dustem_plugin),n_elements(scopes)) ; FOR i=0L,n_elements(scopes)-1 DO BEGIN ;test if there is a '*' string. This means that all scopes containing that string will be either avoided or chosen. ;this means that we need the index of the scope corresponding to this string. ;first we need to extract the string between the two '*' and then match it with the existing plugin ; this enables us to get the indices of the plugins to run... IF strmid(scopes(i),0,1) EQ '*' THEN BEGIN ;Now I need to locate the second '*': ind_ast = STRPOS(scopes(i), '*', /REVERSE_SEARCH) scope_to_match = strmid(scopes(i),1,ind_ast-1) ; ENDIF ELSE BEGIN ;assuming the string is identical to the plugin's scope scope_to_match = strupcase(scopes(i)) ENDELSE ;Now locating the indices of this scope in !dustem_plugin. FOR j=0L,n_tags(*!dustem_plugin)-1 DO BEGIN test = strcmp((*((*!dustem_plugin).(j).scope)),scope_to_match,strlen(scope_to_match),/FOLD_CASE) EQ 1 IF test and keyword_set(avoid) THEN BEGIN if (size(indices))[0] GT 1 then indices[j,i] *= 0 ELSE indices[j] *= 0 ENDIF ;ELSE BEGIN ; if (size(indices))[0] GT 1 then indices[j,i] = 1 ELSE indices[j] = 1 ; ENDELSE ; IF test and not keyword_set(avoid)then BEGIN if (size(indices))[0] GT 1 then indices[j,i] += 1 ELSE indices[j] += 1 ENDIF if ~test and keyword_set(avoid) then begin if (size(indices))[0] GT 1 then indices[j,i] += 1 ELSE indices[j] += 1 endif ENDFOR ENDFOR ;stop ;RUNNING THE PLUGINS ACCORDING TO THEIR SCOPES f=1 ; Initializing the index that is associated to each plugin FOR i=0L,n_elements(param_descs)-1 DO BEGIN parameter_type=dustem_parameter_description2type(param_descs[i],string_name=string_name) ; Looping over the parameter description vector (here saved in a system variable) to sort out its different elements with repect to their types IF parameter_type EQ 'PLUGIN' THEN BEGIN ; Selecting the plugins ;stop ftn = strmid(param_descs(i),0) ; String containing the name of the plugin and the keyword used (ie: dustem_create_continuum_2) ii = strsplit(ftn,'_',count=countx) & ii = ii(countx-1)-1 ; Locating the last underscore to automate the extraction of the plugin's keyword ftn = strmid(ftn,0,ii) ; String containing the name of the plugin without the associated keyword this_plugin = strupcase(strmid(ftn,strlen(dew_prefix))) ; end_char=strlen(this_plugin) ; k=where(tag_names(*!dustem_plugin) eq strmid(this_plugin,0,end_char),counte) ; Selecting a plugin through matching the string name of the plugin from the scope system variable with the one read from the parameter description vector k=where(tag_names(*!dustem_plugin) eq this_plugin,counte) ; Selecting a plugin through matching the string name of the plugin from the scope system variable with the o ; k=where(strmid(tag_names(*!dustem_plugin),0,8) eq strmid(strupcase(strmid(ftn,14)),0,8),counte) ; Selecting a plugin through matching the string name of the plugin from the scope system variable with the one read from the parameter description vector ;help,f ;stop tmp = where(param_func eq f, count) ; Array of indices helping with the sorting out of the different plugins (as a plugin can be called several times with different keywords so as to fit several quantities) PDO_tmp = param_descs(tmp) ; Array containing the different calls of each plugin at a time p_dim_tmp = param_values(tmp) ; Array containing the values of the parameters used for the run vari = strmid(PDO_tmp(0),ii+1) ; Test variable to be used to alter the initialization of the index and value arrays below for a keyword can be a scalar or a string ; Initialization of the index and value arrays. These arrays will be used (below) when calling the plugins (with the default key=key and val=val keywords) ;print,ftn,ii,ftn,counte ;print,PDO_tmp,p_dim_tmp,vari if strlen(strtrim(string(vari),2)) gt 1 then begin index = strarr(count) & value = fltarr(count) endif else begin index = fltarr(count) & value = fltarr(count) endelse ;help,'<<<<',index,count ; Filling the index and value arrays for one plugin at a time FOR l=0, count-1 DO BEGIN index[l] = strmid(PDO_tmp[l],ii+1) value[l] = p_dim_tmp[l] ENDFOR ;==============Running of the plugins according to their scopes============== ;stop if keyword_set(avoid) then idd=product(indices[k,*]) else idd=total(indices[k,*]) IF idd GE 1 THEN BEGIN str='((*!dustem_plugin).('+strtrim(k,2)+').spec)=ptr_new('+ftn+'(key=index,val=value)'+')' & str=str(0) ;print,'>>>>>',str ;stop ;plante fawlty quand ;((*!dustem_plugin).(1).spec)=ptr_new(dustem_plugin_stellar_population(key=index,val=value)) toto=execute(str) ;stop IF !dustem_verbose NE 0 THEN message,strupcase(strmid(ftn,7)),/info ENDIF f=f+1 & i=i+count-1 ; Incrementing the parameter and same-type plugin indices ENDIF ENDFOR the_end: if keyword_set(dustem_run) then begin ;and idd gt 0. : remved this so that dustem is always run when the keyword is mentioned. st = dustem_run(p_dim) !dustem_current = ptr_new(st) endif END