1extern int dbcg_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
2 long* ja,
double* a,
const long* isym,
int ( *matvec )(),
int ( *mttvec )(),
3 int ( *msolve )(),
int ( *mtsolv )(),
const long* itol,
const double* tol,
4 const long* itmax,
long*
iter,
double* err,
long* ierr,
const long* iunit,
5 double* r,
double* z,
double* p,
double* rr,
double* zz,
double* pp,
6 double* dz,
double* rwork,
long* iwork );
7extern int dsdbcg_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
8 long* ja,
double* a,
const long* isym,
const long* itol,
const double* tol,
9 const long* itmax,
long*
iter,
double* err,
long* ierr,
const long* iunit,
10 double* rwork,
const long* lenw,
long* iwork,
const long* leniw );
11extern int dslubc_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
12 long* ja,
double* a,
const long* isym,
const long* itol,
const double* tol,
13 const long* itmax,
long*
iter,
double* err,
long* ierr,
const long* iunit,
14 double* rwork,
const long* lenw,
long* iwork,
const long* leniw );
15extern long isdbcg_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
16 long* ja,
double* a,
const long* isym,
int ( *msolve )(),
17 const long* itol,
const double* tol,
const long* itmax,
long*
iter,
18 double* err,
long* ierr,
const long* iunit,
double* r,
double* z,
19 double* p,
double* rr,
double* zz,
double* pp,
double* dz,
double* rwork,
20 long* iwork,
double* ak,
double* bk,
double* bnrm,
double* solnrm );
36extern int dcg_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
37 long* ja,
double* a,
const long* isym,
int ( *matvec )(),
int ( *msolve )(),
38 const long* itol,
const double* tol,
const long* itmax,
long*
iter,
39 double* err,
long* ierr,
const long* iunit,
double* r,
double* z,
double* p,
40 double* dz,
double* rwork,
long* iwork );
41extern int dsdcg_(
const int*
n,
const double* b,
double* x,
const long* nelt,
int* ia,
42 int* ja,
double* a,
const long* isym,
const long* itol,
const double* tol,
43 const long* itmax,
long*
iter,
double* err,
long* ierr,
const long* iunit,
44 double* rwork,
const long* lenw,
long* iwork,
const long* leniw );
45extern int dsiccg_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
46 long* ja,
double* a,
const long* isym,
const long* itol,
const double* tol,
47 const long* itmax,
long*
iter,
double* err,
long* ierr,
const long* iunit,
48 double* rwork,
const long* lenw,
long* iwork,
const long* leniw );
49extern long isdcg_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
50 long* ja,
double* a,
const long* isym,
int ( *msolve )(),
const long* itol,
51 const double* tol,
const long* itmax,
long*
iter,
double* err,
long* ierr,
52 const long* iunit,
double* r,
double* z,
double* p,
double* dz,
53 double* rwork,
long* iwork,
double* ak,
double* bk,
double* bnrm,
69extern int dcgn_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
70 long* ja,
double* a,
const long* isym,
int ( *matvec )(),
int ( *mttvec )(),
71 int ( *msolve )(),
const long* itol,
const double* tol,
const long* itmax,
72 long*
iter,
double* err,
long* ierr,
const long* iunit,
double* r,
double* z,
73 double* p,
double* atp,
double* atz,
double* dz,
double* atdz,
double* rwork,
75extern int dsdcgn_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
76 long* ja,
double* a,
const long* isym,
const long* itol,
const double* tol,
77 const long* itmax,
long*
iter,
double* err,
long* ierr,
const long* iunit,
78 double* rwork,
const long* lenw,
long* iwork,
const long* leniw );
79extern int dslucn_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
80 long* ja,
double* a,
const long* isym,
const long* itol,
const double* tol,
81 const long* itmax,
long*
iter,
double* err,
long* ierr,
const long* iunit,
82 double* rwork,
const long* lenw,
long* iwork,
const long* leniw );
83extern long isdcgn_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
84 long* ja,
double* a,
const long* isym,
int ( *matvec )(),
85 int ( *mttvec )(),
int ( *msolve )(),
const long* itol,
const double* tol,
86 const long* itmax,
long*
iter,
double* err,
long* ierr,
const long* iunit,
87 double* r,
double* z,
double* p,
double* atp,
double* atz,
double* dz,
88 double* atdz,
double* rwork,
long* iwork,
double* ak,
double* bk,
89 double* bnrm,
double* solnrm );
104extern int dcgs_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
105 long* ja,
double* a,
const long* isym,
int ( *matvec )(),
int ( *msolve )(),
106 const long* itol,
const double* tol,
const long* itmax,
long*
iter,
107 double* err,
long* ierr,
const long* iunit,
double* r,
double* r0,
double* p,
108 double*
q,
double* u,
double* v1,
double* v2,
double* rwork,
long* iwork );
109extern int dsdcgs_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
110 long* ja,
double* a,
const long* isym,
const long* itol,
const double* tol,
111 const long* itmax,
long*
iter,
double* err,
long* ierr,
const long* iunit,
112 double* rwork,
const long* lenw,
long* iwork,
const long* leniw );
113extern int dslucs_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
114 long* ja,
double* a,
const long* isym,
const long* itol,
const double* tol,
115 const long* itmax,
long*
iter,
double* err,
long* ierr,
const long* iunit,
116 double* rwork,
const long* lenw,
long* iwork,
const long* leniw );
117extern long isdcgs_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
118 long* ja,
double* a,
const long* isym,
int ( *matvec )(),
119 int ( *msolve )(),
const long* itol,
const double* tol,
const long* itmax,
120 long*
iter,
double* err,
long* ierr,
const long* iunit,
double* r,
121 double* r0,
double* p,
double*
q,
double* u,
double* v1,
double* v2,
122 double* rwork,
long* iwork,
double* ak,
double* bk,
double* bnrm,
136extern int dgmres_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
137 long* ja,
double* a,
const long* isym,
int ( *matvec )(),
138 int ( *msolve )(),
const long* itol,
const double* tol,
const long* itmax,
139 long*
iter,
double* err,
long* ierr,
const long* iunit,
double* sb,
140 double* sx,
double* rgwk,
long* lrgw,
long* igwk,
long* ligw,
141 double* rwork,
long* iwork );
142extern int dsdgmr_(
const int*
n,
const double* b,
double* x,
const long* nelt,
int* ia,
143 int* ja,
double* a,
const long* isym,
long* nsave,
const long* itol,
144 const double* tol,
const long* itmax,
long*
iter,
double* err,
long* ierr,
145 const long* iunit,
double* rwork,
const long* lenw,
long* iwork,
147extern int dslugm_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
148 long* ja,
double* a,
const long* isym,
long* nsave,
const long* itol,
149 const double* tol,
const long* itmax,
long*
iter,
double* err,
long* ierr,
150 const long* iunit,
double* rwork,
const long* lenw,
long* iwork,
152extern int dhels_(
double* a,
long* lda,
const long*
n,
double*
q,
double* b );
153extern int dheqr_(
double* a,
long* lda,
const long*
n,
double*
q,
long* info,
long* ijob );
154extern int dorth_(
double* vnew,
double*
v,
double* hes,
const long*
n,
long* ll,
long* ldhes,
155 long* kmp,
double* snormw );
156extern int dpigmr_(
const long*
n,
double* r0,
double* sr,
double* sz,
long* jscal,
long* maxl,
157 long* maxlp1,
long* kmp,
long* nrsts,
long* jpre,
int ( *matvec )(),
158 int ( *msolve )(),
long* nmsl,
double* z,
double*
v,
double* hes,
159 double*
q,
long* lgmr,
double* rpar,
long* ipar,
double* wk,
double* dl,
160 double* rhol,
long* nrmax,
const double* b,
double* bnrm,
double* x,
161 double* xl,
const long* itol,
const double* tol,
const long* nelt,
162 long* ia,
long* ja,
double* a,
const long* isym,
const long* iunit,
163 long* iflag,
double* err );
164extern int drlcal_(
const long*
n,
long* kmp,
long* ll,
long* maxl,
double*
v,
double*
q,
165 double* rl,
double* snormw,
double* prod,
double* r0nrm );
166extern int dxlcal_(
const long*
n,
long* lgmr,
double* x,
double* xl,
double* zl,
double* hes,
167 long* maxlp1,
double*
q,
double*
v,
double* r0nrm,
double* wk,
double* sz,
168 long* jscal,
long* jpre,
int ( *msolve )(),
long* nmsl,
double* rpar,
169 long* ipar,
const long* nelt,
long* ia,
long* ja,
double* a,
long* isym );
170extern long isdgmr_(
const long*
n,
const double* b,
double* x,
double* xl,
const long* nelt,
171 double* ia,
double* ja,
double* a,
const long* isym,
int ( *msolve )(),
172 long* nmsl,
const long* itol,
const double* tol,
const long* itmax,
173 long*
iter,
double* err,
const long* iunit,
double* r,
double* z,
174 double* dz,
double* rwork,
double* iwork,
double* rnrm,
double* bnrm,
175 double* sb,
double* sx,
long* jscal,
long* kmp,
long* lgmr,
long* maxl,
176 long* maxlp1,
double*
v,
double*
q,
double* snormw,
double* prod,
177 double* r0nrm,
double* hes,
long* jpre );
193extern int dir_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
194 long* ja,
double* a,
const long* isym,
int ( *matvec )(),
int ( *msolve )(),
195 const long* itol,
const double* tol,
const long* itmax,
long*
iter,
196 double* err,
long* ierr,
const long* iunit,
double* r,
double* z,
double* dz,
197 double* rwork,
long* iwork );
198extern int dsjac_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
199 long* ja,
double* a,
const long* isym,
const long* itol,
const double* tol,
200 const long* itmax,
long*
iter,
double* err,
long* ierr,
const long* iunit,
201 double* rwork,
const long* lenw,
long* iwork,
const long* leniw );
202extern int dsgs_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
203 long* ja,
double* a,
const long* isym,
const long* itol,
const double* tol,
204 const long* itmax,
long*
iter,
double* err,
long* ierr,
const long* iunit,
205 double* rwork,
const long* lenw,
long* iwork,
const long* leniw );
206extern int dsilur_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
207 long* ja,
double* a,
const long* isym,
const long* itol,
const double* tol,
208 const long* itmax,
long*
iter,
double* err,
long* ierr,
const long* iunit,
209 double* rwork,
const long* lenw,
long* iwork,
const long* leniw );
210extern long isdir_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
211 long* ja,
double* a,
const long* isym,
int ( *msolve )(),
const long* itol,
212 const double* tol,
const long* itmax,
long*
iter,
double* err,
long* ierr,
213 const long* iunit,
double* r,
double* z,
double* dz,
double* rwork,
214 long* iwork,
double* bnrm,
double* solnrm );
227extern int dbhin_(
const long*
n,
const long* nelt,
long* ia,
long* ja,
double* a,
228 const long* isym,
double* soln,
double* rhs,
const long* iunit,
long* job );
229extern int dchkw_(
char* name,
long* lociw,
const long* leniw,
long* locw,
const long* lenw,
230 long* ierr,
long*
iter,
double* err,
int name_len );
231extern int qs2i1d_(
long* ia,
long* ja,
double* a,
const long*
n,
long* kflag );
232extern int ds2y_(
const long*
n,
const long* nelt,
long* ia,
long* ja,
double* a,
long* isym );
233extern int dcpplt_(
const long*
n,
const long* nelt,
long* ia,
long* ja,
double* a,
234 const long* isym,
const long* iunit );
235extern int dtout_(
const long*
n,
const long* nelt,
long* ia,
long* ja,
double* a,
236 const long* isym,
double* soln,
double* rhs,
const long* iunit,
long* job );
237extern int dtin_(
const long*
n,
const long* nelt,
long* ia,
long* ja,
double* a,
238 const long* isym,
double* soln,
double* rhs,
const long* iunit,
long* job );
242extern int dsds_(
const long*
n,
const long* nelt,
long* ia,
long* ja,
double* a,
243 const long* isym,
double* dinv );
244extern int dsdscl_(
const long*
n,
const long* nelt,
long* ia,
long* ja,
double* a,
245 const long* isym,
double* x,
const double* b,
double* dinv,
long* job,
247extern int dsd2s_(
const long*
n,
const long* nelt,
long* ia,
long* ja,
double* a,
248 const long* isym,
double* dinv );
249extern int ds2lt_(
const long*
n,
const long* nelt,
long* ia,
long* ja,
double* a,
250 const long* isym,
long* nel,
long* iel,
long* jel,
double* el );
251extern int dsics_(
const long*
n,
const long* nelt,
long* ia,
long* ja,
double* a,
252 const long* isym,
long* nel,
long* iel,
long* jel,
double* el,
double* d,
253 double* r,
long* iwarn );
254extern int dsilus_(
const long*
n,
const long* nelt,
long* ia,
long* ja,
double* a,
255 const long* isym,
long* nl,
long* il,
long* jl,
double* l,
double* dinv,
256 long* nu,
long* iu,
long* ju,
double* u,
long* nrow,
long* ncol );
259extern int dsmv_(
const long*
n,
double* x,
double* y,
const long* nelt,
long* ia,
long* ja,
260 double* a,
long* isym );
261extern int dsmtv_(
const long*
n,
double* x,
double* y,
const long* nelt,
long* ia,
long* ja,
262 double* a,
long* isym );
263extern int dsdi_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
264 long* ja,
double* a,
const long* isym,
double* rwork,
long* iwork );
265extern int dsli_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
266 long* ja,
double* a,
const long* isym,
double* rwork,
long* iwork );
267extern int dsli2_(
const long*
n,
const double* b,
double* x,
long* nel,
long* iel,
long* jel,
269extern int dsllti_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
270 long* ja,
double* a,
const long* isym,
double* rwork,
long* iwork );
271extern int dllti2_(
const long*
n,
const double* b,
double* x,
long* nel,
long* iel,
long* jel,
272 double* el,
double* dinv );
273extern int dslui_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
274 long* ja,
double* a,
const long* isym,
double* rwork,
long* iwork );
275extern int dslui2_(
const long*
n,
const double* b,
double* x,
long* il,
long* jl,
double* l,
276 double* dinv,
long* iu,
long* ju,
double* u );
277extern int dsluti_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
278 long* ja,
double* a,
const long* isym,
double* rwork,
long* iwork );
279extern int dslui4_(
const long*
n,
const double* b,
double* x,
long* il,
long* jl,
double* l,
280 double* dinv,
long* iu,
long* ju,
double* u );
281extern int dsmmti_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
282 long* ja,
double* a,
const long* isym,
double* rwork,
long* iwork );
283extern int dsmmi2_(
const long*
n,
const double* b,
double* x,
long* il,
long* jl,
double* l,
284 double* dinv,
long* iu,
long* ju,
double* u );
285extern int domn_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
286 long* ja,
double* a,
const long* isym,
int ( *matvec )(),
int ( *msolve )(),
287 long* nsave,
const long* itol,
const double* tol,
const long* itmax,
288 long*
iter,
double* err,
long* ierr,
const long* iunit,
double* r,
double* z,
289 double* p,
double* ap,
double* emap,
double* dz,
double* csav,
double* rwork,
291extern int dsdomn_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
292 long* ja,
double* a,
const long* isym,
long* nsave,
const long* itol,
293 const double* tol,
const long* itmax,
long*
iter,
double* err,
long* ierr,
294 const long* iunit,
double* rwork,
const long* lenw,
long* iwork,
296extern int dsluom_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
297 long* ja,
double* a,
const long* isym,
long* nsave,
const long* itol,
298 const double* tol,
const long* itmax,
long*
iter,
double* err,
long* ierr,
299 const long* iunit,
double* rwork,
const long* lenw,
long* iwork,
301extern long isdomn_(
const long*
n,
const double* b,
double* x,
const long* nelt,
long* ia,
302 long* ja,
double* a,
const long* isym,
int ( *msolve )(),
long* nsave,
303 const long* itol,
const double* tol,
const long* itmax,
long*
iter,
304 double* err,
long* ierr,
const long* iunit,
double* r,
double* z,
305 double* p,
double* ap,
double* emap,
double* dz,
double* csav,
306 double* rwork,
long* iwork,
double* ak,
double* bnrm,
double* solnrm );
****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
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)