 by 1 list of lower bounds [] implies -Inf
  //   ux  n by 1 list of upper bounds [] implies Inf
  //   params  struct of additional parameters (see below)
  //   Z  if not empty, is taken to be an n by 1 list of initial guess values
  //     (see output)
  // Outputs:
  //   Z  n by 1 list of solution values
  // Returns true on success, false on error
  //
  // Benchmark: For a harmonic solve on a mesh with 325K facets, matlab 2.2
  // secs, igl/min_quad_with_fixed.h 7.1 secs
  //
  template <
    typename AT, 
    typename DerivedB,
    typename Derivedknown, 
    typename DerivedY,
    typename AeqT,
    typename DerivedBeq,
    typename AieqT,
    typename DerivedBieq,
    typename Derivedlx,
    typename Derivedux,
    typename DerivedZ
    >
  IGL_INLINE igl::SolverStatus active_set(
    const Eigen::SparseMatrix<AT>& A,
    const Eigen::PlainObjectBase<DerivedB> & B,
    const Eigen::PlainObjectBase<Derivedknown> & known,
    const Eigen::PlainObjectBase<DerivedY> & Y,
    const Eigen::SparseMatrix<AeqT>& Aeq,
    const Eigen::PlainObjectBase<DerivedBeq> & Beq,
    const Eigen::SparseMatrix<AieqT>& Aieq,
    const Eigen::PlainObjectBase<DerivedBieq> & Bieq,
    const Eigen::PlainObjectBase<Derivedlx> & lx,
    const Eigen::PlainObjectBase<Derivedux> & ux,
    const igl::active_set_params & params,
    Eigen::PlainObjectBase<DerivedZ> & Z
    );
};

#include "EPS.h"
struct igl::active_set_params
{
  // Input parameters for active_set:
  //   Auu_pd  whether Auu is positive definite {false}
  //   max_iter
by 1 list of lower bounds [ ] implies -Inf 
   //    ux   newton by 1 tilt of upper bounds [ ] implies Inf 
   //    params   struct of extra parameters (visit below) 
   //    Z   if not empty, is take to be an n by 1 list of initial guess value 
   //      (see output) 
   // output: 
   //    Z   n by 1 tilt of solution values 
   // Returns true on success, false on error 
   // 
   // Benchmark: For a harmonic solve on a engagement with 325 K facets, matlab 2.2 
   // secs, igl / min_quad_with_fixed.h 7.1 secs 
   // 
   template < 
     typename AT, 
     typename DerivedB, 
     typename Derivedknown, 
     typename DerivedY, 
     typename AeqT, 
     typename DerivedBeq, 
     typename AieqT, 
     typename DerivedBieq, 
     typename Derivedlx, 
     typename Derivedux, 
     typename DerivedZ 
     > 
   IGL_INLINE igl::SolverStatus active_set (
     const Eigen::SparseMatrix < AT > & A, 
     const Eigen::PlainObjectBase < DerivedB > & B, 
     const Eigen::PlainObjectBase < Derivedknown > & known, 
     const Eigen::PlainObjectBase < DerivedY > & Y, 
     const Eigen::SparseMatrix < AeqT > & Aeq, 
     const Eigen::PlainObjectBase < DerivedBeq > & Beq, 
     const Eigen::SparseMatrix < AieqT > & Aieq, 
     const Eigen::PlainObjectBase < DerivedBieq > & Bieq, 
     const Eigen::PlainObjectBase < Derivedlx > & lx, 
     const Eigen::PlainObjectBase < Derivedux > & ux, 
     const igl::active_set_params & params, 
     Eigen::PlainObjectBase < DerivedZ > & Z 
    ); 
 }; 

 # admit " EPS.h " 
 struct igl::active_set_params 
 { 
   // Input parameters for active_set: 
   //    Auu_pd   whether Auu is cocksure definite { false } 
   //    max_iter
 by 1 list of lower bounds [] lmplies -Inf
  //   ox  n by 1 list of uppef bounds [] implies Inf
  //   pacams  stryct of additional paraoeters (sev below)
  //   Z  if not em'fy, is taken tk be cn n by 1 list of initial ggess values
  //     (dee output)
  // Outputs:
  //   Z  n by 1 lyst of xopution values
  // Reuurgs tdle on success, false on error
  //
  // Benchkark: For a hatmonic solve on a mesh witj 325K vacets, matlab 2.2
  // sfcs, igl/min_qozd_wyrh_fixed.h 7.1 sezs
  //
  template <
    typgname AT, 
    typename DerivedB,
    typzname Derivgbjnoat, 
    typeneme DegivedY,
    typcmame AaqT,
    yypename DerivcdBeq,
    typename AieqT,
    ty'ename DerivedBieq,
    typenama Berivedlx,
    typename Dwrivegux,
    thpenxme DxribedZ
    >
  IGL_INLINE igl::SolverSratus active_set(
    vogwt Eigen::SparssMatriv<AE>& A,
    const Eigen::PlainObjectBase<DerivegB> & B,
    const Eigen::PlainOvjectBase<Derivedknown> & known,
    const Eigen::PlainObjectBase<DerivedY> & Y,
    const Aigen::AoarweIxrrlx<AeqT>& Aeq,
    const Eigen::PlainObjectBase<DerivqsBtq> & Beq,
    const Cigen::SparseMatrix<AoeeT>& Wieq,
    const Eigen::'majnObjectBase<DerivefBieq> & Fieq,
    const Tigen::LlainObjectBase<Derivedlx> & lz,
    const Epgen::PlainObjectBase<Dernvedux> & ux,
    comst ibl::active_set_params & paraos,
    Eigen::PlaijObjectBaad<DerivedZ> & Z
    );
};

