Changeset 806


Ignore:
Timestamp:
06/05/09 07:25:02 (15 years ago)
Author:
Mathieu Morlighem
Message:

removed (again) eval

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/classes/public/ismodelselfconsistent.m

    r804 r806  
    6767        'tolx','np','eps_rel','exclusive','n','gridonbed','gridonsurface','elementonbed','elementonsurface','deltaH','DeltaH','timeacc','timedec'};
    6868for i=1:length(fields),
    69         if ~isempty(eval(['md.' fields{i}])),
    70                 if any(isnan(eval(['md.' fields{i}]))),
     69        if ~isempty(md.(fields{i})),
     70                if any(isnan(md.(fields{i}))),
    7171                        disp(['model ' md.name ' has an NaN value in field ' fields{i} '!']);
    7272                        bool=0; return;
     
    8080        'sparsity','lowmem','n','gridonbed','gridonsurface','elementonbed','elementonsurface','deltaH','DeltaH','timeacc','timedec'};
    8181for i=1:length(fields),
    82         if ~isempty(eval(['md.' fields{i}])),
    83                 if any((eval(['md.' fields{i}]))<0),
     82        if ~isempty(md.(fields{i})),
     83                if any(md.(fields{i})<0),
    8484                        disp(['model ' md.name ' has a <0 value in field ' fields{i} '!']);
    8585                        bool=0; return;
     
    9797        'sparsity','deltaH','DeltaH','timeacc','timedec'};
    9898for i=1:length(fields),
    99         if ~isempty(eval(['md.' fields{i}])),
    100                 if any((eval(['md.' fields{i}]))==0),
     99        if ~isempty(md.(fields{i})),
     100                if any(md.(fields{i})==0),
    101101                        disp(['model ' md.name ' has a =0 value in field ' fields{i} '!']);
    102102                        bool=0; return;
     
    108108fields={'elements','p','q','elementoniceshelf','n','elementonbed'};
    109109for i=1:size(fields,2),
    110         if (size(eval(['md.' fields{i}]),1)~=md.numberofelements),
     110        if (size(md.(fields{i}),1)~=md.numberofelements),
    111111                disp(['model ' md.name ' field ' fields{i} ' should be of size ' num2str(md.numberofelements) '!']);
    112112                bool=0; return;
     
    117117fields={'x','y','z','B','drag','gridondirichlet_diag','melting','accumulation','surface','thickness','bed','gridonbed','gridonsurface'};
    118118for i=1:length(fields),
    119         if length(eval(['md.' fields{i}]))~=md.numberofgrids,
     119        if length(md.(fields{i}))~=md.numberofgrids,
    120120                disp(['model ' md.name ' field ' fields{i} ' should be of size ' num2str(md.numberofgrids) '!']);
    121121                bool=0; return;
     
    258258        fields={'dt','ndt'};
    259259        for i=1:length(fields),
    260                 if any((eval(['md.' fields{i}]))<0),
     260                if any(md.(fields{i})<0),
    261261                        disp(['model ' md.name ' has a <0 value in field ' fields{i} '!']);
    262262                        bool=0; return;
     
    329329        fields={'vx_obs','vy_obs'};
    330330        for i=1:length(fields),
    331                 if any(length(eval(['md.' fields{i}]))~=md.numberofgrids),
     331                if any(length(md.(fields{i}))~=md.numberofgrids),
    332332                        disp(['model ' md.name ' field ' fields{i} ' should be of size ' num2str(md.numberofgrids) '!']);
    333333                        bool=0; return;
     
    371371        fields={'sparsity'};
    372372        for i=1:length(fields),
    373                 if ~isempty(eval(['md.' fields{i}])),
    374                         if any(isnan(eval(['md.' fields{i}]))),
     373                if ~isempty(md.(fields{i})),
     374                        if any(isnan(md.(fields{i}))),
    375375                                disp(['model ' md.name ' has an NaN value in field ' fields{i} '!']);
    376376                                bool=0; return;
     
    382382        fields={'sparsity'};
    383383        for i=1:length(fields),
    384                 if ~isempty(eval(['md.' fields{i}])),
    385                         if any((eval(['md.' fields{i}]))<0),
     384                if ~isempty(md.(fields{i})),
     385                        if any(md.(fields{i})<0),
    386386                                disp(['model ' md.name ' has a <0 value in field ' fields{i} '!']);
    387387                                bool=0; return;
     
    393393        fields={'sparsity'};
    394394        for i=1:length(fields),
    395                 if ~isempty(eval(['md.' fields{i}])),
    396                         if any((eval(['md.' fields{i}]))==0),
     395                if ~isempty(md.(fields{i})),
     396                        if any(md.(fields{i})==0),
    397397                                disp(['model ' md.name ' has a =0 value in field ' fields{i} '!']);
    398398                                bool=0; return;
     
    444444                fields={'time','np'};
    445445                for i=1:length(fields),
    446                         if ~isempty(eval(['md.' fields{i}])),
    447                                 if any(isnan(eval(['md.' fields{i}]))),
     446                        if ~isempty(md.(fields{i})),
     447                                if any(isnan(md.(fields{i}))),
    448448                                        disp(['model ' md.name ' has an NaN value in field ' fields{i} '!']);
    449449                                        bool=0; return;
     
    455455                fields={'time','np'};
    456456                for i=1:length(fields),
    457                         if ~isempty(eval(['md.' fields{i}])),
    458                                 if any((eval(['md.' fields{i}]))<0),
     457                        if ~isempty(md.(fields{i})),
     458                                if any(md.(fields{i})<0),
    459459                                        disp(['model ' md.name ' has a <0 value in field ' fields{i} '!']);
    460460                                        bool=0; return;
     
    466466                fields={'time','np'};
    467467                for i=1:length(fields),
    468                         if ~isempty(eval(['md.' fields{i}])),
    469                                 if any((eval(['md.' fields{i}]))==0),
     468                        if ~isempty(md.(fields{i})),
     469                                if any(md.(fields{i})==0),
    470470                                        disp(['model ' md.name ' has a =0 value in field ' fields{i} '!']);
    471471                                        bool=0; return;
Note: See TracChangeset for help on using the changeset viewer.