#include <iostream>

// メモリリークを検出した位置を表示するためのマクロ
#define _CRTDBG_MAP_ALLOC

// メモリリークの検出
#include <crtdbg.h>

// メモリリークしたメモリのnewした場所のファイル名と行数を表示する
// #define _CRTDBG_MAP_ALLOC で動く機能のはずなのだがなぜか機能しないので自分で定義する
#define new ::new(_NORMAL_BLOCK, __FILE__, __LINE__)

class Dummy
{
public:

	char* myleak() {
        char* p = new char[3] {0, 0, 0}; // メモリリークしている (19行目)
        return p;
	}

	void mysafety() {
        int* q = new int[3] {1, 1, 1}; // メモリリークしていない (24行目)
        printf("%d %d %d\n", q[0], q[1], q[2]);
        delete [] q;
	}

    Dummy() {}
    ~Dummy() {}
};

int main()
{
    // アプリケーション終了時に_CrtDumpMemoryLeaks()が呼ばれるように指示
    _CrtSetDbgFlag(_CRTDBG_ALLOC_MEM_DF | _CRTDBG_LEAK_CHECK_DF);
    
    Dummy d;
    std::cout << "Hello World!\n";
    char* c = d.myleak();
    d.mysafety();
    
    // 出力内容と出力先
    // https://learn.microsoft.com/ja-jp/cpp/c-runtime-library/reference/crtsetreportmode?view=msvc-170
    // _CRT_WARN          ... 警告、メッセージ、およびすぐに注意する必要のない情報。
    // _CRTDBG_MODE_DEBUG ... デバッガーの出力ウィンドウにメッセージを書き込みます。
    _CrtSetReportMode(_CRT_WARN, _CRTDBG_MODE_DEBUG);
    
    // メモリリークの情報を出力
    // この関数は、呼び出された時点で解放されていないメモリを全て検知する。
    // アプリケーション終了時でのリークを検出したいならここで書かないで
    // _CrtSetDbgFlagでプログラムが終了したときに自動的に呼び出されるように指示したほうがいい。
    // _CrtDumpMemoryLeaks();
    
    return 0;
}

Detected memory leaks!
Dumping objects ->
C:\myproj\ConsoleApplication5\ConsoleApplication5.cpp(19) : {161} normal block at 0x0000029DA8896EA0, 3 bytes long.
 Data: <   > 00 00 00
Object dump complete.
プログラム '[26332] ConsoleApplication5.exe' はコード 0 (0x0) で終了しました。

<?php
/****************************************/
/* 投稿・固定ページ一覧のカラム幅指定 */
function column_title_width(){

echo '

<style>
  .posts .column-cb { width: 40px; }
  .pages .column-cb { width: 40px;}

  .posts .column-title { width: 300px; }
  .pages .column-title { width: 300px; }
  
  .posts .column-author { width: 50px; }
  .pages .column-author { width: 50px; }
  
  .posts .column-categories { width: 100px; }
  .pages .column-categories { width: 100px; }

  .posts .column-tags { width: 150px; }
  .pages .column-tags { width: 150px; }

</style>

';

}

add_action('admin_head','column_title_width');
/****************************************/

?>

#include <windows.h>
#include <memory>

#include "mycuda.cuh"

#pragma warning(disable:4996)

void pnmP3_Write(const char* const fname, const int width, const int height, const unsigned char* const p);

static const int N = 500;//作成画像 N × N


LRESULT CALLBACK WndProc(HWND hwnd, UINT msg, WPARAM wp, LPARAM lp) {

  static std::unique_ptr<unsigned char[]> p_cpu;

  switch (msg) {
  case WM_DESTROY:
    PostQuitMessage(0);
    return 0;
  case WM_CREATE:
    p_cpu.reset(new unsigned char[N * N * 3]);
    call_mycuda(N, p_cpu.get());

    return 0;
  case WM_PAINT:
  {
    // dc設定
    PAINTSTRUCT ps;
    HDC hdc = BeginPaint(hwnd, &ps);
    for(int x=0;x<N;x++)
      for (int y = 0; y < N; y++) {
        size_t pos = N * y + x;//二次元座標を一次元座標に
        SetPixel(hdc, x, y,
          RGB(
            p_cpu[pos * 3 + 0],
            p_cpu[pos * 3 + 1],
            p_cpu[pos * 3 + 2]
          )
        );
      }
    EndPaint(hwnd, &ps);
  }
  }
  return DefWindowProc(hwnd, msg, wp, lp);
}

int WINAPI WinMain(HINSTANCE hInstance, HINSTANCE hPrevInstance,
  PSTR lpCmdLine, int nCmdShow) {
  HWND hwnd;
  MSG msg;
  WNDCLASS winc;

  winc.style = CS_HREDRAW | CS_VREDRAW;
  winc.lpfnWndProc = WndProc;
  winc.cbClsExtra = winc.cbWndExtra = 0;
  winc.hInstance = hInstance;
  winc.hIcon = LoadIcon(NULL, IDI_APPLICATION);
  winc.hCursor = LoadCursor(NULL, IDC_ARROW);
  winc.hbrBackground = (HBRUSH)GetStockObject(WHITE_BRUSH);
  winc.lpszMenuName = NULL;
  winc.lpszClassName = TEXT("SZL-WINDOW");

  if (!RegisterClass(&winc)) return -1;

  hwnd = CreateWindow(
    TEXT("SZL-WINDOW"), TEXT("mywindow"),
    WS_OVERLAPPEDWINDOW | WS_VISIBLE,
    CW_USEDEFAULT, CW_USEDEFAULT,
    N+100, N+100,
    NULL, NULL, hInstance, NULL
  );

  if (hwnd == NULL) return -1;

  while (GetMessage(&msg, NULL, 0, 0)) DispatchMessage(&msg);
  return msg.wParam;
}