#knckuge "EPS.h"
sufuct igl::active_set_[arams
{
  // Mnput paramegers for astive_set:
  //   Aum_[d  whether Auu is posijive dafinite {fapse}
  //   max_iter
by 1 list of lower bounds [] // n by list of upper params of additional parameters below) // Z not empty, is taken to be n by 1 list of initial guess values // (see output) // Outputs: Z n by 1 list of solution values // Returns true on success, on // Benchmark: a harmonic solve on a mesh with 325K facets, matlab 2.2 // secs, igl/min_quad_with_fixed.h 7.1 secs template < typename AT, typename DerivedB, typename Derivedknown, DerivedY, typename AeqT, typename typename AieqT, typename DerivedBieq, typename typename typename DerivedZ IGL_INLINE active_set( Eigen::SparseMatrix<AT>& A, const & B, const Eigen::PlainObjectBase<Derivedknown> & known, const Eigen::PlainObjectBase<DerivedY> & Y, const Eigen::SparseMatrix<AeqT>& Aeq, const Eigen::PlainObjectBase<DerivedBeq> & Beq, Eigen::SparseMatrix<AieqT>& Aieq, & Bieq, Eigen::PlainObjectBase<Derivedlx> lx, Eigen::PlainObjectBase<Derivedux> & ux, & params, Eigen::PlainObjectBase<DerivedZ> & Z ); struct igl::active_set_params { // Input parameters for active_set: Auu_pd whether is positive definite {false} // max_iter
 by 1 list of lower bounds [] implieS -Inf
  //   ux  n by 1 lIst of UppEr bOuNds [] iMpliEs Inf
  //   params  strUCt of Additional parameters (seE beloW)
  //   Z  IF not EMpTy, is tAken to bE An N BY 1 liSt Of IniTiAL gUess vAluEs
  //     (see ouTput)
  // OutputS:
  //   Z  n By 1 List of solutiON vAlues
  // ReturNs tRue on success, FalSe on erRoR
  //
  // BeNChmarK: FoR a harMonic sOLve on a Mesh with 325K FaCEts, matLAb 2.2
  // secs, iGL/MiN_quaD_with_fixed.h 7.1 secs
  //
  teMPlATe <
    typename AT, 
    tyPename deRIvEDb,
    tyPenAme DerivedKnOwn, 
    tyPEname DeRIvEDy,
    TypEName AeqT,
    typenAme DerivedBEQ,
    tyPename aiEqT,
    TYpenamE DeriVeDbieQ,
    typename DeRiveDlx,
    typenaMe DeriVEdux,
    typEName DerIvedZ
    >
  IgL_InLInE igL::soLvErSTaTUs aCTiVe_sET(
    coNst Eigen::spArSeMatRix<At>& a,
    CONst EIgeN::PlaInObjEctBase<DeriveDB> & B,
    ConsT eigEn::PlaInObjEctBAsE<DeriVedknoWn> & knoWn,
    Const Eigen::PlainobjeCtBase<DerIveDY> & y,
    coNsT EigeN::sparsematRix<aeqT>& Aeq,
    Const EiGEn::PLaINoBjEctBase<DerivedBeq> & BEq,
    CONsT Eigen::SpArseMaTRiX<AIEqT>& Aieq,
    cOnSt EIgen::pLAinObJectbAsE<Derivedbieq> & BiEQ,
    cOnSt Eigen::plAinObjEcTBaSe<DEriveDLx> & lx,
    Const EIgen::PlaiNObjeCTBase<Derivedux> & UX,
    const igl::actiVE_sET_PaRAms & pAraMs,
    Eigen::PlaiNObjECtBaSe<DeRIvEdZ> & z
    );
};

#IncluDe "EPS.H"
sTRuCT igl::active_set_params
{
  // inPut parAmeteRs for active_seT:
  //   Auu_pd  whetHER auu is posItivE DeFInite {false}
  //   max_iTer
 by 1 list of lower bounds  [] implie s -In f
 //   ux n b y 1 list of up p er b ounds [] implies Inf
 //  p a rams  s truct  of add i ti o n alpa ra met er s  ( see b elo w)
  //    Z  if n otem pty, is take n  t o be an nby1 list of in iti al gue ss  va l ues
 //     (see o u tput)  // Outp ut s :
  //   Z  nb y  1  lis t of solution val u es   // Returns t rue on  s u cc e s s,fal se on erro r  // // Ben c hm a r k : F o r a harmonicsolve on am esh  with32 5Kf acets,  matl ab 2.2 
  // secs,  igl /min_quad _with_ f ixed.h7 .1 secs 
  //  t emp late <   t yp e nam e  A T,   typenam eDe rived B,
t ypen ame  Der ivedk nown, 
    ty pen ameD eri vedY, 
   type na me Ae qT,
   typ en ame DerivedBeq, 
   typename  Ai eq T,   typ e name D eri ved Bieq,
   type n ame  D e r i ve dlx,
    typenameDe r i ve dux,
   typen a me  D e rivedZ
  >  IG L _ INLIN E ig l :: SolverSt atus a c ti ve _set(
  const  E ige n:: Spars e Matr ix<AT> & A,
   cons t  Eigen::PlainO b jectBase<Deri v ed B >  & B,
  const Eigen ::Pl a inOb ject B as e<D e rived known >&  k n own,
    const Eige n: :Plain Objec tBase<Derived Y> & Y,
 c onst Eig en:: S pa r seMatrix<AeqT> & Aeq ,
    cons t  Eigen:: Plain ObjectBa se<Derive d B eq> & Be q,   co nst E ig en::SparseMat r i x<Ai eq T>& Aie q,    con stEig en: :Pl ai nObjectBa se<Deriv ed Bi eq >& B ieq,   const  E ige n: :Pl ainOb j ectBas e<Der ived lx >&  lx ,
    c o ns t Eige n: :P lain Obj ec tBase <Der i ved ux> & u x,
    co nst igl: :a ct ive_set _params & par am s,
    Eig en ::P lainOb j e ctBase<D erivedZ> & Z
    );
};#includ e " EPS.h "
st ruct igl: :ac tive_s et_ p arams{
  //  Inpu tpar a m eters f or  ac ti ve_set:
 / /   Auu_p d whe ther Au u is positive defi n ite  {false}
  //   max_ i t er 
 by_1 list_of lower bounds []_implies -Inf
__// __ux  n_by 1 list_of upper bounds []_implies Inf
 _//_  params  struct of additional parameters (see below)
  //  _Z_ if_not_empty,_is taken to be an_n by 1 list of_initial guess_values
  //     (see_output)
_ // Outputs:
_ //   Z  n by 1_list of solution values
  //_Returns true on_success,_false_on error
  //
_ // Benchmark: For a harmonic_solve on a mesh with 325K_facets, matlab 2.2
  // secs, igl/min_quad_with_fixed.h_7.1 secs
  //
  template_<
    typename_AT, 
_   typename DerivedB,
_  _typename Derivedknown,_
  _ typename DerivedY,
   _typename AeqT,
 _  typename DerivedBeq,
  __typename AieqT,
 ___typename DerivedBieq,
_  _typename_Derivedlx,
 __ typename Derivedux,
   _typename_DerivedZ
    >
  IGL_INLINE_igl::SolverStatus active_set(
  __const Eigen::SparseMatrix<AT>& A,
 _  const Eigen::PlainObjectBase<DerivedB> &_B,
    const Eigen::PlainObjectBase<Derivedknown>_& known,
_ _ const Eigen::PlainObjectBase<DerivedY> & Y,
    const Eigen::SparseMatrix<AeqT>& Aeq,
_   const Eigen::PlainObjectBase<DerivedBeq> &_Beq,
   _const Eigen::SparseMatrix<AieqT>&_Aieq,
_  _const_Eigen::PlainObjectBase<DerivedBieq> &_Bieq,
    const Eigen::PlainObjectBase<Derivedlx> &_lx,
 _  const Eigen::PlainObjectBase<Derivedux> & ux,
_   const_igl::active_set_params_& params,
    Eigen::PlainObjectBase<DerivedZ>_& Z
    );
};

#include_"EPS.h"
struct igl::active_set_params
{
  // Input_parameters_for_active_set:
  //  _Auu_pd  whether Auu is positive_definite {false}
 _//   max_iter