#include <qop_internal.h>
#include <string.h>

#define OPTNUM (QOP_wilson_style+4*QOP_wilson_nsvec)
int QOP_wilson_inited = 0;
int QOP_wilson_style = 0;
int QOP_wilson_nsvec = 4;
int QOP_wilson_nvec = 4;
int QOP_wilson_cgtype = 1;
int QOP_wilson_optnum = 16;
int QOP_wilson_eigcg_numax = 20;
int QOP_wilson_eigcg_m = 10;
int QOP_wilson_eigcg_nev = 3;

QOP_status_t
QOP_wilson_invert_set_opts(QOP_opt_t opts[], int nopts)
{
  int i;

  WILSON_INVERT_BEGIN;

  for(i=0; i<nopts; i++) {
    char *tag = opts[i].tag;
    double value = opts[i].value;

    if(!strcmp(tag,"ns")) {
      if( (value==1)||(value==2)||(value==4)||
	  (((QOP_wilson_style&1)==1)&&(value==8)) ) {
	QOP_wilson_nsvec = (int) value;
      } else {
	return QOP_FAIL;
      }
    } else if(!strcmp(tag,"nm")) {
      if( (value==1)||(value==2)||(value==4)||
	  (((QOP_wilson_style&1)==1)&&(value==8)) ) {
	QOP_wilson_nvec = (int) value;
      } else {
	return QOP_FAIL;
      }
    } else if(!strcmp(tag,"st")) {
      if((value==0)||(value==1)||(value==2)||(value==3)) {
	QOP_wilson_style = (int) value;
	if((QOP_wilson_style&1)==0) {
	  if(QOP_wilson_nsvec>4) QOP_wilson_nsvec = 4;
	  if(QOP_wilson_nvec>4) QOP_wilson_nvec = 4;
	}
      } else {
	return QOP_FAIL;
      }
    } else if(!strcmp(tag,"cg")) {
      if((value==0)||(value==1)||(value==2)) {
	QOP_wilson_cgtype = (int) value;
      } else {
	return QOP_FAIL;
      }
    } else if(!strcmp(tag,"eigcg_l")) {
      QOP_wilson_eigcg_numax = (int) value;
    } else if(!strcmp(tag,"eigcg_m")) {
      QOP_wilson_eigcg_m = (int) value;
    } else if(!strcmp(tag,"eigcg_nev")) {
      QOP_wilson_eigcg_nev = (int) value;
    }
  }
  QOP_wilson_optnum = OPTNUM;

  WILSON_INVERT_END;
  return QOP_SUCCESS;
}