/////////////////////////////////////////////
//画像ファイル書き出し/////////////////////////
//! @brief PPM(RGB各1byte,カラー,テキスト)を書き込む
//! @param [in] fname ファイル名
//! @param [in] width 画像の幅
//! @param [in] height 画像の高さ
//! @param [in] p 画像のメモリへのアドレス
//! @details RGBRGBRGB....のメモリを渡すと、RGBテキストでファイル名fnameで書き込む
void pnmP3_Write(const char* const fname, const int width, const int height, const unsigned char* const p) { // PPM ASCII

  FILE* fp = fopen(fname, "wb");
  fprintf(fp, "P3\n%d %d\n%d\n", width, height, 255);

  size_t k = 0;
  for (size_t i = 0; i < (size_t)height; i++) {
    for (size_t j = 0; j < (size_t)width; j++) {
      fprintf(fp, "%d %d %d ", p[k * 3 + 0], p[k * 3 + 1], p[k * 3 + 2]);
      k++;
    }
    fprintf(fp, "\n");
  }

  fclose(fp);
}

#include <stdio.h>
#include <cstring>

#include <cuda_runtime.h>

void call_mycuda(const int N, unsigned char* p_cpu);


///////////////////////////////////////////////
// GPU側 ////////////////////////////////////// 
__global__ void gradation(unsigned char* c, const int N);

#include "mycuda.cuh"

#include <device_launch_parameters.h>

#pragma comment(lib,"cudart.lib")
#pragma comment(lib,"cuda.lib")


void call_mycuda(const int N, unsigned char* p_cpu) {
  //作成画像 N × N × 3

  //16×16の領域に分けて計算する
  dim3 block(16, 16);

  //グリッド数
  dim3 grid(N / 16 + 1, N / 16 + 1);

  unsigned char* p_gpu;//GPU側メモリ確保
  cudaMalloc((void**)&p_gpu, N * N * 3);
  gradation << <grid, block >> > (p_gpu, N);

  //GPU→CPU側へメモリコピー
  cudaMemcpy(p_cpu, p_gpu, N * N * 3, cudaMemcpyDeviceToHost);
  cudaFree(p_gpu);

}

__global__ void gradation(unsigned char* c, const int N) {

    //アクセス法
    //このスレッドが担当する画素の位置を二次元座標で求める
    size_t xpos = blockIdx.x * blockDim.x + threadIdx.x;
    size_t ypos = blockIdx.y * blockDim.y + threadIdx.y;

    if (xpos < N && ypos < N) {
        size_t pos = N * ypos + xpos;//二次元座標を一次元座標に
        float R = ypos / (float)N;
        float G = (N - ypos) / (float)N;
        float B = xpos / (float)N;
        c[pos * 3 + 0] = R * 255;//0～255にして色書き込み
        c[pos * 3 + 1] = G * 255;
        c[pos * 3 + 2] = B * 255;
    }
}

// 進捗管理用
#include <vtkProgressObserver.h>

// PLYファイル読み込み用
#include <vtkPLYReader.h>

// 読み込めているか確認用
#include <vtkPointData.h>

//////////////////////////////////////////////////
// 進捗管理のためにvtkProgressObserverを継承したクラスを作成
class MyProgressObserver : public vtkProgressObserver {
public:

  // UpdateProgressをオーバーライド
  virtual void UpdateProgress(double amount)override {
    // 進捗を表示
    printf("== %d\n", int(amount*100) );// amountは進捗状況を0.0～1.0で表現
  }
};

void myload() {

  const char* filename = "C:\\test\\mill.ply";

  ////////////////////////////////////////////////////////////

  // 進捗管理用のオブザーバを作成
  vtkSmartPointer<MyProgressObserver> observer = vtkSmartPointer<MyProgressObserver>(new MyProgressObserver);

  ////////////////////////////////////////////////////////////
  // PLYファイルを読み込むオブジェクト作成
  vtkSmartPointer<vtkPLYReader> plyReader = vtkSmartPointer<vtkPLYReader>::New();
  plyReader->SetProgressObserver(observer);// オブザーバを設定
  plyReader->SetFileName(filename);
  plyReader->Update();
  ////////////////////////////////////////////////////////////

  // 読み込んだデータを取得
  vtkSmartPointer<vtkPolyData> polyData = plyReader->GetOutput();
  vtkSmartPointer<vtkPoints> points = polyData->GetPoints();

  if (points == nullptr) {
    return;
  }
  for (vtkIdType i = 0; i < points->GetNumberOfPoints(); i++)
  {
    double p[3];
    points->GetPoint(i, p);
    printf("%lf %lf %lf\n", p[0], p[1], p[2]);
  }

}