int APIENTRY WinMain (
HINSTANCE hInstance, HINSTANCE hPrevInstance, LPSTR lpCmdLine, int nCmdShow
) {
Tobj co;
Psrc_t src;
hinstance = hInstance;
hprevinstance = hPrevInstance;
idlerunmode = 0;
exprstr = NULL;
fp = NULL;
init (NULL);
Ginit ();
#ifndef FEATURE_MS
FD_SET (Gxfd, &inputfds);
#endif
Eerrlevel = 1;
Estackdepth = 2;
Eshowbody = 1;
Eshowcalls = 1;
processstr (leftyoptions);
__argv++, __argc--;
processargs (__argc, __argv);
if (setjmp (exitljbuf))
goto eop;
Cinit ();
IOinit ();
Minit (GFXprune);
Tinit ();
Pinit ();
Einit ();
Sinit ();
Dinit ();
Iinit ();
TXTinit (txtcoords);
GFXinit ();
if (exprstr) {
src.flag = CHARSRC, src.s = exprstr, src.fp = NULL;
src.tok = -1, src.lnum = 1;
while ((co = Punit (&src)))
Eunit (co);
}
if (fp) {
src.flag = FILESRC, src.s = NULL, src.fp = fp;
src.tok = -1, src.lnum = 1;
while ((co = Punit (&src)))
Eunit (co);
}
if (endflag)
goto eop;
TXTupdate ();
Gneedredraw = FALSE;
for (;;) {
if (Gneedredraw)
GFXredraw (), Gneedredraw = FALSE;
if (Gbuttonsdown > 0)
GFXmove (), Gprocessevents (FALSE, G_ONEEVENT);
else {
if (Mcouldgc) {
if (!processinput (FALSE))
Mdogc (M_GCINCR);
} else if (idlerunmode) {
if (!processinput (FALSE))
GFXidle ();
} else
processinput (TRUE);
}
if (Erun)
TXTupdate (), Erun = FALSE;
}
eop:
#ifdef PARANOID
GFXterm ();
TXTterm ();
Iterm ();
Dterm ();
Sterm ();
Eterm ();
Pterm ();
Tterm ();
Mterm ();
IOterm ();
Cterm ();
Gterm ();
term ();
#endif
printusage ();
exit (0);
}
int main (int argc, char **argv) {
Tobj co;
Psrc_t src;
#ifdef MTRACE
extern int Mt_certify;
Mt_certify = 1;
#endif
#ifdef STATS
stime = time (NULL);
#endif
idlerunmode = 0;
exprstr = NULL;
fp = NULL;
init (argv[0]);
Minit (GFXprune);
Ginit ();
FD_SET (Gxfd, &inputfds);
Eerrlevel = 1;
Estackdepth = 2;
Eshowbody = 1;
Eshowcalls = 1;
processstr (leftyoptions);
argv++, argc--;
processargs (argc, argv);
if (setjmp (exitljbuf))
goto eop;
Cinit ();
IOinit ();
Tinit ();
Pinit ();
Einit ();
Sinit ();
Dinit ();
Iinit ();
TXTinit (txtcoords);
GFXinit ();
#ifdef FEATURE_GMAP
gmapon = TRUE, G2Linit ();
#endif
if (exprstr) {
src.flag = CHARSRC, src.s = exprstr, src.fp = NULL;
src.tok = -1, src.lnum = 1;
while ((co = Punit (&src)))
Eunit (co);
}
if (fp) {
src.flag = FILESRC, src.s = NULL, src.fp = fp;
src.tok = -1, src.lnum = 1;
while ((co = Punit (&src)))
Eunit (co);
}
if (endflag)
goto eop;
TXTupdate ();
Gneedredraw = FALSE;
for (;;) {
if (Gneedredraw)
GFXredraw (), Gneedredraw = FALSE;
if (Gbuttonsdown > 0) {
GFXmove (), Gprocessevents (FALSE, G_ONEEVENT);
processinput (FALSE);
} else {
if (Mcouldgc) {
if (!processinput (FALSE))
Mdogc (M_GCINCR);
}
if (idlerunmode) {
if (!processinput (FALSE))
GFXidle ();
#ifdef FEATURE_GMAP
} else if (GMAPneedupdate) {
processinput (FALSE);
#endif
} else
processinput (TRUE);
}
#ifdef FEATURE_GMAP
if (gmapon)
GMAPupdate ();
#endif
if (Erun)
TXTupdate (), Erun = FALSE;
}
eop:
#ifdef PARANOID
#ifdef FEATURE_GMAP
if (gmapon)
G2Lterm ();
#endif
GFXterm ();
TXTterm ();
Iterm ();
Dterm ();
Sterm ();
Eterm ();
Pterm ();
Tterm ();
IOterm ();
Cterm ();
Gterm ();
Mterm ();
FD_CLR (Gxfd, &inputfds);
term ();
#endif
printusage ();
return 0;
}
void thermicSolver::assemble(linearSystem<double> *lsys)
{
if(pAssembler) delete pAssembler;
pAssembler = new dofManager<double>(lsys);
// we first do all fixations. the behavior of the dofManager is to
// give priority to fixations : when a dof is fixed, it cannot be
// numbered afterwards
// Dirichlet conditions
for(std::size_t i = 0; i < allDirichlet.size(); i++) {
FilterDofTrivial filter;
FixNodalDofs(*LagSpace, allDirichlet[i].g->begin(),
allDirichlet[i].g->end(), *pAssembler, *allDirichlet[i]._f,
filter);
}
// LagrangeMultipliers
for(std::size_t i = 0; i < LagrangeMultiplierFields.size(); ++i) {
NumberDofs(*LagrangeMultiplierSpace, LagrangeMultiplierFields[i].g->begin(),
LagrangeMultiplierFields[i].g->end(), *pAssembler);
}
// Thermic Fields
for(std::size_t i = 0; i < thermicFields.size(); ++i) {
NumberDofs(*LagSpace, thermicFields[i].g->begin(),
thermicFields[i].g->end(), *pAssembler);
}
// Neumann conditions
GaussQuadrature Integ_Boundary(GaussQuadrature::Val);
for(std::size_t i = 0; i < allNeumann.size(); i++) {
std::cout << "Neumann BC" << std::endl;
LoadTerm<double> Lterm(*LagSpace, allNeumann[i]._f);
Assemble(Lterm, *LagSpace, allNeumann[i].g->begin(), allNeumann[i].g->end(),
Integ_Boundary, *pAssembler);
}
// Assemble cross term, laplace term and rhs term for LM
GaussQuadrature Integ_LagrangeMult(GaussQuadrature::ValVal);
GaussQuadrature Integ_Laplace(GaussQuadrature::GradGrad);
for(std::size_t i = 0; i < LagrangeMultiplierFields.size(); i++) {
printf("Lagrange Mult Lag\n");
LagrangeMultiplierTerm<double> LagTerm(*LagSpace, *LagrangeMultiplierSpace,
1.);
Assemble(LagTerm, *LagSpace, *LagrangeMultiplierSpace,
LagrangeMultiplierFields[i].g->begin(),
LagrangeMultiplierFields[i].g->end(), Integ_LagrangeMult,
*pAssembler);
printf("Lagrange Mult Lap\n");
LaplaceTerm<double, double> LapTerm(*LagrangeMultiplierSpace,
-LagrangeMultiplierFields[i]._tau);
Assemble(LapTerm, *LagrangeMultiplierSpace,
LagrangeMultiplierFields[i].g->begin(),
LagrangeMultiplierFields[i].g->end(), Integ_Laplace, *pAssembler);
printf("Lagrange Mult Load\n");
LoadTermOnBorder<double> Lterm(*LagrangeMultiplierSpace,
LagrangeMultiplierFields[i]._f);
Assemble(Lterm, *LagrangeMultiplierSpace,
LagrangeMultiplierFields[i].g->begin(),
LagrangeMultiplierFields[i].g->end(), Integ_Boundary, *pAssembler);
}
// Assemble thermic term
GaussQuadrature Integ_Bulk(GaussQuadrature::ValVal);
for(std::size_t i = 0; i < thermicFields.size(); i++) {
printf("Thermic Term\n");
LaplaceTerm<double, double> Tterm(*LagSpace, thermicFields[i]._k);
Assemble(Tterm, *LagSpace, thermicFields[i].g->begin(),
thermicFields[i].g->end(), Integ_Bulk, *pAssembler);
}
/*for (int i = 0;i<pAssembler->sizeOfR();i++){
for (int j = 0;j<pAssembler->sizeOfR();j++){
double d; lsys->getFromMatrix(i, j, d);
printf("%g ", d);
}
double d; lsys->getFromRightHandSide(i, d);
printf(" | %g\n", d);
}*/
printf("nDofs=%d\n", pAssembler->sizeOfR());
printf("nFixed=%d\n", pAssembler->sizeOfF());
}