/* -*- C++ -*- */

//=============================================================================
/**
 *  @文件    TinCore.h
 *
 *  TinCore.h, 版本 0.10 2007/09/27
 *
 *  @作者 鄧雪清 <xueqingdeng@sohu.com>
 *  基于Incremental Construction算法的單線程實現
 */
//=============================================================================

#ifndef _TIN_CORE_H_
#define _TIN_CORE_H_

#include <windows.h>
#include "..\Tss_Inc\Point2D.h"
#include "..\Tss_Inc\Point3D.h"
#include "..\Tss_Inc\TssTset.h"
#include "..\Tss_API\TssHash.h"
#include "..\Tss_API\TssList.h"
typedef	struct _timeb	TSS_TIMEB;

// 網格索引結構
typedef struct tagIDX_CELL {
	TSS_SNT32	 num;	// 容量
	TSS_SNT32	 cnt;	// 計數
	TSS_SNT32	*pid;	// 數組
	TSS_SNT32	 sum;	// 總數
} IDX_CELL;

// 統一網格結構
typedef struct tagIDX_GRID {
	TSS_FLT64	 bnd[4];	// 坐標邊界
	TSS_FLT64	 spn[2];	// 網格間距
	TSS_SNT32	 num[2];	// 矩陣大小
	TSS_FLT32	 dup;		// 坐標精度
	IDX_CELL	*idx;		// 單元數組
} IDX_GRID;

// 邊(面)結構
typedef struct tagTIN_FACE {
	TSS_SNT32	 sid[2];	// 起點(0:pid, 1:cid)
	TSS_SNT32	 eid[2];	// 末點(0:pid, 1:cid)
	TSS_SNT32	 pid;		// 對面點
	TSS_SNT32	 tid;		// 鄰接三角形
} TIN_FACE;

// 三角形結構
typedef struct tagTIN_TGON {
	TSS_SNT32	 pid[3];
	TSS_SNT32	 tid[3];
} TIN_TGON;

long GridToTin2D(POINT2D *ppts, TSS_SNT32 *pnum, TSS_SNT32 *stop,
			   TSS_SNT32 *vnum, TSS_SNT32 **vidx, TSS_TIMEB *betm, TSS_TIMEB *entm);
long DelaunayTriangulation2D(POINT3D *ppts, TSS_SNT32 *pnum, TSS_FLT32 ddup, TSS_SNT32 *stop,
							 TSS_SNT32 *vnum, TSS_SNT32 **vidx, TSS_TIMEB *betm, TSS_TIMEB *entm);

// 哈希函數
TSS_SNT32 FACE_EQLH(void const *f1, void const *f2);
TSS_UNT32 FACE_HASH(void const *f1);

// 計算網格矩陣大小
void CalcGridSize(IDX_GRID *uidx, TSS_SNT32 num, TSS_SNT32 avg);
// 計算點所處單元-數組偏移
long CalcGridCell(IDX_GRID *uidx, TSS_FLT64 *pnt);
// 計算點所處單元-坐標位置
void CalcGridCell(IDX_GRID *uidx, TSS_FLT64 *pnt, TSS_SNT32 *pos);
// 增加點到單元
void AddPid(IDX_CELL *cell, TSS_SNT32 pid);
// 從單元刪除點
void DelPid(IDX_CELL *cell, TSS_SNT32 pos);
void DelPidEx(IDX_CELL *cell, TSS_SNT32 pid);

// 判斷單元中是否有重合點
long IsHaveDup(IDX_CELL *cell, TSS_FLT64 *ppts, TSS_SNT32 ddim, TSS_FLT32 dup,
			   TSS_FLT64 *pt, TSS_SNT32 s);

// 消除重合點
void RemoveDupA(IDX_GRID *uidx, TSS_SNT32 *pnum, TSS_SNT32 ddim, TSS_FLT64 *ppts,
				TSS_SNT32 *ppid, TSS_SNT32 *stop);
// 按列處理重合點
void RemoveDupC(IDX_GRID *uidx, TSS_SNT32 *pnum, TSS_SNT32 ddim, TSS_FLT64 *ppts,
				TSS_SNT32 *ppid, TSS_SNT32 c, TSS_SNT32 r1, TSS_SNT32 r2);
// 按行處理重合點
void RemoveDupR(IDX_GRID *uidx, TSS_SNT32 *pnum, TSS_SNT32 ddim, TSS_FLT64 *ppts,
				TSS_SNT32 *ppid, TSS_SNT32 r, TSS_SNT32 c1, TSS_SNT32 c2);
// 處理邊界重合點
void RemoveDupS(IDX_GRID *uidx, TSS_SNT32 *pnum, TSS_SNT32 ddim, TSS_FLT64 *ppts,
				TSS_SNT32 *ppid);
// 按矩陣處理重合點
void RemoveDupM(IDX_GRID *uidx, TSS_SNT32 *pnum, TSS_SNT32 ddim, TSS_FLT64 *ppts,
				TSS_SNT32 *ppid, TSS_SNT32 *stop, TSS_SNT32 c1, TSS_SNT32 c2,
				TSS_SNT32 r1, TSS_SNT32 r2);
// 設置邊參數
void SetFace(TIN_FACE *face, TSS_SNT32 sid, TSS_SNT32 eid, TSS_SNT32 pid,
			 TSS_SNT32 tid, TSS_SNT32 sic, TSS_SNT32 eic);
