0001 function status= call_netgen(geo_file, vol_file, msz_file, finelevel)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023 if nargin<3
0024 msz_file= '';
0025 end
0026
0027 if ~isempty(msz_file)
0028 eidors_msg('call_netgen: Warning. Using an *.msz file. This often fails.');
0029 end
0030
0031 if nargin<4
0032
0033
0034 finelevel= eidors_default('get','call_netgen_finelevel');
0035 end
0036
0037 islinux = isunix;
0038
0039
0040 cache_path = eidors_cache('cache_path');
0041 if islinux
0042 ng_name = 'netgen';
0043 else
0044 ng_name = [cache_path,'/ng'];
0045 end
0046
0047 while( 1 )
0048
0049 fid= fopen('ng.opt','a');
0050 if fid==-1
0051 error(['Netgen requires writing files in the current directory(%s). ', ...
0052 'Unfortunately, you don''t have permission. ' ...
0053 'Your options are: 1) change your working directory to one in which you have write permission, or ' ...
0054 '2) change the permissions on the current working directory.'], pwd);
0055 end
0056 if ~isempty(msz_file)
0057
0058
0059
0060 fprintf(fid,'options.meshsizefilename= %s\n',msz_file);
0061 end
0062 fclose(fid);
0063
0064 if ~islinux
0065
0066
0067 oldpath = getenv('NETGEN_USER_DIR');
0068 setenv('NETGEN_USER_DIR', cd);
0069 end
0070
0071 if eidors_debug('query','call_netgen')
0072 sys_cmd = sprintf('"%s" %s -geofile=%s -meshfile=%s ', ...
0073 ng_name, finelevel,geo_file,vol_file);
0074 else
0075 sys_cmd = sprintf('"%s" %s -batchmode -geofile="%s" -meshfile="%s" ', ...
0076 ng_name, finelevel,geo_file,vol_file);
0077 end
0078 status= system_cmd( sys_cmd );
0079
0080 if status==0; break; end
0081 try
0082 if islinux
0083 global eidors_objects;
0084 ignore_errors = 0;
0085 if isfield(eidors_objects,'ng_ignore_errors');
0086 ignore_errors = eidors_objects.ng_ignore_errors;
0087 end
0088 fprintf(['It seems that netgen has not worked (on linux). ' ...
0089 'Check that netgen is installed and on the path. ' ...
0090 'Check the environment variable NETGENDIR is set. ' ...
0091 ' This can be set via the following commands:\n' ...
0092 ' setenv(''NETGENDIR'',''/path/to/netgen/bin'');\n' ...
0093 ' setenv(''PATH'',[''/path/to/netgen/bin:'',getenv(''PATH'')]);\n' ...
0094 'Please enter a new netgen file name\n' ]);
0095 if ignore_errors == 1;
0096 eidors_msg('Ignoring Netgen Error (as requested)',1);
0097 break;
0098 end
0099 ng_name = input( ...
0100 'netgen file name (with path)? [or i=ignore, e=error, a=always ignore]','s');
0101 if strcmp(ng_name,'i'); break;end
0102 if strcmp(ng_name,'a');
0103 eidors_objects.ng_ignore_errors = 1;
0104 break;
0105 end
0106 if strcmp(ng_name,'e'); error('user requested'),end;
0107 else
0108
0109
0110 netgen_path = [eidors_cache('eidors_path'),'/../Netgen-5.3_x64/bin'];
0111 if exist(netgen_path)
0112 fprintf('Creating batch file to access Netgen.\n')
0113 fid= fopen([cache_path, '/ng.bat'],'w');
0114 if fid<0; error('Unable to write to %s',cache_path); end
0115 fprintf(fid,'set TCL_LIBRARY=%s/lib/tcl8.3\n', netgen_path);
0116 fprintf(fid,'set TIX_LIBRARY=%s/lib/tix8.1\n', netgen_path);
0117 fprintf(fid,'set NETGENDIR=%s\n', netgen_path);
0118 fprintf(fid,'"%s/netgen.exe" %%*\n', netgen_path);
0119 fclose(fid);
0120 else
0121 fprintf([ ...
0122 'Netgen call failed. Is netgen installed and on the search path?\n' ...
0123 'If you are running under windows, I can attempt to create\n' ...
0124 'a batch file to access netgen.\n' ...
0125 'Please enter the directory in which to find netgen.exe.\n' ...
0126 'A typical path is "C:\\Program Files (x86)\\Netgen-5.0_x64\\bin"\n' ...
0127 'If you don''t have a copy, download it from' ...
0128 'http://sourceforge.net/projects/netgen-mesher/ \n\n']);
0129 netgen_path = input('netgen_path? [or i=ignore, e=error] ','s');
0130 if strcmp(ng_name,'i'); break;end
0131 if strcmp(ng_name,'e'); error('user requested'),end;
0132 if exist( sprintf('%s/netgen.exe',netgen_path) , 'file' ) || ...
0133 exist( sprintf('%s/bin/netgen.exe',netgen_path) , 'file' )
0134 disp('Found netgen version 4.4 or higher');
0135 netgen_exe = netgen_path;
0136 if exist( sprintf('%s/bin/netgen.exe',netgen_path) , 'file' )
0137 netgen_exe = [netgen_path '/bin'];
0138 end
0139
0140 fid= fopen([cache_path, '/ng.bat'],'w');
0141 if fid<0; error('Unable to write to %s',cache_path); end
0142 fprintf(fid,'set TCL_LIBRARY=%s/lib/tcl8.3\n', netgen_path);
0143 fprintf(fid,'set TIX_LIBRARY=%s/lib/tix8.1\n', netgen_path);
0144 fprintf(fid,'set NETGENDIR=%s\n', netgen_path);
0145 fprintf(fid,'"%s/netgen.exe" %%*\n', netgen_exe);
0146 fclose(fid);
0147 elseif exist( sprintf('%s/ng431.exe',netgen_path) , 'file' )
0148 disp('Found netgen version 4.3.1');
0149 fid= fopen([cache_path, '/ng.bat'],'w');
0150 if fid<0; error('Unable to write to %s',cache_path); end
0151 fprintf(fid,'set TCL_LIBRARY=%s/lib/tcl8.3\n', netgen_path);
0152 fprintf(fid,'set TIX_LIBRARY=%s/lib/tcl8.2\n', netgen_path);
0153 fprintf(fid,'"%s/ng431.exe" %%*\n', netgen_path);
0154 fclose(fid);
0155 else
0156 warning(['cannot find a version of netgen that I know about\n' ...
0157 'Install netgen or check the path\n']);
0158 end
0159 end
0160 end
0161 catch e
0162 if strncmp(computer,'PC',2)
0163
0164 setenv('NETGEN_USER_DIR', oldpath);
0165 end
0166 rethrow(e)
0167 end
0168 end