三维坐标下的平面点阵（射线点源）

#include <math.h>
#include <stdlib.h>
#include <iostream>
#include <fstream>
using namespace std;

#include "Vector.h"

#define  PI 3.14159265
#define  CC 3E8

void main()
{
    Vector3d dir(-1,-1,-1);//光线照射方向
    
    dir = dir.normalize();//归一化入射方向

    Vector3d cen(0,0,0);//模型的包围球球心
    double R = 2;//半径
    
    Vector3d ori = cen - 2*R*dir;//入射光线原点在物体之外

    Vector3d x(dir.y, -dir.x, 0);

    if (x==Vector3d(0,0,0))
    {
        x = Vector3d(1,0,0);
    }

    x = x.normalize();

    Vector3d y = Cross(dir,x);
    y = y.normalize();
    
    double freq = 30e9;
    double lumda = CC / freq;

    double dx = 0.2;//lumda / 10 ;
    double dy = dx;

    int NX= (int) (2*R/dx);
    int NY= (int) (2*R/dy);
    
    Vector3d** points = new Vector3d*[NX];
    
    for(int i=0;i<NX;i++)
    {
        points[i] = new Vector3d[NY];
    }
        
    fstream file("outdata2.txt",ios::out);

    file<<"x"<<"\t"<<"y"<<"\t"<<"z"<<endl;

    for(int nx = 0;nx < NX;nx++)
    {
        for(int ny=0;ny<NY;ny++)
        {
            points[nx][ny] = ori+dx*(nx-NX/2)*x+dy*(ny-NY/2)*y;
        
            file<<points[nx][ny].x<<"\t"<<points[nx][ny].y<<"\t"<<points[nx][ny].z<<endl;
        }
    }
    
    for(int i=0;i<NX;i++)
    {
        delete []points[i];
    }
    
    file.close();

    cout<<"done!"<<endl;

    getchar();
}