EGwsStatus CGwsFdoInsertCommand::Init (const wchar_t* pFDOCommandClass /*NULL*/)
{
EGwsStatus fdoes = eGwsOk;
try {
eGwsOkThrow (CGwsFdoCommand::Init ());
m_pCommand = m_connection->CreateCommand(FdoCommandType_Insert);
assert (m_pCommand != NULL);
// Set class name
WSTR name;
if(pFDOCommandClass)
name = pFDOCommandClass;
else
name = QualifiedClassName();
((FdoIInsert *)m_pCommand.p)->SetFeatureClassName(name.c_str());
} catch(FdoException *e) {
PushFdoException (eGwsFailedToPrepareCommand, e);
e->Release();
fdoes = eGwsFailedToPrepareCommand;
} catch(EGwsStatus es) {
PushStatus (es);
fdoes = es;
}
return fdoes;
}
EGwsStatus CGwsNestedLoopsJoinQueryResults::InitializeReader (
IGWSQuery * query,
CGwsPreparedJoinQuery * prepquery,
bool bScrollable
)
{
EGwsStatus stat = eGwsOk;
CGwsPreparedQuery * leftquery = prepquery->LeftQuery ();
CGwsPreparedQuery * rightquery = prepquery->RightQuery ();
FdoPtr<FdoStringCollection> leftcols = prepquery->LeftProperties ();
FdoPtr<FdoStringCollection> rightcols = prepquery->RightProperties ();
stat = CGwsJoinQueryResults::InitializeReader (leftcols, query, leftquery, bScrollable);
if (IGWSException::IsError (stat)) {
PushStatus (stat);
return stat;
}
m_prepquery = prepquery;
CGwsRightJoinQueryResults * results =
(CGwsRightJoinQueryResults *) rightquery->CreateFeatureIterator (eGwsRightNestedLoopsIterator);
stat = results->InitializeReader (query, rightquery, rightcols, bScrollable);
if (IGWSException::IsError (stat)) {
delete results;
} else {
m_right = results;
m_right->AddRef ();
}
return stat;
}
void CGwsObject::PushFdoException (
const CGwsStatus & stat,
FdoException * fdoEx
) const
{
FdoException * e = fdoEx;
e->AddRef ();
// do I need to reverse fdo exception chain?
while (e != NULL) {
PushStatus (CGwsStatus (eGwsFdoProviderError, e->GetExceptionMessage ()));
FdoException * e1 = e->GetCause ();
e->Release ();
e = e1;
}
PushStatus (stat);
}
void CPrivateSendServer::ProcessMessage(CNode* pfrom, std::string& strCommand, CDataStream& vRecv, CConnman& connman)
{
if(!fMasterNode) return;
if(fLiteMode) return; // ignore all Monoeci related functionality
if(!masternodeSync.IsBlockchainSynced()) return;
if(strCommand == NetMsgType::DSACCEPT) {
if(pfrom->nVersion < MIN_PRIVATESEND_PEER_PROTO_VERSION) {
LogPrintf("DSACCEPT -- incompatible version! nVersion: %d\n", pfrom->nVersion);
PushStatus(pfrom, STATUS_REJECTED, ERR_VERSION, connman);
return;
}
if(IsSessionReady()) {
// too many users in this session already, reject new ones
LogPrintf("DSACCEPT -- queue is already full!\n");
PushStatus(pfrom, STATUS_ACCEPTED, ERR_QUEUE_FULL, connman);
return;
}
int nDenom;
CTransaction txCollateral;
vRecv >> nDenom >> txCollateral;
LogPrint("privatesend", "DSACCEPT -- nDenom %d (%s) txCollateral %s", nDenom, CPrivateSend::GetDenominationsToString(nDenom), txCollateral.ToString());
masternode_info_t mnInfo;
if(!mnodeman.GetMasternodeInfo(activeMasternode.outpoint, mnInfo)) {
PushStatus(pfrom, STATUS_REJECTED, ERR_MN_LIST, connman);
return;
}
if(vecSessionCollaterals.size() == 0 && mnInfo.nLastDsq != 0 &&
mnInfo.nLastDsq + mnodeman.CountEnabled(MIN_PRIVATESEND_PEER_PROTO_VERSION)/5 > mnodeman.nDsqCount)
{
LogPrintf("DSACCEPT -- last dsq too recent, must wait: addr=%s\n", pfrom->addr.ToString());
PushStatus(pfrom, STATUS_REJECTED, ERR_RECENT, connman);
return;
}
PoolMessage nMessageID = MSG_NOERR;
bool fResult = nSessionID == 0 ? CreateNewSession(nDenom, txCollateral, nMessageID, connman)
: AddUserToExistingSession(nDenom, txCollateral, nMessageID);
if(fResult) {
LogPrintf("DSACCEPT -- is compatible, please submit!\n");
PushStatus(pfrom, STATUS_ACCEPTED, nMessageID, connman);
return;
} else {
LogPrintf("DSACCEPT -- not compatible with existing transactions!\n");
PushStatus(pfrom, STATUS_REJECTED, nMessageID, connman);
return;
}
} else if(strCommand == NetMsgType::DSQUEUE) {
EGwsStatus CGwsPreparedJoinQuery::Execute (
FdoFilter * filter,
IGWSFeatureIterator ** results,
bool bScrollable
)
{
EGwsStatus stat = eGwsOk;
try {
eGwsOkThrow (SetFilter (filter));
return Execute (results,bScrollable);
} catch (EGwsStatus es) {
PushStatus (es);
stat = es;
}
return stat;
}
EGwsStatus CGwsPreparedFeatureQuery::Execute (
const GWSFeatureId & featid,
IGWSFeatureIterator ** results
)
{
FdoPtr<FdoFilter> filter;
EGwsStatus stat = eGwsOk;
try {
eGwsOkThrow (BuildFilter (featid, filter.p));
return Execute (filter, results);
} catch (EGwsStatus es) {
PushStatus (es);
stat = es;
}
return stat;
}
EGwsStatus CGwsPreparedFeatureQuery::Execute (
const GwsFeaturesIdVector & featids,
int lBound,
int uBound,
IGWSFeatureIterator ** results,
bool bScrollable
)
{
FdoPtr<FdoFilter> filter;
EGwsStatus stat = eGwsOk;
try {
eGwsOkThrow (BuildFilter (featids, lBound, uBound, filter.p));
return Execute (filter, results, bScrollable);
} catch (EGwsStatus es) {
PushStatus (es);
stat = es;
}
return stat;
}
EGwsStatus CGwsPreparedJoinQuery::Execute (
const GWSFeatureId & featid,
IGWSFeatureIterator ** results
)
{
CGwsPreparedFeatureQuery * pquery = NULL;
FdoPtr<FdoFilter> filter;
EGwsStatus stat = eGwsOk;
try {
pquery = GetPrimaryQuery ();
eGwsOkThrow (pquery->BuildFilter (featid, filter.p));
return Execute (filter, results);
} catch (EGwsStatus es) {
PushStatus (es);
stat = es;
}
return stat;
}
EGwsStatus CGwsPreparedJoinQuery::Execute (
const GwsFeaturesIdVector & featids,
int lbound,
int ubound,
IGWSFeatureIterator ** results,
bool bScrollable
)
{
CGwsPreparedFeatureQuery * pquery = NULL;
FdoPtr<FdoFilter> filter;
EGwsStatus stat = eGwsOk;
try {
pquery = GetPrimaryQuery ();
eGwsOkThrow (pquery->BuildFilter (featids, lbound, ubound, filter.p));
return Execute (filter, results, bScrollable);
} catch (EGwsStatus es) {
PushStatus (es);
stat = es;
}
return stat;
}
void RepairBisection(Mesh & mesh, Array<ElementIndex> & bad_elements,
const BitArray & isnewpoint, const Refinement & refinement,
const Array<double> & pure_badness,
double max_worsening, const bool uselocalworsening,
const Array< Array<int,PointIndex::BASE>* > & idmaps)
{
ostringstream ostrstr;
const int maxtrials = 100;
//bool doit;
//cout << "DOIT: " << flush;
//cin >> doit;
int ne = mesh.GetNE();
int np = mesh.GetNP();
int numbadneighbours = 3;
const int numtopimprove = 3;
PrintMessage(1,"repairing");
PushStatus("Repair Bisection");
Array<Point<3>* > should(np);
Array<Point<3>* > can(np);
Array<Vec<3>* > nv(np);
for(int i=0; i<np; i++)
{
nv[i] = new Vec<3>;
should[i] = new Point<3>;
can[i] = new Point<3>;
}
BitArray isboundarypoint(np),isedgepoint(np);
isboundarypoint.Clear();
isedgepoint.Clear();
for(int i = 1; i <= mesh.GetNSeg(); i++)
{
const Segment & seg = mesh.LineSegment(i);
isedgepoint.Set(seg[0]);
isedgepoint.Set(seg[1]);
}
Array<int> surfaceindex(np);
surfaceindex = -1;
for (int i = 1; i <= mesh.GetNSE(); i++)
{
const Element2d & sel = mesh.SurfaceElement(i);
for (int j = 1; j <= sel.GetNP(); j++)
if(!isedgepoint.Test(sel.PNum(j)))
{
isboundarypoint.Set(sel.PNum(j));
surfaceindex[sel.PNum(j) - PointIndex::BASE] =
mesh.GetFaceDescriptor(sel.GetIndex()).SurfNr();
}
}
Validate(mesh,bad_elements,pure_badness,
((uselocalworsening) ? (0.8*(max_worsening-1.) + 1.) : (0.1*(max_worsening-1.) + 1.)),
uselocalworsening); // -> larger working area
BitArray working_elements(ne);
BitArray working_points(np);
GetWorkingArea(working_elements,working_points,mesh,bad_elements,numbadneighbours);
//working_elements.Set();
//working_points.Set();
ostrstr.str("");
ostrstr << "worsening: " <<
Validate(mesh,bad_elements,pure_badness,max_worsening,uselocalworsening);
PrintMessage(4,ostrstr.str());
int auxnum=0;
for(int i=1; i<=np; i++)
if(working_points.Test(i))
auxnum++;
ostrstr.str("");
ostrstr << "Percentage working points: " << 100.*double(auxnum)/np;
PrintMessage(5,ostrstr.str());
BitArray isworkingboundary(np);
for(int i=1; i<=np; i++)
if(working_points.Test(i) && isboundarypoint.Test(i))
isworkingboundary.Set(i);
else
isworkingboundary.Clear(i);
for(int i=0; i<np; i++)
*should[i] = mesh.Point(i+1);
for(int i=0; i<np; i++)
{
if(isnewpoint.Test(i+PointIndex::BASE) &&
//working_points.Test(i+PointIndex::BASE) &&
mesh.mlbetweennodes[i+PointIndex::BASE][0] > 0)
*can[i] = Center(*can[mesh.mlbetweennodes[i+PointIndex::BASE][0]-PointIndex::BASE],
*can[mesh.mlbetweennodes[i+PointIndex::BASE][1]-PointIndex::BASE]);
else
*can[i] = mesh.Point(i+1);
}
int cnttrials = 1;
double lamedge = 0.5;
double lamface = 0.5;
double facokedge = 0;
double facokface = 0;
double factryedge;
double factryface = 0;
double oldlamedge,oldlamface;
MeshOptimize2d * optimizer2d = refinement.Get2dOptimizer();
if(!optimizer2d)
{
cerr << "No 2D Optimizer!" << endl;
return;
}
while ((facokedge < 1.-1e-8 || facokface < 1.-1e-8) &&
cnttrials < maxtrials &&
multithread.terminate != 1)
{
(*testout) << " facokedge " << facokedge << " facokface " << facokface << " cnttrials " << cnttrials << endl
<< " perc. " << 95. * max2( min2(facokedge,facokface),
double(cnttrials)/double(maxtrials)) << endl;
SetThreadPercent(95. * max2( min2(facokedge,facokface),
double(cnttrials)/double(maxtrials)));
ostrstr.str("");
ostrstr << "max. worsening " << max_worsening;
PrintMessage(5,ostrstr.str());
oldlamedge = lamedge;
lamedge *= 6;
if (lamedge > 2)
lamedge = 2;
if(1==1 || facokedge < 1.-1e-8)
{
for(int i=0; i<nv.Size(); i++)
*nv[i] = Vec<3>(0,0,0);
for (int i = 1; i <= mesh.GetNSE(); i++)
{
const Element2d & sel = mesh.SurfaceElement(i);
Vec<3> auxvec = Cross(mesh.Point(sel.PNum(2))-mesh.Point(sel.PNum(1)),
mesh.Point(sel.PNum(3))-mesh.Point(sel.PNum(1)));
auxvec.Normalize();
for (int j = 1; j <= sel.GetNP(); j++)
if(!isedgepoint.Test(sel.PNum(j)))
*nv[sel.PNum(j) - PointIndex::BASE] += auxvec;
}
for(int i=0; i<nv.Size(); i++)
nv[i]->Normalize();
do // move edges
{
lamedge *= 0.5;
cnttrials++;
if(cnttrials % 10 == 0)
max_worsening *= 1.1;
factryedge = lamedge + (1.-lamedge) * facokedge;
ostrstr.str("");
ostrstr << "lamedge = " << lamedge << ", trying: " << factryedge;
PrintMessage(5,ostrstr.str());
for (int i = 1; i <= np; i++)
{
if (isedgepoint.Test(i))
{
for (int j = 0; j < 3; j++)
mesh.Point(i)(j) =
lamedge * (*should.Get(i))(j) +
(1.-lamedge) * (*can.Get(i))(j);
}
else
mesh.Point(i) = *can.Get(i);
}
if(facokedge < 1.-1e-8)
{
ostrstr.str("");
ostrstr << "worsening: " <<
Validate(mesh,bad_elements,pure_badness,max_worsening,uselocalworsening);
PrintMessage(5,ostrstr.str());
}
else
Validate(mesh,bad_elements,pure_badness,-1,uselocalworsening);
ostrstr.str("");
ostrstr << bad_elements.Size() << " bad elements";
PrintMessage(5,ostrstr.str());
}
while (bad_elements.Size() > 0 &&
cnttrials < maxtrials &&
multithread.terminate != 1);
}
if(cnttrials < maxtrials &&
multithread.terminate != 1)
{
facokedge = factryedge;
// smooth faces
mesh.CalcSurfacesOfNode();
MeshingParameters dummymp;
mesh.ImproveMeshJacobianOnSurface(dummymp,isworkingboundary,nv,OPT_QUALITY, &idmaps);
for (int i = 1; i <= np; i++)
*can.Elem(i) = mesh.Point(i);
if(optimizer2d)
optimizer2d->ProjectBoundaryPoints(surfaceindex,can,should);
}
oldlamface = lamface;
lamface *= 6;
if (lamface > 2)
lamface = 2;
if(cnttrials < maxtrials &&
multithread.terminate != 1)
{
do // move faces
{
lamface *= 0.5;
cnttrials++;
if(cnttrials % 10 == 0)
max_worsening *= 1.1;
factryface = lamface + (1.-lamface) * facokface;
ostrstr.str("");
ostrstr << "lamface = " << lamface << ", trying: " << factryface;
PrintMessage(5,ostrstr.str());
for (int i = 1; i <= np; i++)
{
if (isboundarypoint.Test(i))
{
for (int j = 0; j < 3; j++)
mesh.Point(i)(j) =
lamface * (*should.Get(i))(j) +
(1.-lamface) * (*can.Get(i))(j);
}
else
mesh.Point(i) = *can.Get(i);
}
ostrstr.str("");
ostrstr << "worsening: " <<
Validate(mesh,bad_elements,pure_badness,max_worsening,uselocalworsening);
PrintMessage(5,ostrstr.str());
ostrstr.str("");
ostrstr << bad_elements.Size() << " bad elements";
PrintMessage(5,ostrstr.str());
}
while (bad_elements.Size() > 0 &&
cnttrials < maxtrials &&
multithread.terminate != 1);
}
if(cnttrials < maxtrials &&
multithread.terminate != 1)
{
facokface = factryface;
// smooth interior
mesh.CalcSurfacesOfNode();
MeshingParameters dummymp;
mesh.ImproveMeshJacobian (dummymp, OPT_QUALITY,&working_points);
//mesh.ImproveMeshJacobian (OPT_WORSTCASE,&working_points);
for (int i = 1; i <= np; i++)
*can.Elem(i) = mesh.Point(i);
}
//!
if((facokedge < 1.-1e-8 || facokface < 1.-1e-8) &&
cnttrials < maxtrials &&
multithread.terminate != 1)
{
MeshingParameters dummymp;
MeshOptimize3d optmesh(dummymp);
for(int i=0; i<numtopimprove; i++)
{
optmesh.SwapImproveSurface(mesh,OPT_QUALITY,&working_elements,&idmaps);
optmesh.SwapImprove(mesh,OPT_QUALITY,&working_elements);
}
// mesh.mglevels = 1;
ne = mesh.GetNE();
working_elements.SetSize(ne);
for (int i = 1; i <= np; i++)
mesh.Point(i) = *should.Elem(i);
Validate(mesh,bad_elements,pure_badness,
((uselocalworsening) ? (0.8*(max_worsening-1.) + 1.) : (0.1*(max_worsening-1.) + 1.)),
uselocalworsening);
if(lamedge < oldlamedge || lamface < oldlamface)
numbadneighbours++;
GetWorkingArea(working_elements,working_points,mesh,bad_elements,numbadneighbours);
for(int i=1; i<=np; i++)
if(working_points.Test(i) && isboundarypoint.Test(i))
isworkingboundary.Set(i);
else
isworkingboundary.Clear(i);
auxnum=0;
for(int i=1; i<=np; i++)
if(working_points.Test(i))
auxnum++;
ostrstr.str("");
ostrstr << "Percentage working points: " << 100.*double(auxnum)/np;
PrintMessage(5,ostrstr.str());
for (int i = 1; i <= np; i++)
mesh.Point(i) = *can.Elem(i);
}
//!
}
MeshingParameters dummymp;
MeshOptimize3d optmesh(dummymp);
for(int i=0; i<numtopimprove && multithread.terminate != 1; i++)
{
optmesh.SwapImproveSurface(mesh,OPT_QUALITY,NULL,&idmaps);
optmesh.SwapImprove(mesh,OPT_QUALITY);
//mesh.UpdateTopology();
}
mesh.UpdateTopology();
/*
if(cnttrials < 100)
{
nv = Vec3d(0,0,0);
for (int i = 1; i <= mesh.GetNSE(); i++)
{
const Element2d & sel = mesh.SurfaceElement(i);
Vec3d auxvec = Cross(mesh.Point(sel.PNum(2))-mesh.Point(sel.PNum(1)),
mesh.Point(sel.PNum(3))-mesh.Point(sel.PNum(1)));
auxvec.Normalize();
for (int j = 1; j <= sel.GetNP(); j++)
if(!isedgepoint.Test(sel.PNum(j)))
nv[sel.PNum(j) - PointIndex::BASE] += auxvec;
}
for(int i=0; i<nv.Size(); i++)
nv[i].Normalize();
mesh.ImproveMeshJacobianOnSurface(isboundarypoint,nv,OPT_QUALITY);
mesh.CalcSurfacesOfNode();
// smooth interior
for (int i = 1; i <= np; i++)
if(isboundarypoint.Test(i))
can.Elem(i) = mesh.Point(i);
if(optimizer2d)
optimizer2d->ProjectBoundaryPoints(surfaceindex,can,should);
for (int i = 1; i <= np; i++)
if(isboundarypoint.Test(i))
for(int j=1; j<=3; j++)
mesh.Point(i).X(j) = should.Get(i).X(j);
}
*/
if(cnttrials == maxtrials)
{
for (int i = 1; i <= np; i++)
mesh.Point(i) = *should.Get(i);
Validate(mesh,bad_elements,pure_badness,max_worsening,uselocalworsening);
for(int i=0; i<bad_elements.Size(); i++)
{
ostrstr.str("");
ostrstr << "bad element:" << endl
<< mesh[bad_elements[i]][0] << ": " << mesh.Point(mesh[bad_elements[i]][0]) << endl
<< mesh[bad_elements[i]][1] << ": " << mesh.Point(mesh[bad_elements[i]][1]) << endl
<< mesh[bad_elements[i]][2] << ": " << mesh.Point(mesh[bad_elements[i]][2]) << endl
<< mesh[bad_elements[i]][3] << ": " << mesh.Point(mesh[bad_elements[i]][3]);
PrintMessage(5,ostrstr.str());
}
for (int i = 1; i <= np; i++)
mesh.Point(i) = *can.Get(i);
}
for(int i=0; i<np; i++)
{
delete nv[i];
delete can[i];
delete should[i];
}
PopStatus();
}
EGwsStatus CGwsFdoInsertCommand::Execute (CGwsMutableFeature & feat)
{
Clear ();
EGwsStatus fdoes = eGwsOk;
if(!m_bIsPrepared) {
PushStatus (eGwsFdoCommandIsNotPrepared);
return eGwsFdoCommandIsNotPrepared;
}
try
{
GetPropertyValues (); // initialize property value collection
eGwsOkThrow(SetProperties (feat, true)) ;
FdoPtr<FdoIFeatureReader> pReader = ((FdoIInsert*)m_pCommand.p)->Execute();
//!FDO may return null pointer even if the command succeed. So remove the following assert.
//assert (pReader != NULL);
// Read assigned Id
if( pReader != NULL )
{
if (pReader->ReadNext()) {
// Read the key of just inserted record
//
CGwsDataValueCollection * keyvals = NULL;
for (FdoInt32 idx = 0; m_identity != NULL && idx < m_identity->GetCount (); idx ++) {
FdoPtr<FdoDataPropertyDefinition> keyprop =
m_identity->GetItem (idx);
FdoPtr<FdoDataValue> val =
GwsQueryUtils::GetDataPropertyValue (pReader,
keyprop->GetDataType (),
keyprop->GetName ());
if (val == NULL)
continue;
if (keyvals == NULL)
keyvals = (CGwsDataValueCollection *) CGwsDataValueCollection::Create ();
keyvals->Add (val);
feat.SetValue (keyprop->GetName (),val);
}
#ifdef OBJECT_PROPERTIES
//insert object property values that were cached in SetProperties()
if(!m_objPropNames.empty())
{
//Get object property class instance array that was cached in SetProperties()
int numProps = m_objPropNames.size();
for( int i = 0; i<numProps; i++)
{
const std::wstring strPropName = m_objPropNames[i];
const std::wstring objectClassName =
acmapFdoUtils::constructObjectClassName(aTOw(m_classname.Class()),strPropName.c_str());
//get the object property values
AcMapFdoFeatureIterator* featIter = NULL;
eGwsOkThrow(feature.GetFeature(wTOa(strPropName.c_str()),featIter));
assert(featIter != NULL);
//insert the primary keys of the parent feature and the object property values
eGwsOkThrow(InsertObjectProperties(objectClassName.c_str(), valList, featIter));
}
}
#endif
// long ver;
// Read sequence number
if (!m_revisionprop.empty()) {
FdoInt32 revnum = pReader->GetInt32 (m_revisionprop.c_str ());
feat.SetRevisionNumber (revnum);
}
feat.SetFeatureId (GWSFeatureId (keyvals));
}
pReader->Close();
}
} catch(FdoException *e) {
PushFdoException (eGwsFailedToExecuteCommand, e);
e->Release();
fdoes = eGwsFailedToExecuteCommand;
} catch(EGwsStatus es) {
PushStatus (es);
fdoes = es;
}
ReleasePropertyValues (); // releases property value collection
return fdoes;
}
