/* Update the vertices and texture coordinates. Orientation must be in {0,90,180,270} */
void QSGVideoNode::setTexturedRectGeometry(const QRectF &rect, const QRectF &textureRect, int orientation)
{
if (rect == m_rect && textureRect == m_textureRect && orientation == m_orientation)
return;
m_rect = rect;
m_textureRect = textureRect;
m_orientation = orientation;
QSGGeometry *g = geometry();
if (g == 0)
g = new QSGGeometry(QSGGeometry::defaultAttributes_TexturedPoint2D(), 4);
QSGGeometry::TexturedPoint2D *v = g->vertexDataAsTexturedPoint2D();
// Set geometry first
qSetGeom(v + 0, rect.topLeft());
qSetGeom(v + 1, rect.bottomLeft());
qSetGeom(v + 2, rect.topRight());
qSetGeom(v + 3, rect.bottomRight());
// and then texture coordinates
switch (orientation) {
default:
// tl, bl, tr, br
qSetTex(v + 0, textureRect.topLeft());
qSetTex(v + 1, textureRect.bottomLeft());
qSetTex(v + 2, textureRect.topRight());
qSetTex(v + 3, textureRect.bottomRight());
break;
case 90:
// tr, tl, br, bl
qSetTex(v + 0, textureRect.topRight());
qSetTex(v + 1, textureRect.topLeft());
qSetTex(v + 2, textureRect.bottomRight());
qSetTex(v + 3, textureRect.bottomLeft());
break;
case 180:
// br, tr, bl, tl
qSetTex(v + 0, textureRect.bottomRight());
qSetTex(v + 1, textureRect.topRight());
qSetTex(v + 2, textureRect.bottomLeft());
qSetTex(v + 3, textureRect.topLeft());
break;
case 270:
// bl, br, tl, tr
qSetTex(v + 0, textureRect.bottomLeft());
qSetTex(v + 1, textureRect.bottomRight());
qSetTex(v + 2, textureRect.topLeft());
qSetTex(v + 3, textureRect.topRight());
break;
}
if (!geometry())
setGeometry(g);
markDirty(DirtyGeometry);
}
void SpanLayoutView::layout_subviews(Canvas &canvas)
{
View::layout_subviews(canvas);
impl->layout_views(canvas, geometry().content.get_width());
}
void MainWindow::saveGeometry() {
// Did not use geometry() on purpose,
// see http://doc.qt.io/qt-5/qsettings.html#restoring-the-state-of-a-gui-application
QRect geometry(pos(), size());
SettingHandles::windowGeometry.set(geometry);
}
PxRigidActor* createBoxActor( const osg::Vec3& dim, double density, PxMaterial* mtl )
{
PxBoxGeometry geometry( PxVec3(dim[0]*0.5f, dim[1]*0.5f, dim[2]*0.5f) );
return createActor( geometry, density, mtl );
}
PxRigidActor* createCapsuleActor( double radius, double height, double density, PxMaterial* mtl )
{
PxCapsuleGeometry geometry( radius, height*0.5 );
return createActor( geometry, density, mtl );
}
bool Model::Load(Deserializer& source)
{
PROFILE(LoadModel);
// Check ID
if (source.ReadFileID() != "UMDL")
{
LOGERROR(source.GetName() + " is not a valid model file");
return false;
}
geometries_.Clear();
geometryBoneMappings_.Clear();
geometryCenters_.Clear();
morphs_.Clear();
vertexBuffers_.Clear();
indexBuffers_.Clear();
unsigned memoryUse = sizeof(Model);
// Read vertex buffers
unsigned numVertexBuffers = source.ReadUInt();
vertexBuffers_.Reserve(numVertexBuffers);
morphRangeStarts_.Resize(numVertexBuffers);
morphRangeCounts_.Resize(numVertexBuffers);
for (unsigned i = 0; i < numVertexBuffers; ++i)
{
unsigned vertexCount = source.ReadUInt();
unsigned elementMask = source.ReadUInt();
morphRangeStarts_[i] = source.ReadUInt();
morphRangeCounts_[i] = source.ReadUInt();
SharedPtr<VertexBuffer> buffer(new VertexBuffer(context_));
buffer->SetShadowed(true);
buffer->SetSize(vertexCount, elementMask);
void* dest = buffer->Lock(0, vertexCount);
unsigned vertexSize = buffer->GetVertexSize();
source.Read(dest, vertexCount * vertexSize);
buffer->Unlock();
memoryUse += sizeof(VertexBuffer) + vertexCount * vertexSize;
vertexBuffers_.Push(buffer);
}
// Read index buffers
unsigned numIndexBuffers = source.ReadUInt();
indexBuffers_.Reserve(numIndexBuffers);
for (unsigned i = 0; i < numIndexBuffers; ++i)
{
unsigned indexCount = source.ReadUInt();
unsigned indexSize = source.ReadUInt();
SharedPtr<IndexBuffer> buffer(new IndexBuffer(context_));
buffer->SetShadowed(true);
buffer->SetSize(indexCount, indexSize > sizeof(unsigned short));
void* dest = buffer->Lock(0, indexCount);
source.Read(dest, indexCount * indexSize);
buffer->Unlock();
memoryUse += sizeof(IndexBuffer) + indexCount * indexSize;
indexBuffers_.Push(buffer);
}
// Read geometries
unsigned numGeometries = source.ReadUInt();
geometries_.Reserve(numGeometries);
geometryBoneMappings_.Reserve(numGeometries);
geometryCenters_.Reserve(numGeometries);
for (unsigned i = 0; i < numGeometries; ++i)
{
// Read bone mappings
unsigned boneMappingCount = source.ReadUInt();
PODVector<unsigned> boneMapping(boneMappingCount);
for (unsigned j = 0; j < boneMappingCount; ++j)
boneMapping[j] = source.ReadUInt();
geometryBoneMappings_.Push(boneMapping);
unsigned numLodLevels = source.ReadUInt();
Vector<SharedPtr<Geometry> > geometryLodLevels;
geometryLodLevels.Reserve(numLodLevels);
for (unsigned j = 0; j < numLodLevels; ++j)
{
float distance = source.ReadFloat();
PrimitiveType type = (PrimitiveType)source.ReadUInt();
unsigned vertexBufferRef = source.ReadUInt();
unsigned indexBufferRef = source.ReadUInt();
unsigned indexStart = source.ReadUInt();
unsigned indexCount = source.ReadUInt();
if (vertexBufferRef >= vertexBuffers_.Size())
{
LOGERROR("Vertex buffer index out of bounds");
return false;
}
if (indexBufferRef >= indexBuffers_.Size())
{
LOGERROR("Index buffer index out of bounds");
return false;
}
SharedPtr<Geometry> geometry(new Geometry(context_));
geometry->SetVertexBuffer(0, vertexBuffers_[vertexBufferRef]);
geometry->SetIndexBuffer(indexBuffers_[indexBufferRef]);
geometry->SetDrawRange(type, indexStart, indexCount);
geometry->SetLodDistance(distance);
geometryLodLevels.Push(geometry);
memoryUse += sizeof(Geometry);
}
geometries_.Push(geometryLodLevels);
}
// Read morphs
unsigned numMorphs = source.ReadUInt();
morphs_.Reserve(numMorphs);
for (unsigned i = 0; i < numMorphs; ++i)
{
ModelMorph newMorph;
newMorph.name_ = source.ReadString();
newMorph.nameHash_ = newMorph.name_;
newMorph.weight_ = 0.0f;
unsigned nubuffers_ = source.ReadUInt();
for (unsigned j = 0; j < nubuffers_; ++j)
{
VertexBufferMorph newBuffer;
unsigned bufferIndex = source.ReadUInt();
newBuffer.elementMask_ = source.ReadUInt();
newBuffer.vertexCount_ = source.ReadUInt();
// Base size: size of each vertex index
unsigned vertexSize = sizeof(unsigned);
// Add size of individual elements
if (newBuffer.elementMask_ & MASK_POSITION)
vertexSize += sizeof(Vector3);
if (newBuffer.elementMask_ & MASK_NORMAL)
vertexSize += sizeof(Vector3);
if (newBuffer.elementMask_ & MASK_TANGENT)
vertexSize += sizeof(Vector3);
newBuffer.morphData_ = new unsigned char[newBuffer.vertexCount_ * vertexSize];
source.Read(&newBuffer.morphData_[0], newBuffer.vertexCount_ * vertexSize);
newMorph.buffers_[bufferIndex] = newBuffer;
memoryUse += sizeof(VertexBufferMorph) + newBuffer.vertexCount_ * vertexSize;
}
morphs_.Push(newMorph);
memoryUse += sizeof(ModelMorph);
}
// Read skeleton
skeleton_.Load(source);
memoryUse += skeleton_.GetNumBones() * sizeof(Bone);
// Read bounding box
boundingBox_ = source.ReadBoundingBox();
// Read geometry centers
for (unsigned i = 0; i < geometries_.Size() && !source.IsEof(); ++i)
geometryCenters_.Push(source.ReadVector3());
while (geometryCenters_.Size() < geometries_.Size())
geometryCenters_.Push(Vector3::ZERO);
memoryUse += sizeof(Vector3) * geometries_.Size();
SetMemoryUse(memoryUse);
return true;
}
void ossimPlanetAnnotationPlacemark::traverse(osg::NodeVisitor& nv)
{
if(!enableFlag()) return;
if(!isStaged())
{
return;
}
bool needRedraw = false;
OpenThreads::ScopedLock<OpenThreads::Mutex> lock(theUpdateMutex);
switch(nv.getVisitorType())
{
case osg::NodeVisitor::UPDATE_VISITOR:
{
if(isDirtyBitSet(LABEL_DIRTY)&&theLabel.valid())
{
needRedraw = true;
if(theLabel.valid())
{
theLabel->setText(name().c_str());
}
clearDirtyBit(LABEL_DIRTY);
}
if(isDirtyBitSet(COLOR_DIRTY))
{
if(theLabelStyle.valid()&&theLabel.valid())
{
needRedraw = true;
if(theLabelStyle->colorMode() == ossimPlanetAnnotationColorMode_NORMAL)
{
theLabel->setColor(theLabelStyle->color());
}
else
{
theLabel->setColor(osg::Vec4d((ossim_float64)rand()/(ossim_float64)RAND_MAX,
(ossim_float64)rand()/(ossim_float64)RAND_MAX,
(ossim_float64)rand()/(ossim_float64)RAND_MAX,
theLabel->getColor()[3]));
}
}
clearDirtyBit(COLOR_DIRTY);
}
if(theExpireTime.valid())
{
if(theExpireTime->hasExpired())
{
if(theLayer)
{
theLayer->needsRemoving(this);
}
}
needRedraw = true;
}
break;
}
default:
{
break;
}
}
if(theExpireTime.valid())
{
setRedrawFlag(true);
}
osg::ref_ptr<ossimPlanetAnnotationGeometry> geom = geometry();
if(geom.valid())
{
geom->traverse(nv);
}
ossimPlanetAnnotationLayerNode::traverse(nv);
if(needRedraw) setRedrawFlag(true);
}
int main(int argc, char **argv) {
int c, i, j, mu, nu;
int count = 0;
int filename_set = 0;
int dims[4] = {0,0,0,0};
int l_LX_at, l_LXstart_at;
int x0, x1, x2, x3, ix, iix, it;
int sid, status, gid;
double **corr=NULL, **corr2=NULL;
double *tcorr=NULL, *tcorr2=NULL;
double *work = (double*)NULL;
double q[4], fnorm;
int verbose = 0;
int do_gt = 0;
int nsource=0;
char filename[100], contype[200];
double ratime, retime;
double plaq;
double spinor1[24], spinor2[24], U_[18];
double *gauge_trafo=(double*)NULL;
double mom2, mom4;
complex w, w1, *cp1, *cp2, *cp3;
FILE *ofs;
#ifdef MPI
// MPI_Init(&argc, &argv);
fprintf(stderr, "[jc_ud_x] Error, only non-mpi version implemented\n");
exit(1);
#endif
while ((c = getopt(argc, argv, "h?f:")) != -1) {
switch (c) {
case 'f':
strcpy(filename, optarg);
filename_set=1;
break;
case 'h':
case '?':
default:
usage();
break;
}
}
/* set the default values */
if(filename_set==0) strcpy(filename, "cvc.input");
fprintf(stdout, "# Reading input from file %s\n", filename);
read_input_parser(filename);
/* some checks on the input data */
if((T_global == 0) || (LX==0) || (LY==0) || (LZ==0)) {
if(g_proc_id==0) fprintf(stdout, "T and L's must be set\n");
usage();
}
if(g_kappa == 0.) {
if(g_proc_id==0) fprintf(stdout, "kappa should be > 0.n");
usage();
}
fprintf(stdout, "\n**************************************************\n");
fprintf(stdout, "* jc_ud_x\n");
fprintf(stdout, "**************************************************\n\n");
/*********************************
* initialize MPI parameters
*********************************/
// mpi_init(argc, argv);
/* initialize */
dims[0]=T_global; dims[1]=LX; dims[2]=LY; dims[3]=LZ;
T = T_global;
Tstart = 0;
l_LX_at = LX;
l_LXstart_at = 0;
FFTW_LOC_VOLUME = T*LX*LY*LZ;
fprintf(stdout, "# [%2d] parameters:\n"\
"# [%2d] T = %3d\n"\
"# [%2d] Tstart = %3d\n"\
"# [%2d] l_LX_at = %3d\n"\
"# [%2d] l_LXstart_at = %3d\n"\
"# [%2d] FFTW_LOC_VOLUME = %3d\n",
g_cart_id, g_cart_id, T, g_cart_id, Tstart, g_cart_id, l_LX_at,
g_cart_id, l_LXstart_at, g_cart_id, FFTW_LOC_VOLUME);
if(init_geometry() != 0) {
fprintf(stderr, "ERROR from init_geometry\n");
exit(1);
}
geometry();
/*************************************************
* allocate mem for gauge field and spinor fields
*************************************************/
alloc_gauge_field(&g_gauge_field, VOLUMEPLUSRAND);
no_fields = 2;
g_spinor_field = (double**)calloc(no_fields, sizeof(double*));
for(i=0; i<no_fields; i++) alloc_spinor_field(&g_spinor_field[i], VOLUMEPLUSRAND);
/****************************************
* allocate memory for the contractions
****************************************/
nsource = (g_sourceid2 - g_sourceid + 1) / g_sourceid_step;
if(g_cart_id==0) fprintf(stdout, "# nsource = %d\n", nsource);
corr = (double**)calloc( nsource, sizeof(double*));
corr[0] = (double*)calloc( nsource*T*8, sizeof(double));
for(i=1;i<nsource;i++) corr[i] = corr[i-1] + 8*T;
corr2 = (double**)calloc( nsource, sizeof(double*));
corr2[0] = (double*)calloc( nsource*8*T, sizeof(double));
for(i=1;i<nsource;i++) corr2[i] = corr2[i-1] + 8*T;
tcorr = (double*)calloc(T*8, sizeof(double));
tcorr2 = (double*)calloc(T*8, sizeof(double));
/***********************************************
* start loop on gauge id.s
***********************************************/
for(gid=g_gaugeid; gid<=g_gaugeid2; gid++) {
sprintf(filename, "%s.%.4d", gaugefilename_prefix, gid);
if(g_cart_id==0) fprintf(stdout, "# reading gauge field from file %s\n", filename);
read_lime_gauge_field_doubleprec(filename);
xchange_gauge();
plaquette(&plaq);
if(g_cart_id==0) fprintf(stdout, "# measured plaquette value: %25.16e\n", plaq);
/* reset disc to zero */
for(ix=0; ix<nsource*8*T; ix++) corr[0][ix] = 0.;
for(ix=0; ix<nsource*8*T; ix++) corr2[0][ix] = 0.;
count=0;
/***********************************************
* start loop on source id.s
***********************************************/
for(sid=g_sourceid; sid<=g_sourceid2; sid+=g_sourceid_step) {
/* read the new propagator to g_spinor_field[0] */
ratime = (double)clock() / CLOCKS_PER_SEC;
if(format==0) {
sprintf(filename, "%s.%.4d.%.2d.inverted", filename_prefix, gid, sid);
if(read_lime_spinor(g_spinor_field[0], filename, 0) != 0) break;
}
else if(format==1) {
sprintf(filename, "%s.%.4d.%.5d.inverted", filename_prefix, gid, sid);
if(read_cmi(g_spinor_field[0], filename) != 0) {
fprintf(stderr, "\nError from read_cmi\n");
break;
}
}
xchange_field(g_spinor_field[0]);
retime = (double)clock() / CLOCKS_PER_SEC;
if(g_cart_id==0) fprintf(stdout, "# time to read prop.: %e seconds\n", retime-ratime);
ratime = (double)clock() / CLOCKS_PER_SEC;
/* apply [1] = D_tm [0] */
Q_phi_tbc(g_spinor_field[1], g_spinor_field[0]);
xchange_field(g_spinor_field[1]);
retime = (double)clock() / CLOCKS_PER_SEC;
if(g_cart_id==0) fprintf(stdout, "# time to apply D_W: %e seconds\n", retime-ratime);
ratime = (double)clock() / CLOCKS_PER_SEC;
/* calculate real and imaginary part */
for(mu=0; mu<4; mu++) {
for(x0=0; x0<T; x0++) {
for(x1=0; x1<LX; x1++) {
for(x2=0; x2<LY; x2++) {
for(x3=0; x3<LZ; x3++) {
ix = g_ipt[x0][x1][x2][x3];
_cm_eq_cm_ti_co(U_, g_gauge_field+_GGI(ix,mu), &(co_phase_up[mu]));
_fv_eq_cm_ti_fv(spinor1, U_, &g_spinor_field[0][_GSI(g_iup[ix][mu])]);
_fv_eq_gamma_ti_fv(spinor2, mu, spinor1);
_fv_mi_eq_fv(spinor2, spinor1);
_co_eq_fv_dag_ti_fv(&w, &g_spinor_field[1][_GSI(ix)], spinor2);
corr[count][2*(mu*T+x0) ] -= 0.5*w.re;
corr[count][2*(mu*T+x0)+1] -= 0.5*w.im;
_fv_eq_cm_dag_ti_fv(spinor1, U_, &g_spinor_field[0][_GSI(ix)]);
_fv_eq_gamma_ti_fv(spinor2, mu, spinor1);
_fv_pl_eq_fv(spinor2, spinor1);
_co_eq_fv_dag_ti_fv(&w, &g_spinor_field[1][_GSI(g_iup[ix][mu])], spinor2);
corr[count][2*(mu*T+x0) ] -= 0.5*w.re;
corr[count][2*(mu*T+x0)+1] -= 0.5*w.im;
_fv_eq_gamma_ti_fv(spinor1, mu, &g_spinor_field[0][_GSI(ix)]);
_co_eq_fv_dag_ti_fv(&w, &g_spinor_field[1][_GSI(ix)], spinor1);
corr2[count][2*(mu*T+x0) ] -= w.re;
corr2[count][2*(mu*T+x0)+1] -= w.im;
}}}
}
} // of mu
count++;
} // of sid
retime = (double)clock() / CLOCKS_PER_SEC;
if(g_cart_id==0) fprintf(stdout, "# time to calculate contractions: %e seconds\n", retime-ratime);
for(ix=0;ix<8*T;ix++) tcorr[ix] = 0.;
for(ix=0;ix<8*T;ix++) tcorr2[ix] = 0.;
for(i=0;i<nsource-1;i++) {
for(j=i+1;j<nsource;j++) {
for(mu=0;mu<4;mu++) {
for(x0=0;x0<T;x0++) { // times at source
for(x1=0;x1<T;x1++) { // times at sink
it = (x1 - x0 + T) % T;
// conserved current
tcorr[2*(mu*T+it) ] += corr[i][2*(mu*T+x1)] * corr[j][2*(mu*T+x0) ] - corr[i][2*(mu*T+x1)+1] * corr[j][2*(mu*T+x0)+1];
tcorr[2*(mu*T+it)+1] += corr[i][2*(mu*T+x1)] * corr[j][2*(mu*T+x0)+1] + corr[i][2*(mu*T+x1)+1] * corr[j][2*(mu*T+x0) ];
tcorr[2*(mu*T+it) ] += corr[j][2*(mu*T+x1)] * corr[i][2*(mu*T+x0) ] - corr[j][2*(mu*T+x1)+1] * corr[i][2*(mu*T+x0)+1];
tcorr[2*(mu*T+it)+1] += corr[j][2*(mu*T+x1)] * corr[i][2*(mu*T+x0)+1] + corr[j][2*(mu*T+x1)+1] * corr[i][2*(mu*T+x0) ];
// local current
tcorr2[2*(mu*T+it) ] += corr2[i][2*(mu*T+x1)] * corr2[j][2*(mu*T+x0) ] - corr2[i][2*(mu*T+x1)+1] * corr2[j][2*(mu*T+x0)+1];
tcorr2[2*(mu*T+it)+1] += corr2[i][2*(mu*T+x1)] * corr2[j][2*(mu*T+x0)+1] + corr2[i][2*(mu*T+x1)+1] * corr2[j][2*(mu*T+x0) ];
tcorr2[2*(mu*T+it) ] += corr2[j][2*(mu*T+x1)] * corr2[i][2*(mu*T+x0) ] - corr2[j][2*(mu*T+x1)+1] * corr2[i][2*(mu*T+x0)+1];
tcorr2[2*(mu*T+it)+1] += corr2[j][2*(mu*T+x1)] * corr2[i][2*(mu*T+x0)+1] + corr2[j][2*(mu*T+x1)+1] * corr2[i][2*(mu*T+x0) ];
}}
}
}}
fnorm = 1. / ( g_prop_normsqr * g_prop_normsqr * (double)(LX*LY*LZ) * (double)(LX*LY*LZ) * nsource * (nsource-1));
if(g_cart_id==0) fprintf(stdout, "X-fnorm = %e\n", fnorm);
for(ix=0;ix<8*T;ix++) tcorr[ix] *= fnorm;
for(ix=0;ix<8*T;ix++) tcorr2[ix] *= fnorm;
/************************************************
* save results
************************************************/
if(g_cart_id == 0) fprintf(stdout, "# save results for gauge id %d and sid %d\n", gid, sid);
/* save the result in position space */
sprintf(filename, "jc_u_tp0.%.4d.%.4d", gid, sid);
ofs = fopen(filename, "w");
for(x0=0;x0<T;x0++) fprintf(ofs, "%d%25.16e%25.16e%25.16e%25.16e%25.16e%25.16e%25.16e%25.16e\n", x0,
tcorr[2*(0*T+x0)], tcorr[2*(0*T+x0)+1],
tcorr[2*(1*T+x0)], tcorr[2*(1*T+x0)+1],
tcorr[2*(2*T+x0)], tcorr[2*(2*T+x0)+1],
tcorr[2*(3*T+x0)], tcorr[2*(3*T+x0)+1]);
fclose(ofs);
} /* of loop on gid */
/***********************************************
* free the allocated memory, finalize
***********************************************/
free(g_gauge_field);
for(i=0; i<no_fields; i++) free(g_spinor_field[i]);
free(g_spinor_field);
free_geometry();
free(corr);
free(corr2);
free(tcorr);
free(tcorr2);
return(0);
}
void PrinterEmulator::saveSettings()
{
settings.beginGroup(WINDOW_SETTINGS_GROUP);
settings.setValue(WINDOW_GEOMETRY_SETTING, geometry());
settings.endGroup();
}
/**
* @brief encoding the diff scheme, i.e. the output value for input intensity I is equal to the sum
* of the neighborhood divided by the number of pixels traversed in the neighborhood
*
* @param _in input stack of type raw_type
* @param _out output stack of type compressed_type but same extent than the input
* @param _shape dimensions of the input _in
* @return sqeazy::error_code
*/
compressed_type* encode( const raw_type* _raw,
compressed_type* _compressed,
const std::vector<std::size_t>& _shape) override final {
typedef std::size_t size_type;
if(_shape.size()!=3){
std::cerr << "[diff_scheme] unable to process input data that is not 3D\n";
return _compressed;
}
std::size_t length = std::accumulate(_shape.begin(), _shape.end(),1,std::multiplies<std::size_t>());
std::copy(_raw, _raw + length, _compressed);//crossing fingers due to possible type mismatch
std::vector<size_type> offsets;
sqeazy::halo<Neighborhood, size_type> geometry(_shape[row_major::w],
_shape[row_major::h],
_shape[row_major::d]);
geometry.compute_offsets_in_x(offsets);
size_type halo_size_x = geometry.non_halo_end(0)-geometry.non_halo_begin(0);
if(offsets.size()==1)//no offsets in other dimensions than x
{
halo_size_x = length - offsets.front();
}
const auto n_traversed_pixels = sqeazy::num_traversed_pixels<Neighborhood>();
out_type* signed_compressed = reinterpret_cast<out_type*>(_compressed);
auto offsets_begin = offsets.begin();
const omp_size_type offset_size = offsets.size();
const int nthreads = this->n_threads();
//auto shape_begin = _shape.data();
const auto pshape = _shape.data();
#pragma omp parallel for \
shared(signed_compressed) \
firstprivate( offset_size, halo_size_x, offsets_begin, _raw, pshape,n_traversed_pixels) \
num_threads(nthreads)
for(omp_size_type offset = 0; offset < offset_size; ++offset) {
for(std::size_t index = 0; index < (std::size_t)halo_size_x; ++index) {
auto local_index = index + *(offsets_begin + offset);
sum_type local_sum = naive_sum<Neighborhood>(_raw,local_index,
pshape[row_major::w],
pshape[row_major::h],
pshape[row_major::d]
);
*(signed_compressed + local_index) = _raw[local_index] - local_sum/n_traversed_pixels;
}
}
return _compressed+length;
}
void QSGBasicInternalImageNode::updateGeometry()
{
Q_ASSERT(!m_targetRect.isEmpty());
const QSGTexture *t = materialTexture();
if (!t) {
QSGGeometry *g = geometry();
g->allocate(4);
g->setDrawingMode(QSGGeometry::DrawTriangleStrip);
memset(g->vertexData(), 0, g->sizeOfVertex() * 4);
} else {
QRectF sourceRect = t->normalizedTextureSubRect();
QRectF innerSourceRect(sourceRect.x() + m_innerSourceRect.x() * sourceRect.width(),
sourceRect.y() + m_innerSourceRect.y() * sourceRect.height(),
m_innerSourceRect.width() * sourceRect.width(),
m_innerSourceRect.height() * sourceRect.height());
bool hasMargins = m_targetRect != m_innerTargetRect;
int floorLeft = qFloor(m_subSourceRect.left());
int ceilRight = qCeil(m_subSourceRect.right());
int floorTop = qFloor(m_subSourceRect.top());
int ceilBottom = qCeil(m_subSourceRect.bottom());
int hTiles = ceilRight - floorLeft;
int vTiles = ceilBottom - floorTop;
bool hasTiles = hTiles != 1 || vTiles != 1;
bool fullTexture = innerSourceRect == QRectF(0, 0, 1, 1);
// An image can be rendered as a single quad if:
// - There are no margins, and either:
// - the image isn't repeated
// - the source rectangle fills the entire texture so that texture wrapping can be used,
// and NPOT is supported
if (!hasMargins && (!hasTiles || (fullTexture && supportsWrap(t->textureSize())))) {
QRectF sr;
if (!fullTexture) {
sr = QRectF(innerSourceRect.x() + (m_subSourceRect.left() - floorLeft) * innerSourceRect.width(),
innerSourceRect.y() + (m_subSourceRect.top() - floorTop) * innerSourceRect.height(),
m_subSourceRect.width() * innerSourceRect.width(),
m_subSourceRect.height() * innerSourceRect.height());
} else {
sr = QRectF(m_subSourceRect.left() - floorLeft, m_subSourceRect.top() - floorTop,
m_subSourceRect.width(), m_subSourceRect.height());
}
if (m_mirror) {
qreal oldLeft = sr.left();
sr.setLeft(sr.right());
sr.setRight(oldLeft);
}
if (m_antialiasing) {
QSGGeometry *g = geometry();
Q_ASSERT(g != &m_geometry);
g->allocate(8, 14);
g->setDrawingMode(QSGGeometry::DrawTriangleStrip);
SmoothVertex *vertices = reinterpret_cast<SmoothVertex *>(g->vertexData());
float delta = float(qAbs(m_targetRect.width()) < qAbs(m_targetRect.height())
? m_targetRect.width() : m_targetRect.height()) * 0.5f;
float sx = float(sr.width() / m_targetRect.width());
float sy = float(sr.height() / m_targetRect.height());
for (int d = -1; d <= 1; d += 2) {
for (int j = 0; j < 2; ++j) {
for (int i = 0; i < 2; ++i, ++vertices) {
vertices->x = m_targetRect.x() + i * m_targetRect.width();
vertices->y = m_targetRect.y() + j * m_targetRect.height();
vertices->u = sr.x() + i * sr.width();
vertices->v = sr.y() + j * sr.height();
vertices->dx = (i == 0 ? delta : -delta) * d;
vertices->dy = (j == 0 ? delta : -delta) * d;
vertices->du = (d < 0 ? 0 : vertices->dx * sx);
vertices->dv = (d < 0 ? 0 : vertices->dy * sy);
}
}
}
Q_ASSERT(vertices - g->vertexCount() == g->vertexData());
static const quint16 indices[] = {
0, 4, 1, 5, 3, 7, 2, 6, 0, 4,
4, 6, 5, 7
};
Q_ASSERT(g->sizeOfIndex() * g->indexCount() == sizeof(indices));
memcpy(g->indexDataAsUShort(), indices, sizeof(indices));
} else {
m_geometry.allocate(4);
m_geometry.setDrawingMode(QSGGeometry::DrawTriangleStrip);
QSGGeometry::updateTexturedRectGeometry(&m_geometry, m_targetRect, sr);
}
} else {
QSGGeometry *g = m_antialiasing ? geometry() : &m_geometry;
updateGeometry(m_targetRect, m_innerTargetRect,
sourceRect, innerSourceRect, m_subSourceRect,
g, m_mirror, m_antialiasing);
}
}
markDirty(DirtyGeometry);
m_dirtyGeometry = false;
}
void Window::synchronizeLayout() {
if(visible() && applicationState.quit == false) setGeometry(geometry());
}
/*
Purpose:
RUN is the main program for FEM1D.
Discussion:
FEM1D solves a one dimensional ODE using the finite element method.
The differential equation solved is
- d/dX (P dU/dX) + Q U = F
The finite-element method uses piecewise linear basis functions.
Here U is an unknown scalar function of X defined on the
interval [XL,XR], and P, Q and F are given functions of X.
The values of U or U' at XL and XR are also specified.
The interval [XL,XR] is "meshed" with NSUB+1 points,
XN(0) = XL, XN(1)=XL+H, XN(2)=XL+2*H, ..., XN(NSUB)=XR.
This creates NSUB subintervals, with interval number 1
having endpoints XN(0) and XN(1), and so on up to interval
NSUB, which has endpoints XN(NSUB-1) and XN(NSUB).
Licensing:
This code is distributed under the GNU LGPL license.
Modified:
29 May 2009
Author:
C version by John Burkardt
Parameters:
double ADIAG(NU), the "diagonal" coefficients. That is, ADIAG(I) is
the coefficient of the I-th unknown in the I-th equation.
double ALEFT(NU), the "left hand" coefficients. That is, ALEFT(I)
is the coefficient of the (I-1)-th unknown in the I-th equation.
There is no value in ALEFT(1), since the first equation
does not refer to a "0-th" unknown.
double ARITE(NU).
ARITE(I) is the "right hand" coefficient of the I-th
equation in the linear system. ARITE(I) is the coefficient
of the (I+1)-th unknown in the I-th equation. There is
no value in ARITE(NU) because the NU-th equation does not
refer to an "NU+1"-th unknown.
double F(NU).
ASSEMBLE stores into F the right hand side of the linear
equations.
SOLVE replaces those values of F by the solution of the
linear equations.
double H(NSUB)
H(I) is the length of subinterval I. This code uses
equal spacing for all the subintervals.
int IBC.
IBC declares what the boundary conditions are.
1, at the left endpoint, U has the value UL,
at the right endpoint, U' has the value UR.
2, at the left endpoint, U' has the value UL,
at the right endpoint, U has the value UR.
3, at the left endpoint, U has the value UL,
and at the right endpoint, U has the value UR.
4, at the left endpoint, U' has the value UL,
at the right endpoint U' has the value UR.
int INDX[NSUB+1].
For a node I, INDX(I) is the index of the unknown
associated with node I.
If INDX(I) is equal to -1, then no unknown is associated
with the node, because a boundary condition fixing the
value of U has been applied at the node instead.
Unknowns are numbered beginning with 1.
If IBC is 2 or 4, then there is an unknown value of U
at node 0, which will be unknown number 1. Otherwise,
unknown number 1 will be associated with node 1.
If IBC is 1 or 4, then there is an unknown value of U
at node NSUB, which will be unknown NSUB or NSUB+1,
depending on whether there was an unknown at node 0.
int NL.
The number of basis functions used in a single
subinterval. (NL-1) is the degree of the polynomials
used. For this code, NL is fixed at 2, meaning that
piecewise linear functions are used as the basis.
int NODE[NL*NSUB].
For each subinterval I:
NODE[0+I*2] is the number of the left node, and
NODE[1+I*2] is the number of the right node.
int NQUAD.
The number of quadrature points used in a subinterval.
This code uses NQUAD = 1.
int NSUB.
The number of subintervals into which the interval [XL,XR] is broken.
int NU.
NU is the number of unknowns in the linear system.
Depending on the value of IBC, there will be NSUB-1,
NSUB, or NSUB+1 unknown values, which are the coefficients
of basis functions.
double UL.
If IBC is 1 or 3, UL is the value that U is required
to have at X = XL.
If IBC is 2 or 4, UL is the value that U' is required
to have at X = XL.
double UR.
If IBC is 2 or 3, UR is the value that U is required
to have at X = XR.
If IBC is 1 or 4, UR is the value that U' is required
to have at X = XR.
double XL.
XL is the left endpoint of the interval over which the
differential equation is being solved.
double XN(0:NSUB).
XN(I) is the location of the I-th node. XN(0) is XL,
and XN(NSUB) is XR.
double XQUAD(NSUB)
XQUAD(I) is the location of the single quadrature point
in interval I.
double XR.
XR is the right endpoint of the interval over which the
differential equation is being solved.
*/
double run(int version){
if(version == 0){
if((fp_out = fopen("old_out.txt", "w+")) == NULL ||
(fp_sol = fopen("old_sol.txt", "a")) == NULL){
printf("Old Version files not found.\n");
exit(EXIT_FAILURE);
}
}
else if (version == 1){
if((fp_out = fopen("new_out.txt", "w+")) == NULL ||
(fp_sol = fopen("new_sol.txt", "a")) == NULL){
printf("New Version files not found.\n");
exit(EXIT_FAILURE);
}
}
//Allocate array memory
adiag = (double *)malloc(sizeof(double)*(double)(NSUB+1));
aleft = (double *)malloc(sizeof(double)*(double)(NSUB+1));
arite = (double *)malloc(sizeof(double)*(double)(NSUB+1));
f = (double *)malloc(sizeof(double)*(double)(NSUB+1));
h = (double *)malloc(sizeof(double)*(double)(NSUB));
indx = (int *)malloc(sizeof(int)*(int)(NSUB+1));
node = (int *)malloc(sizeof(int)*((int)NL*(int)NSUB));
xn = (double *)malloc(sizeof(double)*(double)(NSUB+1));
xquad = (double *)malloc(sizeof(double)*(double)(NSUB));
//START TIMER//
double begin, end, time_spent;
begin = omp_get_wtime();
timestamp ();
fprintf (fp_out, "\n" );
fprintf (fp_out, "FEM1D\n" );
fprintf (fp_out, " C version\n" );
fprintf (fp_out, "\n" );
fprintf (fp_out, " Solve the two-point boundary value problem\n" );
fprintf (fp_out, "\n" );
fprintf (fp_out, " - d/dX (P dU/dX) + Q U = F\n" );
fprintf (fp_out, "\n" );
fprintf (fp_out, " on the interval [XL,XR], specifying\n" );
fprintf (fp_out," the value of U or U' at each end.\n" );
fprintf (fp_out, "\n" );
fprintf (fp_out," The interval [XL,XR] is broken into NSUB = %ld subintervals\n", NSUB );
fprintf (fp_out, " Number of basis functions per element is NL = %ld\n", NL );
/*
Initialize the data.
*/
init ();
/*
Compute the geometric quantities.
*/
geometry (version);
/*
Assemble the linear system.
*/
assemble (version);
/*
Print out the linear system.
*/
prsys ();
/*
Solve the linear system.
*/
solve ();
/*
Print out the solution.
*/
output (version);
/*
Terminate.
*/
fprintf (fp_out, "\n" );
fprintf (fp_out,"FEM1D:\n" );
fprintf (fp_out, " Normal end of execution.\n" );
fprintf ( fp_out,"\n" );
//END TIMER//
end = omp_get_wtime();
time_spent = end - begin;
timestamp ( );
//CLOSE STREAMS
fclose(fp_out);
fclose(fp_sol);
//FREE MEMORY
free(adiag);
free(aleft);
free(arite);
free(f);
free(h);
free(indx);
free(node);
free(xn);
free(xquad);
return time_spent;
}
// gpSummaryTab
bool gpSummaryTab::Init(gpTraceDataContainer* pDataContainer, gpTraceView* pSessionView, quint64 timelineStart, quint64 timelineEnd)
{
m_pTraceView = pSessionView;
m_pSessionDataContainer = pDataContainer;
if (m_callType == API_CALL || m_callType == GPU_CALL)
{
m_pSummaryTable = new gpTraceSummaryTable(pDataContainer, pSessionView, m_callType);
}
else
{
m_pSummaryTable = new gpCommandListSummaryTable(pDataContainer, pSessionView);
}
m_pTop20Table = new acListCtrl(this);
m_pTop20Caption = new QLabel(this);
QHBoxLayout* pHLayout = new QHBoxLayout;
QVBoxLayout* pLeftVLayout = new QVBoxLayout;
QHBoxLayout* pLeftHLayout = new QHBoxLayout;
QVBoxLayout* pRightVLayout = new QVBoxLayout;
QSplitter* pVMainSplitter = new QSplitter(Qt::Horizontal);
// set the margins to all layout to zero
pHLayout->setContentsMargins(0, 0, 0, 0);
pLeftVLayout->setContentsMargins(0, 0, 0, 0);
pLeftHLayout->setContentsMargins(0, 0, 0, 0);
pRightVLayout->setContentsMargins(0, 0, 0, 0);
QWidget* pWidgetLeft = new QWidget;
QWidget* pWidgetRight = new QWidget;
QString caption;
if (m_callType == eCallType::API_CALL || m_callType == eCallType::GPU_CALL)
{
caption = table_captions[m_callType];
}
else
{
int commandCount = m_pSessionDataContainer->CommandListsData().size();
if (pDataContainer->SessionAPIType() == ProfileSessionDataItem::ProfileItemAPIType::DX12_API_PROFILE_ITEM)
{
caption = QString(table_captions[m_callType]).arg(commandCount);
}
else
{
caption = QString(table_captions[m_callType + 1]).arg(commandCount);
}
}
QLabel* pLblCaption = new QLabel(caption);
m_pChkboxUseScope = new QCheckBox(GPU_STR_Use_Scope_Summary);
pLeftHLayout->addWidget(pLblCaption, 0, Qt::AlignTop);
QString styleSheetStr = QString("QDialog{border:1px solid gray;}");
pLeftHLayout->addWidget(m_pChkboxUseScope, 0, Qt::AlignHCenter);
pLeftVLayout->addLayout(pLeftHLayout);
pWidgetLeft->setLayout(pLeftVLayout);
pWidgetLeft->setStyleSheet(styleSheetStr);
m_pSummaryTable->resizeColumnsToContents();
pLeftVLayout->addWidget(m_pSummaryTable);
m_pSummaryTable->SetSelectionBackgroundColor(acQAMD_CYAN_SELECTION_BKG_COLOUR);
pLeftVLayout->setSpacing(0);
pLeftVLayout->setMargin(1);
QFont qFont = m_pTop20Caption->font();
qFont.setBold(true);
qFont.setPointSize(10);
m_pTop20Caption->setFont(qFont);
pRightVLayout->addWidget(m_pTop20Caption);
QStringList columnCaptions;
columnCaptions << "#";
if (m_callType == API_CALL)
{
columnCaptions << GP_STR_SummaryTop20TableColumnThreadId;
}
columnCaptions << GP_STR_SummaryTop20TableColumnCallIndex;
columnCaptions << GP_STR_SummaryTop20TableColumnTime;
m_pTop20Table->initHeaders(columnCaptions, false);
m_pTop20Table->setContextMenuPolicy(Qt::NoContextMenu);
m_pTop20Table->resizeColumnsToContents();
m_pTop20Table->setShowGrid(true);
pRightVLayout->addWidget(m_pTop20Table);
pRightVLayout->setSpacing(1);
pRightVLayout->setMargin(1);
pWidgetRight->setLayout(pRightVLayout);
pWidgetRight->setStyleSheet(styleSheetStr);
pVMainSplitter->addWidget(pWidgetLeft);
pVMainSplitter->addWidget(pWidgetRight);
int widgetWidth = geometry().right() - geometry().left();
QList<int> sizes;
sizes << (int)(widgetWidth * 0.75);
sizes << (int)(widgetWidth * 0.25);
pVMainSplitter->setSizes(sizes);
pHLayout->addWidget(pVMainSplitter);
setLayout(pHLayout);
SetTimelineScope(false, timelineStart, timelineEnd);
bool rc = connect(m_pChkboxUseScope, SIGNAL(toggled(bool)), this, SLOT(OnUseTimelineSelectionScopeChanged(bool)));
GT_ASSERT(rc);
rc = connect(m_pSummaryTable, SIGNAL(itemSelectionChanged()), this, SLOT(OnSummaryTableSelectionChanged()));
GT_ASSERT(rc);
rc = connect(m_pSummaryTable, SIGNAL(cellClicked(int, int)), this, SLOT(OnSummaryTableCellClicked(int, int)));
GT_ASSERT(rc);
rc = connect(m_pSummaryTable, SIGNAL(cellDoubleClicked(int, int)), this, SLOT(OnSummaryTableCellDoubleClicked(int, int)));
GT_ASSERT(rc);
rc = connect(m_pTop20Table, SIGNAL(cellClicked(int, int)), this, SLOT(OnTop20TableCellClicked(int, int)));
GT_ASSERT(rc);
rc = connect(m_pTop20Table, SIGNAL(cellDoubleClicked(int, int)), this, SLOT(OnTop20TableCellDoubleClicked(int, int)));
GT_ASSERT(rc);
m_pSummaryTable->selectRow(0);
return true;
}
AboutDlg::AboutDlg(QWidget *parent) : QDialog(parent) {
this->setWindowTitle(tr("About QtRptDesigner"));
QSize MaxSize(420, 450);
QSize MinSize(420, 450);
setMaximumSize(MaxSize);
setMinimumSize(MinSize);
QRect gry = geometry();
gry.moveCenter(qApp->desktop()->availableGeometry().center());
setGeometry(gry);
QPushButton *btnOK = new QPushButton(this);
btnOK->setText("OK");
QObject::connect(btnOK, SIGNAL(clicked()), this, SLOT(close()));
QLabel *labImg = new QLabel(this);
labImg->setPixmap(QPixmap(":/new/prefix1/images/Logo128.png"));
QString lbl1 = "<h2><b><p style='color:#0033FF'>"+QApplication::applicationName()+"</p></b></h2>"+
tr("Version: ")+QApplication::applicationVersion()+"<br>"+
tr("Programmer: Aleksey Osipov")+"<br>"+
"WebSite: <a href='http://www.aliks-os.tk'>http://www.aliks-os.tk</a>"+"<br>"+
"E-mail: [email protected]"+"<br>"+
"<a href='http://www.qtrpt.tk'>http://www.qtrpt.tk</a>"+"<br>"+
tr("2012-2015 years")+"<br><br>";
QString lbl2 = "<b>"+tr("Thanks for donation:")+"</b>"+
"<ul>"+
"<li>"+tr("Sailendram")+"</li>"+
"</ul>"+
"<b>"+tr("Thanks for project developing:")+"</b>"+
"<ul>"+
"<li>"+tr("Lukas Lalinsky for DBmodel")+"</li>"+
"<li>"+tr("Norbert Schlia for help in developing")+"</li>"+
"<li>"+tr("Muhamad Bashir Al-Noimi for Arabic translation")+"</li>"+
"<li>"+tr("Luis Brochado for Portuguese translation")+"</li>"+
"<li>"+tr("Li Wei for Chinese translation")+"</li>"+
"<li>"+tr("Laurent Guilbert for French translation")+"</li>"+
"<li>"+tr("David Heremans for Dutch translation")+"</li>"+
"<li>"+tr("Mirko Marx for German translation")+"</li>"+
"<li>"+tr("Manuel Soriano for Spanish translation")+"</li>"+
"</ul>";
QLabel *lab1 = new QLabel(lbl1, this);
QObject::connect(lab1, SIGNAL(linkActivated(const QString)), this, SLOT(openLink(const QString)));
QLabel *lab2 = new QLabel(lbl2, this);
QHBoxLayout *hLayout2 = new QHBoxLayout;
hLayout2->addWidget(labImg);
hLayout2->addWidget(lab1);
hLayout2->addStretch();
QHBoxLayout *hLayout1 = new QHBoxLayout;
//hLayout->addSpacerItem(spacer1);
hLayout1->addStretch();
hLayout1->addWidget(btnOK);
//hLayout->addSpacerItem(spacer2);
hLayout1->addStretch();
QVBoxLayout *vLayout = new QVBoxLayout;
//vLayout->addSpacerItem(spacer3);
vLayout->addLayout(hLayout2);
vLayout->addWidget(lab2);
vLayout->addStretch();
vLayout->addLayout(hLayout1);
this->setLayout(vLayout);
}
/// @brief
/// @todo Doc me!
/// @return
GeometryType type() const
{
return geometry().type();
}
int main(int argc,char *argv[])
{
int irw,isp,ispp[2],status[6],mnkv;
int bs[4],Ns,nmx,nkv,nmr,ncy,ninv;
double kappa,m0,dm,mu0,mu,res,mres;
double sqne,sqnp[2];
complex_dble lnw1[2],lnr,dr,drmx;
solver_parms_t sp;
mrw_masses_t ms;
MPI_Init(&argc,&argv);
MPI_Comm_rank(MPI_COMM_WORLD,&my_rank);
if (my_rank==0)
{
flog=freopen("check2.log","w",stdout);
fin=freopen("check2.in","r",stdin);
printf("\n");
printf("Direct check of mrw2\n");
printf("----------------------\n\n");
printf("%dx%dx%dx%d lattice, ",NPROC0*L0,NPROC1*L1,NPROC2*L2,NPROC3*L3);
printf("%dx%dx%dx%d process grid, ",NPROC0,NPROC1,NPROC2,NPROC3);
printf("%dx%dx%dx%d local lattice\n\n",L0,L1,L2,L3);
}
mnkv=0;
mres=0.0;
for (isp=0;isp<3;isp++)
{
read_solver_parms(isp);
sp=solver_parms(isp);
if (sp.res>mres)
mres=sp.res;
if (sp.nkv>mnkv)
mnkv=sp.nkv;
}
read_bc_parms();
if (my_rank==0)
{
find_section("SAP");
read_line("bs","%d %d %d %d",bs,bs+1,bs+2,bs+3);
}
MPI_Bcast(bs,4,MPI_INT,0,MPI_COMM_WORLD);
set_sap_parms(bs,0,1,1);
if (my_rank==0)
{
find_section("Deflation subspace");
read_line("bs","%d %d %d %d",bs,bs+1,bs+2,bs+3);
read_line("Ns","%d",&Ns);
}
MPI_Bcast(bs,4,MPI_INT,0,MPI_COMM_WORLD);
MPI_Bcast(&Ns,1,MPI_INT,0,MPI_COMM_WORLD);
set_dfl_parms(bs,Ns);
if (my_rank==0)
{
find_section("Deflation subspace generation");
read_line("kappa","%lf",&kappa);
read_line("mu","%lf",&mu);
read_line("ninv","%d",&ninv);
read_line("nmr","%d",&nmr);
read_line("ncy","%d",&ncy);
}
MPI_Bcast(&kappa,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
MPI_Bcast(&mu,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
MPI_Bcast(&ninv,1,MPI_INT,0,MPI_COMM_WORLD);
MPI_Bcast(&nmr,1,MPI_INT,0,MPI_COMM_WORLD);
MPI_Bcast(&ncy,1,MPI_INT,0,MPI_COMM_WORLD);
set_dfl_gen_parms(kappa,mu,ninv,nmr,ncy);
if (my_rank==0)
{
find_section("Deflation projection");
read_line("nkv","%d",&nkv);
read_line("nmx","%d",&nmx);
read_line("res","%lf",&res);
fclose(fin);
}
MPI_Bcast(&nkv,1,MPI_INT,0,MPI_COMM_WORLD);
MPI_Bcast(&nmx,1,MPI_INT,0,MPI_COMM_WORLD);
MPI_Bcast(&res,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
set_dfl_pro_parms(nkv,nmx,res);
set_lat_parms(6.0,1.0,0,NULL,1.234);
print_solver_parms(status,status+1);
print_sap_parms(0);
print_dfl_parms(0);
start_ranlux(0,1245);
geometry();
mnkv=2*mnkv+2;
if (mnkv<(Ns+2))
mnkv=Ns+2;
if (mnkv<5)
mnkv=5;
alloc_ws(mnkv);
alloc_wsd(7);
alloc_wv(2*nkv+2);
alloc_wvd(4);
drmx.re=0.0;
drmx.im=0.0;
for (irw=0;irw<3;irw++)
{
dm=1.0e-2;
for (isp=0;isp<3;isp++)
{
ispp[0]=isp;
ispp[1]=isp;
if (isp==0)
{
m0=1.0877;
mu0=0.1;
}
else if (isp==1)
{
m0=0.0877;
mu0=0.01;
}
else
{
m0=-0.0123;
mu0=0.001;
}
random_ud();
if (isp==2)
{
dfl_modes(status);
error_root(status[0]<0,1,"main [check2.c]",
"dfl_modes failed");
}
if (irw==0)
{
ms.m1=m0;
ms.d1=dm;
ms.mu1=mu0;
ms.m2=m0;
ms.d2=dm;
ms.mu2=mu0;
lnr=mrw2(ms,0,ispp,lnw1,sqnp,&sqne,status);
dr.re=fabs(lnw1[0].re-lnw1[1].re);
dr.im=fabs(lnw1[0].im-lnw1[1].im);
lnr=mrw2(ms,1,ispp,lnw1,sqnp,&sqne,status);
dr.re+=fabs(lnr.re-(2.0*mu0*dm+dm*dm)*sqnp[0]);
dr.re+=fabs(lnw1[0].re-lnw1[1].re);
dr.re+=fabs(sqnp[0]-sqnp[1]);
dr.im+=fabs(lnr.im);
dr.im+=fabs(lnw1[0].im-lnw1[1].im);
}
else if (irw==1)
{
ms.m1=m0;
ms.d1=dm;
ms.mu1=mu0;
ms.m2=m0;
ms.d2=-dm;
ms.mu2=mu0;
lnr=mrw2(ms,0,ispp,lnw1,sqnp,&sqne,status);
dr.re=fabs(lnw1[0].re+lnw1[1].re);
dr.im=fabs(lnw1[0].im+lnw1[1].im);
lnr=mrw2(ms,1,ispp,lnw1,sqnp,&sqne,status);
dr.re+=fabs(lnr.re+dm*dm*sqnp[0]);
dr.re+=fabs(lnw1[0].re+lnw1[1].re);
dr.re+=fabs(sqnp[0]-sqnp[1]);
dr.im+=fabs(lnr.im-2.0*lnw1[0].im);
dr.im+=fabs(lnw1[0].im+lnw1[1].im);
}
else
{
ms.m1=m0;
ms.d1=dm;
ms.mu1=mu0;
ms.m2=m0+dm;
ms.d2=-dm;
ms.mu2=mu0;
lnr=mrw2(ms,0,ispp,lnw1,sqnp,&sqne,status);
dr.re=fabs(lnr.re);
dr.im=fabs(lnr.im);
}
if (dr.re>drmx.re)
drmx.re=dr.re;
if (dr.im>drmx.im)
drmx.im=dr.im;
if (my_rank==0)
{
if (irw==0)
printf("mrw2(d2=d1): ");
else if (irw==1)
printf("mrw2(d2=-d1): ");
else
printf("mrw2(m2=m1+d1,d2=-d1): ");
if ((isp==0)||(isp==1))
printf("status = %d\n",status[0]);
else if (isp==2)
printf("status = (%d,%d,%d)\n",
status[0],status[1],status[2]);
printf("diff = %.1e + i%.1e\n\n",dr.re,dr.im);
}
error_chk();
}
}
if (my_rank==0)
{
printf("\n");
printf("max diff = %.1e + i%.1e\n",drmx.re,drmx.im);
printf("(should be smaller than %.1e)\n\n",mres*sqrt((double)(VOLUME*NPROC*24)));
fclose(flog);
}
MPI_Finalize();
exit(0);
}
/// @brief
/// @todo Doc me!
/// @param
/// @return
FieldVector<ctype, 3> integrationOuterNormal(const FieldVector<ctype, 2>& unused) const
{
FieldVector<ctype, 3> n = pgrid_->face_normals_[faces_of_cell_[subindex_]];
return n*=geometry().integrationElement(unused);
}
int main(int argc, char **argv) {
int c, mu, nu, status;
int i, j, ncon=-1, ir, is, ic, id;
int filename_set = 0;
int x0, x1, x2, x3, ix, iix;
int y0, y1, y2, y3, iy, iiy;
int start_valuet=0, start_valuex=0, start_valuey=0;
int num_threads=1, threadid, nthreads;
int seed, seed_set=0;
double diff1, diff2;
/* double *chi=NULL, *psi=NULL; */
double plaq=0., pl_ts, pl_xs, pl_global;
double *gauge_field_smeared = NULL;
double s[18], t[18], u[18], pl_loc;
double spinor1[24], spinor2[24];
double *pl_gather=NULL;
double dtmp;
complex prod, w, w2;
int verbose = 0;
char filename[200];
char file1[200];
char file2[200];
FILE *ofs=NULL;
double norm, norm2;
fermion_propagator_type *prop=NULL, prop2=NULL, seq_prop=NULL, seq_prop2=NULL, prop_aux=NULL, prop_aux2=NULL;
int idx, eoflag, shift;
float *buffer = NULL;
unsigned int VOL3;
size_t items, bytes;
#ifdef MPI
MPI_Init(&argc, &argv);
#endif
while ((c = getopt(argc, argv, "h?vf:N:c:C:t:s:")) != -1) {
switch (c) {
case 'v':
verbose = 1;
break;
case 'f':
strcpy(filename, optarg);
filename_set=1;
break;
case 'N':
ncon = atoi(optarg);
break;
case 'c':
strcpy(file1, optarg);
break;
case 'C':
strcpy(file2, optarg);
break;
case 't':
num_threads = atoi(optarg);
break;
case 's':
seed = atoi(optarg);
fprintf(stdout, "# [] use seed value %d\n", seed);
seed_set = 1;
break;
case 'h':
case '?':
default:
usage();
break;
}
}
/* set the default values */
if(filename_set==0) strcpy(filename, "cvc.input");
if(g_cart_id==0) fprintf(stdout, "# Reading input from file %s\n", filename);
read_input_parser(filename);
/* some checks on the input data */
if((T_global == 0) || (LX==0) || (LY==0) || (LZ==0)) {
if(g_proc_id==0) fprintf(stdout, "T and L's must be set\n");
usage();
}
if(g_kappa == 0.) {
if(g_proc_id==0) fprintf(stdout, "kappa should be > 0.n");
usage();
}
/* initialize MPI parameters */
mpi_init(argc, argv);
/* initialize T etc. */
fprintf(stdout, "# [%2d] parameters:\n"\
"# [%2d] T_global = %3d\n"\
"# [%2d] T = %3d\n"\
"# [%2d] Tstart = %3d\n"\
"# [%2d] LX_global = %3d\n"\
"# [%2d] LX = %3d\n"\
"# [%2d] LXstart = %3d\n"\
"# [%2d] LY_global = %3d\n"\
"# [%2d] LY = %3d\n"\
"# [%2d] LYstart = %3d\n",\
g_cart_id, g_cart_id, T_global, g_cart_id, T, g_cart_id, Tstart,
g_cart_id, LX_global, g_cart_id, LX, g_cart_id, LXstart,
g_cart_id, LY_global, g_cart_id, LY, g_cart_id, LYstart);
if(init_geometry() != 0) {
fprintf(stderr, "ERROR from init_geometry\n");
exit(101);
}
geometry();
if(init_geometry_5d() != 0) {
fprintf(stderr, "ERROR from init_geometry_5d\n");
exit(102);
}
geometry_5d();
VOL3 = LX*LY*LZ;
/* read the gauge field */
alloc_gauge_field(&g_gauge_field, VOLUMEPLUSRAND);
sprintf(filename, "%s.%.4d", gaugefilename_prefix, Nconf);
if(g_cart_id==0) fprintf(stdout, "# reading gauge field from file %s\n", filename);
if(strcmp(gaugefilename_prefix, "identity")==0) {
status = unit_gauge_field(g_gauge_field, VOLUME);
} else {
// status = read_nersc_gauge_field_3x3(g_gauge_field, filename, &plaq);
// status = read_ildg_nersc_gauge_field(g_gauge_field, filename);
status = read_lime_gauge_field_doubleprec(filename);
// status = read_nersc_gauge_field(g_gauge_field, filename, &plaq);
// status = 0;
}
if(status != 0) {
fprintf(stderr, "[apply_Dtm] Error, could not read gauge field\n");
exit(11);
}
xchange_gauge();
// measure the plaquette
if(g_cart_id==0) fprintf(stdout, "# read plaquette value 1st field: %25.16e\n", plaq);
plaquette(&plaq);
if(g_cart_id==0) fprintf(stdout, "# measured plaquette value 1st field: %25.16e\n", plaq);
g_kappa5d = 0.5 / (5. + g_m0);
fprintf(stdout, "# [] g_kappa5d = %e\n", g_kappa5d);
no_fields=4;
g_spinor_field = (double**)calloc(no_fields, sizeof(double*));
for(i=0; i<no_fields; i++) alloc_spinor_field(&g_spinor_field[i], L5*VOLUMEPLUSRAND);
/*
items = VOL3 * 288;
bytes = items * sizeof(float);
if( (buffer = (float*)malloc( bytes ) ) == NULL ) {
fprintf(stderr, "[] Error, could not allocate buffer\n");
exit(20);
}
*/
/****************************************
* read read the spinor fields
****************************************/
/*
prop = create_fp_field(VOL3);
create_fp(&prop2);
create_fp(&prop_aux);
create_fp(&prop_aux2);
create_fp(&seq_prop);
create_fp(&seq_prop2);
*/
#ifdef MPI
if(!seed_set) { seed = g_seed; }
srand(seed+g_cart_id);
for(ix=0;ix<VOLUME*L5;ix++) {
for(i=0;i<24;i++) {
spinor1[i] = 2* (double)rand() / (double)RAND_MAX - 1.;
}
_fv_eq_fv(g_spinor_field[0]+_GSI(ix), spinor1 );
}
for(i=0;i<g_nproc;i++) {
if(g_cart_id==i) {
if(i==0) ofs = fopen("source", "w");
else ofs = fopen("source", "a");
for(is=0;is<L5;is++) {
for(x0=0;x0<T; x0++) {
for(x1=0;x1<LX; x1++) {
for(x2=0;x2<LX; x2++) {
for(x3=0;x3<LX; x3++) {
iix = is*VOLUME*g_nproc + (((x0+g_proc_coords[0]*T)*LX*g_nproc_x+ x1+g_proc_coords[1]*LX )*LY*g_nproc_y + x2+g_proc_coords[2]*LY )*LZ*g_nproc_z + x3+g_proc_coords[3]*LZ;
ix = g_ipt_5d[is][x0][x1][x2][x3];
for(c=0;c<24;c++) {
fprintf(ofs, "%8d%8d%3d%25.16e\n", iix, ix, c, g_spinor_field[0][_GSI(ix)+c]);
}
}}}}
}
fclose(ofs);
}
#ifdef MPI
MPI_Barrier(g_cart_grid);
#endif
}
#else
ofs = fopen("source", "r");
for(ix=0;ix<24*VOLUME*L5;ix++) {
fscanf(ofs, "%d%d%d%lf", &x1,&x2,&x3, &dtmp);
g_spinor_field[0][_GSI(x1)+x3] = dtmp;
}
fclose(ofs);
#endif
xchange_field_5d(g_spinor_field[0]);
Q_DW_Wilson_dag_phi(g_spinor_field[1], g_spinor_field[0]);
xchange_field_5d(g_spinor_field[1]);
Q_DW_Wilson_phi(g_spinor_field[2], g_spinor_field[1]);
sprintf(filename, "prop_%.2d.%.2d", g_nproc, g_cart_id);
ofs = fopen(filename, "w");
printf_spinor_field_5d(g_spinor_field[2], ofs);
fclose(ofs);
// for(ix=0;ix<VOLUME*L5;ix++) {
// for(i=0;i<24;i++) {
// spinor1[i] = 2* (double)rand() / (double)RAND_MAX - 1.;
// }
// _fv_eq_fv(g_spinor_field[1]+_GSI(ix), spinor1 );
// }
/*
xchange_field_5d(g_spinor_field[0]);
sprintf(filename, "spinor.%.2d", g_cart_id);
ofs = fopen(filename, "w");
printf_spinor_field_5d(g_spinor_field[0], ofs);
fclose(ofs);
*/
/*
// 2 = D 0
Q_DW_Wilson_phi(g_spinor_field[2], g_spinor_field[0]);
// 3 = D^dagger 1
Q_DW_Wilson_dag_phi(g_spinor_field[3], g_spinor_field[1]);
// <1, 2> = <1, D 0 >
spinor_scalar_product_co(&w, g_spinor_field[1], g_spinor_field[2], VOLUME*L5);
// <3, 0> = < D^dagger 1, 0 >
spinor_scalar_product_co(&w2, g_spinor_field[3], g_spinor_field[0], VOLUME*L5);
fprintf(stdout, "# [] w = %e + %e*1.i\n", w.re, w.im);
fprintf(stdout, "# [] w2 = %e + %e*1.i\n", w2.re, w2.im);
fprintf(stdout, "# [] abs difference = %e \n", sqrt(_SQR(w2.re-w.re)+_SQR(w2.im-w.im)) );
*/
/*
for(i=0;i<12;i++) {
fprintf(stdout, "s1[%2d] <- %25.16e + %25.16e*1.i\n", i+1, spinor1[2*i], spinor1[2*i+1]);
}
for(i=0;i<24;i++) {
spinor2[i] = 2* (double)rand() / (double)RAND_MAX - 1.;
}
for(i=0;i<12;i++) {
fprintf(stdout, "s2[%2d] <- %25.16e + %25.16e*1.i\n", i+1, spinor2[2*i], spinor2[2*i+1]);
}
_fv_mi_eq_PRe_fv(spinor2, spinor1);
for(i=0;i<12;i++) {
fprintf(stdout, "s3[%2d] <- %25.16e + %25.16e*1.i\n", i+1, spinor2[2*i], spinor2[2*i+1]);
}
*/
/*
ofs = fopen("dw_spinor", "w");
Q_DW_Wilson_phi(g_spinor_field[1], g_spinor_field[0]);
printf_spinor_field(g_spinor_field[1], ofs);
fclose(ofs);
g_kappa = g_kappa5d;
ofs = fopen("wilson_spinor", "w");
Q_Wilson_phi(g_spinor_field[2], g_spinor_field[0]);
printf_spinor_field(g_spinor_field[2], ofs);
fclose(ofs);
*/
#ifdef _UNDEF
/*******************************************************************
* propagators
*******************************************************************/
// for(i=0; i<12;i++)
for(i=0; i<1;i++)
{
//sprintf(file1, "source.%.4d.t00x00y00z00.%.2d.inverted", Nconf, i);
sprintf(file1, "/home/mpetschlies/quda-0.3.2/tests/prop");
if(g_cart_id==0) fprintf(stdout, "# Reading prop. from file %s\n", file1);
fflush(stdout);
//if( read_lime_spinor(g_spinor_field[0], file1, 0) != 0 ) {
ofs = fopen(file1, "rb");
if( fread(g_spinor_field[0], sizeof(double), 24*L5*VOLUME, ofs) != 24*L5*VOLUME) {
fprintf(stderr, "Error, could not read proper amount of data from file %s\n", file1);
exit(100);
}
fclose(ofs);
for(ix=0;ix<VOLUME*L5;ix++) {
_fv_ti_eq_re(g_spinor_field[0]+_GSI(ix), 2.*g_kappa5d);
}
/*
if( (ofs = fopen("prop_full", "w")) == NULL ) exit(22);
for(ix=0;ix<L5;ix++) {
fprintf(ofs, "# [] s = %d\n", ix);
printf_spinor_field(g_spinor_field[0]+_GSI(ix*VOLUME), ofs);
}
fclose(ofs);
*/
// reorder, multiply with g2
for(is=0,iix=0; is<L5; is++) {
for(ix=0; ix<VOLUME; ix++) {
iiy = lexic2eot_5d (is, ix);
_fv_eq_fv(spinor1, g_spinor_field[0]+_GSI(iiy));
_fv_eq_gamma_ti_fv(g_spinor_field[1]+_GSI(iix), 2, spinor1 );
iix++;
}}
Q_DW_Wilson_phi(g_spinor_field[2], g_spinor_field[1]);
// Q_DW_Wilson_dag_phi(g_spinor_field[2], g_spinor_field[1]);
fprintf(stdout, "# [] finished application of Dirac operator\n");
fflush(stdout);
// reorder, multiply with g2
for(is=0, iix=0;is<L5;is++) {
for(ix=0; ix<VOLUME; ix++) {
iiy = lexic2eot_5d(is, ix);
_fv_eq_fv(spinor1, g_spinor_field[2]+_GSI(iix));
_fv_eq_gamma_ti_fv(g_spinor_field[1]+_GSI(iiy), 2, spinor1 );
iix++;
}}
if( (ofs = fopen("my_out", "w")) == NULL ) exit(23);
for(ix=0;ix<L5;ix++) {
fprintf(ofs, "# [] s = %d\n", ix);
printf_spinor_field(g_spinor_field[1]+_GSI(ix*VOLUME), ofs);
}
fclose(ofs);
sprintf(file1, "/home/mpetschlies/quda-0.3.2/tests/source");
if(g_cart_id==0) fprintf(stdout, "# Reading prop. from file %s\n", file1);
fflush(stdout);
//if( read_lime_spinor(g_spinor_field[0], file1, 0) != 0 ) {
ofs = fopen(file1, "rb");
if( fread(g_spinor_field[2], sizeof(double), 24*L5*VOLUME, ofs) != 24*L5*VOLUME) {
fprintf(stderr, "Error, could not read proper amount of data from file %s\n", file1);
exit(100);
}
fclose(ofs);
/*
if( (ofs = fopen("v_out", "w")) == NULL ) exit(23);
for(ix=0;ix<L5;ix++) {
fprintf(ofs, "# [] s = %d\n", ix);
printf_spinor_field(g_spinor_field[2]+_GSI(ix*VOLUME), ofs);
}
fclose(ofs);
*/
spinor_scalar_product_re(&norm2, g_spinor_field[2], g_spinor_field[2], VOLUME*L5);
for(ix=0;ix<VOLUME*L5;ix++) {
_fv_mi_eq_fv(g_spinor_field[1]+_GSI(ix), g_spinor_field[2]+_GSI(ix));
}
spinor_scalar_product_re(&norm, g_spinor_field[1], g_spinor_field[1], VOLUME*L5);
fprintf(stdout, "\n# [] absolut residuum squared: %e; relative residuum %e\n", norm, sqrt(norm/norm2) );
} // of loop on spin color indices
#endif
/***********************************************
* free the allocated memory, finalize
***********************************************/
free(g_gauge_field);
free_geometry();
if(gauge_field_smeared != NULL) free(gauge_field_smeared);
if(g_spinor_field != NULL) {
for(i=0; i<no_fields; i++) free(g_spinor_field[i]);
free(g_spinor_field);
}
free(buffer);
free_fp_field(&prop);
free_fp(&prop2);
free_fp(&prop_aux);
free_fp(&prop_aux2);
free_fp(&seq_prop);
free_fp(&seq_prop2);
g_the_time = time(NULL);
fprintf(stdout, "# [] %s# [] end fo run\n", ctime(&g_the_time));
fflush(stdout);
fprintf(stderr, "# [] %s# [] end fo run\n", ctime(&g_the_time));
fflush(stderr);
#ifdef MPI
MPI_Finalize();
#endif
return(0);
}
int main(int argc,char *argv[])
{
int j,j_max,k,k_max = 1;
#ifdef HAVE_LIBLEMON
paramsXlfInfo *xlfInfo;
#endif
int status = 0;
static double t1,t2,dt,sdt,dts,qdt,sqdt;
double antioptaway=0.0;
#ifdef MPI
static double dt2;
DUM_DERI = 6;
DUM_SOLVER = DUM_DERI+2;
DUM_MATRIX = DUM_SOLVER+6;
NO_OF_SPINORFIELDS = DUM_MATRIX+2;
# ifdef OMP
int mpi_thread_provided;
MPI_Init_thread(&argc, &argv, MPI_THREAD_SERIALIZED, &mpi_thread_provided);
# else
MPI_Init(&argc, &argv);
# endif
MPI_Comm_rank(MPI_COMM_WORLD, &g_proc_id);
#else
g_proc_id = 0;
#endif
g_rgi_C1 = 1.;
/* Read the input file */
if((status = read_input("benchmark.input")) != 0) {
fprintf(stderr, "Could not find input file: benchmark.input\nAborting...\n");
exit(-1);
}
#ifdef OMP
if(omp_num_threads > 0)
{
omp_set_num_threads(omp_num_threads);
}
else {
if( g_proc_id == 0 )
printf("# No value provided for OmpNumThreads, running in single-threaded mode!\n");
omp_num_threads = 1;
omp_set_num_threads(omp_num_threads);
}
init_omp_accumulators(omp_num_threads);
#endif
tmlqcd_mpi_init(argc, argv);
if(g_proc_id==0) {
#ifdef SSE
printf("# The code was compiled with SSE instructions\n");
#endif
#ifdef SSE2
printf("# The code was compiled with SSE2 instructions\n");
#endif
#ifdef SSE3
printf("# The code was compiled with SSE3 instructions\n");
#endif
#ifdef P4
printf("# The code was compiled for Pentium4\n");
#endif
#ifdef OPTERON
printf("# The code was compiled for AMD Opteron\n");
#endif
#ifdef _GAUGE_COPY
printf("# The code was compiled with -D_GAUGE_COPY\n");
#endif
#ifdef BGL
printf("# The code was compiled for Blue Gene/L\n");
#endif
#ifdef BGP
printf("# The code was compiled for Blue Gene/P\n");
#endif
#ifdef _USE_HALFSPINOR
printf("# The code was compiled with -D_USE_HALFSPINOR\n");
#endif
#ifdef _USE_SHMEM
printf("# The code was compiled with -D_USE_SHMEM\n");
# ifdef _PERSISTENT
printf("# The code was compiled for persistent MPI calls (halfspinor only)\n");
# endif
#endif
#ifdef MPI
# ifdef _NON_BLOCKING
printf("# The code was compiled for non-blocking MPI calls (spinor and gauge)\n");
# endif
#endif
printf("\n");
fflush(stdout);
}
#ifdef _GAUGE_COPY
init_gauge_field(VOLUMEPLUSRAND + g_dbw2rand, 1);
#else
init_gauge_field(VOLUMEPLUSRAND + g_dbw2rand, 0);
#endif
init_geometry_indices(VOLUMEPLUSRAND + g_dbw2rand);
if(even_odd_flag) {
j = init_spinor_field(VOLUMEPLUSRAND/2, 2*k_max+1);
}
else {
j = init_spinor_field(VOLUMEPLUSRAND, 2*k_max);
}
if ( j!= 0) {
fprintf(stderr, "Not enough memory for spinor fields! Aborting...\n");
exit(0);
}
j = init_moment_field(VOLUME, VOLUMEPLUSRAND + g_dbw2rand);
if ( j!= 0) {
fprintf(stderr, "Not enough memory for moment fields! Aborting...\n");
exit(0);
}
if(g_proc_id == 0) {
fprintf(stdout,"# The number of processes is %d \n",g_nproc);
printf("# The lattice size is %d x %d x %d x %d\n",
(int)(T*g_nproc_t), (int)(LX*g_nproc_x), (int)(LY*g_nproc_y), (int)(g_nproc_z*LZ));
printf("# The local lattice size is %d x %d x %d x %d\n",
(int)(T), (int)(LX), (int)(LY),(int) LZ);
if(even_odd_flag) {
printf("# benchmarking the even/odd preconditioned Dirac operator\n");
}
else {
printf("# benchmarking the standard Dirac operator\n");
}
fflush(stdout);
}
/* define the geometry */
geometry();
/* define the boundary conditions for the fermion fields */
boundary(g_kappa);
#ifdef _USE_HALFSPINOR
j = init_dirac_halfspinor();
if ( j!= 0) {
fprintf(stderr, "Not enough memory for halfspinor fields! Aborting...\n");
exit(0);
}
if(g_sloppy_precision_flag == 1) {
g_sloppy_precision = 1;
j = init_dirac_halfspinor32();
if ( j!= 0) {
fprintf(stderr, "Not enough memory for 32-Bit halfspinor fields! Aborting...\n");
exit(0);
}
}
# if (defined _PERSISTENT)
init_xchange_halffield();
# endif
#endif
status = check_geometry();
if (status != 0) {
fprintf(stderr, "Checking of geometry failed. Unable to proceed.\nAborting....\n");
exit(1);
}
#if (defined MPI && !(defined _USE_SHMEM))
check_xchange();
#endif
start_ranlux(1, 123456);
random_gauge_field(reproduce_randomnumber_flag);
#ifdef MPI
/*For parallelization: exchange the gaugefield */
xchange_gauge(g_gauge_field);
#endif
if(even_odd_flag) {
/*initialize the pseudo-fermion fields*/
j_max=2048;
sdt=0.;
for (k = 0; k < k_max; k++) {
random_spinor_field(g_spinor_field[k], VOLUME/2, 0);
}
while(sdt < 30.) {
#ifdef MPI
MPI_Barrier(MPI_COMM_WORLD);
#endif
t1 = gettime();
antioptaway=0.0;
for (j=0;j<j_max;j++) {
for (k=0;k<k_max;k++) {
Hopping_Matrix(0, g_spinor_field[k+k_max], g_spinor_field[k]);
Hopping_Matrix(1, g_spinor_field[2*k_max], g_spinor_field[k+k_max]);
antioptaway+=creal(g_spinor_field[2*k_max][0].s0.c0);
}
}
t2 = gettime();
dt = t2-t1;
#ifdef MPI
MPI_Allreduce (&dt, &sdt, 1, MPI_DOUBLE, MPI_SUM, MPI_COMM_WORLD);
#else
sdt = dt;
#endif
qdt=dt*dt;
#ifdef MPI
MPI_Allreduce (&qdt, &sqdt, 1, MPI_DOUBLE, MPI_SUM, MPI_COMM_WORLD);
#else
sqdt = qdt;
#endif
sdt=sdt/((double)g_nproc);
sqdt=sqrt(sqdt/g_nproc-sdt*sdt);
j_max*=2;
}
j_max=j_max/2;
dts=dt;
sdt=1.0e6f*sdt/((double)(k_max*j_max*(VOLUME)));
sqdt=1.0e6f*sqdt/((double)(k_max*j_max*(VOLUME)));
if(g_proc_id==0) {
printf("# The following result is just to make sure that the calculation is not optimized away: %e\n", antioptaway);
printf("# Total compute time %e sec, variance of the time %e sec. (%d iterations).\n", sdt, sqdt, j_max);
printf("# Communication switched on:\n# (%d Mflops [%d bit arithmetic])\n", (int)(1608.0f/sdt),(int)sizeof(spinor)/3);
#ifdef OMP
printf("# Mflops per OpenMP thread ~ %d\n",(int)(1608.0f/(omp_num_threads*sdt)));
#endif
printf("\n");
fflush(stdout);
}
#ifdef MPI
/* isolated computation */
t1 = gettime();
antioptaway=0.0;
for (j=0;j<j_max;j++) {
for (k=0;k<k_max;k++) {
Hopping_Matrix_nocom(0, g_spinor_field[k+k_max], g_spinor_field[k]);
Hopping_Matrix_nocom(1, g_spinor_field[2*k_max], g_spinor_field[k+k_max]);
antioptaway += creal(g_spinor_field[2*k_max][0].s0.c0);
}
}
t2 = gettime();
dt2 = t2-t1;
/* compute the bandwidth */
dt=dts-dt2;
MPI_Allreduce (&dt, &sdt, 1, MPI_DOUBLE, MPI_SUM, MPI_COMM_WORLD);
sdt=sdt/((double)g_nproc);
MPI_Allreduce (&dt2, &dt, 1, MPI_DOUBLE, MPI_SUM, MPI_COMM_WORLD);
dt=dt/((double)g_nproc);
dt=1.0e6f*dt/((double)(k_max*j_max*(VOLUME)));
if(g_proc_id==0) {
printf("# The following result is printed just to make sure that the calculation is not optimized away: %e\n",antioptaway);
printf("# Communication switched off: \n# (%d Mflops [%d bit arithmetic])\n", (int)(1608.0f/dt),(int)sizeof(spinor)/3);
#ifdef OMP
printf("# Mflops per OpenMP thread ~ %d\n",(int)(1608.0f/(omp_num_threads*dt)));
#endif
printf("\n");
fflush(stdout);
}
sdt=sdt/((double)k_max);
sdt=sdt/((double)j_max);
sdt=sdt/((double)(2*SLICE));
if(g_proc_id==0) {
printf("# The size of the package is %d bytes.\n",(SLICE)*192);
#ifdef _USE_HALFSPINOR
printf("# The bandwidth is %5.2f + %5.2f MB/sec\n", 192./sdt/1024/1024, 192./sdt/1024./1024);
#else
printf("# The bandwidth is %5.2f + %5.2f MB/sec\n", 2.*192./sdt/1024/1024, 2.*192./sdt/1024./1024);
#endif
}
#endif
fflush(stdout);
}
else {
/* the non even/odd case now */
/*initialize the pseudo-fermion fields*/
j_max=1;
sdt=0.;
for (k=0;k<k_max;k++) {
random_spinor_field(g_spinor_field[k], VOLUME, 0);
}
while(sdt < 3.) {
#ifdef MPI
MPI_Barrier(MPI_COMM_WORLD);
#endif
t1 = gettime();
for (j=0;j<j_max;j++) {
for (k=0;k<k_max;k++) {
D_psi(g_spinor_field[k+k_max], g_spinor_field[k]);
antioptaway+=creal(g_spinor_field[k+k_max][0].s0.c0);
}
}
t2 = gettime();
dt=t2-t1;
#ifdef MPI
MPI_Allreduce (&dt, &sdt, 1, MPI_DOUBLE, MPI_SUM, MPI_COMM_WORLD);
#else
sdt = dt;
#endif
qdt=dt*dt;
#ifdef MPI
MPI_Allreduce (&qdt, &sqdt, 1, MPI_DOUBLE, MPI_SUM, MPI_COMM_WORLD);
#else
sqdt = qdt;
#endif
sdt=sdt/((double)g_nproc);
sqdt=sqrt(sqdt/g_nproc-sdt*sdt);
j_max*=2;
}
j_max=j_max/2;
dts=dt;
sdt=1.0e6f*sdt/((double)(k_max*j_max*(VOLUME)));
sqdt=1.0e6f*sqdt/((double)(k_max*j_max*(VOLUME)));
if(g_proc_id==0) {
printf("# The following result is just to make sure that the calculation is not optimized away: %e\n", antioptaway);
printf("# Total compute time %e sec, variance of the time %e sec. (%d iterations).\n", sdt, sqdt, j_max);
printf("\n# (%d Mflops [%d bit arithmetic])\n", (int)(1680.0f/sdt),(int)sizeof(spinor)/3);
#ifdef OMP
printf("# Mflops per OpenMP thread ~ %d\n",(int)(1680.0f/(omp_num_threads*sdt)));
#endif
printf("\n");
fflush(stdout);
}
}
#ifdef HAVE_LIBLEMON
if(g_proc_id==0) {
printf("# Performing parallel IO test ...\n");
}
xlfInfo = construct_paramsXlfInfo(0.5, 0);
write_gauge_field( "conf.test", 64, xlfInfo);
free(xlfInfo);
if(g_proc_id==0) {
printf("# done ...\n");
}
#endif
#ifdef MPI
MPI_Finalize();
#endif
#ifdef OMP
free_omp_accumulators();
#endif
free_gauge_field();
free_geometry_indices();
free_spinor_field();
free_moment_field();
return(0);
}
void UBWebPluginWidget::resizeEvent(QResizeEvent *event)
{
Q_UNUSED(event);
mLoadingProgressBar.move(geometry().center() - mLoadingProgressBar.geometry().center());
}
void Message::setSendLines(int lines)
{
msg_send_lines=lines;
setGeometry(geometry());
}
PxRigidActor* createSphereActor( double radius, double density, PxMaterial* mtl )
{
PxSphereGeometry geometry( radius );
return createActor( geometry, density, mtl );
}
int main(int argc,char *argv[])
{
int opt;
PsimagLite::String file("");
SizeType total=0;
RealType offset = 0;
RealType step = 0;
SizeType site3 = 0;
while ((opt = getopt(argc, argv, "f:t:o:i:p:")) != -1) {
switch (opt) {
case 'f':
file=optarg;
break;
case 't':
total = atoi(optarg);
break;
case 'i':
step = atof(optarg);
break;
case 'o':
offset = atof(optarg);
break;
case 'p':
site3 = atoi(optarg);
break;
default: /* '?' */
throw std::runtime_error("Wrong usage\n");
}
}
if (file=="") {
throw std::runtime_error("Wrong usage\n");
}
FreeFermions::InputCheck inputCheck;
InputNgType::Writeable ioWriteable(file,inputCheck);
InputNgType::Readable io(ioWriteable);
GeometryParamsType geometryParams(io);
SizeType electronsUp = GeometryParamsType::readElectrons(io,geometryParams.sites);
PsimagLite::Vector<SizeType>::Type sites;
GeometryParamsType::readVector(sites,file,"TSPSites");
sites.resize(3);
sites[2] = site3;
SizeType dof = 1; // spinless
GeometryLibraryType geometry(geometryParams);
std::cerr<<geometry;
SizeType npthreads = 1;
ConcurrencyType concurrency(&argc,&argv,npthreads);
EngineType engine(geometry.matrix(),dof,true);
PsimagLite::Vector<SizeType>::Type ne(dof,electronsUp); // 8 up and 8 down
bool debug = false;
HilbertStateType gs(engine,ne,debug);
SizeType sigma =0;
OpNormalFactoryType opNormalFactory(engine);
OperatorType& myOp = opNormalFactory(OperatorType::DESTRUCTION,sites[0],sigma);
HilbertStateType phi2 = gs;
myOp.applyTo(phi2);
// FieldType density = scalarProduct(phi2,phi2);
// std::cerr<<"density="<<density<<"\n";
std::cout<<"#site="<<sites[0]<<"\n";
std::cout<<"#site2="<<sites[1]<<"\n";
for (SizeType it = 0; it<total; it++) {
OpDiagonalFactoryType opDiagonalFactory(engine);
RealType time = it * step + offset;
EtoTheIhTimeType eih(time,engine,0);
DiagonalOperatorType& eihOp = opDiagonalFactory(eih);
HilbertStateType phi = gs;
myOp.applyTo(phi);
eihOp.applyTo(phi);
OperatorType& myOp2 = opNormalFactory(OperatorType::DESTRUCTION,sites[1],sigma);
myOp2.applyTo(phi);
std::cout<<time<<" "<<scalarProduct(phi,phi)<<"\n";
}
}
PxRigidActor* createConvexMeshActor( PxConvexMesh* mesh, double density, PxMaterial* mtl )
{
PxConvexMeshGeometry geometry( mesh );
return createActor( geometry, density, mtl );
}
void PreviewWidget::paintEvent(QPaintEvent *event)
{
//#define DEBUG_LAYOUT
#ifdef DEBUG_LAYOUT
{
if (mImage.isNull())
return;
if (mSheetNum < 0)
return;
Sheet *sheet = project->previewSheet(mSheetNum);
if (!sheet)
return;
QPainter painter(this);
QRectF rect = project->printer()->paperRect();
painter.fillRect(rect, Qt::white);
for (int i=0; i< sheet->count(); ++i)
{
QPen pen = painter.pen();
pen.setStyle(Qt::DotLine);
pen.setColor(Qt::darkGray);
painter.setPen(pen);
TransformSpec spec = project->layout()->transformSpec(sheet, i, project->rotation());
painter.drawRect(spec.rect);
QFont font = painter.font();
font.setPixelSize(spec.rect.height() / 2 );
painter.setFont(font);
painter.drawText(spec.rect, Qt::AlignCenter, QString("%1").arg(i+1));
pen.setStyle(Qt::SolidLine);
pen.setColor(Qt::red);
painter.setPen(pen);
switch (spec.rotation)
{
case NoRotate:
painter.drawLine(spec.rect.topLeft(), spec.rect.topRight());
break;
case Rotate90:
painter.drawLine(spec.rect.topRight(), spec.rect.bottomRight());
break;
case Rotate180:
painter.drawLine(spec.rect.bottomLeft(), spec.rect.bottomRight());
break;
case Rotate270:
painter.drawLine(spec.rect.topLeft(), spec.rect.bottomLeft());
break;
}
}
return;
}
#endif
if (mImage.isNull())
return;
QSizeF printerSize = project->printer()->paperRect().size();
Rotation rotation = project->rotation();
if (isLandscape(rotation))
printerSize.transpose();
mScaleFactor = qMin((this->geometry().width() - 2.0 * MARGIN_H) * 1.0 / printerSize.width(),
(this->geometry().height() - 2.0 * MARGIN_V) * 1.0 / printerSize.height());
if (mScaleFactor == 0)
{
mDrawRect = QRect();
return;
}
QSize size = QSize(printerSize.width() * mScaleFactor,
printerSize.height() * mScaleFactor);
mDrawRect = QRect(QPoint(0, 0), size);
mDrawRect.moveCenter(QPoint(0, 0));
QRectF clipRect = mDrawRect;
if (mHints.testFlag(Sheet::HintOnlyLeft))
{
switch (rotation)
{
case NoRotate: clipRect.setBottom(0); break;
case Rotate90: clipRect.setRight(0); break;
case Rotate180: clipRect.setBottom(0); break;
case Rotate270: clipRect.setRight(0); break;
}
}
if (mHints.testFlag(Sheet::HintOnlyRight))
{
switch (rotation)
{
case NoRotate: clipRect.setTop(0); break;
case Rotate90: clipRect.setLeft(0); break;
case Rotate180: clipRect.setTop(0); break;
case Rotate270: clipRect.setLeft(0); break;
}
}
QImage img = mImage.scaled(mDrawRect.size(), Qt::IgnoreAspectRatio, Qt::SmoothTransformation);
// Draw .....................................
QPainter painter(this);
painter.save();
QPoint center = QRect(0, 0, geometry().width(), geometry().height()).center();
painter.translate(center);
if (mHints.testFlag(Sheet::HintSubBooklet))
{
painter.save();
QPoint center = mDrawRect.center();
QRectF imgRect;
clipRect = mDrawRect;
clipRect.setRight(center.x());
clipRect.adjust(-MARGIN_BOOKLET, 0, -MARGIN_BOOKLET, 0);
painter.setClipRect(clipRect);
imgRect = img.rect();
imgRect.setRight(img.rect().center().x());
painter.drawImage(clipRect, img, imgRect);
drawShadow(painter, clipRect);
clipRect = mDrawRect;
clipRect.setLeft(center.x());
clipRect.adjust(MARGIN_BOOKLET, 0, MARGIN_BOOKLET, 0);
painter.setClipRect(clipRect);
imgRect = img.rect();
imgRect.setLeft(img.rect().center().x());
painter.drawImage(clipRect, img, imgRect);
drawShadow(painter, clipRect);
painter.restore();
}
else
{
painter.save();
painter.setClipRect(clipRect);
painter.drawImage(mDrawRect, img);
drawShadow(painter, clipRect);
painter.restore();
}
if (mHints.testFlag(Sheet::HintDrawFold) && !mHints.testFlag(Sheet::HintSubBooklet))
{
QPen pen = painter.pen();
pen.setStyle(Qt::SolidLine);
pen.setColor(Qt::lightGray);
painter.setPen(pen);
if (isLandscape(rotation))
painter.drawLine(0, mDrawRect.top(), 0, mDrawRect.bottom());
else
painter.drawLine(mDrawRect.left(), 0, mDrawRect.right(), 0);
}
painter.restore();
mDrawRect.moveCenter(center);
// Draw current page rect ...................
Sheet *sheet = project->currentSheet();
if (sheet)
{
ProjectPage *curPage = project->currentPage();
if (curPage)
{
painter.save();
QPen pen = painter.pen();
pen.setStyle(Qt::DashLine);
//pen.setColor(QColor(142, 188, 226, 128));
pen.setColor(QColor(105, 101, 98, 70));
painter.setPen(pen);
painter.drawRect(this->pageRect(sheet->indexOfPage(curPage)));
painter.restore();
}
}
//#define DEBUG_CLICK_RECT
#ifdef DEBUG_CLICK_RECT
{
Sheet *sheet = project->currentSheet();
if (sheet)
{
ProjectPage *curPage = project->currentPage();
painter.save();
for (int i=0; i< sheet->count(); ++i)
{
QPen pen = painter.pen();
pen.setStyle(Qt::DotLine);
if (sheet->page(i) == curPage)
pen.setColor(Qt::red);
else
pen.setColor(QColor(142, 188, 226));
painter.setPen(pen);
painter.drawRect(this->pageRect(i));
painter.drawText(this->pageRect(i).translated(10, 10), QString("%1").arg(i));
}
painter.restore();
}
}
#endif
}
void SpanLayoutView::render_content(Canvas &canvas)
{
return impl->render_content(canvas, geometry().content.get_width());
}
int QWidget::qt_metacall(QMetaObject::Call _c, int _id, void **_a)
{
_id = QObject::qt_metacall(_c, _id, _a);
if (_id < 0)
return _id;
if (_c == QMetaObject::InvokeMetaMethod) {
if (_id < 23)
qt_static_metacall(this, _c, _id, _a);
_id -= 23;
}
#ifndef QT_NO_PROPERTIES
else if (_c == QMetaObject::ReadProperty) {
void *_v = _a[0];
switch (_id) {
case 0: *reinterpret_cast< bool*>(_v) = isModal(); break;
case 1: *reinterpret_cast< Qt::WindowModality*>(_v) = windowModality(); break;
case 2: *reinterpret_cast< bool*>(_v) = isEnabled(); break;
case 3: *reinterpret_cast< QRect*>(_v) = geometry(); break;
case 4: *reinterpret_cast< QRect*>(_v) = frameGeometry(); break;
case 5: *reinterpret_cast< QRect*>(_v) = normalGeometry(); break;
case 6: *reinterpret_cast< int*>(_v) = x(); break;
case 7: *reinterpret_cast< int*>(_v) = y(); break;
case 8: *reinterpret_cast< QPoint*>(_v) = pos(); break;
case 9: *reinterpret_cast< QSize*>(_v) = frameSize(); break;
case 10: *reinterpret_cast< QSize*>(_v) = size(); break;
case 11: *reinterpret_cast< int*>(_v) = width(); break;
case 12: *reinterpret_cast< int*>(_v) = height(); break;
case 13: *reinterpret_cast< QRect*>(_v) = rect(); break;
case 14: *reinterpret_cast< QRect*>(_v) = childrenRect(); break;
case 15: *reinterpret_cast< QRegion*>(_v) = childrenRegion(); break;
case 16: *reinterpret_cast< QSizePolicy*>(_v) = sizePolicy(); break;
case 17: *reinterpret_cast< QSize*>(_v) = minimumSize(); break;
case 18: *reinterpret_cast< QSize*>(_v) = maximumSize(); break;
case 19: *reinterpret_cast< int*>(_v) = minimumWidth(); break;
case 20: *reinterpret_cast< int*>(_v) = minimumHeight(); break;
case 21: *reinterpret_cast< int*>(_v) = maximumWidth(); break;
case 22: *reinterpret_cast< int*>(_v) = maximumHeight(); break;
case 23: *reinterpret_cast< QSize*>(_v) = sizeIncrement(); break;
case 24: *reinterpret_cast< QSize*>(_v) = baseSize(); break;
case 25: *reinterpret_cast< QPalette*>(_v) = palette(); break;
case 26: *reinterpret_cast< QFont*>(_v) = font(); break;
case 27: *reinterpret_cast< QCursor*>(_v) = cursor(); break;
case 28: *reinterpret_cast< bool*>(_v) = hasMouseTracking(); break;
case 29: *reinterpret_cast< bool*>(_v) = isActiveWindow(); break;
case 30: *reinterpret_cast< Qt::FocusPolicy*>(_v) = focusPolicy(); break;
case 31: *reinterpret_cast< bool*>(_v) = hasFocus(); break;
case 32: *reinterpret_cast< Qt::ContextMenuPolicy*>(_v) = contextMenuPolicy(); break;
case 33: *reinterpret_cast< bool*>(_v) = updatesEnabled(); break;
case 34: *reinterpret_cast< bool*>(_v) = isVisible(); break;
case 35: *reinterpret_cast< bool*>(_v) = isMinimized(); break;
case 36: *reinterpret_cast< bool*>(_v) = isMaximized(); break;
case 37: *reinterpret_cast< bool*>(_v) = isFullScreen(); break;
case 38: *reinterpret_cast< QSize*>(_v) = sizeHint(); break;
case 39: *reinterpret_cast< QSize*>(_v) = minimumSizeHint(); break;
case 40: *reinterpret_cast< bool*>(_v) = acceptDrops(); break;
case 41: *reinterpret_cast< QString*>(_v) = windowTitle(); break;
case 42: *reinterpret_cast< QIcon*>(_v) = windowIcon(); break;
case 43: *reinterpret_cast< QString*>(_v) = windowIconText(); break;
case 44: *reinterpret_cast< double*>(_v) = windowOpacity(); break;
case 45: *reinterpret_cast< bool*>(_v) = isWindowModified(); break;
case 46: *reinterpret_cast< QString*>(_v) = toolTip(); break;
case 47: *reinterpret_cast< QString*>(_v) = statusTip(); break;
case 48: *reinterpret_cast< QString*>(_v) = whatsThis(); break;
case 49: *reinterpret_cast< QString*>(_v) = accessibleName(); break;
case 50: *reinterpret_cast< QString*>(_v) = accessibleDescription(); break;
case 51: *reinterpret_cast< Qt::LayoutDirection*>(_v) = layoutDirection(); break;
case 52: *reinterpret_cast< bool*>(_v) = autoFillBackground(); break;
case 53: *reinterpret_cast< QString*>(_v) = styleSheet(); break;
case 54: *reinterpret_cast< QLocale*>(_v) = locale(); break;
case 55: *reinterpret_cast< QString*>(_v) = windowFilePath(); break;
case 56: *reinterpret_cast< Qt::InputMethodHints*>(_v) = inputMethodHints(); break;
}
_id -= 57;
} else if (_c == QMetaObject::WriteProperty) {
void *_v = _a[0];
switch (_id) {
case 1: setWindowModality(*reinterpret_cast< Qt::WindowModality*>(_v)); break;
case 2: setEnabled(*reinterpret_cast< bool*>(_v)); break;
case 3: setGeometry(*reinterpret_cast< QRect*>(_v)); break;
case 8: move(*reinterpret_cast< QPoint*>(_v)); break;
case 10: resize(*reinterpret_cast< QSize*>(_v)); break;
case 16: setSizePolicy(*reinterpret_cast< QSizePolicy*>(_v)); break;
case 17: setMinimumSize(*reinterpret_cast< QSize*>(_v)); break;
case 18: setMaximumSize(*reinterpret_cast< QSize*>(_v)); break;
case 19: setMinimumWidth(*reinterpret_cast< int*>(_v)); break;
case 20: setMinimumHeight(*reinterpret_cast< int*>(_v)); break;
case 21: setMaximumWidth(*reinterpret_cast< int*>(_v)); break;
case 22: setMaximumHeight(*reinterpret_cast< int*>(_v)); break;
case 23: setSizeIncrement(*reinterpret_cast< QSize*>(_v)); break;
case 24: setBaseSize(*reinterpret_cast< QSize*>(_v)); break;
case 25: setPalette(*reinterpret_cast< QPalette*>(_v)); break;
case 26: setFont(*reinterpret_cast< QFont*>(_v)); break;
case 27: setCursor(*reinterpret_cast< QCursor*>(_v)); break;
case 28: setMouseTracking(*reinterpret_cast< bool*>(_v)); break;
case 30: setFocusPolicy(*reinterpret_cast< Qt::FocusPolicy*>(_v)); break;
case 32: setContextMenuPolicy(*reinterpret_cast< Qt::ContextMenuPolicy*>(_v)); break;
case 33: setUpdatesEnabled(*reinterpret_cast< bool*>(_v)); break;
case 34: setVisible(*reinterpret_cast< bool*>(_v)); break;
case 40: setAcceptDrops(*reinterpret_cast< bool*>(_v)); break;
case 41: setWindowTitle(*reinterpret_cast< QString*>(_v)); break;
case 42: setWindowIcon(*reinterpret_cast< QIcon*>(_v)); break;
case 43: setWindowIconText(*reinterpret_cast< QString*>(_v)); break;
case 44: setWindowOpacity(*reinterpret_cast< double*>(_v)); break;
case 45: setWindowModified(*reinterpret_cast< bool*>(_v)); break;
case 46: setToolTip(*reinterpret_cast< QString*>(_v)); break;
case 47: setStatusTip(*reinterpret_cast< QString*>(_v)); break;
case 48: setWhatsThis(*reinterpret_cast< QString*>(_v)); break;
case 49: setAccessibleName(*reinterpret_cast< QString*>(_v)); break;
case 50: setAccessibleDescription(*reinterpret_cast< QString*>(_v)); break;
case 51: setLayoutDirection(*reinterpret_cast< Qt::LayoutDirection*>(_v)); break;
case 52: setAutoFillBackground(*reinterpret_cast< bool*>(_v)); break;
case 53: setStyleSheet(*reinterpret_cast< QString*>(_v)); break;
case 54: setLocale(*reinterpret_cast< QLocale*>(_v)); break;
case 55: setWindowFilePath(*reinterpret_cast< QString*>(_v)); break;
case 56: setInputMethodHints(*reinterpret_cast< Qt::InputMethodHints*>(_v)); break;
}
_id -= 57;
} else if (_c == QMetaObject::ResetProperty) {
switch (_id) {
case 27: unsetCursor(); break;
case 51: unsetLayoutDirection(); break;
case 54: unsetLocale(); break;
}
_id -= 57;
} else if (_c == QMetaObject::QueryPropertyDesignable) {
bool *_b = reinterpret_cast<bool*>(_a[0]);
switch (_id) {
case 41: *_b = isWindow(); break;
case 42: *_b = isWindow(); break;
case 43: *_b = isWindow(); break;
case 44: *_b = isWindow(); break;
case 45: *_b = isWindow(); break;
case 55: *_b = isWindow(); break;
}
_id -= 57;
} else if (_c == QMetaObject::QueryPropertyScriptable) {
_id -= 57;
} else if (_c == QMetaObject::QueryPropertyStored) {
_id -= 57;
} else if (_c == QMetaObject::QueryPropertyEditable) {
_id -= 57;
} else if (_c == QMetaObject::QueryPropertyUser) {
_id -= 57;
}
#endif // QT_NO_PROPERTIES
return _id;
}
void AnalogClock::mousePressEvent( QMouseEvent *e )
{
if(isTopLevel())
clickPos = e->pos() + QPoint(geometry().topLeft() - frameGeometry().topLeft());
}
void FancyTabWidget::setBackgroundBrush(const QBrush &brush)
{
QPalette pal;
pal.setBrush(QPalette::Mid, brush);
m_tabBar->setPalette(pal);
m_cornerWidgetContainer->setPalette(pal);
}
void FancyTabWidget::paintEvent(QPaintEvent *event)
{
Q_UNUSED(event)
if (m_selectionWidget->isVisible()) {
QPainter painter(this);
QRect rect = m_selectionWidget->rect().adjusted(0, 0, 1, 0);
rect = style()->visualRect(layoutDirection(), geometry(), rect);
StyleHelper::verticalGradient(&painter, rect, rect);
painter.setPen(StyleHelper::borderColor());
painter.drawLine(rect.topRight(), rect.bottomRight());
QColor light = StyleHelper::sidebarHighlight();
painter.setPen(light);
painter.drawLine(rect.bottomLeft(), rect.bottomRight());
}
}
void FancyTabWidget::insertCornerWidget(int pos, QWidget *widget)
{
QVBoxLayout *layout = static_cast<QVBoxLayout *>(m_cornerWidgetContainer->layout());
layout->insertWidget(pos, widget);
}
