/// <summary>
/// 問題を設定する
/// </summary>
/// <returns></returns>
private bool setNewProblem()
{
bool success = false;
try
{
if( ProbNo == 0 )
{
CMesher2D mesh2d = null;
using(CCadObj2D cad2d = new CCadObj2D())
{
IList<CVector2D> pts = new List<CVector2D>();
pts.Add( new CVector2D(0.0,0.0) );
pts.Add( new CVector2D(1.0,0.0) );
pts.Add( new CVector2D(1.0,0.2) );
pts.Add( new CVector2D(0.5,0.2) );
pts.Add( new CVector2D(0.5,0.8) );
pts.Add( new CVector2D(1.0,0.8) );
pts.Add( new CVector2D(1.0,1.0) );
pts.Add( new CVector2D(0.0,1.0) );
cad2d.AddPolygon( pts );
cad2d.ConnectVertex_Line(6,3);
mesh2d = new CMesher2D(cad2d,0.05);
////
using(CMesher2D msh_tmp = (CMesher2D)mesh2d.Clone())
{
mesh2d.Clear();
mesh2d.Dispose();
mesh2d = (CMesher2D)msh_tmp.Clone();
}
}
using(CSerializer fout = new CSerializer("hoge.txt", false))
{
// write file
mesh2d.Serialize(fout);
}
using(CSerializer fin = new CSerializer("hoge.txt", true))
{ // load file
mesh2d.Serialize(fin);
}
////
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh2D(mesh2d) );
DrawerAry.InitTrans( Camera );
mesh2d.Clear();
mesh2d.Dispose();
}
else if( ProbNo == 1 )
{
string mshfn = "../../../input_file/hexa_tri.msh";
if (File.Exists(mshfn))
{
using(CMesher2D mesh2d = new CMesher2D())
{
mesh2d.ReadFromFile_GiDMsh(mshfn);
using(CSerializer fout = new CSerializer("hoge.txt", false))
{
// write file
mesh2d.Serialize(fout);
}
using(CSerializer fin = new CSerializer("hoge.txt", true))
{ // load file
mesh2d.Serialize(fin);
}
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh2D(mesh2d) );
DrawerAry.InitTrans( Camera );
}
}
else
{
Console.WriteLine("not exist:{0}", mshfn);
MessageBox.Show(string.Format("メッシュファイル:{0}がありません", mshfn));
}
}
else if( ProbNo == 2 )
{
string mshfn = "../../../input_file/rect_quad.msh";
if (File.Exists(mshfn))
{
using(CMesher2D mesh2d = new CMesher2D())
{
mesh2d.ReadFromFile_GiDMsh(mshfn);
using(CSerializer fout = new CSerializer("hoge.txt", false))
{
// write file
mesh2d.Serialize(fout);
}
using(CSerializer fin = new CSerializer("hoge.txt", true))
{ // load file
mesh2d.Serialize(fin);
}
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh2D(mesh2d) );
DrawerAry.InitTrans( Camera );
}
}
else
{
Console.WriteLine("not exist:{0}", mshfn);
MessageBox.Show(string.Format("メッシュファイル:{0}がありません", mshfn));
}
}
else if( ProbNo == 3 )
{
using(CCadObj2D cad2d = new CCadObj2D())
{
IList<CVector2D> pts = new List<CVector2D>();
pts.Add( new CVector2D(0.0,0.0) );
pts.Add( new CVector2D(1.0,0.0) );
pts.Add( new CVector2D(1.0,0.5) );
pts.Add( new CVector2D(0.5,0.5) );
pts.Add( new CVector2D(0.5,1.0) );
pts.Add( new CVector2D(0.0,1.0) );
cad2d.AddPolygon( pts );
using(CMesher2D mesh2d = new CMesher2D(cad2d, 0.1))
using(CMesh3D_Extrude mesh3d = new CMesh3D_Extrude())
{
mesh3d.Extrude(mesh2d, 0.5, 0.1 );
using(CSerializer fout = new CSerializer("hoge.txt", false))
{
// write file
mesh3d.Serialize(fout);
}
using(CSerializer fin = new CSerializer("hoge.txt", true))
{ // load file
mesh3d.Serialize(fin);
}
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh3D(mesh3d) );
DrawerAry.InitTrans( Camera );
}
}
}
else if( ProbNo == 4 ) // load mesh of GiD
{
string mshfn = "../../../input_file/hexa_tri.msh";
if (File.Exists(mshfn))
{
using(CMesher2D mesh2d = new CMesher2D())
using(CMesh3D_Extrude mesh3d = new CMesh3D_Extrude())
{
mesh2d.ReadFromFile_GiDMsh(mshfn);
mesh3d.Extrude(mesh2d, 5.0, 0.5 );
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh3D(mesh3d) );
DrawerAry.InitTrans( Camera );
}
}
else
{
Console.WriteLine("not exist:{0}", mshfn);
MessageBox.Show(string.Format("メッシュファイル:{0}がありません", mshfn));
}
}
else if( ProbNo == 5 )
{
string mshfn = "../../../input_file/cylinder_hex.msh";
if (File.Exists(mshfn))
{
using(CMesher3D mesh3d = new CMesher3D())
{
mesh3d.ReadFromFile_GiDMsh(mshfn);
using(CSerializer fout = new CSerializer("hoge.txt", false))
{
// write file
mesh3d.Serialize(fout);
}
using(CSerializer fin = new CSerializer("hoge.txt", true))
{ // load file
mesh3d.Serialize(fin);
}
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh3D(mesh3d) );
DrawerAry.InitTrans( Camera );
}
}
else
{
Console.WriteLine("not exist:{0}", mshfn);
MessageBox.Show(string.Format("メッシュファイル:{0}がありません", mshfn));
}
}
else if( ProbNo == 6 )
{
string mshfn = "../../../input_file/cylinder_tet.msh";
if (File.Exists(mshfn))
{
using(CMesher3D mesh3d = new CMesher3D())
{
mesh3d.ReadFromFile_GiDMsh(mshfn);
using(CSerializer fout = new CSerializer("hoge.txt", false))
{
// write file
mesh3d.Serialize(fout);
}
using(CSerializer fin = new CSerializer("hoge.txt", true))
{ // load file
mesh3d.Serialize(fin);
}
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh3D(mesh3d) );
DrawerAry.InitTrans( Camera );
}
}
else
{
Console.WriteLine("not exist:{0}", mshfn);
MessageBox.Show(string.Format("メッシュファイル:{0}がありません", mshfn));
}
}
else if( ProbNo == 7 )
{
string mshfn = "../../../input_file/rect_quad.msh";
if (File.Exists(mshfn))
{
using(CMesher2D mesh2d = new CMesher2D())
using(CMesh3D_Extrude mesh3d = new CMesh3D_Extrude())
{
mesh2d.ReadFromFile_GiDMsh(mshfn);
mesh3d.Extrude(mesh2d, 5.0, 0.5 );
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh3D(mesh3d) );
DrawerAry.InitTrans( Camera );
}
}
else
{
Console.WriteLine("not exist:{0}", mshfn);
MessageBox.Show(string.Format("メッシュファイル:{0}がありません", mshfn));
}
}
else if( ProbNo == 8 )
{
using(CCadObj2D cad2d = new CCadObj2D())
{
uint id_l = 0;
{
IList<CVector2D> pts = new List<CVector2D>();
pts.Add( new CVector2D(0.0,0.0) );
pts.Add( new CVector2D(1.0,0.0) );
pts.Add( new CVector2D(1.0,1.0) );
pts.Add( new CVector2D(0.0,1.0) );
id_l = cad2d.AddPolygon( pts ).id_l_add;
}
cad2d.AddVertex(CAD_ELEM_TYPE.LOOP,id_l, new CVector2D(0.8,0.6) );
cad2d.AddVertex(CAD_ELEM_TYPE.LOOP,id_l, new CVector2D(0.6,0.6) );
cad2d.AddVertex(CAD_ELEM_TYPE.LOOP,id_l, new CVector2D(0.4,0.6) );
cad2d.AddVertex(CAD_ELEM_TYPE.LOOP,id_l, new CVector2D(0.2,0.6) );
cad2d.AddVertex(CAD_ELEM_TYPE.LOOP,id_l, new CVector2D(0.8,0.4) );
cad2d.AddVertex(CAD_ELEM_TYPE.LOOP,id_l, new CVector2D(0.6,0.4) );
cad2d.AddVertex(CAD_ELEM_TYPE.LOOP,id_l, new CVector2D(0.4,0.4) );
cad2d.AddVertex(CAD_ELEM_TYPE.LOOP,id_l, new CVector2D(0.2,0.4) );
using(CMesher2D mesh2d = new CMesher2D(cad2d, 0.02))
{
using(CSerializer fout = new CSerializer("hoge.txt", false))
{
// write file
mesh2d.Serialize(fout);
}
using(CSerializer fin = new CSerializer("hoge.txt", true))
{ // load file
mesh2d.Serialize(fin);
}
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh2D(mesh2d) );
DrawerAry.InitTrans( Camera );
}
}
}
else if( ProbNo == 9 )
{
using(CCadObj2D cad2d = new CCadObj2D())
{
{
IList<CVector2D> pts = new List<CVector2D>();
pts.Add( new CVector2D(0.0,0.0) );
pts.Add( new CVector2D(1.0,0.0) );
pts.Add( new CVector2D(1.0,1.0) );
pts.Add( new CVector2D(0.0,1.0) );
cad2d.AddPolygon( pts );
}
cad2d.SetCurve_Arc(1,true, -0.5);
cad2d.SetCurve_Arc(2,false,-0.5);
cad2d.SetCurve_Arc(3,true, -0.5);
cad2d.SetCurve_Arc(4,false,-0.5);
using(CMesher2D mesh2d = new CMesher2D(cad2d, 0.05))
{
using(CSerializer fout = new CSerializer("hoge.txt", false))
{
// write file
mesh2d.Serialize(fout);
}
using(CSerializer fin = new CSerializer("hoge.txt", true))
{ // load file
mesh2d.Serialize(fin);
}
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh2D(mesh2d) );
DrawerAry.InitTrans( Camera );
}
}
}
else if( ProbNo == 10 )
{
using(CCadObj2D cad2d = new CCadObj2D())
{
{
IList<CVector2D> pts = new List<CVector2D>();
pts.Add( new CVector2D(0.0,0.0) );
pts.Add( new CVector2D(1.0,0.0) );
pts.Add( new CVector2D(1.0,1.0) );
pts.Add( new CVector2D(0.0,1.0) );
uint id_l = cad2d.AddPolygon( pts ).id_l_add;
uint id_v1 = cad2d.AddVertex(CAD_ELEM_TYPE.LOOP,id_l, new CVector2D(0.3,0.2)).id_v_add;
uint id_v2 = cad2d.AddVertex(CAD_ELEM_TYPE.LOOP,id_l, new CVector2D(0.7,0.2)).id_v_add;
uint id_v3 = cad2d.AddVertex(CAD_ELEM_TYPE.LOOP,id_l, new CVector2D(0.7,0.8)).id_v_add;
uint id_v4 = cad2d.AddVertex(CAD_ELEM_TYPE.LOOP,id_l, new CVector2D(0.3,0.8)).id_v_add;
cad2d.ConnectVertex_Line(id_v1,id_v2);
cad2d.ConnectVertex_Line(id_v2,id_v3);
cad2d.ConnectVertex_Line(id_v3,id_v4);
cad2d.ConnectVertex_Line(id_v4,id_v1);
}
using(CMesher2D mesh2d = new CMesher2D())
{
mesh2d.AddIdLCad_CutMesh(1); // cut mesh to loop whitch have id 1
mesh2d.SetMeshingMode_ElemLength(0.05);
mesh2d.Meshing(cad2d);
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh2D(mesh2d) );
DrawerAry.InitTrans( Camera );
}
}
}
else if( ProbNo == 11 ) // mesh with cut
{
using(CCadObj2D cad2d = new CCadObj2D())
{
uint id_l;
uint id_e1, id_e2, id_e3, id_e4, id_e5;
{
IList<CVector2D> pts = new List<CVector2D>();
pts.Add( new CVector2D(0.0,0.0) );
pts.Add( new CVector2D(0.3,0.0) );
pts.Add( new CVector2D(1.0,0.0) );
pts.Add( new CVector2D(1.0,1.0) );
pts.Add( new CVector2D(0.0,1.0) );
id_l = cad2d.AddPolygon( pts ).id_l_add;
uint id_v1 = cad2d.AddVertex(CAD_ELEM_TYPE.LOOP, id_l, new CVector2D(0.3,0.5) ).id_v_add;
id_e1 = cad2d.ConnectVertex_Line(2,id_v1).id_e_add;
uint id_v2 = cad2d.AddVertex(CAD_ELEM_TYPE.LOOP, id_l, new CVector2D(0.7,0.5) ).id_v_add;
uint id_v3 = cad2d.AddVertex(CAD_ELEM_TYPE.LOOP, id_l, new CVector2D(0.7,0.2) ).id_v_add;
uint id_v4 = cad2d.AddVertex(CAD_ELEM_TYPE.LOOP, id_l, new CVector2D(0.7,0.8) ).id_v_add;
uint id_v5 = cad2d.AddVertex(CAD_ELEM_TYPE.LOOP, id_l, new CVector2D(0.5,0.5) ).id_v_add;
uint id_v6 = cad2d.AddVertex(CAD_ELEM_TYPE.LOOP, id_l, new CVector2D(0.9,0.5) ).id_v_add;
id_e2 = cad2d.ConnectVertex_Line(id_v2,id_v3).id_e_add;
id_e3 = cad2d.ConnectVertex_Line(id_v2,id_v4).id_e_add;
id_e4 = cad2d.ConnectVertex_Line(id_v2,id_v5).id_e_add;
id_e5 = cad2d.ConnectVertex_Line(id_v2,id_v6).id_e_add;
}
using(CMesher2D mesh2d = new CMesher2D(cad2d, 0.2))
{
IList<uint> aIdMsh_Inc = new List<uint>();
aIdMsh_Inc.Add( mesh2d.GetElemID_FromCadID(id_l,CAD_ELEM_TYPE.LOOP) );
IList<uint> aIdMshBar_Cut = new List<uint>();
aIdMshBar_Cut.Add( mesh2d.GetElemID_FromCadID(id_e1,CAD_ELEM_TYPE.EDGE) );
aIdMshBar_Cut.Add( mesh2d.GetElemID_FromCadID(id_e2,CAD_ELEM_TYPE.EDGE) );
aIdMshBar_Cut.Add( mesh2d.GetElemID_FromCadID(id_e3,CAD_ELEM_TYPE.EDGE) );
aIdMshBar_Cut.Add( mesh2d.GetElemID_FromCadID(id_e4,CAD_ELEM_TYPE.EDGE) );
aIdMshBar_Cut.Add( mesh2d.GetElemID_FromCadID(id_e5,CAD_ELEM_TYPE.EDGE) );
////////////////
IList< IList<int> > aLnods = new List<IList<int>>();
IList<uint> mapVal2Co = new List<uint>();
mesh2d.GetClipedMesh(out aLnods, out mapVal2Co, aIdMsh_Inc, aIdMshBar_Cut);
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh2D(mesh2d) );
DrawerAry.InitTrans( Camera );
}
}
}
else if( ProbNo == 12 ) // mesh with cut
{
using(CCadObj2D cad2d = new CCadObj2D())
{
uint id_l1, id_l2;
uint id_e3, id_e4;
{
IList<CVector2D> pts = new List<CVector2D>();
pts.Add( new CVector2D(0.0, 0.0) );
pts.Add( new CVector2D(1.0, 0.0) );
pts.Add( new CVector2D(1.5, 0.0) );
pts.Add( new CVector2D(2.0, 0.0) );
pts.Add( new CVector2D(2.0, 1.0) );
pts.Add( new CVector2D(1.5, 1.0) );
pts.Add( new CVector2D(1.0, 1.0) );
pts.Add( new CVector2D(0.0, 1.0) );
uint id_l0 = cad2d.AddPolygon( pts ).id_l_add;
uint id_v1 = cad2d.AddVertex(CAD_ELEM_TYPE.LOOP,id_l0, new CVector2D(0.5,0.5) ).id_v_add;
uint id_e1 = cad2d.ConnectVertex_Line(2,7).id_e_add;
uint id_e2 = cad2d.ConnectVertex_Line(3,6).id_e_add;
uint id_v2 = cad2d.AddVertex(CAD_ELEM_TYPE.EDGE,id_e1, new CVector2D(1.0,0.5) ).id_v_add;
uint id_v3 = cad2d.AddVertex(CAD_ELEM_TYPE.EDGE,1, new CVector2D(0.5,0.0) ).id_v_add;
id_e3 = cad2d.ConnectVertex_Line(id_v1,id_v2).id_e_add;
id_e4 = cad2d.ConnectVertex_Line(id_v1,id_v3).id_e_add;
id_l1 = 1;
id_l2 = 2;
}
using(CMesher2D mesh2d = new CMesher2D(cad2d, 0.2))
{
IList<uint> aIdMsh_Inc = new List<uint>();
aIdMsh_Inc.Add( mesh2d.GetElemID_FromCadID(id_l1,CAD_ELEM_TYPE.LOOP) );
aIdMsh_Inc.Add( mesh2d.GetElemID_FromCadID(id_l2,CAD_ELEM_TYPE.LOOP) );
IList<uint> aIdMshBar_Cut = new List<uint>();
//aIdMshBar_Cut.Add( mesh2d.GetElemID_FromCadID(id_e3,CAD_ELEM_TYPE.EDGE) );
aIdMshBar_Cut.Add( mesh2d.GetElemID_FromCadID(id_e4,CAD_ELEM_TYPE.EDGE) );
////////////////
IList< IList<int> > aLnods = new List<IList<int>>();
IList<uint> mapVal2Co = new List<uint>();
mesh2d.GetClipedMesh(out aLnods, out mapVal2Co, aIdMsh_Inc, aIdMshBar_Cut);
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh2D(mesh2d) );
DrawerAry.InitTrans( Camera );
}
}
}
else if( ProbNo == 13 )
{
using(CCadObj2D cad2d = new CCadObj2D())
{
{
IList<CVector2D> pts = new List<CVector2D>();
pts.Add( new CVector2D(0.0, 0.0) ); // 1
pts.Add( new CVector2D(1.5, 0.0) ); // 2
pts.Add( new CVector2D(1.5, 0.4) ); // 3
pts.Add( new CVector2D(1.0, 0.4) ); // 4
pts.Add( new CVector2D(1.0, 0.5) ); // 5
pts.Add( new CVector2D(2.0, 0.5) ); // 6
pts.Add( new CVector2D(2.0, 1.0) ); // 7
pts.Add( new CVector2D(0.0, 1.0) ); // 8
pts.Add( new CVector2D(0.0, 0.5) ); // 9
uint id_l0 = cad2d.AddPolygon( pts ).id_l_add;
uint id_e1 = cad2d.ConnectVertex_Line(5,9).id_e_add;
cad2d.ShiftLayer_Loop(id_l0,true);
double[] col = new double[3] { 0.9, 0.4, 0.4 };
cad2d.SetColor_Loop(id_l0, col);
cad2d.AddVertex(CAD_ELEM_TYPE.EDGE,3, new CVector2D(1.3,0.5));
}
using(CMesher2D mesh2d = new CMesher2D(cad2d, 0.05))
{
DrawerAry.Clear();
DrawerAry.PushBack( new CDrawerMsh2D(mesh2d) );
DrawerAry.InitTrans( Camera );
}
}
}
else if( ProbNo == 14 )
{
string svgfn = "../../../input_file/shape2d_0.svg";
if (File.Exists(svgfn))
{
CCadObj2D cad2d = null;
CCadSVG.ReadSVG_AddLoopCad(svgfn, out cad2d);
using (CMesher2D mesh2d = new CMesher2D())
{
mesh2d.SetMeshingMode_ElemSize(400);
mesh2d.AddIdLCad_CutMesh(1);
mesh2d.Meshing(cad2d);
DrawerAry.Clear();
DrawerAry.PushBack(new CDrawerMsh2D(mesh2d));
DrawerAry.InitTrans(Camera);
}
cad2d.Clear();
cad2d.Dispose();
}
else
{
Console.WriteLine("not exist:{0}", svgfn);
MessageBox.Show(string.Format("SVGファイル:{0}がありません", svgfn));
}
}
success = true;
}
catch (Exception exception)
{
Console.WriteLine(exception.Message + " " + exception.StackTrace);
}
Gl.glMatrixMode(Gl.GL_PROJECTION);
Gl.glLoadIdentity();
DrawerGlUtility.SetProjectionTransform(Camera);
Glut.glutPostRedisplay();
ProbNo++;
if (ProbNo == ProbCnt) ProbNo = 0;
return success;
}