static void MrEdSchemeMessages(char *msg, ...)
{
GC_CAN_IGNORE va_list args;
scheme_start_atomic();
XFORM_HIDE_EXPR(va_start(args, msg));
if (!console_out) {
AllocConsole();
console_out = GetStdHandle(STD_OUTPUT_HANDLE);
if (!wx_in_terminal) {
has_stdio = 1;
waiting_sema = CreateSemaphore(NULL, 0, 1, NULL);
orig_signal_handle = scheme_get_signal_handle();
orig_break_handle = scheme_get_main_thread_break_handle();
SetConsoleCtrlHandler(ConsoleHandler, TRUE);
{
HMODULE hm;
gcw_proc gcw;
hm = LoadLibrary("kernel32.dll");
if (hm)
gcw = (gcw_proc)GetProcAddress(hm, "GetConsoleWindow");
else
gcw = NULL;
if (gcw)
console_hwnd = gcw();
}
if (console_hwnd) {
EnableMenuItem(GetSystemMenu(console_hwnd, FALSE), SC_CLOSE,
MF_BYCOMMAND | MF_GRAYED);
RemoveMenu(GetSystemMenu(console_hwnd, FALSE), SC_CLOSE, MF_BYCOMMAND);
}
}
}
if (!msg) {
char *s;
intptr_t l, d;
DWORD wrote;
s = va_arg(args, char*);
d = va_arg(args, intptr_t);
l = va_arg(args, intptr_t);
WriteConsole(console_out, s XFORM_OK_PLUS d, l, &wrote, NULL);
} else {
void
AdvDiffPhysicalBoundaryUtilities::setPhysicalBoundaryConditions(
Pointer<CellData<NDIM,double> > Q_data,
Pointer<FaceData<NDIM,double> > u_ADV_data,
Pointer<Patch<NDIM> > patch,
const std::vector<RobinBcCoefStrategy<NDIM>*>& bc_coefs,
const double fill_time,
const bool inflow_boundaries_only,
const bool homogeneous_bc)
{
Pointer<CartesianPatchGeometry<NDIM> > pgeom = patch->getPatchGeometry();
if (!pgeom->getTouchesRegularBoundary()) return;
const Array<BoundaryBox<NDIM> > physical_codim1_boxes = PhysicalBoundaryUtilities::getPhysicalBoundaryCodim1Boxes(*patch);
if (physical_codim1_boxes.size() == 0) return;
// Loop over the boundary fill boxes and set boundary conditions.
const Box<NDIM>& patch_box = patch->getBox();
const double* const dx = pgeom->getDx();
// Setup any extended Robin BC coef objects.
for (int depth = 0; depth < Q_data->getDepth(); ++depth)
{
ExtendedRobinBcCoefStrategy* extended_bc_coef = dynamic_cast<ExtendedRobinBcCoefStrategy*>(bc_coefs[depth]);
if (extended_bc_coef)
{
extended_bc_coef->clearTargetPatchDataIndex();
extended_bc_coef->setHomogeneousBc(homogeneous_bc);
}
}
// Set the boundary conditions.
const IntVector<NDIM>& gcw = Q_data->getGhostCellWidth();
for (int n = 0; n < physical_codim1_boxes.size(); ++n)
{
const BoundaryBox<NDIM>& bdry_box = physical_codim1_boxes[n];
const unsigned int location_index = bdry_box.getLocationIndex();
const unsigned int bdry_normal_axis = location_index/2;
const bool is_lower = location_index%2 == 0;
static const IntVector<NDIM> gcw_to_fill = 1;
const Box<NDIM> bc_fill_box = pgeom->getBoundaryFillBox(bdry_box, patch_box, gcw_to_fill);
const BoundaryBox<NDIM> trimmed_bdry_box(bdry_box.getBox() * bc_fill_box, bdry_box.getBoundaryType(), bdry_box.getLocationIndex());
Box<NDIM> bc_coef_box = PhysicalBoundaryUtilities::makeSideBoundaryCodim1Box(trimmed_bdry_box);
for (unsigned int d = 0; d < NDIM; ++d)
{
if (d != bdry_normal_axis)
{
bc_coef_box.lower(d) = std::max(bc_coef_box.lower(d), patch_box.lower(d));
bc_coef_box.upper(d) = std::min(bc_coef_box.upper(d), patch_box.upper(d));
}
}
Pointer<ArrayData<NDIM,double> > acoef_data = new ArrayData<NDIM,double>(bc_coef_box, 1);
Pointer<ArrayData<NDIM,double> > bcoef_data = new ArrayData<NDIM,double>(bc_coef_box, 1);
Pointer<ArrayData<NDIM,double> > gcoef_data = new ArrayData<NDIM,double>(bc_coef_box, 1);
for (int depth = 0; depth < Q_data->getDepth(); ++depth)
{
bc_coefs[depth]->setBcCoefs(acoef_data, bcoef_data, gcoef_data, NULL, *patch, trimmed_bdry_box, fill_time);
ExtendedRobinBcCoefStrategy* extended_bc_coef = dynamic_cast<ExtendedRobinBcCoefStrategy*>(bc_coefs[depth]);
if (homogeneous_bc && !extended_bc_coef) gcoef_data->fillAll(0.0);
for (Box<NDIM>::Iterator bc(bc_coef_box); bc; bc++)
{
const Index<NDIM>& i = bc();
const FaceIndex<NDIM> i_f(i, bdry_normal_axis, FaceIndex<NDIM>::Lower);
bool is_inflow_bdry = (is_lower && (*u_ADV_data)(i_f) > 0.0) || (!is_lower && (*u_ADV_data)(i_f) < 0.0);
if (!inflow_boundaries_only || is_inflow_bdry)
{
const double& a = (*acoef_data)(i,0);
const double& b = (*bcoef_data)(i,0);
const double& g = (*gcoef_data)(i,0);
const double& h = dx[bdry_normal_axis];
int sgn;
Index<NDIM> i_intr(i);
if (is_lower)
{
sgn = -1;
}
else
{
sgn = +1;
i_intr(bdry_normal_axis) -= 1;
}
Index<NDIM> i_true(i_intr), i_ghost(i_intr);
for (int k = 1; k <= gcw(bdry_normal_axis); ++k)
{
i_ghost(bdry_normal_axis) = i_intr(bdry_normal_axis) + sgn*k;
i_true (bdry_normal_axis) = i_intr(bdry_normal_axis) - sgn*(k-1);
const double Q_i = (*Q_data)(i_true,depth);
(*Q_data)(i_ghost,depth) = -(-4.0*g*h*k-2.0*g*h+2.0*a*Q_i*h*k+a*Q_i*h-2.0*b*Q_i)/(2.0*a*h*k+a*h+2.0*b);
}
}
}
}
}
return;
}// setPhysicalBoundaryConditions
