public CMobPool(CField parentField) : base(parentField)
{
//Spawns = new List<CMob>();
aMobGen = new List <MOBGEN>();
aMobGenExcept = new List <int>();
m_nSubMobCount = -1;
}
private static void ProcessRows(RowSet rows)
{
_logger.LogInformation("Processing rows...");
_ = Directory.CreateDirectory(Path.GetDirectoryName(Settings.Values.Files.Extraction.Path));
using var fileWriter = new StreamWriter(Settings.Values.Files.Extraction.Path);
if (Settings.Values.Files.Extraction.HasHeader)
{
string columnNames = string.Join(Settings.Values.Files.Extraction.Delimiter, rows.Columns.Select(c => c.Name));
fileWriter.WriteLine(columnNames);
}
foreach (Row row in rows)
{
CField[] rowFields = new CField[row.Length];
for (int i = 0; i < row.Length; i++)
{
rowFields[i] = rows.Columns[i].Type.IsAssignableFrom(typeof(DateTimeOffset))
? new CField(((DateTimeOffset)row[i]).ToUnixTimeMilliseconds(), rows.Columns[i].Name, typeof(long))
: new CField(row[i], rows.Columns[i].Name, rows.Columns[i].Type);
}
string rowToWrite = PrepareRowToBeWritten(rowFields);
fileWriter.WriteLine(rowToWrite);
}
}
public Field()
{
FieldData = new CField();
x_size = 7;
y_size = 7;
NotDoingStuff = true;
}
public void Visit(CField field)
{
if (field.FieldType == null)
{
throw new ArgumentException("FieldType not set");
}
var readOnly = field.IsReadonly ? "readonly " : string.Empty;
var isConst = field.IsConst ? "const " : string.Empty;
var staticTest = field.IsStatic ? "static " : string.Empty;
if (!string.IsNullOrEmpty(field.Comment))
{
_codeWriter.WriteLine("/// <summary>");
_codeWriter.WriteLine($"///{field.Comment}");
_codeWriter.WriteLine("/// </summary>");
}
_codeWriter.Write(
$"{field.AccessModifier.GetString()} {staticTest}{readOnly}{isConst}{field.FieldType} {field.FieldName}");
if (field.DefaultValue != null)
{
_codeWriter.Write($" = {field.DefaultValue}");
}
_codeWriter.WriteLine(";");
}
public static void BroadcastDojoEffect(CField field, DojoEffect nEffectType)
{
string sName;
switch (nEffectType)
{
case DojoEffect.Clear:
sName = "dojang/end/clear";
break;
case DojoEffect.TimeOver:
sName = "dojang/timeOver";
break;
default:
throw new InvalidOperationException($"Unhandled dojo map effect ({nameof(nEffectType)}).");
}
field.Broadcast(
new FieldEffectPacket(FieldEffect.Screen)
{
sName = sName
}
.GenerateForBroadcast());
}
public static void BroadcastCakeVsPieEffect(CField field, CakeEvent nEffectType)
{
string sName;
switch (nEffectType)
{
case CakeEvent.CakeWin:
sName = "event/5th/cakewin";
break;
case CakeEvent.PieWin:
sName = "event/5th/piewin";
break;
case CakeEvent.Start:
sName = "event/5th/start";
break;
case CakeEvent.TimeOver:
sName = "event/5th/timeover";
break;
default:
throw new InvalidOperationException($"Unhandled cake vs pie map effect ({nameof(nEffectType)}).");
}
field.Broadcast(
new FieldEffectPacket(FieldEffect.Screen)
{
sName = sName
}
.GenerateForBroadcast());
}
//-----------------------------------------------------------------------------
public CField GetField(int id)
{
if (!CFieldMan.ContainsKey(id))
{
var field = CField.Load(id, ParentServer);
CFieldMan.Add(id, field);
}
return(CFieldMan[id]);
}
public void ParseCfg(CIniFile config, string[] separators, CIniFile iniFileInstance)
{
string[] separeted_string = null;
CView lastView = null;
CField lastField = null;
CSignal lastSignal = null;
CSigMultitext lastSigMultitext = null;
string[] lines = System.IO.File.ReadAllLines(CfgPath, System.Text.Encoding.Default);
for (int i = 0; i < lines.Length; i++)
{
separeted_string = lines[i].Split(separators, StringSplitOptions.RemoveEmptyEntries);
try
{
if (!(lines[i].StartsWith("#")) && (lines[i].Length != 0))
{
switch (separeted_string[0])
{
case "View":
separeted_string = lines[i].Split(Const.separators_view, StringSplitOptions.RemoveEmptyEntries);
lastView = parseView(config, separeted_string);
iniFileInstance.AddView(lastView);
break;
case "Field":
lastField = parseField(config, separeted_string);
lastView.AddField(lastField);
break;
case "Signal":
separeted_string = lines[i].Split(Const.separators_signal, StringSplitOptions.RemoveEmptyEntries);
lastSignal = parseSignal(config, separeted_string);
lastField.AddSignal(lastSignal);
break;
case "SigMultitext":
separeted_string = lines[i].Split(Const.separators_signal, StringSplitOptions.RemoveEmptyEntries);
lastSigMultitext = parseSigMultitext(config, separeted_string);
lastField.AddSignalMultitext(lastSigMultitext);
break;
}
}
}
catch (Exception e) {
Error.toFile(e.Message.ToString(), this.GetType().Name.ToString());
}
}
}
public Field(bool isPuzzle,
string seed)
{
NotDoingStuff = true;
if (isPuzzle)
{
var puzzle = Puzzle.Decode(seed);
FieldData = puzzle.Field;
AuthorName = puzzle.AuthorName;
PuzzleName = puzzle.Name;
ClearConditions = puzzle.ClearConditions;
AvailableMoves = puzzle.AvailableMoves;
PuzzleVersion = puzzle.Version;
x_size = 7;
y_size = 7;
}
}
/// <summary>
/// 境界上の節点番号と座標の取得
/// </summary>
/// <param name="world">ワールド座標系</param>
/// <param name="fieldValId">フィールド値ID</param>
/// <param name="no_c_all">全体節点番号配列</param>
/// <param name="to_no_boundary">全体節点番号→境界上の節点番号のマップ</param>
/// <param name="coord_c_all">座標リスト</param>
/// <returns></returns>
public static bool GetBoundaryCoordList(
CFieldWorld world, uint fieldValId,
double rotAngle, double[] rotOrigin,
out uint[] no_c_all, out Dictionary <uint, uint> to_no_boundary, out double[][] coord_c_all)
{
no_c_all = null;
to_no_boundary = null;
coord_c_all = null;
// フィールドを取得
CField valField = world.GetField(fieldValId);
if (valField.GetFieldType() != FIELD_TYPE.ZSCALAR)
{
return(false);
}
IList <uint> aIdEA = valField.GetAryIdEA();
// 全体節点番号→境界節点番号変換テーブル(to_no_boundary)作成
to_no_boundary = new Dictionary <uint, uint>();
IList <double[]> coord_c_list = new List <double[]>();
foreach (uint eaId in aIdEA)
{
bool res = getBoundaryCoordList_EachElementAry(
world, fieldValId, eaId,
rotAngle, rotOrigin,
ref to_no_boundary, ref coord_c_list);
if (!res)
{
return(false);
}
}
//境界節点番号→全体節点番号変換テーブル(no_c_all)作成
no_c_all = new uint[to_no_boundary.Count];
foreach (KeyValuePair <uint, uint> kvp in to_no_boundary)
{
uint ino_boundary = kvp.Value;
uint ino = kvp.Key;
no_c_all[ino_boundary] = ino;
}
coord_c_all = coord_c_list.ToArray();
return(true);
}
private void btnClick(object sender, EventArgs e)
{
int iPlayer = 0;
MouseButtons p = (e as MouseEventArgs).Button;
if ((p & MouseButtons.Middle) != 0)
{
iPlayer = 1;
}
CField pTag = (sender as Control).Tag as CField;
if (pTag == null)
{
return;
}
m_pGame.Click(pTag.X, pTag.Y, -1);
h_Refresh();
}
/// <summary>
/// 境界上の節点番号の取得(要素アレイ単位)
/// </summary>
/// <param name="world">ワールド座標系</param>
/// <param name="fieldValId">フィールド値ID</param>
/// <param name="eaId">要素アレイID</param>
/// <param name="to_no_boundary">全体節点番号→境界上節点番号マップ</param>
/// <returns></returns>
private static bool getBoundaryNodeList_EachElementAry(CFieldWorld world, uint fieldValId, uint eaId, Dictionary <uint, uint> to_no_boundary)
{
// 要素アレイを取得する
CElemAry ea = world.GetEA(eaId);
System.Diagnostics.Debug.Assert(ea.ElemType() == ELEM_TYPE.LINE);
if (ea.ElemType() != ELEM_TYPE.LINE)
{
return(false);
}
// フィールドを取得
CField valField = world.GetField(fieldValId);
// 座標セグメントを取得
CElemAry.CElemSeg es_c_co = valField.GetElemSeg(eaId, ELSEG_TYPE.CORNER, false, world);
//境界節点番号→全体節点番号変換テーブル(no_c_all)作成
// 全体節点番号→境界節点番号変換テーブル(to_no_boundary)作成
uint node_cnt = ea.Size() + 1; // 全節点数
// 線要素の節点数
uint nno = 2;
uint[] no_c = new uint[nno]; // 要素節点の全体節点番号
for (uint ielem = 0; ielem < ea.Size(); ielem++)
{
// 要素配列から要素セグメントの節点番号を取り出す
es_c_co.GetNodes(ielem, no_c);
for (uint ino = 0; ino < nno; ino++)
{
if (!to_no_boundary.ContainsKey(no_c[ino]))
{
uint ino_boundary_tmp = (uint)to_no_boundary.Count;
to_no_boundary[no_c[ino]] = ino_boundary_tmp;
}
}
}
return(true);
}
private void h_DrawField(CField pField)
{
Button pButton = null;
foreach (object pC in tabPage2.Controls)
{
if (!(pC is Button))
{
continue;
}
CField pTag = (pC as Button).Tag as CField;
if (pTag == null)
{
continue;
}
if (pField.X == pTag.X && pField.Y == pTag.Y)
{
pButton = pC as Button;
}
}
if (pButton == null)
{
const int Size = 30;
const int Padding = 10;
pButton = new Button();
pButton.Parent = tabPage2;
pButton.Size = new Size(Size, Size);
pButton.Left = Padding + pField.X * (Size + Padding);
pButton.Top = Padding + pField.Y * (Size + Padding);
pButton.Tag = pField;
pButton.Click += btnClick;
}
pButton.Text =
pField.Owner == null
? " "
: pField.Owner.Symbol.ToString();
}
/// <summary>
/// Sets the target Y pos to the foothold below the Position.X and StartPosY values
/// </summary>
/// <param name="field">Field that the drop will calculated the foothold in</param>
public void CalculateY(CField field, short nDefaultY)
{
Foothold fh = null;
for (var i = 100; i <= 1000; i += 100)
{
fh = field.Footholds.FindBelow(Position.X, (short)(StartPosY - i));
if (fh != null)
{
break;
}
}
if (fh is null)
{
Position.Y = nDefaultY;
}
else
{
Position.Y = fh.Y1;
}
}
public const double eps0 = 8.85418782e-12;//1.0 / (mu0 * c0 * c0);
/// <summary>
/// 境界上の節点番号の取得
/// </summary>
/// <param name="world">ワールド座標系</param>
/// <param name="fieldValId">フィールド値ID</param>
/// <param name="no_c_all">全体節点番号配列</param>
/// <param name="to_no_boundary">全体節点番号→境界上の節点番号のマップ</param>
/// <returns></returns>
public static bool GetBoundaryNodeList(CFieldWorld world, uint fieldValId, out uint[] no_c_all, out Dictionary <uint, uint> to_no_boundary)
{
no_c_all = null;
to_no_boundary = null;
// フィールドを取得
CField valField = world.GetField(fieldValId);
if (valField.GetFieldType() != FIELD_TYPE.ZSCALAR)
{
return(false);
}
IList <uint> aIdEA = valField.GetAryIdEA();
// 全体節点番号→境界節点番号変換テーブル(to_no_boundary)作成
to_no_boundary = new Dictionary <uint, uint>();
foreach (uint eaId in aIdEA)
{
bool res = getBoundaryNodeList_EachElementAry(world, fieldValId, eaId, to_no_boundary);
if (!res)
{
return(false);
}
}
//境界節点番号→全体節点番号変換テーブル(no_c_all)作成
no_c_all = new uint[to_no_boundary.Count];
foreach (KeyValuePair <uint, uint> kvp in to_no_boundary)
{
uint ino_boundary = kvp.Value;
uint ino = kvp.Key;
no_c_all[ino_boundary] = ino;
}
return(true);
}
public CSwipe(CField field, CIcon selected, CIcon targeted)
{
this.field = field;
this.selected = selected;
this.targeted = targeted;
}
public CUserPool(CField parentField)
: base(parentField)
{
}
public CNpcPool(CField parentField)
: base(parentField)
{
}
public CReactorPool(CField parentField)
: base(parentField)
{
}
public CDropPool(CField parentField)
: base(parentField)
{
}
public BlowWeather(CField parent)
{
Field = parent;
StartTime = DateTime.Now.AddMinutes(-3);
}
private CField LoadNew(int nInstanceId)
{
var template = MasterManager.MapTemplates[TemplateID];
CField entry;
switch (template.FieldType)
{
case FieldType.WAITINGROOM:
entry = new CField_WaitingRoom(Parent, TemplateID, nInstanceId);
break;
case FieldType.CUSTOM_AREA_BOSS:
entry = new CField_AreaBoss(Parent, TemplateID, nInstanceId);
break;
case FieldType.SNOWBALL:
if (nInstanceId == 0)
{
goto default;
}
entry = new CField_Snowball(Parent, TemplateID, nInstanceId);
break;
case FieldType.CUSTOM_PINKBEAN:
if (nInstanceId == 0)
{
goto default;
}
entry = new CField_PinkBean(Parent, TemplateID, nInstanceId);
break;
case FieldType.HONTAIL:
if (nInstanceId == 0)
{
goto default;
}
entry = new CField_Hontail(Parent, TemplateID, nInstanceId);
break;
case FieldType.ZAKUM:
case FieldType.CHAOSZAKUM:
if (nInstanceId == 0)
{
goto default;
}
entry = new CField_Zakum(Parent, TemplateID, nInstanceId);
break;
case FieldType.DOJANG:
if (nInstanceId == 0)
{
goto default;
}
entry = new CField_Dojang(Parent, TemplateID, nInstanceId);
break;
case FieldType.CUSTOM_DUNGEONRAID_GOLEMTEMPLE:
if (nInstanceId == 0)
{
goto default;
}
entry = new CField_DungeonRaid_GolemTemple(Parent, TemplateID, nInstanceId);
break;
case FieldType.CUSTOM_COMMUNITYEVENT_RUSSIANROULETTE:
if (nInstanceId == 0)
{
goto default;
}
entry = new CField_RussianRoulette(Parent, TemplateID, nInstanceId);
break;
case FieldType.PERSONALTIMELIMIT:
if (nInstanceId == 0)
{
goto default;
}
entry = new CField_PersonalTimeLimit(Parent, TemplateID, nInstanceId);
break;
case FieldType.BATTLEFIELD:
if (nInstanceId == 0)
{
goto default;
}
entry = new CField_Battlefield(Parent, TemplateID, nInstanceId);
break;
default:
entry = new CField(Parent, TemplateID, nInstanceId);
break;
}
entry.Portals.Load(template);
entry.Footholds.Load(template);
entry.Mobs.Load(TemplateID);
entry.Npcs.Load(TemplateID);
entry.Reactors.Load(TemplateID);
return(entry);
}
/// <summary>
/// ループ内の節点番号と座標の取得(要素アレイ単位)
/// </summary>
/// <param name="world">ワールド座標系</param>
/// <param name="fieldValId">フィールド値ID</param>
/// <param name="eaId">要素アレイID</param>
/// <param name="to_no_boundary">全体節点番号→境界上節点番号マップ</param>
/// <returns></returns>
private static bool getLoopCoordList_EachElementAry(
CFieldWorld world, uint fieldValId, uint eaId,
double rotAngle, double[] rotOrigin,
ref Dictionary <uint, uint> to_no_loop, ref IList <double[]> coord_c_list, ref IList <uint[]> elem_no_c_list, ref IList <uint> elem_loopId_list)
{
// 要素アレイを取得する
CElemAry ea = world.GetEA(eaId);
System.Diagnostics.Debug.Assert(ea.ElemType() == ELEM_TYPE.TRI);
if (ea.ElemType() != ELEM_TYPE.TRI)
{
return(false);
}
// フィールドを取得
CField valField = world.GetField(fieldValId);
// 座標セグメントを取得
CElemAry.CElemSeg es_c_co = valField.GetElemSeg(eaId, ELSEG_TYPE.CORNER, false, world);
// 座標のノードセグメントを取得
CNodeAry.CNodeSeg ns_c_co = valField.GetNodeSeg(ELSEG_TYPE.CORNER, false, world);
uint node_cnt = ea.Size() + 1; // 全節点数
// 三角形要素の節点数
uint nno = 3;
// 座標の次元
uint ndim = 2;
// 要素節点の全体節点番号
uint[] no_c = new uint[nno];
// 要素節点の座標
double[][] coord_c = new double[nno][];
for (int inoes = 0; inoes < nno; inoes++)
{
coord_c[inoes] = new double[ndim];
}
for (uint ielem = 0; ielem < ea.Size(); ielem++)
{
// 要素配列から要素セグメントの節点番号を取り出す
es_c_co.GetNodes(ielem, no_c);
// ループID
elem_loopId_list.Add(eaId);
// 要素内節点番号
{
uint[] work_no_c = new uint[nno];
no_c.CopyTo(work_no_c, 0);
elem_no_c_list.Add(work_no_c);
}
// 座標を取得
for (uint inoes = 0; inoes < nno; inoes++)
{
double[] tmpval = null;
ns_c_co.GetValue(no_c[inoes], out tmpval);
System.Diagnostics.Debug.Assert(tmpval.Length == ndim);
for (int i = 0; i < tmpval.Length; i++)
{
coord_c[inoes][i] = tmpval[i];
}
}
if (Math.Abs(rotAngle) >= Constants.PrecisionLowerLimit)
{
// 座標を回転移動する
for (uint inoes = 0; inoes < nno; inoes++)
{
double[] srcPt = coord_c[inoes];
double[] destPt = GetRotCoord(srcPt, rotAngle, rotOrigin);
for (int i = 0; i < ndim; i++)
{
coord_c[inoes][i] = destPt[i];
}
}
}
for (uint ino = 0; ino < nno; ino++)
{
if (!to_no_loop.ContainsKey(no_c[ino]))
{
uint ino_loop_tmp = (uint)to_no_loop.Count;
to_no_loop[no_c[ino]] = ino_loop_tmp;
coord_c_list.Add(new double[] { coord_c[ino][0], coord_c[ino][1] });
System.Diagnostics.Debug.Assert(coord_c_list.Count == (ino_loop_tmp + 1));
}
}
}
return(true);
}
public CTownPortalPool(CField parent)
: base(parent)
{
}
public CSwipeBack(CField field, CIcon selected, CIcon targeted) : base(field, selected, targeted)
{
}
////////////////////////////////////////////////////////////////////////////////////////////////
// DekFEMの処理を隠ぺいした簡易版
////////////////////////////////////////////////////////////////////////////////////////////////
/// <summary>
/// PCOCGで線形方程式を解く
/// </summary>
/// <param name="mat"></param>
/// <param name="resVec"></param>
/// <param name="value_c_all"></param>
/// <param name="isConverged"></param>
/// <returns></returns>
//public static bool SolvePCOCG(ref MyUtilLib.Matrix.MyComplexMatrix mat, ref System.Numerics.Complex[] resVec, out System.Numerics.Complex[] value_c_all, out bool isConverged)
public static bool SolvePCOCG(ref MyUtilLib.Matrix.MyComplexMatrix mat, ref KrdLab.clapack.Complex[] resVec, out KrdLab.clapack.Complex[] value_c_all, out bool isConverged)
{
bool success = false;
int matLen = mat.RowSize;
System.Diagnostics.Debug.WriteLine("SolvePCOCG 1");
//------------------------------------------------------------------
// ワールド座標系を生成
//------------------------------------------------------------------
uint baseId = 0;
CFieldWorld World = new CFieldWorld();
setupWorld((uint)matLen, ref World, out baseId);
// 界の値を扱うバッファ？を生成する。
// フィールド値IDが返却される。
// 要素の次元: 2次元-->ポイントにしたので0次元？ 界: 複素数スカラー 微分タイプ: 値 要素セグメント: 角節点
uint FieldValId = World.MakeField_FieldElemDim(baseId, 0,
FIELD_TYPE.ZSCALAR, FIELD_DERIVATION_TYPE.VALUE, ELSEG_TYPE.CORNER);
System.Diagnostics.Debug.WriteLine("SolvePCOCG 2");
//------------------------------------------------------------------
// リニアシステム
//------------------------------------------------------------------
CZLinearSystem Ls = new CZLinearSystem();
CZPreconditioner_ILU Prec = new CZPreconditioner_ILU();
//------------------------------------------------------------------
// 界パターン追加
//------------------------------------------------------------------
Ls.AddPattern_Field(FieldValId, World);
// POINTのフィールドなので対角要素のパターンのみ追加される。非対角要素のパターンは追加されない
// 非対角要素のパターンの追加
{
uint node_cnt = (uint)matLen;
// 要素剛性行列(コーナ-コーナー)
CZMatDia_BlkCrs_Ptr mat_cc = Ls.GetMatrixPtr(FieldValId, ELSEG_TYPE.CORNER, World);
// 要素残差ベクトル(コーナー)
CZVector_Blk_Ptr res_c = Ls.GetResidualPtr(FieldValId, ELSEG_TYPE.CORNER, World);
// パターンを追加
using (CIndexedArray crs = new CIndexedArray())
{
crs.InitializeSize(mat_cc.NBlkMatCol());
//crs.Fill(mat_cc.NBlkMatCol(), mat_cc.NBlkMatRow());
using (UIntVectorIndexer index = crs.index)
using (UIntVectorIndexer ary = crs.array)
{
for (int iblk = 0; iblk < (int)crs.Size(); iblk++)
{
index[iblk] = 0;
}
for (int iblk = 0; iblk < mat_cc.NBlkMatCol(); iblk++)
{
// 現在のマトリクスのインデックス設定をコピーする
uint npsup = 0;
ConstUIntArrayIndexer cur_rows = mat_cc.GetPtrIndPSuP((uint)iblk, out npsup);
foreach (uint row_index in cur_rows)
{
ary.Add(row_index);
}
index[iblk + 1] = (uint)ary.Count;
//System.Diagnostics.Debug.WriteLine("chk:{0} {1}", iblk, npsup);
// Note: インデックス設定はされていないので常にnpsup == 0 (cur_rows.Count == 0)
// すべての節点に関して列を追加
int ino_global = iblk;
int col = iblk;
if (col != -1 && ino_global != -1)
{
// 関連付けられていない節点をその行の列データの最後に追加
int last_index = (int)index[col + 1] - 1;
int add_cnt = 0;
for (int jno_global = 0; jno_global < node_cnt; jno_global++)
{
// 行列の０要素を除く
//if (System.Numerics.Complex.Abs(mat[(int)ino_global, (int)jno_global]) < Constants.PrecisionLowerLimit)
//if (System.Numerics.Complex.Abs(mat[(int)ino_global, (int)jno_global]) < 1.0e-15)
//if (mat._body[ino_global + jno_global * mat.RowSize].Magnitude < Constants.PrecisionLowerLimit)
KrdLab.clapack.Complex cvalue = mat._body[ino_global + jno_global * mat.RowSize];
if (cvalue.Real == 0 && cvalue.Imaginary == 0)
{
continue;
}
uint row = (uint)jno_global;
if (ino_global != jno_global) // 対角要素は除く
{
if (!cur_rows.Contains(row))
{
ary.Insert(last_index + 1 + add_cnt, row);
add_cnt++;
//System.Diagnostics.Debug.WriteLine("added:" + col + " " + row);
}
}
}
if (add_cnt > 0)
{
index[col + 1] = (uint)ary.Count;
}
}
}
}
crs.Sort();
System.Diagnostics.Debug.Assert(crs.CheckValid());
// パターンを削除する
mat_cc.DeletePattern();
// パターンを追加する
mat_cc.AddPattern(crs);
}
}
System.Diagnostics.Debug.WriteLine("SolvePCOCG 3");
//------------------------------------------------------------------
// プリコンディショナ―
//------------------------------------------------------------------
// set Preconditioner
//Prec.SetFillInLevel(1);
// フィルインなしで不完全LU分解した方が収束が早く、また解もそれなりの結果だったので０に設定した
//Prec.SetFillInLevel(0);
// フィルインなしでは、収束しない場合がある。フィルインレベルを最大にする
//Prec.SetFillInLevel((uint)mat.RowSize);
// 完全LU分解（だと思われる) Prec.SetLinearSystemの処理が遅い
//Prec.SetFillInLevel(3);
// 誘電体スラブ導波路だと最低このレベル。これでも収束しない場合がある
//以上から
//メインの導波管の計算時間を短くしたいのでフィルインなしにする
//Prec.SetFillInLevel(0);
// 導波管だとこれで十分
// フィルインレベルの設定(既定値)
Prec.SetFillInLevel(DefFillInLevel);
// ILU(0)のパターン初期化
Prec.SetLinearSystem(Ls);
System.Diagnostics.Debug.WriteLine("SolvePCOCG 4");
//------------------------------------------------------------------
// 剛性行列、残差ベクトルのマージ
//------------------------------------------------------------------
Ls.InitializeMarge();
{
uint ntmp = Ls.GetTmpBufferSize();
int[] tmpBuffer = new int[ntmp];
for (int i = 0; i < ntmp; i++)
{
tmpBuffer[i] = -1;
}
uint node_cnt = (uint)matLen;
// 要素剛性行列(コーナ-コーナー)
CZMatDia_BlkCrs_Ptr mat_cc = Ls.GetMatrixPtr(FieldValId, ELSEG_TYPE.CORNER, World);
// 要素残差ベクトル(コーナー)
CZVector_Blk_Ptr res_c = Ls.GetResidualPtr(FieldValId, ELSEG_TYPE.CORNER, World);
bool[,] add_flg = new bool[node_cnt, node_cnt];
for (int i = 0; i < node_cnt; i++)
{
for (int j = 0; j < node_cnt; j++)
{
add_flg[i, j] = false;
// ０要素はマージ対象でないので最初から除外する
//if (System.Numerics.Complex.Abs(mat[i, j]) < Constants.PrecisionLowerLimit)
//if (mat._body[i + j * mat.RowSize].Magnitude < Constants.PrecisionLowerLimit)
KrdLab.clapack.Complex cvalue = mat._body[i + j * mat.RowSize];
if (cvalue.Real == 0 && cvalue.Imaginary == 0)
{
add_flg[i, j] = true;
}
}
}
for (int iblk = 0; iblk < mat_cc.NBlkMatCol(); iblk++)
{
// 自分及び自分自身と関連のある"row"の要素をマージする( 1 x rowcntのベクトル : 行列の横ベクトルに対応)
uint npsup = 0;
ConstUIntArrayIndexer cur_rows = mat_cc.GetPtrIndPSuP((uint)iblk, out npsup);
uint colcnt = 1;
uint rowcnt = (uint)cur_rows.Count + 1; // cur_rowsには自分自身は含まれないので+1する
Complex[] emattmp = new Complex[colcnt * rowcnt];
uint[] no_c_tmp_col = new uint[colcnt];
uint[] no_c_tmp_row = new uint[rowcnt];
no_c_tmp_col[0] = (uint)iblk;
no_c_tmp_row[0] = (uint)iblk;
for (int irow = 0; irow < cur_rows.Count; irow++)
{
no_c_tmp_row[irow + 1] = cur_rows[irow];
}
{
uint ino = 0;
uint ino_global = no_c_tmp_col[ino];
for (int jno = 0; jno < rowcnt; jno++)
{
uint jno_global = no_c_tmp_row[jno];
//emat[ino, jno]
if (!add_flg[ino_global, jno_global])
{
//System.Numerics.Complex cvalue = (System.Numerics.Complex)mat[ino_global, jno_global];
//System.Numerics.Complex cvalue = (System.Numerics.Complex)mat._body[ino_global + jno_global * mat.RowSize];
KrdLab.clapack.Complex cvalue = mat._body[ino_global + jno_global * mat.RowSize];
emattmp[ino * rowcnt + jno] = new Complex(cvalue.Real, cvalue.Imaginary);
add_flg[ino_global, jno_global] = true;
}
else
{
// ここにはこない
System.Diagnostics.Debug.Assert(false);
emattmp[ino * rowcnt + jno] = new Complex(0, 0);
}
}
}
// マージ!
mat_cc.Mearge(1, no_c_tmp_col, rowcnt, no_c_tmp_row, 1, emattmp, ref tmpBuffer);
}
for (int i = 0; i < node_cnt; i++)
{
for (int j = 0; j < node_cnt; j++)
{
//System.Diagnostics.Debug.WriteLine( i + " " + j + " " + add_flg[i, j] );
System.Diagnostics.Debug.Assert(add_flg[i, j]);
}
}
// 残差ベクトルにマージ
for (uint ino_global = 0; ino_global < node_cnt; ino_global++)
{
// 残差ベクトルにマージする
uint no_tmp = ino_global;
//System.Numerics.Complex cvalue = resVec[ino_global];
KrdLab.clapack.Complex cvalue = resVec[ino_global];
Complex val = new Complex(cvalue.Real, cvalue.Imaginary);
res_c.AddValue(no_tmp, 0, val);
}
}
double res = Ls.FinalizeMarge();
//------------------------------------------------------------------
// リニアシステムを解く
//------------------------------------------------------------------
System.Diagnostics.Debug.WriteLine("Solve_PCOCG 5: CZSolverLsIter.Solve_PCOCG");
// プリコンディショナに値を設定してILU分解を行う
Prec.SetValue(Ls);
//double tol = 1.0e-1;
//double tol = 1.0e-2;
//double tol = 0.5e-3;
//double tol = 1.0e-4;
double tol = 1.0e-6;
//uint maxIter = 500;
//uint maxIter = 1000;
//uint maxIter = 2000;
//uint maxIter = 5000;
//uint maxIter = 10000;
uint maxIter = 20000;
//uint maxIter = uint.MaxValue;
uint iter = maxIter;
bool ret = CZSolverLsIter.Solve_PCOCG(ref tol, ref iter, Ls, Prec);
System.Diagnostics.Debug.WriteLine("iter : " + iter + " Res : " + tol + " ret : " + ret);
if (iter == maxIter)
{
System.Diagnostics.Debug.WriteLine("Not converged");
isConverged = false;
}
else
{
isConverged = true;
}
System.Diagnostics.Debug.WriteLine("Solve_PCOCG 6");
//------------------------------------------------------------------
// 計算結果の後処理
//------------------------------------------------------------------
// 計算結果をワールド座標系に反映する
Ls.UpdateValueOfField(FieldValId, World, FIELD_DERIVATION_TYPE.VALUE);
{
uint node_cnt = (uint)matLen;
//value_c_all = new System.Numerics.Complex[node_cnt]; // コーナーの値
value_c_all = new KrdLab.clapack.Complex[node_cnt]; // コーナーの値
CField valField = World.GetField(FieldValId);
// 要素アレイIDのリストを取得
IList <uint> aIdEA = valField.GetAryIdEA();
// 要素アレイを走査
foreach (uint eaId in aIdEA)
{
CElemAry ea = World.GetEA(eaId);
CElemAry.CElemSeg es_c_va = valField.GetElemSeg(eaId, ELSEG_TYPE.CORNER, true, World);
CElemAry.CElemSeg es_c_co = valField.GetElemSeg(eaId, ELSEG_TYPE.CORNER, false, World);
// 要素の節点数
uint nno = 1; //点
// 要素節点の全体節点番号
uint[] no_c = new uint[nno];
// 要素節点の値
Complex[] value_c = new Complex[nno];
// 節点の値の節点セグメント
CNodeAry.CNodeSeg ns_c_val = valField.GetNodeSeg(ELSEG_TYPE.CORNER, true, World);
for (uint ielem = 0; ielem < ea.Size(); ielem++)
{
// 要素配列から要素セグメントの節点番号を取り出す
es_c_co.GetNodes(ielem, no_c);
// 節点の値を取って来る
es_c_va.GetNodes(ielem, no_c);
for (uint inoes = 0; inoes < nno; inoes++)
{
Complex[] tmpval = null;
ns_c_val.GetValue(no_c[inoes], out tmpval);
System.Diagnostics.Debug.Assert(tmpval.Length == 1);
value_c[inoes] = tmpval[0];
//System.Diagnostics.Debug.Write( "(" + value_c[inoes].Real + "," + value_c[inoes].Imag + ")" );
}
for (uint ino = 0; ino < nno; ino++)
{
//DelFEMの節点番号は０開始
//uint ino_global = no_c[ino] - 1;
uint ino_global = no_c[ino];
Complex cvalue = value_c[ino];
//value_c_all[ino_global] = new KrdLab.clapack.Complex(cvalue.Real, cvalue.Imag);
value_c_all[ino_global].Real = cvalue.Real;
value_c_all[ino_global].Imaginary = cvalue.Imag;
}
}
}
}
Ls.Clear();
Ls.Dispose();
Prec.Clear();
Prec.Dispose();
World.Clear();
World.Dispose();
System.Diagnostics.Debug.WriteLine("Solve_PCOCG End");
return(success);
}
public void Visit(CField field)
{
_fieldVisitor.Visit(field);
}
public CFieldItem(CField field)
{
Field = field;
LastUpdate = DateTime.Now;
}
private void GenerateDef(CDBM db, CSharpProject p, string namespaceName, TextWriter mensajes)
{
string className = DefName(db.Name);
p.AddCSharpFile(className + ".generated", namespaceName);
CClass c = p.AddClass(className, isStatic: true);
c.lUsing.Add(CSharpStandarUsing.System);
c.lUsing.Add(new CUsing("IvanCruz.Util.UtilBD"));
foreach (CTableM t in db.lTable)
{
c.AddField(new CField(t.Name, new TextualType(ClassDefName(t.Name)), CSharpVisibility.cvPublic, isStatic: true, InitialValue: " new " + ClassDefName(t.Name) + "()"));
//CConstructor consVacio = new CConstructor(c, CSharpVisibility.cvPublic);
//CConstructor cons = new CConstructor(c, CSharpVisibility.cvPublic);
////Se añaden al final por estética de la clase generada
//foreach (CFieldM f in t.lFieldAll) {
// CType aux = CGenerator.SqlServerTD2CSharp(f.Td, f.IsNullable);
// if (aux.Equals(CSharpPrimitiveType.cDateTime)) {
// c.lUsing.Add(CSharpStandarUsing.System);
// }
// c.AddField(new CField(f.Name, aux, CSharpVisibility.cvPublic));
// cons.AddParam(SUtil.LowerFirst(f.Name), CGenerator.SqlServerTD2CSharp(f.Td, f.IsNullable));
// cons.AddSentence(new CTextualSentence($"this.{f.Name} = {SUtil.LowerFirst(f.Name)};"));
//}
//c.AddConstructor(consVacio);
//c.AddConstructor(cons);
CClass cTable = p.AddClass(ClassDefName(t.Name), isStatic: false);
cTable.ParentClassName = "CTableDef";
CProperty prop;
prop = new CProperty("TableName", CSharpVisibility.cvPublic, new TextualType("string"), hasGet: true, hasSet: false);
prop.Override = true;
prop.lSentenceGet.Add(new CTextualSentence("return " + SUtil.DoubleQuote(t.Name) + ";"));
cTable.AddProperty(prop);
prop = new CProperty("SingularTitle", CSharpVisibility.cvPublic, new TextualType("string"), hasGet: true, hasSet: false);
prop.Override = true;
prop.lSentenceGet.Add(new CTextualSentence("return " + SUtil.DoubleQuote(t.SingularTitle) + ";"));
cTable.AddProperty(prop);
CField auxField;
StringBuilder sb = new StringBuilder();
prop = new CProperty("Fields", CSharpVisibility.cvPublic, new TextualType("System.Collections.Generic.IEnumerable<CFieldDef>"), hasGet: true, hasSet: false);
prop.Override = true;
foreach (CFieldM f in t.lFieldAll)
{
auxField = new CField(f.Name, new TextualType("CFieldDef"), CSharpVisibility.cvPublic, isStatic: false);
sb.Clear();
sb.Append("new CFieldDef");
sb.Append("(");
sb.Append(SUtil.DoubleQuote(f.Name));
if (!string.IsNullOrWhiteSpace(f.Title))
{
sb.Append(",");
sb.Append(SUtil.DoubleQuote(f.Title));
}
sb.Append(")");
auxField.InitialValue = sb.ToString();
cTable.AddField(auxField);
prop.lSentenceGet.Add(new CTextualSentence("yield return " + f.Name + ";"));
}
cTable.AddProperty(prop);
//public override System.Collections.Generic.IEnumerable<CFieldDef> Fields {
// get {
// yield return Id;
// yield return CodigoMedico;
// yield return Nomb;
// }
//}
}
mensajes.WriteLine("Definición generada en: " + p.RootDirectory.ActualDirectory.FullName);
}
/// <summary>
///
/// </summary>
/// <param name="World"></param>
/// <param name="FieldValId"></param>
/// <param name="Ls"></param>
/// <param name="Prec"></param>
/// <param name="tmpBuffer"></param>
/// <param name="resVec"></param>
/// <param name="value_c_all"></param>
/// <param name="isConverged"></param>
/// <returns></returns>
public static bool SolvePCOCG(
ref CFieldWorld World, ref uint FieldValId,
ref CZLinearSystem Ls, ref CZPreconditioner_ILU Prec, ref int[] tmpBuffer,
ref KrdLab.clapack.Complex[] resVec, out KrdLab.clapack.Complex[] value_c_all, out bool isConverged)
{
bool success = false;
//------------------------------------------------------------------
// マージの終了
//------------------------------------------------------------------
double res = Ls.FinalizeMarge();
tmpBuffer = null;
// マトリクスサイズの取得
// 要素剛性行列(コーナ-コーナー)
CZMatDia_BlkCrs_Ptr mat_cc = Ls.GetMatrixPtr(FieldValId, ELSEG_TYPE.CORNER, World);
// マトリクスのサイズ
int matLen = (int)mat_cc.NBlkMatCol();
//------------------------------------------------------------------
// リニアシステムを解く
//------------------------------------------------------------------
System.Diagnostics.Debug.WriteLine("Solve_PCOCG : CZSolverLsIter.Solve_PCOCG");
// プリコンディショナに値を設定してILU分解を行う
Prec.SetValue(Ls);
//double tol = 1.0e-1;
//double tol = 1.0e-2;
//double tol = 0.5e-3;
//double tol = 1.0e-4;
double tol = 1.0e-6;
//uint maxIter = 500;
//uint maxIter = 1000;
//uint maxIter = 2000;
//uint maxIter = 5000;
//uint maxIter = 10000;
uint maxIter = 20000;
//uint maxIter = uint.MaxValue;
uint iter = maxIter;
bool ret = CZSolverLsIter.Solve_PCOCG(ref tol, ref iter, Ls, Prec);
System.Diagnostics.Debug.WriteLine("iter : " + iter + " Res : " + tol + " ret : " + ret);
if (iter == maxIter)
{
System.Diagnostics.Debug.WriteLine("Not converged");
isConverged = false;
}
else
{
isConverged = true;
}
System.Diagnostics.Debug.WriteLine("Solve_PCOCG solved");
//------------------------------------------------------------------
// 計算結果の後処理
//------------------------------------------------------------------
// 計算結果をワールド座標系に反映する
Ls.UpdateValueOfField(FieldValId, World, FIELD_DERIVATION_TYPE.VALUE);
{
uint node_cnt = (uint)matLen;
//value_c_all = new System.Numerics.Complex[node_cnt]; // コーナーの値
value_c_all = new KrdLab.clapack.Complex[node_cnt]; // コーナーの値
CField valField = World.GetField(FieldValId);
// 要素アレイIDのリストを取得
IList <uint> aIdEA = valField.GetAryIdEA();
// 要素アレイを走査
foreach (uint eaId in aIdEA)
{
CElemAry ea = World.GetEA(eaId);
CElemAry.CElemSeg es_c_va = valField.GetElemSeg(eaId, ELSEG_TYPE.CORNER, true, World);
CElemAry.CElemSeg es_c_co = valField.GetElemSeg(eaId, ELSEG_TYPE.CORNER, false, World);
// 要素の節点数
uint nno = 1; //点
// 要素節点の全体節点番号
uint[] no_c = new uint[nno];
// 要素節点の値
Complex[] value_c = new Complex[nno];
// 節点の値の節点セグメント
CNodeAry.CNodeSeg ns_c_val = valField.GetNodeSeg(ELSEG_TYPE.CORNER, true, World);
for (uint ielem = 0; ielem < ea.Size(); ielem++)
{
// 要素配列から要素セグメントの節点番号を取り出す
es_c_co.GetNodes(ielem, no_c);
// 節点の値を取って来る
es_c_va.GetNodes(ielem, no_c);
for (uint inoes = 0; inoes < nno; inoes++)
{
Complex[] tmpval = null;
ns_c_val.GetValue(no_c[inoes], out tmpval);
System.Diagnostics.Debug.Assert(tmpval.Length == 1);
value_c[inoes] = tmpval[0];
//System.Diagnostics.Debug.Write( "(" + value_c[inoes].Real + "," + value_c[inoes].Imag + ")" );
}
for (uint ino = 0; ino < nno; ino++)
{
//DelFEMの節点番号は０開始
//uint ino_global = no_c[ino] - 1;
uint ino_global = no_c[ino];
Complex cvalue = value_c[ino];
//value_c_all[ino_global] = new KrdLab.clapack.Complex(cvalue.Real, cvalue.Imag);
value_c_all[ino_global].Real = cvalue.Real;
value_c_all[ino_global].Imaginary = cvalue.Imag;
}
}
}
}
Ls.Clear();
Ls.Dispose();
Ls = null;
Prec.Clear();
Prec.Dispose();
Prec = null;
World.Clear();
World.Dispose();
World = null;
System.Diagnostics.Debug.WriteLine("Solve_PCOCG End");
return(success);
}
