1extern int dbcg_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
int (*matvec)(),
int (*mttvec)(),
int (*msolve)(),
int (*mtsolv)(),
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *r,
double *z,
double *p,
double *rr,
double *zz,
double *pp,
double *dz,
double *rwork,
long *iwork);
2extern int dsdbcg_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
3extern int dslubc_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
4extern long isdbcg_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
int (*msolve)(),
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *r,
double *z,
double *p,
double *rr,
double *zz,
double *pp,
double *dz,
double *rwork,
long *iwork,
double *ak,
double *bk,
double *bnrm,
double *solnrm);
20extern int dcg_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
int (*matvec)(),
int (*msolve)(),
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *r,
double *z,
double *p,
double *dz,
double *rwork,
long *iwork);
21extern int dsdcg_(
const int *n,
const double *
b,
double *x,
const long *nelt,
int *ia,
int *ja,
double *a,
const long *isym,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
22extern int dsiccg_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
23extern long isdcg_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
int (*msolve)(),
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *r,
double *z,
double *p,
double *dz,
double *rwork,
long *iwork,
double *ak,
double *bk,
double *bnrm,
double *solnrm);
38extern int dcgn_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
int (*matvec)(),
int (*mttvec)(),
int (*msolve)(),
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *r,
double *z,
double *p,
double *atp,
double *atz,
double *dz,
double *atdz,
double *rwork,
long *iwork);
39extern int dsdcgn_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
40extern int dslucn_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
41extern long isdcgn_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
int (*matvec)(),
int (*mttvec)(),
int (*msolve)(),
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *r,
double *z,
double *p,
double *atp,
double *atz,
double *dz,
double *atdz,
double *rwork,
long *iwork,
double *ak,
double *bk,
double *bnrm,
double *solnrm);
56extern int dcgs_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
int (*matvec)(),
int (*msolve)(),
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *r,
double *r0,
double *p,
double *
q,
double *u,
double *v1,
double *v2,
double *rwork,
long *iwork);
57extern int dsdcgs_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
58extern int dslucs_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
59extern long isdcgs_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
int (*matvec)(),
int (*msolve)(),
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *r,
double *r0,
double *p,
double *
q,
double *u,
double *v1,
double *v2,
double *rwork,
long *iwork,
double *ak,
double *bk,
double *bnrm,
double *solnrm);
72extern int dgmres_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
int (*matvec)(),
int (*msolve)(),
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *sb,
double *sx,
double *rgwk,
long *lrgw,
long *igwk,
long *ligw,
double *rwork,
long *iwork);
73extern int dsdgmr_(
const int *n,
const double *
b,
double *x,
const long *nelt,
int *ia,
int *ja,
double *a,
const long *isym,
long *nsave,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
74extern int dslugm_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
long *nsave,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
75extern int dhels_(
double *a,
long *lda,
const long *n,
double *
q,
double *
b);
76extern int dheqr_(
double *a,
long *lda,
const long *n,
double *
q,
long *info,
long *ijob);
77extern int dorth_(
double *vnew,
double *
v,
double *hes,
const long *n,
long *ll,
long *ldhes,
long *kmp,
double *snormw);
78extern int dpigmr_(
const long *n,
double *r0,
double *sr,
double *sz,
long *jscal,
long *maxl,
long *maxlp1,
long *kmp,
long *nrsts,
long *jpre,
int (*matvec)(),
int (*msolve)(),
long *nmsl,
double *z,
double *
v,
double *hes,
double *
q,
long *lgmr,
double *rpar,
long *ipar,
double *wk,
double *dl,
double *rhol,
long *nrmax,
const double *
b,
double *bnrm,
double *x,
double *xl,
const long *itol,
const double *tol,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
const long *iunit,
long *iflag,
double *err);
79extern int drlcal_(
const long *n,
long *kmp,
long *ll,
long *maxl,
double *
v,
double *
q,
double *rl,
double *snormw,
double *prod,
double *r0nrm);
80extern int dxlcal_(
const long *n,
long *lgmr,
double *x,
double *xl,
double *zl,
double *hes,
long *maxlp1,
double *
q,
double *
v,
double *r0nrm,
double *wk,
double *sz,
long *jscal,
long *jpre,
int (*msolve)(),
long *nmsl,
double *rpar,
long *ipar,
const long *nelt,
long *ia,
long *ja,
double *a,
long *isym);
81extern long isdgmr_(
const long *n,
const double *
b,
double *x,
double *xl,
const long *nelt,
double *ia,
double *ja,
double *a,
const long *isym,
int (*msolve)(),
long *nmsl,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
const long *iunit,
double *r,
double *z,
double *dz,
double *rwork,
double *iwork,
double *rnrm,
double *bnrm,
double *sb,
double *sx,
long *jscal,
long *kmp,
long *lgmr,
long *maxl,
long *maxlp1,
double *
v,
double *
q,
double *snormw,
double *prod,
double *r0nrm,
double *hes,
long *jpre);
97extern int dir_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
int (*matvec)(),
int (*msolve)(),
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *r,
double *z,
double *dz,
double *rwork,
long *iwork);
98extern int dsjac_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
99extern int dsgs_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
100extern int dsilur_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
101extern long isdir_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
int (*msolve)(),
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *r,
double *z,
double *dz,
double *rwork,
long *iwork,
double *bnrm,
double *solnrm);
114extern int dbhin_(
const long *n,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
double *soln,
double *rhs,
const long *iunit,
long *job);
115extern int dchkw_(
char *name,
long *lociw,
const long *leniw,
long *locw,
const long *lenw,
long *ierr,
long *
iter,
double *err,
int name_len);
116extern int qs2i1d_(
long *ia,
long *ja,
double *a,
const long *n,
long *kflag);
117extern int ds2y_(
const long *n,
const long *nelt,
long *ia,
long *ja,
double *a,
long *isym);
118extern int dcpplt_(
const long *n,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
const long *iunit);
119extern int dtout_(
const long *n,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
double *soln,
double *rhs,
const long *iunit,
long *job);
120extern int dtin_(
const long *n,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
double *soln,
double *rhs,
const long *iunit,
long *job);
124extern int dsds_(
const long *n,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
double *dinv);
125extern int dsdscl_(
const long *n,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
double *x,
const double *
b,
double *dinv,
long *job,
long *itol);
126extern int dsd2s_(
const long *n,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
double *dinv);
127extern int ds2lt_(
const long *n,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
long *nel,
long *iel,
long *jel,
double *el);
128extern int dsics_(
const long *n,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
long *nel,
long *iel,
long *jel,
double *el,
double *d,
double *r,
long *iwarn);
129extern int dsilus_(
const long *n,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
long *nl,
long *il,
long *jl,
double *l,
double *dinv,
long *nu,
long *iu,
long *ju,
double *u,
long *nrow,
long *ncol);
132extern int dsmv_(
const long *n,
double *x,
double *
y,
const long *nelt,
long *ia,
long *ja,
double *a,
long *isym);
133extern int dsmtv_(
const long *n,
double *x,
double *
y,
const long *nelt,
long *ia,
long *ja,
double *a,
long *isym);
134extern int dsdi_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
double *rwork,
long *iwork);
135extern int dsli_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
double *rwork,
long *iwork);
136extern int dsli2_(
const long *n,
const double *
b,
double *x,
long *nel,
long *iel,
long *jel,
double *el);
137extern int dsllti_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
double *rwork,
long *iwork);
138extern int dllti2_(
const long *n,
const double *
b,
double *x,
long *nel,
long *iel,
long *jel,
double *el,
double *dinv);
139extern int dslui_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
double *rwork,
long *iwork);
140extern int dslui2_(
const long *n,
const double *
b,
double *x,
long *il,
long *jl,
double *l,
double *dinv,
long *iu,
long *ju,
double *u);
141extern int dsluti_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
double *rwork,
long *iwork);
142extern int dslui4_(
const long *n,
const double *
b,
double *x,
long *il,
long *jl,
double *l,
double *dinv,
long *iu,
long *ju,
double *u);
143extern int dsmmti_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
double *rwork,
long *iwork);
144extern int dsmmi2_(
const long *n,
const double *
b,
double *x,
long *il,
long *jl,
double *l,
double *dinv,
long *iu,
long *ju,
double *u);
145extern int domn_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
int (*matvec)(),
int (*msolve)(),
long *nsave,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *r,
double *z,
double *p,
double *ap,
double *emap,
double *dz,
double *csav,
double *rwork,
long *iwork);
146extern int dsdomn_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
long *nsave,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
147extern int dsluom_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
long *nsave,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *rwork,
const long *lenw,
long *iwork,
const long *leniw);
148extern long isdomn_(
const long *n,
const double *
b,
double *x,
const long *nelt,
long *ia,
long *ja,
double *a,
const long *isym,
int (*msolve)(),
long *nsave,
const long *itol,
const double *tol,
const long *itmax,
long *
iter,
double *err,
long *ierr,
const long *iunit,
double *r,
double *z,
double *p,
double *ap,
double *emap,
double *dz,
double *csav,
double *rwork,
long *iwork,
double *ak,
double *bnrm,
double *solnrm);
int dcgn_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, int(*matvec)(), int(*mttvec)(), int(*msolve)(), const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *r, double *z, double *p, double *atp, double *atz, double *dz, double *atdz, double *rwork, long *iwork)
int dsdi_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, double *rwork, long *iwork)
int dsmmti_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, double *rwork, long *iwork)
int dsdcg_(const int *n, const double *b, double *x, const long *nelt, int *ia, int *ja, double *a, const long *isym, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int dtin_(const long *n, const long *nelt, long *ia, long *ja, double *a, const long *isym, double *soln, double *rhs, const long *iunit, long *job)
long isdcg_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, int(*msolve)(), const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *r, double *z, double *p, double *dz, double *rwork, long *iwork, double *ak, double *bk, double *bnrm, double *solnrm)
long isdcgs_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, int(*matvec)(), int(*msolve)(), const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *r, double *r0, double *p, double *q, double *u, double *v1, double *v2, double *rwork, long *iwork, double *ak, double *bk, double *bnrm, double *solnrm)
int dheqr_(double *a, long *lda, const long *n, double *q, long *info, long *ijob)
int qs2i1d_(long *ia, long *ja, double *a, const long *n, long *kflag)
int dslucn_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int dir_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, int(*matvec)(), int(*msolve)(), const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *r, double *z, double *dz, double *rwork, long *iwork)
int dsjac_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int dsdscl_(const long *n, const long *nelt, long *ia, long *ja, double *a, const long *isym, double *x, const double *b, double *dinv, long *job, long *itol)
int dslui2_(const long *n, const double *b, double *x, long *il, long *jl, double *l, double *dinv, long *iu, long *ju, double *u)
int dtout_(const long *n, const long *nelt, long *ia, long *ja, double *a, const long *isym, double *soln, double *rhs, const long *iunit, long *job)
int dsdbcg_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int dslui_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, double *rwork, long *iwork)
int dsli_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, double *rwork, long *iwork)
int dsds_(const long *n, const long *nelt, long *ia, long *ja, double *a, const long *isym, double *dinv)
int dbhin_(const long *n, const long *nelt, long *ia, long *ja, double *a, const long *isym, double *soln, double *rhs, const long *iunit, long *job)
int dsdcgs_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
long isdgmr_(const long *n, const double *b, double *x, double *xl, const long *nelt, double *ia, double *ja, double *a, const long *isym, int(*msolve)(), long *nmsl, const long *itol, const double *tol, const long *itmax, long *iter, double *err, const long *iunit, double *r, double *z, double *dz, double *rwork, double *iwork, double *rnrm, double *bnrm, double *sb, double *sx, long *jscal, long *kmp, long *lgmr, long *maxl, long *maxlp1, double *v, double *q, double *snormw, double *prod, double *r0nrm, double *hes, long *jpre)
int drlcal_(const long *n, long *kmp, long *ll, long *maxl, double *v, double *q, double *rl, double *snormw, double *prod, double *r0nrm)
int domn_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, int(*matvec)(), int(*msolve)(), long *nsave, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *r, double *z, double *p, double *ap, double *emap, double *dz, double *csav, double *rwork, long *iwork)
int dorth_(double *vnew, double *v, double *hes, const long *n, long *ll, long *ldhes, long *kmp, double *snormw)
int dsli2_(const long *n, const double *b, double *x, long *nel, long *iel, long *jel, double *el)
int dsdcgn_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int dsmmi2_(const long *n, const double *b, double *x, long *il, long *jl, double *l, double *dinv, long *iu, long *ju, double *u)
int dsmv_(const long *n, double *x, double *y, const long *nelt, long *ia, long *ja, double *a, long *isym)
int dsdgmr_(const int *n, const double *b, double *x, const long *nelt, int *ia, int *ja, double *a, const long *isym, long *nsave, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int dsilur_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int dsmtv_(const long *n, double *x, double *y, const long *nelt, long *ia, long *ja, double *a, long *isym)
long isdbcg_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, int(*msolve)(), const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *r, double *z, double *p, double *rr, double *zz, double *pp, double *dz, double *rwork, long *iwork, double *ak, double *bk, double *bnrm, double *solnrm)
int dpigmr_(const long *n, double *r0, double *sr, double *sz, long *jscal, long *maxl, long *maxlp1, long *kmp, long *nrsts, long *jpre, int(*matvec)(), int(*msolve)(), long *nmsl, double *z, double *v, double *hes, double *q, long *lgmr, double *rpar, long *ipar, double *wk, double *dl, double *rhol, long *nrmax, const double *b, double *bnrm, double *x, double *xl, const long *itol, const double *tol, const long *nelt, long *ia, long *ja, double *a, const long *isym, const long *iunit, long *iflag, double *err)
int dhels_(double *a, long *lda, const long *n, double *q, double *b)
int dsics_(const long *n, const long *nelt, long *ia, long *ja, double *a, const long *isym, long *nel, long *iel, long *jel, double *el, double *d, double *r, long *iwarn)
int dxlcal_(const long *n, long *lgmr, double *x, double *xl, double *zl, double *hes, long *maxlp1, double *q, double *v, double *r0nrm, double *wk, double *sz, long *jscal, long *jpre, int(*msolve)(), long *nmsl, double *rpar, long *ipar, const long *nelt, long *ia, long *ja, double *a, long *isym)
long isdomn_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, int(*msolve)(), long *nsave, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *r, double *z, double *p, double *ap, double *emap, double *dz, double *csav, double *rwork, long *iwork, double *ak, double *bnrm, double *solnrm)
int dllti2_(const long *n, const double *b, double *x, long *nel, long *iel, long *jel, double *el, double *dinv)
long isdcgn_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, int(*matvec)(), int(*mttvec)(), int(*msolve)(), const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *r, double *z, double *p, double *atp, double *atz, double *dz, double *atdz, double *rwork, long *iwork, double *ak, double *bk, double *bnrm, double *solnrm)
int dsluti_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, double *rwork, long *iwork)
int dsilus_(const long *n, const long *nelt, long *ia, long *ja, double *a, const long *isym, long *nl, long *il, long *jl, double *l, double *dinv, long *nu, long *iu, long *ju, double *u, long *nrow, long *ncol)
int dsiccg_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int dslugm_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, long *nsave, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int dslui4_(const long *n, const double *b, double *x, long *il, long *jl, double *l, double *dinv, long *iu, long *ju, double *u)
int dgmres_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, int(*matvec)(), int(*msolve)(), const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *sb, double *sx, double *rgwk, long *lrgw, long *igwk, long *ligw, double *rwork, long *iwork)
int dcgs_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, int(*matvec)(), int(*msolve)(), const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *r, double *r0, double *p, double *q, double *u, double *v1, double *v2, double *rwork, long *iwork)
int dsdomn_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, long *nsave, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int dslubc_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int ds2y_(const long *n, const long *nelt, long *ia, long *ja, double *a, long *isym)
long isdir_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, int(*msolve)(), const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *r, double *z, double *dz, double *rwork, long *iwork, double *bnrm, double *solnrm)
int dsluom_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, long *nsave, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int dcpplt_(const long *n, const long *nelt, long *ia, long *ja, double *a, const long *isym, const long *iunit)
int ds2lt_(const long *n, const long *nelt, long *ia, long *ja, double *a, const long *isym, long *nel, long *iel, long *jel, double *el)
int dsgs_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int dbcg_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, int(*matvec)(), int(*mttvec)(), int(*msolve)(), int(*mtsolv)(), const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *r, double *z, double *p, double *rr, double *zz, double *pp, double *dz, double *rwork, long *iwork)
int dsllti_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, double *rwork, long *iwork)
int dchkw_(char *name, long *lociw, const long *leniw, long *locw, const long *lenw, long *ierr, long *iter, double *err, int name_len)
int dslucs_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *rwork, const long *lenw, long *iwork, const long *leniw)
int dsd2s_(const long *n, const long *nelt, long *ia, long *ja, double *a, const long *isym, double *dinv)
int dcg_(const long *n, const double *b, double *x, const long *nelt, long *ia, long *ja, double *a, const long *isym, int(*matvec)(), int(*msolve)(), const long *itol, const double *tol, const long *itmax, long *iter, double *err, long *ierr, const long *iunit, double *r, double *z, double *p, double *dz, double *rwork, long *iwork)
****INTEGER imax DOUBLE PRECISION m_pi *DOUBLE PRECISION m_amfin DOUBLE PRECISION m_Chfin DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_sinw2 DOUBLE PRECISION m_GFermi DOUBLE PRECISION m_MfinMin DOUBLE PRECISION m_ta2 INTEGER m_out INTEGER m_KeyFSR INTEGER m_KeyQCD *COMMON c_Semalib $ !copy of input $ !CMS energy $ !beam mass $ !final mass $ !beam charge $ !final charge $ !smallest final mass $ !Z mass $ !Z width $ !EW mixing angle $ !Gmu Fermi $ alphaQED at q
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyWtm INTEGER m_idyfs DOUBLE PRECISION m_zini DOUBLE PRECISION m_q2 DOUBLE PRECISION m_Wt_KF DOUBLE PRECISION m_WtCut INTEGER m_KFfin *COMMON c_KarLud $ !Input CMS energy[GeV] $ !CMS energy after beam spread beam strahlung[GeV] $ !Beam energy spread[GeV] $ !z boost due to beam spread $ !electron beam mass *ff pair spectrum $ !minimum v