// 釋放CELL索引
void FreeCell(IDX_CELL *cidx, TSS_SNT32 num);

// 構建初始三角形
long MakeFirstSimplex(IDX_GRID *uidx, TSS_FLT64 *ppts, TSS_SNT32 ddim, TSS_SNT32 *stop,
					  TSS_FLT64 dm2, TSS_SNT32 *pid, TSS_SNT32 *cid);
// 構建三角形網絡
long MakeSimplex(IDX_GRID *uidx, TSS_FLT64 *ppts, TSS_SNT32 ddim, TSS_SNT32 *ppid,
				 TSS_SNT32 *pnum, TSS_SNT32 *vnum, TSS_SNT32 *stop, TSS_SNT32 *pid,
				 TSS_SNT32 *cid, CTssList *ATL);

// 單元搜索最近點
// 同單元搜索
void SearchCellNPS(IDX_CELL *cell, TSS_FLT64 *ppts, TSS_SNT32 ddim, TSS_FLT64 *pnt,
				   TSS_SNT32 cin, TSS_SNT32 &pid, TSS_SNT32 &cid, TSS_FLT64 &dnp);
// 異單元搜索
void SearchCellNPD(IDX_CELL *cell, TSS_FLT64 *ppts, TSS_SNT32 ddim, TSS_FLT64 *pnt,
				   TSS_SNT32 cin, TSS_SNT32 &pid, TSS_SNT32 &cid, TSS_FLT64 &dnp);
// 批量異單元搜索
void SearchCellNPD(IDX_CELL *cidx, TSS_FLT64 *ppts, TSS_SNT32 ddim, TSS_FLT64 *pnt,
				   TSS_SNT32 &pid, TSS_SNT32 &cid, TSS_FLT64 &dnp, TSS_SNT32 col,
				   TSS_SNT32 *c, TSS_SNT32 *r);

// 計算搜索邊界
void CalcSearchBound(TSS_FLT64 *bnd, TSS_FLT64 *spn, TSS_SNT32 col, TSS_SNT32 row,
					 TSS_FLT64 *pnt, TSS_FLT64 rad, TSS_SNT32 &c1, TSS_SNT32 &c2,
					 TSS_SNT32 &r1, TSS_SNT32 &r2);

// 搜索第3點:空圓法則
long SearchCellEMR(IDX_GRID *uidx, IDX_CELL *cell, TSS_FLT64 *ppts, TSS_SNT32 ddim,
				   TSS_FLT64 *pt1, TSS_FLT64 *pt2, TSS_SNT32 cin, TSS_SNT32 &pid,
				   TSS_SNT32 &cid, TSS_FLT64 *cen);
// 檢查圓內是否存在點
void SearchCellEMR(IDX_GRID *uidx, TSS_FLT64 *ppts, TSS_SNT32 ddim, TSS_FLT64 *pt1,
				   TSS_FLT64 *pt2, TSS_FLT64 *pt3, TSS_FLT64 rad, TSS_SNT32 &pid,
				   TSS_SNT32 &cid, TSS_FLT64 *cen);
// 批量異單元搜索
long SearchCellEMR(IDX_GRID *uidx, TSS_FLT64 *ppts, TSS_SNT32 ddim, TSS_FLT64 *pt1,
				   TSS_FLT64 *pt2, TSS_SNT32 &pid, TSS_SNT32 &cid, TSS_FLT64 *cen,
				   TSS_SNT32 col, TSS_SNT32 *c, TSS_SNT32 *r);

// 插入邊
long InsertFace(CTssHash *AFH, CTssHash *BFH, TSS_SNT32 *ppid, TIN_FACE *face);

// 搜索第3點:空圓法則
long SearchCell(CTssHash *BFH, IDX_GRID *uidx, TSS_FLT64 *ppts, TSS_SNT32 ddim, TSS_SNT32 *stop,
				TSS_SNT32 *pid,	TSS_SNT32 &cid);
// 方向約束搜索
long SearchCellEMR(CTssHash *BFH, IDX_GRID *uidx, IDX_CELL *cell, TSS_FLT64 *ppts, TSS_SNT32 ddim,
				   TSS_FLT64 *pt1, TSS_FLT64 *pt2, TSS_FLT64 *opt, TSS_SNT32 orn, TSS_SNT32 cin,
				   TSS_SNT32 *pid, TSS_SNT32 &cid, TSS_FLT64 *cen);
// 批量單元方向約束搜索
long SearchCellEMR(CTssHash *BFH, IDX_GRID *uidx, TSS_FLT64 *ppts, TSS_SNT32 ddim, TSS_FLT64 *pt1,
				   TSS_FLT64 *pt2, TSS_FLT64 *opt, TSS_SNT32 orn, TSS_SNT32 *pid, TSS_SNT32 &cid,
				   TSS_FLT64 *cen, TSS_SNT32 col, TSS_SNT32 *c, TSS_SNT32 *r);

// 判斷是否共線
long IsLine(TSS_FLT64 *pt1, TSS_FLT64 *pt2, TSS_FLT64 *pt3);
// 計算同側符號
long IsSide(TSS_FLT64 *pt1, TSS_FLT64 *pt2, TSS_FLT64 *pt3);

#endif // _TIN_CORE_H_