2 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
6 * A PICTURE RECONSTRUCTION PROGRAM *
8 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
10 quad.h,v 1.7 2009/09/26 14:54:43 jklukowska Exp
18 /// common /quad/ aopt,bopt,copt,dopt,eopt,fopt,gopt,hopt,popt,
19 /// 1 l,bw1,bw2,bw3,bbcon,t,cw1,cw2,cw3,cbcon,sa,sb,sc,xbar,error,
20 /// 2 iterx,period,z,eta,del,order
21 /// integer aopt,bopt,copt,dopt,eopt,fopt,gopt,hopt,popt,l,bbcon,
22 /// 1 t,cbcon,iterx,period,del,order
23 /// real bw1,bw2,bw3,cw1,cw2,cw3,sa,sb,sc,xbar,error,z,eta
29 INTEGER aopt,bopt,copt,dopt,eopt,fopt,gopt,hopt,popt,l,bbcon;
30 INTEGER t,cbcon,iterx,period;
34 REAL bw1,bw2,bw3,cw1,cw2,cw3,sa,sb,sc,xbar,error_1,z,eta;
38 class quad_class: public alg_class
43 INTEGER aopt,bopt,copt,dopt,eopt,fopt,gopt,hopt,popt,l,bbcon;
44 INTEGER t,cbcon,iterx,period;
61 REAL bw1,bw2,bw3,cw1,cw2,cw3,sa,sb,sc,xbar,error_1,z,eta;
71 REAL* delta, //corr. hstau
72 REAL* gma, //corr. hstau
87 REAL* delta, //corr. hstau
88 REAL* gma, //corr. hstau
101 REAL* estmin, //corr.hstau
102 REAL* lambda, //corr.hstau
108 void matrix(REAL* y, REAL* x, REAL* d, INTEGER* list, REAL* weight, REAL* s);
109 void mtamx(REAL* x, REAL* y, INTEGER* list, REAL* weight);
110 void costfn(REAL* recon, REAL* r, REAL* g, REAL* cg, REAL* w, INTEGER* list, REAL* weight, REAL* f);
111 void resedu(REAL* x, REAL* r, REAL* w, REAL* d, INTEGER* list, REAL* weight, REAL* f, REAL* s);
112 void adsmos(REAL* x, INTEGER m, REAL w1, REAL w2, REAL w3, INTEGER bcon);
113 void badsmos(REAL* x, INTEGER m, REAL w1, REAL w2, INTEGER bcon); // blob version of adsmos
114 void mtamxp(REAL* x, REAL* y, INTEGER* list, REAL* weight);
115 void semi(REAL* delta, REAL* gma, BOOLEAN* alg);
116 void prsemi(REAL estmin, REAL lambda, BOOLEAN* alg, REAL* delta);
117 void deset(REAL* d, INTEGER* list, REAL* weight, BOOLEAN* alg);
118 void nxtpct(REAL* x, REAL* cg, REAL delta);
119 REAL uerror(INTEGER np, INTEGER nr, REAL pd, BOOLEAN* alg);
120 void sset(REAL* s, INTEGER* list, REAL* weight, BOOLEAN* alg);
121 void seset(REAL* s, INTEGER* list, REAL* weight, BOOLEAN* alg);
122 void dset(REAL* d, INTEGER* list, REAL* weight, BOOLEAN* alg);
123 void back(REAL* x, INTEGER* list, REAL* weight);
124 void blurck(REAL b1, REAL b2, REAL b3, INTEGER bcon, INTEGER npow, BOOLEAN* alg1);
125 void nextr(REAL* r, REAL* v, REAL delta, INTEGER area);
126 void prpict(INTEGER n, REAL* x, INTEGER it, const CHAR* l, const CHAR* m);
127 void bprpict(INTEGER M, INTEGER H, REAL* x, INTEGER it, const CHAR* l, const CHAR* m); // blob version of prpict
130 BOOLEAN Run(REAL* recon, INTEGER* list, REAL* weight, INTEGER iter);