int main()
{
auto a1 = F1();
auto a2 = F2();
auto a3 = F3();
auto a4 = F4();
auto a5 = F5();
return 0;
}
seconlevel_f4::seconlevel_f4(QWidget *parent) :
QWidget(parent),
ui(new Ui::seconlevel_f4)
{
ui->setupUi(this);
group = new GroupDialog(parent);
group->hide();
connect(group ,SIGNAL(finished(int)) ,this ,SLOT(f1_Done(int)));
autotab = new AutoDialog(parent);
autotab->hide();
connect(autotab ,SIGNAL(finished(int)) ,this ,SLOT(f2_Done(int)));
jump = new JumpDialog(parent);
jump->hide();
connect(jump ,SIGNAL(finished(int)) ,this ,SLOT(f6_Done(int)));
connect(ui->pushButton_F1,SIGNAL(clicked()),this ,SLOT(F1()));
ui->pushButton_F1->setCheckable(true);
connect(ui->pushButton_F2,SIGNAL(clicked()),this ,SLOT(F2()));
ui->pushButton_F2->setCheckable(true);
connect(ui->pushButton_F3,SIGNAL(clicked()),this ,SLOT(F3()));
ui->pushButton_F3->setCheckable(true);
connect(ui->pushButton_F4,SIGNAL(clicked()),this ,SLOT(F4()));
ui->pushButton_F4->setCheckable(true);
connect(ui->pushButton_F5,SIGNAL(clicked()),this ,SLOT(F5()));
ui->pushButton_F5->setCheckable(true);
connect(ui->pushButton_F6,SIGNAL(clicked()),this ,SLOT(F6()));
ui->pushButton_F6->setCheckable(true);
connect(ui->pushButton_F7,SIGNAL(clicked()),this ,SLOT(F7()));
ui->pushButton_F7->setCheckable(true);
connect(this ,SIGNAL(enter(int)),parent ,SLOT(funcbarUpdate(int)));
connect(ui->pushButton_F8,SIGNAL(clicked()),this ,SLOT(F8()));
/*数据表模式改变信号*/
connect(this ,SIGNAL(stateChange(char)) ,parent ,SLOT(tableStateUpdate(char)));
/*段选信号*/
connect(this ,SIGNAL(selectRows(bool)) ,parent ,SLOT(tableSelect(bool)));
setFocusPolicy(Qt::NoFocus);
}
int test_main(int,char *[])
{
bu::quantity<mixed_length> a1(2.0 * mixed_length());
bu::quantity<si_area> a2(a1);
BOOST_CHECK((std::abs(a2.value() - .02) < .0001));
bu::quantity<mixed_length> a3(a2);
BOOST_CHECK((std::abs(a3.value() - 2.0) < .0001));
bu::quantity<mixed_energy_1> e1(2.0 * mixed_energy_1());
bu::quantity<mixed_energy_2> e2(e1);
BOOST_CHECK((std::abs(e2.value() - 20.0) < .0001));
bu::quantity<bu::si::energy> e3(e1);
BOOST_CHECK((std::abs(e3.value() - .0002) < .0001));
bu::quantity<mixed_energy_2> e4(e3);
BOOST_CHECK((std::abs(e4.value() - 20.0) < .0001));
bu::quantity<bu::cgs::force> F0 = 20 * bu::cgs::dyne;
BOOST_CHECK((std::abs(F0.value() - 20.0) < .0001));
bu::quantity<bu::si::force> F3(F0);
BOOST_CHECK((std::abs(F3.value() - 2.0e-4) < .000000001));
bu::quantity<bu::si::force> F5(20 * bu::cgs::dyne);
BOOST_CHECK((std::abs(F5.value() - 2.0e-4) < .000000001));
bu::quantity<bu::si::dimensionless> dimensionless_test1(1.0*bu::cgs::dyne/bu::si::newton);
BOOST_CHECK(dimensionless_test1 == 1e-5);
typedef bu::multiply_typeof_helper<bu::si::length, bu::cgs::length>::type m_cm;
typedef bu::divide_typeof_helper<m_cm, m_cm>::type heterogeneous_dimensionless;
bu::quantity<heterogeneous_dimensionless> dimensionless_test2(1.0*bu::cgs::dyne/bu::si::newton);
BOOST_CHECK(dimensionless_test2.value() == 1e-5);
bu::quantity<bu::divide_typeof_helper<bu::cgs::force, bu::si::force>::type> dimensionless_test3(dimensionless_test2);
BOOST_UNITS_CHECK_CLOSE(dimensionless_test3.value(), 1.0);
//m/cm -> g/kg
bu::quantity<bu::divide_typeof_helper<bu::si::length, bu::cgs::length>::type> dimensionless_test4(2.0 * bu::si::meters / bu::cgs::centimeters);
bu::quantity<bu::divide_typeof_helper<bu::cgs::mass, bu::si::mass>::type> dimensionless_test5(dimensionless_test4);
BOOST_UNITS_CHECK_CLOSE(dimensionless_test5.value(), 2e5);
return(0);
}
/* The RIPEMD160 compression function. */
static inline void rmd160_compress(struct rmd160_ctx *ctx, uint32_t *buf)
{
uint8_t w, round;
uint32_t T;
uint32_t AL, BL, CL, DL, EL; /* left line */
uint32_t AR, BR, CR, DR, ER; /* right line */
uint32_t X[16];
/* Byte-swap the buffer if we're on a big-endian machine */
cpu_to_le32_array(X, buf, 16);
/* Load the left and right lines with the initial state */
AL = AR = ctx->h[0];
BL = BR = ctx->h[1];
CL = CR = ctx->h[2];
DL = DR = ctx->h[3];
EL = ER = ctx->h[4];
/* Round 1 */
round = 0;
for (w = 0; w < 16; w++) { /* left line */
T = rol32(AL + F1(BL, CL, DL) + X[RL[round][w]] + KL[round], SL[round][w]) + EL;
AL = EL; EL = DL; DL = rol32(CL, 10); CL = BL; BL = T;
}
for (w = 0; w < 16; w++) { /* right line */
T = rol32(AR + F5(BR, CR, DR) + X[RR[round][w]] + KR[round], SR[round][w]) + ER;
AR = ER; ER = DR; DR = rol32(CR, 10); CR = BR; BR = T;
}
/* Round 2 */
round++;
for (w = 0; w < 16; w++) { /* left line */
T = rol32(AL + F2(BL, CL, DL) + X[RL[round][w]] + KL[round], SL[round][w]) + EL;
AL = EL; EL = DL; DL = rol32(CL, 10); CL = BL; BL = T;
}
for (w = 0; w < 16; w++) { /* right line */
T = rol32(AR + F4(BR, CR, DR) + X[RR[round][w]] + KR[round], SR[round][w]) + ER;
AR = ER; ER = DR; DR = rol32(CR, 10); CR = BR; BR = T;
}
/* Round 3 */
round++;
for (w = 0; w < 16; w++) { /* left line */
T = rol32(AL + F3(BL, CL, DL) + X[RL[round][w]] + KL[round], SL[round][w]) + EL;
AL = EL; EL = DL; DL = rol32(CL, 10); CL = BL; BL = T;
}
for (w = 0; w < 16; w++) { /* right line */
T = rol32(AR + F3(BR, CR, DR) + X[RR[round][w]] + KR[round], SR[round][w]) + ER;
AR = ER; ER = DR; DR = rol32(CR, 10); CR = BR; BR = T;
}
/* Round 4 */
round++;
for (w = 0; w < 16; w++) { /* left line */
T = rol32(AL + F4(BL, CL, DL) + X[RL[round][w]] + KL[round], SL[round][w]) + EL;
AL = EL; EL = DL; DL = rol32(CL, 10); CL = BL; BL = T;
}
for (w = 0; w < 16; w++) { /* right line */
T = rol32(AR + F2(BR, CR, DR) + X[RR[round][w]] + KR[round], SR[round][w]) + ER;
AR = ER; ER = DR; DR = rol32(CR, 10); CR = BR; BR = T;
}
/* Round 5 */
round++;
for (w = 0; w < 16; w++) { /* left line */
T = rol32(AL + F5(BL, CL, DL) + X[RL[round][w]] + KL[round], SL[round][w]) + EL;
AL = EL; EL = DL; DL = rol32(CL, 10); CL = BL; BL = T;
}
for (w = 0; w < 16; w++) { /* right line */
T = rol32(AR + F1(BR, CR, DR) + X[RR[round][w]] + KR[round], SR[round][w]) + ER;
AR = ER; ER = DR; DR = rol32(CR, 10); CR = BR; BR = T;
}
/* Final mixing stage */
T = ctx->h[1] + CL + DR;
ctx->h[1] = ctx->h[2] + DL + ER;
ctx->h[2] = ctx->h[3] + EL + AR;
ctx->h[3] = ctx->h[4] + AL + BR;
ctx->h[4] = ctx->h[0] + BL + CR;
ctx->h[0] = T;
}
/* The RIPEMD160 compression function. Operates on self->buf */
static void ripemd160_compress(hash_state *self)
{
unsigned w, round;
uint32_t T;
uint32_t AL, BL, CL, DL, EL; /* left line */
uint32_t AR, BR, CR, DR, ER; /* right line */
uint32_t bufw[16];
for (w=0; w<16; w++)
bufw[w] = LOAD_U32_LITTLE(&self->buf[w*4]);
/* Load the left and right lines with the initial state */
AL = AR = self->h[0];
BL = BR = self->h[1];
CL = CR = self->h[2];
DL = DR = self->h[3];
EL = ER = self->h[4];
/* Round 1 */
round = 0;
for (w = 0; w < 16; w++) { /* left line */
T = ROL(SL[round][w], AL + F1(BL, CL, DL) + bufw[RL[round][w]] + KL[round]) + EL;
AL = EL; EL = DL; DL = ROL(10, CL); CL = BL; BL = T;
}
for (w = 0; w < 16; w++) { /* right line */
T = ROL(SR[round][w], AR + F5(BR, CR, DR) + bufw[RR[round][w]] + KR[round]) + ER;
AR = ER; ER = DR; DR = ROL(10, CR); CR = BR; BR = T;
}
/* Round 2 */
round++;
for (w = 0; w < 16; w++) { /* left line */
T = ROL(SL[round][w], AL + F2(BL, CL, DL) + bufw[RL[round][w]] + KL[round]) + EL;
AL = EL; EL = DL; DL = ROL(10, CL); CL = BL; BL = T;
}
for (w = 0; w < 16; w++) { /* right line */
T = ROL(SR[round][w], AR + F4(BR, CR, DR) + bufw[RR[round][w]] + KR[round]) + ER;
AR = ER; ER = DR; DR = ROL(10, CR); CR = BR; BR = T;
}
/* Round 3 */
round++;
for (w = 0; w < 16; w++) { /* left line */
T = ROL(SL[round][w], AL + F3(BL, CL, DL) + bufw[RL[round][w]] + KL[round]) + EL;
AL = EL; EL = DL; DL = ROL(10, CL); CL = BL; BL = T;
}
for (w = 0; w < 16; w++) { /* right line */
T = ROL(SR[round][w], AR + F3(BR, CR, DR) + bufw[RR[round][w]] + KR[round]) + ER;
AR = ER; ER = DR; DR = ROL(10, CR); CR = BR; BR = T;
}
/* Round 4 */
round++;
for (w = 0; w < 16; w++) { /* left line */
T = ROL(SL[round][w], AL + F4(BL, CL, DL) + bufw[RL[round][w]] + KL[round]) + EL;
AL = EL; EL = DL; DL = ROL(10, CL); CL = BL; BL = T;
}
for (w = 0; w < 16; w++) { /* right line */
T = ROL(SR[round][w], AR + F2(BR, CR, DR) + bufw[RR[round][w]] + KR[round]) + ER;
AR = ER; ER = DR; DR = ROL(10, CR); CR = BR; BR = T;
}
/* Round 5 */
round++;
for (w = 0; w < 16; w++) { /* left line */
T = ROL(SL[round][w], AL + F5(BL, CL, DL) + bufw[RL[round][w]] + KL[round]) + EL;
AL = EL; EL = DL; DL = ROL(10, CL); CL = BL; BL = T;
}
for (w = 0; w < 16; w++) { /* right line */
T = ROL(SR[round][w], AR + F1(BR, CR, DR) + bufw[RR[round][w]] + KR[round]) + ER;
AR = ER; ER = DR; DR = ROL(10, CR); CR = BR; BR = T;
}
/* Final mixing stage */
T = self->h[1] + CL + DR;
self->h[1] = self->h[2] + DL + ER;
self->h[2] = self->h[3] + EL + AR;
self->h[3] = self->h[4] + AL + BR;
self->h[4] = self->h[0] + BL + CR;
self->h[0] = T;
/* Clear the buffer and wipe the temporary variables */
T = AL = BL = CL = DL = EL = AR = BR = CR = DR = ER = 0;
memset(&self->buf, 0, sizeof(self->buf));
self->bufpos = 0;
}
/*
* RIPEMD-160 Compression Function
*/
void RIPEMD_160::compress_n(const uint8_t input[], size_t blocks)
{
const uint32_t MAGIC2 = 0x5A827999, MAGIC3 = 0x6ED9EBA1,
MAGIC4 = 0x8F1BBCDC, MAGIC5 = 0xA953FD4E,
MAGIC6 = 0x50A28BE6, MAGIC7 = 0x5C4DD124,
MAGIC8 = 0x6D703EF3, MAGIC9 = 0x7A6D76E9;
for(size_t i = 0; i != blocks; ++i)
{
load_le(m_M.data(), input, m_M.size());
uint32_t A1 = m_digest[0], A2 = A1, B1 = m_digest[1], B2 = B1,
C1 = m_digest[2], C2 = C1, D1 = m_digest[3], D2 = D1,
E1 = m_digest[4], E2 = E1;
F1(A1,B1,C1,D1,E1,m_M[ 0],11 ); F5(A2,B2,C2,D2,E2,m_M[ 5], 8,MAGIC6);
F1(E1,A1,B1,C1,D1,m_M[ 1],14 ); F5(E2,A2,B2,C2,D2,m_M[14], 9,MAGIC6);
F1(D1,E1,A1,B1,C1,m_M[ 2],15 ); F5(D2,E2,A2,B2,C2,m_M[ 7], 9,MAGIC6);
F1(C1,D1,E1,A1,B1,m_M[ 3],12 ); F5(C2,D2,E2,A2,B2,m_M[ 0],11,MAGIC6);
F1(B1,C1,D1,E1,A1,m_M[ 4], 5 ); F5(B2,C2,D2,E2,A2,m_M[ 9],13,MAGIC6);
F1(A1,B1,C1,D1,E1,m_M[ 5], 8 ); F5(A2,B2,C2,D2,E2,m_M[ 2],15,MAGIC6);
F1(E1,A1,B1,C1,D1,m_M[ 6], 7 ); F5(E2,A2,B2,C2,D2,m_M[11],15,MAGIC6);
F1(D1,E1,A1,B1,C1,m_M[ 7], 9 ); F5(D2,E2,A2,B2,C2,m_M[ 4], 5,MAGIC6);
F1(C1,D1,E1,A1,B1,m_M[ 8],11 ); F5(C2,D2,E2,A2,B2,m_M[13], 7,MAGIC6);
F1(B1,C1,D1,E1,A1,m_M[ 9],13 ); F5(B2,C2,D2,E2,A2,m_M[ 6], 7,MAGIC6);
F1(A1,B1,C1,D1,E1,m_M[10],14 ); F5(A2,B2,C2,D2,E2,m_M[15], 8,MAGIC6);
F1(E1,A1,B1,C1,D1,m_M[11],15 ); F5(E2,A2,B2,C2,D2,m_M[ 8],11,MAGIC6);
F1(D1,E1,A1,B1,C1,m_M[12], 6 ); F5(D2,E2,A2,B2,C2,m_M[ 1],14,MAGIC6);
F1(C1,D1,E1,A1,B1,m_M[13], 7 ); F5(C2,D2,E2,A2,B2,m_M[10],14,MAGIC6);
F1(B1,C1,D1,E1,A1,m_M[14], 9 ); F5(B2,C2,D2,E2,A2,m_M[ 3],12,MAGIC6);
F1(A1,B1,C1,D1,E1,m_M[15], 8 ); F5(A2,B2,C2,D2,E2,m_M[12], 6,MAGIC6);
F2(E1,A1,B1,C1,D1,m_M[ 7], 7,MAGIC2); F4(E2,A2,B2,C2,D2,m_M[ 6], 9,MAGIC7);
F2(D1,E1,A1,B1,C1,m_M[ 4], 6,MAGIC2); F4(D2,E2,A2,B2,C2,m_M[11],13,MAGIC7);
F2(C1,D1,E1,A1,B1,m_M[13], 8,MAGIC2); F4(C2,D2,E2,A2,B2,m_M[ 3],15,MAGIC7);
F2(B1,C1,D1,E1,A1,m_M[ 1],13,MAGIC2); F4(B2,C2,D2,E2,A2,m_M[ 7], 7,MAGIC7);
F2(A1,B1,C1,D1,E1,m_M[10],11,MAGIC2); F4(A2,B2,C2,D2,E2,m_M[ 0],12,MAGIC7);
F2(E1,A1,B1,C1,D1,m_M[ 6], 9,MAGIC2); F4(E2,A2,B2,C2,D2,m_M[13], 8,MAGIC7);
F2(D1,E1,A1,B1,C1,m_M[15], 7,MAGIC2); F4(D2,E2,A2,B2,C2,m_M[ 5], 9,MAGIC7);
F2(C1,D1,E1,A1,B1,m_M[ 3],15,MAGIC2); F4(C2,D2,E2,A2,B2,m_M[10],11,MAGIC7);
F2(B1,C1,D1,E1,A1,m_M[12], 7,MAGIC2); F4(B2,C2,D2,E2,A2,m_M[14], 7,MAGIC7);
F2(A1,B1,C1,D1,E1,m_M[ 0],12,MAGIC2); F4(A2,B2,C2,D2,E2,m_M[15], 7,MAGIC7);
F2(E1,A1,B1,C1,D1,m_M[ 9],15,MAGIC2); F4(E2,A2,B2,C2,D2,m_M[ 8],12,MAGIC7);
F2(D1,E1,A1,B1,C1,m_M[ 5], 9,MAGIC2); F4(D2,E2,A2,B2,C2,m_M[12], 7,MAGIC7);
F2(C1,D1,E1,A1,B1,m_M[ 2],11,MAGIC2); F4(C2,D2,E2,A2,B2,m_M[ 4], 6,MAGIC7);
F2(B1,C1,D1,E1,A1,m_M[14], 7,MAGIC2); F4(B2,C2,D2,E2,A2,m_M[ 9],15,MAGIC7);
F2(A1,B1,C1,D1,E1,m_M[11],13,MAGIC2); F4(A2,B2,C2,D2,E2,m_M[ 1],13,MAGIC7);
F2(E1,A1,B1,C1,D1,m_M[ 8],12,MAGIC2); F4(E2,A2,B2,C2,D2,m_M[ 2],11,MAGIC7);
F3(D1,E1,A1,B1,C1,m_M[ 3],11,MAGIC3); F3(D2,E2,A2,B2,C2,m_M[15], 9,MAGIC8);
F3(C1,D1,E1,A1,B1,m_M[10],13,MAGIC3); F3(C2,D2,E2,A2,B2,m_M[ 5], 7,MAGIC8);
F3(B1,C1,D1,E1,A1,m_M[14], 6,MAGIC3); F3(B2,C2,D2,E2,A2,m_M[ 1],15,MAGIC8);
F3(A1,B1,C1,D1,E1,m_M[ 4], 7,MAGIC3); F3(A2,B2,C2,D2,E2,m_M[ 3],11,MAGIC8);
F3(E1,A1,B1,C1,D1,m_M[ 9],14,MAGIC3); F3(E2,A2,B2,C2,D2,m_M[ 7], 8,MAGIC8);
F3(D1,E1,A1,B1,C1,m_M[15], 9,MAGIC3); F3(D2,E2,A2,B2,C2,m_M[14], 6,MAGIC8);
F3(C1,D1,E1,A1,B1,m_M[ 8],13,MAGIC3); F3(C2,D2,E2,A2,B2,m_M[ 6], 6,MAGIC8);
F3(B1,C1,D1,E1,A1,m_M[ 1],15,MAGIC3); F3(B2,C2,D2,E2,A2,m_M[ 9],14,MAGIC8);
F3(A1,B1,C1,D1,E1,m_M[ 2],14,MAGIC3); F3(A2,B2,C2,D2,E2,m_M[11],12,MAGIC8);
F3(E1,A1,B1,C1,D1,m_M[ 7], 8,MAGIC3); F3(E2,A2,B2,C2,D2,m_M[ 8],13,MAGIC8);
F3(D1,E1,A1,B1,C1,m_M[ 0],13,MAGIC3); F3(D2,E2,A2,B2,C2,m_M[12], 5,MAGIC8);
F3(C1,D1,E1,A1,B1,m_M[ 6], 6,MAGIC3); F3(C2,D2,E2,A2,B2,m_M[ 2],14,MAGIC8);
F3(B1,C1,D1,E1,A1,m_M[13], 5,MAGIC3); F3(B2,C2,D2,E2,A2,m_M[10],13,MAGIC8);
F3(A1,B1,C1,D1,E1,m_M[11],12,MAGIC3); F3(A2,B2,C2,D2,E2,m_M[ 0],13,MAGIC8);
F3(E1,A1,B1,C1,D1,m_M[ 5], 7,MAGIC3); F3(E2,A2,B2,C2,D2,m_M[ 4], 7,MAGIC8);
F3(D1,E1,A1,B1,C1,m_M[12], 5,MAGIC3); F3(D2,E2,A2,B2,C2,m_M[13], 5,MAGIC8);
F4(C1,D1,E1,A1,B1,m_M[ 1],11,MAGIC4); F2(C2,D2,E2,A2,B2,m_M[ 8],15,MAGIC9);
F4(B1,C1,D1,E1,A1,m_M[ 9],12,MAGIC4); F2(B2,C2,D2,E2,A2,m_M[ 6], 5,MAGIC9);
F4(A1,B1,C1,D1,E1,m_M[11],14,MAGIC4); F2(A2,B2,C2,D2,E2,m_M[ 4], 8,MAGIC9);
F4(E1,A1,B1,C1,D1,m_M[10],15,MAGIC4); F2(E2,A2,B2,C2,D2,m_M[ 1],11,MAGIC9);
F4(D1,E1,A1,B1,C1,m_M[ 0],14,MAGIC4); F2(D2,E2,A2,B2,C2,m_M[ 3],14,MAGIC9);
F4(C1,D1,E1,A1,B1,m_M[ 8],15,MAGIC4); F2(C2,D2,E2,A2,B2,m_M[11],14,MAGIC9);
F4(B1,C1,D1,E1,A1,m_M[12], 9,MAGIC4); F2(B2,C2,D2,E2,A2,m_M[15], 6,MAGIC9);
F4(A1,B1,C1,D1,E1,m_M[ 4], 8,MAGIC4); F2(A2,B2,C2,D2,E2,m_M[ 0],14,MAGIC9);
F4(E1,A1,B1,C1,D1,m_M[13], 9,MAGIC4); F2(E2,A2,B2,C2,D2,m_M[ 5], 6,MAGIC9);
F4(D1,E1,A1,B1,C1,m_M[ 3],14,MAGIC4); F2(D2,E2,A2,B2,C2,m_M[12], 9,MAGIC9);
F4(C1,D1,E1,A1,B1,m_M[ 7], 5,MAGIC4); F2(C2,D2,E2,A2,B2,m_M[ 2],12,MAGIC9);
F4(B1,C1,D1,E1,A1,m_M[15], 6,MAGIC4); F2(B2,C2,D2,E2,A2,m_M[13], 9,MAGIC9);
F4(A1,B1,C1,D1,E1,m_M[14], 8,MAGIC4); F2(A2,B2,C2,D2,E2,m_M[ 9],12,MAGIC9);
F4(E1,A1,B1,C1,D1,m_M[ 5], 6,MAGIC4); F2(E2,A2,B2,C2,D2,m_M[ 7], 5,MAGIC9);
F4(D1,E1,A1,B1,C1,m_M[ 6], 5,MAGIC4); F2(D2,E2,A2,B2,C2,m_M[10],15,MAGIC9);
F4(C1,D1,E1,A1,B1,m_M[ 2],12,MAGIC4); F2(C2,D2,E2,A2,B2,m_M[14], 8,MAGIC9);
F5(B1,C1,D1,E1,A1,m_M[ 4], 9,MAGIC5); F1(B2,C2,D2,E2,A2,m_M[12], 8 );
F5(A1,B1,C1,D1,E1,m_M[ 0],15,MAGIC5); F1(A2,B2,C2,D2,E2,m_M[15], 5 );
F5(E1,A1,B1,C1,D1,m_M[ 5], 5,MAGIC5); F1(E2,A2,B2,C2,D2,m_M[10],12 );
F5(D1,E1,A1,B1,C1,m_M[ 9],11,MAGIC5); F1(D2,E2,A2,B2,C2,m_M[ 4], 9 );
F5(C1,D1,E1,A1,B1,m_M[ 7], 6,MAGIC5); F1(C2,D2,E2,A2,B2,m_M[ 1],12 );
F5(B1,C1,D1,E1,A1,m_M[12], 8,MAGIC5); F1(B2,C2,D2,E2,A2,m_M[ 5], 5 );
F5(A1,B1,C1,D1,E1,m_M[ 2],13,MAGIC5); F1(A2,B2,C2,D2,E2,m_M[ 8],14 );
F5(E1,A1,B1,C1,D1,m_M[10],12,MAGIC5); F1(E2,A2,B2,C2,D2,m_M[ 7], 6 );
F5(D1,E1,A1,B1,C1,m_M[14], 5,MAGIC5); F1(D2,E2,A2,B2,C2,m_M[ 6], 8 );
F5(C1,D1,E1,A1,B1,m_M[ 1],12,MAGIC5); F1(C2,D2,E2,A2,B2,m_M[ 2],13 );
F5(B1,C1,D1,E1,A1,m_M[ 3],13,MAGIC5); F1(B2,C2,D2,E2,A2,m_M[13], 6 );
F5(A1,B1,C1,D1,E1,m_M[ 8],14,MAGIC5); F1(A2,B2,C2,D2,E2,m_M[14], 5 );
F5(E1,A1,B1,C1,D1,m_M[11],11,MAGIC5); F1(E2,A2,B2,C2,D2,m_M[ 0],15 );
F5(D1,E1,A1,B1,C1,m_M[ 6], 8,MAGIC5); F1(D2,E2,A2,B2,C2,m_M[ 3],13 );
F5(C1,D1,E1,A1,B1,m_M[15], 5,MAGIC5); F1(C2,D2,E2,A2,B2,m_M[ 9],11 );
F5(B1,C1,D1,E1,A1,m_M[13], 6,MAGIC5); F1(B2,C2,D2,E2,A2,m_M[11],11 );
C1 = m_digest[1] + C1 + D2;
m_digest[1] = m_digest[2] + D1 + E2;
m_digest[2] = m_digest[3] + E1 + A2;
m_digest[3] = m_digest[4] + A1 + B2;
m_digest[4] = m_digest[0] + B1 + C2;
m_digest[0] = C1;
input += hash_block_size();
}
}
int main () {
//===============================================================//
//========= PROBLEM SETUP =======================//
//===============================================================//
int d = 4;
int j0 = 2;
cout << "Wavelets: d = " << d << ", j0 = " << j0 << endl << endl;
/// Basis initialization
_Basis basis(d, j0);
Basis2D basis2d(basis, basis);
/*
for (int i = 1; i<= 6; ++i) {
DataType u;
std::string snapshot;
std::string plot;
snapshot = "training_data_toy/snap/snap_";
plot = "plot_";
snapshot += std::to_string(i);
plot += std::to_string(i);
snapshot += ".txt";
std::cout << "Reading snapshot from file " << snapshot << std::endl;
std::cout << "Printing plot to file " << plot << std::endl;
readCoeffVector2D(u, snapshot.c_str());
precon _p;
plot2D(basis2d, u, _p,
zero,
0., 1.,
0., 1.,
1e-02,
plot.c_str());
}
exit(0);
*/
/// Initialization of operators
DenseVectorT no_singPts;
Function<T> zero_Fct(zero_fct,no_singPts);
// Bilinear Forms
Identity1D IdentityBil(basis);
Laplace1D LaplaceBil(basis);
RefIdentity1D RefIdentityBil(basis.refinementbasis);
RefLaplace1D RefLaplaceBil(basis.refinementbasis);
/// Initialization of local operator
LOp_Id1D lOp_Id1D (basis, basis, RefIdentityBil, IdentityBil);
LOp_Lapl1D lOp_Lapl1D (basis, basis, RefLaplaceBil, LaplaceBil);
LOp_Id_Id_2D localIdentityIdentityOp2D (lOp_Id1D, lOp_Id1D);
LOp_Id_Lapl_2D localIdentityLaplaceOp2D (lOp_Id1D, lOp_Lapl1D);
LOp_Lapl_Id_2D localLaplaceIdentityOp2D (lOp_Lapl1D, lOp_Id1D);
localIdentityIdentityOp2D.setJ(9);
localIdentityLaplaceOp2D.setJ(9);
localLaplaceIdentityOp2D.setJ(9);
/// Initialization of preconditioner
Prec2D prec(basis2d);
NoPrec2D noPrec;
/// Initialization of rhs
/// Right Hand Side:
/// No Singular Supports in both dimensions
DenseVectorT sing_support;
/// Forcing Functions
Function2D<T> F1_Fct(f1, sing_support, sing_support);
Function2D<T> F2_Fct(f2, sing_support, sing_support);
Function2D<T> F3_Fct(f3, sing_support, sing_support);
Function2D<T> F4_Fct(f4, sing_support, sing_support);
Function2D<T> F5_Fct(f5, sing_support, sing_support);
Function2D<T> F6_Fct(f6, sing_support, sing_support);
Function2D<T> F7_Fct(f7, sing_support, sing_support);
Function2D<T> F8_Fct(f8, sing_support, sing_support);
Function2D<T> F9_Fct(f9, sing_support, sing_support);
Function2D<T> F10_Fct(f10, sing_support, sing_support);
Function2D<T> F11_Fct(f11, sing_support, sing_support);
Function2D<T> F12_Fct(f12, sing_support, sing_support);
Function2D<T> F13_Fct(f13, sing_support, sing_support);
Function2D<T> F14_Fct(f14, sing_support, sing_support);
Function2D<T> F15_Fct(f15, sing_support, sing_support);
Function2D<T> F16_Fct(f16, sing_support, sing_support);
Function2D<T> F17_Fct(f17, sing_support, sing_support);
Function2D<T> F18_Fct(f18, sing_support, sing_support);
Function2D<T> F19_Fct(f19, sing_support, sing_support);
Function2D<T> F20_Fct(f20, sing_support, sing_support);
Function2D<T> F21_Fct(f21, sing_support, sing_support);
Function2D<T> F22_Fct(f22, sing_support, sing_support);
Function2D<T> F23_Fct(f23, sing_support, sing_support);
Function2D<T> F24_Fct(f24, sing_support, sing_support);
Function2D<T> F25_Fct(f25, sing_support, sing_support);
Function2D<T> F26_Fct(f26, sing_support, sing_support);
Function2D<T> F27_Fct(f27, sing_support, sing_support);
Function2D<T> F28_Fct(f28, sing_support, sing_support);
Function2D<T> F29_Fct(f29, sing_support, sing_support);
Function2D<T> F30_Fct(f30, sing_support, sing_support);
Function2D<T> F31_Fct(f31, sing_support, sing_support);
Function2D<T> F32_Fct(f32, sing_support, sing_support);
Function2D<T> F33_Fct(f33, sing_support, sing_support);
Function2D<T> F34_Fct(f34, sing_support, sing_support);
Function2D<T> F35_Fct(f35, sing_support, sing_support);
Function2D<T> F36_Fct(f36, sing_support, sing_support);
Function2D<T> F37_Fct(f37, sing_support, sing_support);
Function2D<T> F38_Fct(f38, sing_support, sing_support);
Function2D<T> F39_Fct(f39, sing_support, sing_support);
Function2D<T> F40_Fct(f40, sing_support, sing_support);
Function2D<T> F41_Fct(f41, sing_support, sing_support);
Function2D<T> F42_Fct(f42, sing_support, sing_support);
Function2D<T> F43_Fct(f43, sing_support, sing_support);
Function2D<T> F44_Fct(f44, sing_support, sing_support);
Function2D<T> F45_Fct(f45, sing_support, sing_support);
Function2D<T> F46_Fct(f46, sing_support, sing_support);
Function2D<T> F47_Fct(f47, sing_support, sing_support);
Function2D<T> F48_Fct(f48, sing_support, sing_support);
Function2D<T> F49_Fct(f49, sing_support, sing_support);
RhsIntegral2D rhs_1(basis2d, F1_Fct, 100);
RhsIntegral2D rhs_2(basis2d, F2_Fct, 100);
RhsIntegral2D rhs_3(basis2d, F3_Fct, 100);
RhsIntegral2D rhs_4(basis2d, F4_Fct, 100);
RhsIntegral2D rhs_5(basis2d, F5_Fct, 100);
RhsIntegral2D rhs_6(basis2d, F6_Fct, 100);
RhsIntegral2D rhs_7(basis2d, F7_Fct, 100);
RhsIntegral2D rhs_8(basis2d, F8_Fct, 100);
RhsIntegral2D rhs_9(basis2d, F9_Fct, 100);
RhsIntegral2D rhs_10(basis2d, F10_Fct, 100);
RhsIntegral2D rhs_11(basis2d, F11_Fct, 100);
RhsIntegral2D rhs_12(basis2d, F12_Fct, 100);
RhsIntegral2D rhs_13(basis2d, F13_Fct, 100);
RhsIntegral2D rhs_14(basis2d, F14_Fct, 100);
RhsIntegral2D rhs_15(basis2d, F15_Fct, 100);
RhsIntegral2D rhs_16(basis2d, F16_Fct, 100);
RhsIntegral2D rhs_17(basis2d, F17_Fct, 100);
RhsIntegral2D rhs_18(basis2d, F18_Fct, 100);
RhsIntegral2D rhs_19(basis2d, F19_Fct, 100);
RhsIntegral2D rhs_20(basis2d, F20_Fct, 100);
RhsIntegral2D rhs_21(basis2d, F21_Fct, 100);
RhsIntegral2D rhs_22(basis2d, F22_Fct, 100);
RhsIntegral2D rhs_23(basis2d, F23_Fct, 100);
RhsIntegral2D rhs_24(basis2d, F24_Fct, 100);
RhsIntegral2D rhs_25(basis2d, F25_Fct, 100);
RhsIntegral2D rhs_26(basis2d, F26_Fct, 100);
RhsIntegral2D rhs_27(basis2d, F27_Fct, 100);
RhsIntegral2D rhs_28(basis2d, F28_Fct, 100);
RhsIntegral2D rhs_29(basis2d, F29_Fct, 100);
RhsIntegral2D rhs_30(basis2d, F30_Fct, 100);
RhsIntegral2D rhs_31(basis2d, F31_Fct, 100);
RhsIntegral2D rhs_32(basis2d, F32_Fct, 100);
RhsIntegral2D rhs_33(basis2d, F33_Fct, 100);
RhsIntegral2D rhs_34(basis2d, F34_Fct, 100);
RhsIntegral2D rhs_35(basis2d, F35_Fct, 100);
RhsIntegral2D rhs_36(basis2d, F36_Fct, 100);
RhsIntegral2D rhs_37(basis2d, F37_Fct, 100);
RhsIntegral2D rhs_38(basis2d, F38_Fct, 100);
RhsIntegral2D rhs_39(basis2d, F39_Fct, 100);
RhsIntegral2D rhs_40(basis2d, F40_Fct, 100);
RhsIntegral2D rhs_41(basis2d, F41_Fct, 100);
RhsIntegral2D rhs_42(basis2d, F42_Fct, 100);
RhsIntegral2D rhs_43(basis2d, F43_Fct, 100);
RhsIntegral2D rhs_44(basis2d, F44_Fct, 100);
RhsIntegral2D rhs_45(basis2d, F45_Fct, 100);
RhsIntegral2D rhs_46(basis2d, F46_Fct, 100);
RhsIntegral2D rhs_47(basis2d, F47_Fct, 100);
RhsIntegral2D rhs_48(basis2d, F48_Fct, 100);
RhsIntegral2D rhs_49(basis2d, F49_Fct, 100);
Rhs F1(rhs_1,noPrec);
Rhs F2(rhs_2,noPrec);
Rhs F3(rhs_3,noPrec);
Rhs F4(rhs_4,noPrec);
Rhs F5(rhs_5,noPrec);
Rhs F6(rhs_6,noPrec);
Rhs F7(rhs_7,noPrec);
Rhs F8(rhs_8,noPrec);
Rhs F9(rhs_9,noPrec);
Rhs F10(rhs_10,noPrec);
Rhs F11(rhs_11,noPrec);
Rhs F12(rhs_12,noPrec);
Rhs F13(rhs_13,noPrec);
Rhs F14(rhs_14,noPrec);
Rhs F15(rhs_15,noPrec);
Rhs F16(rhs_16,noPrec);
Rhs F17(rhs_17,noPrec);
Rhs F18(rhs_18,noPrec);
Rhs F19(rhs_19,noPrec);
Rhs F20(rhs_20,noPrec);
Rhs F21(rhs_21,noPrec);
Rhs F22(rhs_22,noPrec);
Rhs F23(rhs_23,noPrec);
Rhs F24(rhs_24,noPrec);
Rhs F25(rhs_25,noPrec);
Rhs F26(rhs_26,noPrec);
Rhs F27(rhs_27,noPrec);
Rhs F28(rhs_28,noPrec);
Rhs F29(rhs_29,noPrec);
Rhs F30(rhs_30,noPrec);
Rhs F31(rhs_31,noPrec);
Rhs F32(rhs_32,noPrec);
Rhs F33(rhs_33,noPrec);
Rhs F34(rhs_34,noPrec);
Rhs F35(rhs_35,noPrec);
Rhs F36(rhs_36,noPrec);
Rhs F37(rhs_37,noPrec);
Rhs F38(rhs_38,noPrec);
Rhs F39(rhs_39,noPrec);
Rhs F40(rhs_40,noPrec);
Rhs F41(rhs_41,noPrec);
Rhs F42(rhs_42,noPrec);
Rhs F43(rhs_43,noPrec);
Rhs F44(rhs_44,noPrec);
Rhs F45(rhs_45,noPrec);
Rhs F46(rhs_46,noPrec);
Rhs F47(rhs_47,noPrec);
Rhs F48(rhs_48,noPrec);
Rhs F49(rhs_49,noPrec);
//===============================================================//
//=============== RB SETUP =====================================//
//===============================================================//
// Affine Decompositions:
// Left Hand Side
vector<ThetaStructure<ParamType>::ThetaFct> lhs_theta_fcts;
lhs_theta_fcts.push_back(no_theta);
ThetaStructure<ParamType> lhs_theta(lhs_theta_fcts);
vector<AbstractLocalOperator2D<T>* > lhs_ops;
lhs_ops.push_back(&localLaplaceIdentityOp2D);
lhs_ops.push_back(&localIdentityLaplaceOp2D);
Flex_COp_2D A(lhs_ops);
vector<Flex_COp_2D*> ops_vec;
ops_vec.push_back(&A);
Affine_Op_2D affine_lhs(lhs_theta, ops_vec);
H1_InnProd_2D innprod(localIdentityIdentityOp2D, localLaplaceIdentityOp2D, localIdentityLaplaceOp2D);
// Right Hand Side
vector<ThetaStructure<ParamType>::ThetaFct> rhs_theta_fcts;
rhs_theta_fcts.push_back(theta_chi_1);
rhs_theta_fcts.push_back(theta_chi_2);
rhs_theta_fcts.push_back(theta_chi_3);
rhs_theta_fcts.push_back(theta_chi_4);
rhs_theta_fcts.push_back(theta_chi_5);
rhs_theta_fcts.push_back(theta_chi_6);
rhs_theta_fcts.push_back(theta_chi_7);
rhs_theta_fcts.push_back(theta_chi_8);
rhs_theta_fcts.push_back(theta_chi_9);
rhs_theta_fcts.push_back(theta_chi_10);
rhs_theta_fcts.push_back(theta_chi_11);
rhs_theta_fcts.push_back(theta_chi_12);
rhs_theta_fcts.push_back(theta_chi_13);
rhs_theta_fcts.push_back(theta_chi_14);
rhs_theta_fcts.push_back(theta_chi_15);
rhs_theta_fcts.push_back(theta_chi_16);
rhs_theta_fcts.push_back(theta_chi_17);
rhs_theta_fcts.push_back(theta_chi_18);
rhs_theta_fcts.push_back(theta_chi_19);
rhs_theta_fcts.push_back(theta_chi_20);
rhs_theta_fcts.push_back(theta_chi_21);
rhs_theta_fcts.push_back(theta_chi_22);
rhs_theta_fcts.push_back(theta_chi_23);
rhs_theta_fcts.push_back(theta_chi_24);
rhs_theta_fcts.push_back(theta_chi_25);
rhs_theta_fcts.push_back(theta_chi_26);
rhs_theta_fcts.push_back(theta_chi_27);
rhs_theta_fcts.push_back(theta_chi_28);
rhs_theta_fcts.push_back(theta_chi_29);
rhs_theta_fcts.push_back(theta_chi_30);
rhs_theta_fcts.push_back(theta_chi_31);
rhs_theta_fcts.push_back(theta_chi_32);
rhs_theta_fcts.push_back(theta_chi_33);
rhs_theta_fcts.push_back(theta_chi_34);
rhs_theta_fcts.push_back(theta_chi_35);
rhs_theta_fcts.push_back(theta_chi_36);
rhs_theta_fcts.push_back(theta_chi_37);
rhs_theta_fcts.push_back(theta_chi_38);
rhs_theta_fcts.push_back(theta_chi_39);
rhs_theta_fcts.push_back(theta_chi_40);
rhs_theta_fcts.push_back(theta_chi_41);
rhs_theta_fcts.push_back(theta_chi_42);
rhs_theta_fcts.push_back(theta_chi_43);
rhs_theta_fcts.push_back(theta_chi_44);
rhs_theta_fcts.push_back(theta_chi_45);
rhs_theta_fcts.push_back(theta_chi_46);
rhs_theta_fcts.push_back(theta_chi_47);
rhs_theta_fcts.push_back(theta_chi_48);
rhs_theta_fcts.push_back(theta_chi_49);
ThetaStructure<ParamType> rhs_theta(rhs_theta_fcts);
vector<Rhs*> rhs_fcts;
rhs_fcts.push_back(&F1);
rhs_fcts.push_back(&F2);
rhs_fcts.push_back(&F3);
rhs_fcts.push_back(&F4);
rhs_fcts.push_back(&F5);
rhs_fcts.push_back(&F6);
rhs_fcts.push_back(&F7);
rhs_fcts.push_back(&F8);
rhs_fcts.push_back(&F9);
rhs_fcts.push_back(&F10);
rhs_fcts.push_back(&F11);
rhs_fcts.push_back(&F12);
rhs_fcts.push_back(&F13);
rhs_fcts.push_back(&F14);
rhs_fcts.push_back(&F15);
rhs_fcts.push_back(&F16);
rhs_fcts.push_back(&F17);
rhs_fcts.push_back(&F18);
rhs_fcts.push_back(&F19);
rhs_fcts.push_back(&F20);
rhs_fcts.push_back(&F21);
rhs_fcts.push_back(&F22);
rhs_fcts.push_back(&F23);
rhs_fcts.push_back(&F24);
rhs_fcts.push_back(&F25);
rhs_fcts.push_back(&F26);
rhs_fcts.push_back(&F27);
rhs_fcts.push_back(&F28);
rhs_fcts.push_back(&F29);
rhs_fcts.push_back(&F30);
rhs_fcts.push_back(&F31);
rhs_fcts.push_back(&F32);
rhs_fcts.push_back(&F33);
rhs_fcts.push_back(&F34);
rhs_fcts.push_back(&F35);
rhs_fcts.push_back(&F36);
rhs_fcts.push_back(&F37);
rhs_fcts.push_back(&F38);
rhs_fcts.push_back(&F39);
rhs_fcts.push_back(&F40);
rhs_fcts.push_back(&F41);
rhs_fcts.push_back(&F42);
rhs_fcts.push_back(&F43);
rhs_fcts.push_back(&F44);
rhs_fcts.push_back(&F45);
rhs_fcts.push_back(&F46);
rhs_fcts.push_back(&F47);
rhs_fcts.push_back(&F48);
rhs_fcts.push_back(&F49);
Affine_Rhs_2D affine_rhs(rhs_theta, rhs_fcts);
RieszF_Rhs_2D rieszF_rhs(rhs_fcts);
RieszA_Rhs_2D rieszA_rhs(ops_vec);
// Right Hand Sides for direct Riesz Representor
AffineA_Rhs_2D affineA_rhs(lhs_theta, lhs_ops);
RieszRes_Rhs_2D rieszRes_rhs(affineA_rhs, affine_rhs);
//===============================================================//
//=============== AWGM =========================================//
//===============================================================//
IS_Parameters is_parameters;
AWGM_Parameters awgm_truth_parameters, awgm_riesz_f_parameters, awgm_riesz_a_parameters, awgm_riesz_res_parameters;
awgm_truth_parameters.max_its = 1e+03;
is_parameters.adaptive_tol = true;
is_parameters.init_tol = 0.0001;
is_parameters.res_reduction = 0.01;
//----------- Solver ---------------- //
T gamma = 0.2;
IndexSet<Index2D> Lambda;
getSparseGridIndexSet(basis2d, Lambda, 1, 0, gamma);
MT_AWGM_Truth awgm_u(basis2d, affine_lhs, affine_rhs, prec, awgm_truth_parameters, is_parameters);
awgm_u.set_sol(dummy);
awgm_u.awgm_params.tol = 1e-04;
awgm_u.awgm_params.max_basissize = 1e+06;
awgm_u.set_initial_indexset(Lambda);
MT_AWGM_Riesz_F awgm_rieszF(basis2d, innprod, rieszF_rhs, prec, awgm_riesz_f_parameters, is_parameters);
awgm_rieszF.set_sol(dummy);
awgm_rieszF.set_initial_indexset(Lambda);
awgm_rieszF.awgm_params.tol = 5e-04;
awgm_rieszF.awgm_params.info_filename = "stempel_rieszF_conv_info.txt";
MT_AWGM_Riesz_A awgm_rieszA(basis2d, innprod, rieszA_rhs, prec, awgm_riesz_a_parameters, is_parameters);
awgm_rieszA.set_sol(dummy);
awgm_rieszA.set_initial_indexset(Lambda);
awgm_rieszA.awgm_params.tol = 5e-04;
awgm_rieszA.awgm_params.info_filename = "stempel_rieszA_conv_info.txt";
MT_AWGM_Riesz_Res awgm_rieszRes(basis2d, innprod, rieszRes_rhs, prec, awgm_riesz_res_parameters, is_parameters);
awgm_rieszRes.set_sol(dummy);
awgm_rieszRes.set_initial_indexset(Lambda);
awgm_rieszRes.awgm_params.tol = 5e-04;
awgm_rieszRes.awgm_params.print_info = false;
MTTruthSolver rb_truth(awgm_u, awgm_rieszF, awgm_rieszA, &awgm_rieszRes, innprod, affine_lhs, rieszF_rhs);
//----------- RB System ---------------- //
RB_Model rb_system(lhs_theta, rhs_theta);
rb_system.rb_params.ref_param = {{1.}};
rb_system.rb_params.call = call_cg;
//----------- RB Base ---------------- //
RB_BaseModel rb_base(rb_system, rb_truth, Lambda);
ParamType mu_min = {{0}};
ParamType mu_max = {{1}};
rb_base.greedy_params.training_type = weak;
rb_base.greedy_params.tol = 1e-05;
rb_base.greedy_params.min_param = mu_min;
rb_base.greedy_params.max_param = mu_max;
rb_base.greedy_params.Nmax = 6;
rb_base.greedy_params.nb_training_params = {{nb_stempel}};
rb_base.greedy_params.log_scaling = {{0}};
rb_base.greedy_params.print_file = "toy_greedy_info.txt";
rb_base.greedy_params.trainingdata_folder = "training_data_toy";
rb_base.greedy_params.print_paramset = true;
rb_base.greedy_params.erase_snapshot_params = false;
rb_base.greedy_params.orthonormalize_bfs = true;
rb_base.greedy_params.tighten_tol = true;
rb_base.greedy_params.snapshot_tol_red_crit = repeated_param;
rb_base.greedy_params.tighten_tol_rieszA = false;
rb_base.greedy_params.tighten_tol_rieszF = false;
rb_base.greedy_params.tighten_tol_reduction = 0.1,
rb_base.greedy_params.update_snapshot = true;
rb_base.greedy_params.update_rieszF = false;
rb_base.greedy_params.update_rieszA = false;
rb_base.greedy_params.test_estimator_equivalence = false;
rb_base.greedy_params.write_direct_representors = false;
rb_base.greedy_params.riesz_constant_X = 5.6;
rb_base.greedy_params.riesz_constant_Y = 2;
cout << "Parameters Training: " << std::endl << std::endl;
rb_base.greedy_params.print();
rb_system.rb_params.print();
cout << "Parameters Truth Solver: " << std::endl << std::endl;
awgm_u.awgm_params.print();
awgm_u.is_params.print();
cout << "Parameters Riesz Solver F : " << std::endl << std::endl;
awgm_rieszF.awgm_params.print();
cout << "Parameters Riesz Solver A : " << std::endl << std::endl;
awgm_rieszA.awgm_params.print();
cout << "Parameters Riesz Solver Res : " << std::endl << std::endl;
awgm_rieszRes.awgm_params.print();
//rb_system.read_rb_data("training_data_toy");
//rb_base.read_basisfunctions("training_data_toy/bf");
/*for (int i=0; i<rb_base.n_bf(); ++i) {
std::cout << " i = " << i+1 << std::endl;
std::cout << rb_base.rb_basisfunctions[i].size() << std::endl;
}*/
//rb_base.read_rieszrepresentors("training_data_toy/representors");
//rb_base.calculate_Riesz_RHS_information(true);
//rb_base.calc_rb_data();
//rb_system.write_rb_data("training_data_toy");
rb_base.train_Greedy();
std::cout << "toy_offline exited normally\n";
return 0;
}
/* This is a simple example of multi-dimensional integration
using a simple (not necessarily optimal) spacing of points.
Note that this doesn't perform any error estimation - it
only calculates the value for a given grid size.
*/
double IntegrateExample_dumb(
int functionCode,
int n, // How many points on each dimension
const float *a, // An array of k lower bounds
const float *b, // An array of k upper bounds
const float *params // Parameters to function
){
int k=-1, total=-1, i0, i1, i2, j;
// Accumulate in double, as it avoids floating-point errors when adding large
// numbers of values together. Note that using double in a GPU has implications,
// as some GPUs cannot do doubles, and on others they are much slower than floats
double acc=0;
//float *x=NULL;
int n0=n, n1=n, n2=n; // By default use n points in each dimension
switch(functionCode){
case 0: k=1; break;
case 1: k=2; break;
case 2: k=3; break;
case 3: k=3; break;
case 4: k=3; break;
case 5: k=3; break;
case 6: k=3; break;
default:
fprintf(stderr, "Invalid function code.");
exit(1);
}
// Collapse any dimensions we don't use
if(k<3){
n2=1;
}
if(k<2){
n1=1;
}
float* x = new float[k];
// Loop over highest dimension on outside, as it might be collapsed to zero
for(i2=0;i2<n2;i2++){
if(k>2){
x[2]=a[2]+(b[2]-a[2]) * (i2+0.5f)/n2;
}
for(i1=0;i1<n1;i1++){
if(k>1){
x[1]=a[1]+(b[1]-a[1]) * (i1+0.5f)/n1;
}
// Inner dimension is never collapsed to zero
for(i0=0;i0<n0;i0++){
x[0]=a[0]+(b[0]-a[0]) * (i0+0.5f)/n0;
// Now call the function. Note that it is rather
// inefficient to be choosing the function in the inner loop...
// MAKE A FUNCTION POINTER FURTHER UP THEN USE IT!!!!
switch(functionCode){
case 0: acc+=F0(x,params); break;
case 1: acc+=F1(x,params); break;
case 2: acc+=F2(x,params); break;
case 3: acc+=F3(x,params); break;
case 4: acc+=F4(x,params); break;
case 5: acc+=F5(x,params); break;
case 6: acc+=F6(x,params); break;
}
}
}
}
// Do the final normalisation and return the results
for(j=0;j<k;j++){
acc=acc*(b[j]-a[j]);
}
return acc/(n0*n1*n2);
}
