static int
Qs(q_M_gc)(lua_State *L)
{
Qs(mLatColMat) *b = Qs(qlua_checkLatColMat)(L, 1, NULL, -1);
Qx(QDP_D,_destroy_M)(b->ptr);
b->ptr = 0;
return 0;
}
static int
Qs(q_M_fmt)(lua_State *L)
{
char fmt[72];
Qs(mLatColMat) *b = Qs(qlua_checkLatColMat)(L, 1, NULL, -1);
sprintf(fmt, "QDP:ColorMatrix%d(%p)", QC(b), b->ptr);
lua_pushstring(L, fmt);
return 1;
}
vector<vector<int> > levelOrder(TreeNode *root) {
vector<vector<int> > result;
if(!root)
return result;
vector<queue<TreeNode *> > Qs(2);
Qs[0].push(root);
int flag = 0;
vector<int> item;
while(!Qs[flag].empty()) {
TreeNode *node = Qs[flag].front();
item.push_back(node->val);
Qs[flag].pop();
if(node->left)
Qs[!flag].push(node->left);
if(node->right)
Qs[!flag].push(node->right);
if(Qs[flag].empty()) {
result.push_back(item);
item.clear();
flag = !flag;
}
}
return result;
}
static void
Qs(float_DD_put)(char *buf, size_t index, int count, void *v_arg)
{
Qs(USQCDArgs) *arg = v_arg;
#if QNc == 'N'
typedef QLA_DN_DiracPropagator(arg->nc, Ptype);
#else
typedef Qx(QLA_D,_DiracPropagator) Ptype;
#endif
Ptype *P = arg->P;
Ptype *dst = &P[index];
QLA_F_Complex *src = (void *)buf;
int js = arg->js;
int jc = arg->jc;
int is, ic;
if (count != 1)
luaL_error(arg->L, "qcd.ddpairs.read(): count != 1");
for (is = 0; is < QDP_Ns; is++) {
for (ic = 0; ic < arg->nc; ic++, src++) {
QLA_c_eq_r_plus_ir(QLA_elem_P(*dst, ic, is, jc, js),
QLA_real(*src), QLA_imag(*src));
}
}
}
//Compute subGLCM features.
std::vector<float> ComputeQFeatures(std::vector<cv::Mat> GLCMs)
{
std::vector<float> Qs(8, 0);
for (size_t i = 0; i < GLCMs.size(); i++) {
cv::Mat GLCM = GLCMs[i];
size_t width_t = GLCM.cols;
size_t height_t = GLCM.rows;
size_t wStep = width_t / 2;
size_t hStep = height_t / 2;
int qIndex = 0;
for (size_t u = 0; u < 2; u++) {
for (size_t j = 0; j < 2; j++) {
if (qIndex == 2)
{
qIndex++;
continue;
}
size_t startX = wStep*j, endX = wStep + wStep*j, startY = hStep*u, endY = hStep + hStep*u;
float Q = ComputeQFeature(GLCM, startX, endX, startY, endY);
Qs[qIndex + 4* i] += Q;
qIndex++;
}
}
}
Qs[2] = Qs[3];
Qs[3] = Qs[4];
Qs[4] = Qs[5];
Qs[5] = Qs[7];
Qs.resize(6);
return Qs;
}
static int
Qs(q_g_mul_p_)(lua_State *L)
{
mClifford *m = qlua_checkClifford(L, 1);
Qs(mSeqDirProp) *f = Qs(qlua_checkSeqDirProp)(L, 2, -1);
Qs(mSeqDirProp) *mf = Qs(qlua_newSeqDirProp)(L, QC(f));
Qs(mSeqDirProp) *r = Qs(qlua_newZeroSeqDirProp)(L, QC(f));
int i;
for (i = 0; i < 16; i++) {
switch (m->g[i].t) {
case qG_z: continue;
case qG_p:
Qx(QLA_D,_P_eq_gamma_times_P)(QNC(QC(f)) mf->ptr, f->ptr, i);
Qx(QLA_D,_P_peq_P)(QNC(QC(f)) r->ptr, mf->ptr);
break;
case qG_m:
Qx(QLA_D,_P_eq_gamma_times_P)(QNC(QC(f)) mf->ptr, f->ptr, i);
Qx(QLA_D,_P_meq_P)(QNC(QC(f)) r->ptr, mf->ptr);
break;
case qG_r:
Qx(QLA_D,_P_eq_gamma_times_P)(QNC(QC(f)) mf->ptr, f->ptr, i);
Qx(QLA_D,_P_peq_r_times_P)(QNC(QC(f)) r->ptr, &m->g[i].r, mf->ptr);
break;
case qG_c:
Qx(QLA_D,_P_eq_gamma_times_P)(QNC(QC(f)) mf->ptr, f->ptr, i);
Qx(QLA_D,_P_peq_c_times_P)(QNC(QC(f)) r->ptr, &m->g[i].c, mf->ptr);
break;
}
}
return 1;
}
static int
Qs(q_g_mul_D_)(lua_State *L)
{
mClifford *m = qlua_checkClifford(L, 1);
Qs(mLatDirFerm) *f = Qs(qlua_checkLatDirFerm)(L, 2, NULL, -1);
mLattice *S = qlua_ObjLattice(L, 2);
int Sidx = lua_gettop(L);
Qs(mLatDirFerm) *mf = Qs(qlua_newLatDirFerm)(L, Sidx, QC(f));
Qs(mLatDirFerm) *r = Qs(qlua_newZeroLatDirFerm)(L, Sidx, QC(f));
int i;
CALL_QDP(L);
for (i = 0; i < 16; i++) {
switch (m->g[i].t) {
case qG_z: continue;
case qG_p:
Qx(QDP_D,_D_eq_gamma_times_D)(mf->ptr, f->ptr, i, *S->qss);
Qx(QDP_D,_D_peq_D)(r->ptr, mf->ptr, *S->qss);
break;
case qG_m:
Qx(QDP_D,_D_eq_gamma_times_D)(mf->ptr, f->ptr, i, *S->qss);
Qx(QDP_D,_D_meq_D)(r->ptr, mf->ptr, *S->qss);
break;
case qG_r:
Qx(QDP_D,_D_eq_gamma_times_D)(mf->ptr, f->ptr, i, *S->qss);
Qx(QDP_D,_D_peq_r_times_D)(r->ptr, &m->g[i].r, mf->ptr, *S->qss);
break;
case qG_c:
Qx(QDP_D,_D_eq_gamma_times_D)(mf->ptr, f->ptr, i, *S->qss);
Qx(QDP_D,_D_peq_c_times_D)(r->ptr, &m->g[i].c, mf->ptr, *S->qss);
break;
}
}
return 1;
}
void thermalBaffle1DFvPatchScalarField<solidType>::updateCoeffs()
{
if (updated())
{
return;
}
// Since we're inside initEvaluate/evaluate there might be processor
// comms underway. Change the tag we use.
int oldTag = UPstream::msgType();
UPstream::msgType() = oldTag+1;
const mapDistribute& mapDist = this->mappedPatchBase::map();
const label patchi = patch().index();
const label nbrPatchi = samplePolyPatch().index();
if (baffleActivated_)
{
const fvPatch& nbrPatch = patch().boundaryMesh()[nbrPatchi];
const compressible::turbulenceModel& turbModel =
db().template lookupObject<compressible::turbulenceModel>
(
"turbulenceModel"
);
// local properties
const scalarField kappaw(turbModel.kappaEff(patchi));
const fvPatchScalarField& Tp =
patch().template lookupPatchField<volScalarField, scalar>(TName_);
scalarField Qr(Tp.size(), 0.0);
if (QrName_ != "none")
{
Qr = patch().template lookupPatchField<volScalarField, scalar>
(QrName_);
Qr = QrRelaxation_*Qr + (1.0 - QrRelaxation_)*QrPrevious_;
QrPrevious_ = Qr;
}
tmp<scalarField> Ti = patchInternalField();
scalarField myKDelta(patch().deltaCoeffs()*kappaw);
// nrb properties
scalarField nbrTp =
turbModel.thermo().T().boundaryField()[nbrPatchi];
mapDist.distribute(nbrTp);
// solid properties
scalarField kappas(patch().size(), 0.0);
forAll(kappas, i)
{
kappas[i] = solid().kappa(0.0, (Tp[i] + nbrTp[i])/2.0);
}
const scalarField KDeltaSolid(kappas/baffleThickness());
const scalarField alpha(KDeltaSolid - Qr/Tp);
valueFraction() = alpha/(alpha + myKDelta);
refValue() = (KDeltaSolid*nbrTp + Qs()/2.0)/alpha;
if (debug)
{
scalar Q = gAverage(kappaw*snGrad());
Info<< patch().boundaryMesh().mesh().name() << ':'
<< patch().name() << ':'
<< this->dimensionedInternalField().name() << " <- "
<< nbrPatch.name() << ':'
<< this->dimensionedInternalField().name() << " :"
<< " heat[W]:" << Q
<< " walltemperature "
<< " min:" << gMin(*this)
<< " max:" << gMax(*this)
<< " avg:" << gAverage(*this)
<< endl;
}
}
