function ofc_vfemp1_install(varargin)
  if strversion2num(version)<strversion2num('4.2.0')
    error('Using Octave version %s. Version must be >= 4.2.0. ',version)
  end
  warning('off','Octave:shadowed-function');more off
  v='0.1.0';
  v_fc_simesh='0.2.1';
  p = inputParser; 
  p.KeepUnmatched=true; 
  p.addParamValue('version', v, @ischar ); 
  p.addParamValue('dir', './fc-vfemp1-full', @ischar );
  p.addParamValue('fc_simesh',v_fc_simesh , @ischar );
  p.addParamValue('verbose', 1, @(x) ismember(x,0:2) ); % level of verbosity
  p.parse(varargin{:});
  R=p.Results;
  version=R.version;verbose=R.verbose;
  varargin=deleteCellOptions(varargin,p.Parameters);
  vprintf(verbose,1,'Parts of the GNU Octave <fc-vfemp1> package.\nCopyright (C) 2017 Francois Cuvelier <cuvelier@math.univ-paris13.fr>\n\n');

  vprintf(verbose,1,'**************************************\n');
  vprintf(verbose,1,'Downloading and installing the package\n');
  vprintf(verbose,1,'    <fc-simesh>[%s]\n',p.Results.fc_simesh)
  vprintf(verbose,1,'**************************************\n');
  fc_simesh_dir=get_install('simesh',p.Results.fc_simesh,p.Results.dir,'verbose',verbose,varargin{:});
  vprintf(verbose,1,'    <fc-simesh>[%s]: installed\n\n',p.Results.fc_simesh)
  
  vprintf(verbose,1,'**************************************\n');
  vprintf(verbose,1,'Downloading and installing the package\n');
  vprintf(verbose,1,'    <fc-vfemp1>[%s]\n',p.Results.version)
  vprintf(verbose,1,'**************************************\n');
  fc_vfemp1_dir=get_inst_fc_gen('vfemp1',p.Results.version,p.Results.dir);
  
  vprintf(verbose,1,'*** Setting the <fc-vfemp1> package\n');
  currentdir=pwd;
  cd(fc_vfemp1_dir)
  fc_vfemp1.configure('simesh_dir',fc_simesh_dir);
  cd(currentdir)
  vprintf(verbose,1,'*** The <fc-vfemp1>[%s] is installed\n',p.Results.version)
  vprintf(verbose,1,'*** Using instructions \n');
  vprintf(verbose,1,'   To use the <fc-vfemp1> package:\n');
  vprintf(verbose,1,'   addpath(''%s'')\n',fc_vfemp1_dir);
  vprintf(verbose,1,'   fc_vfemp1.init()\n');
  logfile=[currentdir,filesep,'ofc_vfemp1_set.m'];
  vprintf(verbose,1,'\n   See %s\n',logfile);
  fid=fopen(logfile,'w');
  fprintf(fid,'warning(''off'',''Octave:shadowed-function'');more off\n');
  fprintf(fid,'addpath(''%s'')\n',fc_vfemp1_dir);
  fprintf(fid,'fc_vfemp1.init()\n');
  fclose(fid);
end

function fulldir=get_install(name,ver,dir,varargin)
  mfile=sprintf('ofc_%s_install',name);
  urlfile=sprintf('http://www.math.univ-paris13.fr/~cuvelier/software/codes/Octave/fc-%s/%s/%s.m',name,ver,mfile);
  try
    urlwrite(urlfile,[mfile,'.m']);
  catch
    error('Unable to get :\n  -> %s\n',urlfile);
  end
  eval([mfile,'(''dir'',dir,varargin{:})']) % run installation 
  fulldir=fullfile(pwd,sprintf('%s%sfc_%s-%s',dir,filesep,name,ver));
end

function fulldir=get_inst_fc_gen(name,ver,dir)
  file=sprintf('ofc-%s-%s.tar.gz',name,ver);
  urlfile=sprintf('http://www.math.univ-paris13.fr/~cuvelier/software/codes/Octave/fc-%s/%s/fc_%s-%s.tar.gz',name,ver,name,ver);
  urlwrite(urlfile,file);
  try
    untar(file,dir)
  catch
    error('Unable to get :\n  -> %s\n',urlfile);
  end
  fulldir=fullfile(pwd,sprintf('%s%sfc_%s-%s',dir,filesep,name,ver));
end

function ver=strversion2num(strv)
% convert version string 'x.y.z' in x*10^6+y*10^3+z
S=strsplit(strv,'.');
assert(length(S)==3);
z=str2num(S{3});assert(z<10^3);
ver=z;
y=str2num(S{2});assert(y<10^3);
ver=ver+10^3*y;
x=str2num(S{1});assert(x<10^3);
ver=ver+10^6*x;
end

function NewOptions=deleteCellOptions(Options,delOptions)
  assert(iscell(Options))
  if ischar(delOptions),delOptions={ delOptions};end
  n=length(Options);
  isupp=[];
  delOptions=upper(delOptions);
  for i=1:2:n
    I=find(strcmp(upper(Options{i}),delOptions));
    if ~isempty(I)
      isupp=[isupp,i,i+1];
    end
  end
  NewOptions={Options{setdiff(1:n,isupp)}};
end

function vprintf(verbose,level,varargin)
  if verbose>=level, fprintf(varargin{:});end
end
