void e161122(int dim, int trials, char* path)
{
LTQ ltq;
ltq.SetDimension(dim);
Results r;
r.Add("SR");
r.Add("SR-C");
r.Add("SR-C2");
for (int i = 1; i <= trials; i++)
{
printf_s("\r%5d/%5d", i, trials);
for (int faultRatio = 0; faultRatio < 10; faultRatio++)
{
ltq.GenerateFaults(faultRatio * 10); // 故障を発生させる
uint32_t node1, node2; // 出発ノードと目的ノード
do
{
node1 = ltq.GetNodeRandom();
node2 = ltq.GetConnectedNodeRandom(node1);
} while (node2 == node1); // 連結な候補が見つかるまでループ
r.Update(0, faultRatio, ltq.Routing_Simple(node1, node2));
Score* c = ltq.CalcCapability();
r.Update(1, faultRatio, ltq.Routing_SimpleCapability(node1, node2, c));
delete c;
c = ltq.CalcCapability2();
r.Update(2, faultRatio, ltq.Routing_SimpleCapability2(node1, node2, c));
delete c;
}
}
r.Save(path, trials);
}
string PropGeom::BuildBEMResults()
{
// Calculate prop center and normal vector
vec3d cen = m_ModelMatrix.xform( vec3d( 0, 0, 0 ) );
vec3d norm = m_ModelMatrix.xform( vec3d( -1.0, 0, 0 ) ) - cen;
int n = m_TessU();
//==== Create Results ====//
Results* res = ResultsMgr.CreateResults( "PropBEM" );
res->Add( NameValData( "Num_Sections", n ) );
res->Add( NameValData( "Num_Blade", m_Nblade() ) );
res->Add( NameValData( "Diameter", m_Diameter() ) );
res->Add( NameValData( "Beta34", m_Beta34() ) );
res->Add( NameValData( "Feather", m_Feather() ) );
res->Add( NameValData( "Center", cen ) );
res->Add( NameValData( "Normal", norm ) );
double rfirst = m_ChordCurve.GetRFirst();
double rlast = m_ChordCurve.GetRLast();
// Establish points to evaluate sections at.
vector < double > vtess;
m_MainSurfVec[0].MakeVTess( m_TessW(), vtess, m_CapUMinTess(), false );
vector < double > r_vec(n);
vector < double > chord_vec(n);
vector < double > twist_vec(n);
vector < double > rake_vec(n);
vector < double > skew_vec(n);
double rspan = rlast - rfirst;
for ( int i = 0; i < n; i++ )
{
double t = static_cast < double > ( i ) / ( n - 1 );
double r = rfirst + rspan * Cluster( t, m_RootCluster(), m_TipCluster() );
double u = m_rtou.CompPnt( r );
VspCurve c;
m_FoilSurf.GetUConstCurve( c, u );
vec3d v = c.CompPnt( 0 );
c.OffsetZ( -v.z() );
vector < vec3d > pts;
c.Tesselate( vtess, pts );
vector < double > xpts( pts.size() );
vector < double > ypts( pts.size() );
for ( int j = 0; j < pts.size(); j++ )
{
xpts[j] = pts[j].x();
ypts[j] = pts[j].y();
}
char str[255];
sprintf( str, "%03d", i );
res->Add( NameValData( "XSection_" + string( str ), xpts ) );
res->Add( NameValData( "YSection_" + string( str ), ypts ) );
r_vec[i] = r;
chord_vec[i] = m_ChordCurve.Comp( r );
twist_vec[i] = m_TwistCurve.Comp( r );
rake_vec[i] = m_RakeCurve.Comp( r );
skew_vec[i] = m_SkewCurve.Comp( r );
}
res->Add( NameValData( "Radius", r_vec ) );
res->Add( NameValData( "Chord", chord_vec ) );
res->Add( NameValData( "Twist", twist_vec ) );
res->Add( NameValData( "Rake", rake_vec ) );
res->Add( NameValData( "Skew", skew_vec ) );
return res->GetID();
}