CudaWaveField

Size: px

Start display at page:

Download "CudaWaveField"

すずりくぬぎ
5 years ago
Views:

1 CudaWaveField

2 2 CudaWaveField Rel Rel 1.0 CudaWaveField ( cwfl) / cwfl cwfl Note Acrobat Reader

3 3 I CudaWaveField CudaWaveField Windows x NVIDIA GPU CUDA WaveFieldLib MS Visual Studio.NET (Visual Studio 2003) MS Visual Studio MS Visual Studio Intel C++ Compiler CudaWaveFieldLib FFT

4 ExactAsmProp void StartCWFL(iconst char* fname = NULL,int msglevel = 1,bool device- Property = true,int devicenumber = 0) CUDA int GetCUDADeviceCount(void) cudadeviceprop GetCUDADeviceInformation(int deviceid) void CudaMallocHost(wfl::Complex* plpointer,long nx,long ny) CudaWaveField CudaWaveField(void) CudaWaveField(long nxy) CudaWaveField(long nx, long ny) CudaWaveField(long nx, long ny, double pxy) CudaWaveField(long nx, long ny, double px, double py) CudaWaveField(long nx,long ny,double px,double py, double wavelength) CudaWaveField(const wfl::wavefield& wf) CudaWaveField(const wfl::fieldparam& param) CudaWaveField(const CudaWaveField& temp) CudaWaveField(void) void Init(float2* d = NULL) Dispose(void) void Clear(void) void ThreadsOptimization(void) void ThreadsOptimization(long nx,long ny) void CopyParam(const wfl::fieldparam& param) void CopyParam(const CudaWaveField& param) void CopyParam(const wfl::wavefield& param) void CopyParamAll(const wfl::fieldparam& param) void CopyParamAll(const CudaWaveField& param) void CopyParamAll(const wfl::wavefield& param)

5 5 void CopyWinAreaData(CudaWaveField& dest) CudaWaveField& Embed(int nxy) CudaWaveField& Embed(int nx,int ny) CudaWaveField& Extract(int nxy) CudaWaveField& Extract(int nx,int ny) void AdjustNumPixel(void) CudaWaveField& SetRect(double wx,double wy,float amp = 1.0) CudaWaveField& SetGaussian(double r,double nn = 2.0,double amp = 1.0) 28 CudaWaveField& SetPlaneWave(wfl::Vector v,wfl::phase phs = 0.0) CudaWaveField& SetSphericalWave(wfl::Point sphwcenter, double amp = 1.0) CudaWaveField& MultiplyPlaneWave(Vector v,wfl::phase phs = 0.0) CudaWaveField& MultiplyPlaneWave(double CosA,double CosB, wfl::phase phs = 0.0) CudaWaveField& AddSphericalWave(wfl::Point p) CudaWaveField& AddSphericalWave(wfl::Point p, wfl::phase phs,double a) CudaWaveField& SetRandomPhase(void) CudaWaveField& SetQuadraticPhase(double f) float2 GetDataPointer(void) const size t GetDataSize(void) const void SetParam(wfl::FieldParam& param) void SetBlocks(int x,int y,int z) dim3 GetBlocks(void) void SetThreads(int x,int y,int z) dim3 GetThreads(void) FFT void FakeFft(void) void Fft(int s, bool beforeshift = true,bool aftershift = true) void Fft(int s,cudafftplan cuplan,bool beforeshift = true,bool aftershift = true) void SwitchFs (void) CudaWaveField& RawScaledFft(double s,double t)

6 6 CudaWaveField& Rotate(const CudaWaveField& source, wfl::rmatrix& crmat, wfl::sfrequency c) CudaWaveField& RotateFs(wfl::Rmatrix& crmat,cudawavefield& cpfb,double2 offset) void AsmProp(double d) void AsmPropFs(double d) void RawAsmProp(double d) void RawAsmPropFs(double d) CudaWaveField& ExactAsmProp(double d) CudaWaveField& ShiftedAsmProp(const CudaWaveField& source,int precision = 1) CudaWaveField& FourierProp(double f) CudaWaveField& BackFourierProp(double f) CudaWaveField& ShiftedFresnelProp(const CudaWaveField& source) CudaWaveField& ShiftedFresnelProp(CudaWaveField& source, CudaShiftedFresnelPropDesc& csfpd) CudaWaveField& ShiftedFresnelPropAdd(CudaWaveField& source, CudaShiftedFresnelPropDesc& csfpd) CudaWaveField& ShiftedFresnelPropEx(CudaWaveField& source) CudaWaveField& ShiftedFresnelPropAddEx(CudaWaveField& source) void PaintTriangle(const wfl::pointarray& p, float2 amp,bool memset) PaintPolygonShape(const wfl::pointarray& p,float2 amp) double MaxReduction(void) double2 MaxMinRedution(void) void Normalize(double amp = 1.0) void AddFrom(CudaWaveField& target)

7 7 CudaWaveField& operator=(const CudaWaveField& cwf) CudaWaveField& operator+(const CudaWaveField& cwf) CudaWaveField& operator+=(const CudaWaveField& cwf) CudaWaveField& operator-(const CudaWaveField& cwf) CudaWaveField& operator-=(const CudaWaveField& cwf) CudaWaveField& operator (const CudaWaveField& cwf) CudaWaveField& operator =(const CudaWaveField& cwf) CudaWaveField& operator/(const CudaWaveField& cwf) CudaWaveField& operator/=(const CudaWaveField& cwf) CudaWaveField& operator = (double val) CudaWaveField& opeartor = (const wfl::complexdouble& val) CudaWaveField& operator/= (double val) CudaWaveField& operator/= (const wfl::complexdouble& val) CudaWaveField& operator+= (const wfl::complexdouble& val) CudaWaveField& operator-=(const wfl::complexdouble& val) void operator>> (wfl::wavefield& target) CudaWaveField& operator<< (const wfl::wavefield& wf)

9 I CudaWaveField

11 CudaWaveField NVIDIA GPU WaveField Windows x64 64bit 32bit Windows NVIDIA GPU Compute Capability 2.0 GPU 1.3 GPU CUDA NVIDIA GPU CUDA CUDA CUDA CUDA CUDA Toolkit GPU Computing SDK WaveFieldLib3 WaveFieldLib3(WFL3) WFL WFL3 CUDA WFL CUDA

12 WaveField MS Visual Studio.NET (Visual Studio 2003) MS Visual Studio 2005 MS Visual Studio 2008 IntelliSense 64 CPU (64 Visual Studio ) Intel C++ Compiler WaveField < >=c:\wavefieldtools < >\include\ < >\include\cwfl\ *.h WaveField < >\lib\ cwfl.lib WaveField PATH (.exe ) < >\bin\

13 1.1 CudaWaveField 13 cwfl.dll

15 15 2 CudaWaveFieldLib 2.1 CudaWaveField WaveField GPU WaveField WaveField CudaWaveField WaveField WaveField 2.2 CudaWaveField Example # include <cwfl.h> using namespace wfl ; using namespace cwfl ; int main ( void ) { StartCWFL (); /* */ } return 0; cwfl.h CudaWaveField cwfl.h cwfl.h wfl3.h wfl3.h

16 16 2 CudaWaveFieldLib CudaWaveField cwfl cwfl CudaWaveField StartCWFL() Note StartCWFL() wfl::start() StartCWFL() Start() cwfl.lib ( wfl3.lib ) PATH (*.EXE) cwfl.dll ( wfl3.dll ) 2.3 FFT 1[cm]. Example FFT (: ExRectApertureFFT.cpp) # include <cwfl.h> using namespace std ; using namespace wfl ; using namespace cwfl ; int main ( void ) { // C W F L StartCWFL ("./ FFT / FFT_test. log ",1,true,0); // EventTimer timer ; timer. Start (); // 512 x [ mm] G P U CudaWaveField cwf (512, 512, 1e -3); // 1[ cm] x 1[ cm] cwf. SetRect (0.01,0.01); // cwf. Fft ( -1); // cwf. Normalize (); // timer. Stop (); Printf (" GPU Elapsed time = %f[ms ]\n",timer. GetTime ()); // G P U C P U WaveField wf; cwf >> wf; // b m p wf. SaveAsBmp ("./ FFT /[ GPU ] RectFFT_AMPLITUDE. bmp ",AMPLITUDE, COLOR );

2.4 17 } return 0; WaveField WaveField CudaWaveField SetRect() 1cm 1cm Fft() Normalize(). wfl::wavefield.. GPU CPU GPU! 2.1 2.4 2.4.1 CudaWaveField Example (: ExShortProp.

17 } return 0; WaveField WaveField CudaWaveField SetRect() 1cm 1cm Fft() Normalize(). wfl::wavefield.. GPU CPU GPU! CudaWaveField Example (: ExShortProp.cpp) # include <cwfl.h> using namespace std ; using namespace wfl ; using namespace cwfl ; int main ( void ) { // C W F L StartCWFL ("./ AsmProp / AsmProp. log ",1,true,0); //

$/ AsmProp / Diffraction. bmp ",AMPLITUDE ); return 0; SetGaussian() AsmProp() 2.2! "$#%'&)(+*-,./0132 546 7 2.$

18 18 2 CudaWaveFieldLib CudaWaveField cwf (256,256,2e -6); // cwf. SetGaussian (0.1e -3,50); // 1 m m cwf. AsmProp (1e -3); // cwf. Normalize (); // WaveField saver ; cwf >> saver ; } // saver. SaveAsBmp ("./ AsmProp / Diffraction. bmp ",AMPLITUDE ); return 0; SetGaussian() AsmProp() 2.2! "$#%'&)(+*-,./ mm 10mm 30mm Example (: ExLongProp.cpp) # include <cwfl.h> using namespace std ; using namespace wfl ; using namespace cwfl ;

2.5 19 int main ( void ) { StartCWFL ("./ AsmProp / AsmProp2. log ",1,true,0); CudaWaveField aperture (256,256,2e -6); aperture. SetGaussian (0.

$AsmProp (d); char filename [128]; sprintf ( filename," Diffraction -%d. bmp ",( int )(d / 1e -3)); dest. Normalize (); } WaveField saver ; dest >> saver ; saver.$

19 int main ( void ) { StartCWFL ("./ AsmProp / AsmProp2. log ",1,true,0); CudaWaveField aperture (256,256,2e -6); aperture. SetGaussian (0.1e -3,50); for ( int i = 1; i <= 3; i ++) { CudaWaveField dest ( aperture ); // W a v e F i e l d double d = 10e -3 * i; dest. AsmProp (d); char filename [128]; sprintf ( filename," Diffraction -%d. bmp ",( int )(d / 1e -3)); dest. Normalize (); } WaveField saver ; dest >> saver ; saver. SaveAsBmp ( filename, AMPLITUDE ); } return 0; 2.3 AsmProp() ExactAsmProp() 2.3 AsmProp() (a)10mm (b)20mm (3)30mm ExactAsmProp AsmProp() ExactAsmProp() 2.4 ExactAsmProp() CudaWaveField. Example (: ExImageFormationByLens.cpp)

WaveField saver ; // W F // /////// ///////// WaveField temp ; // W F // b m p temp. LoadBmp ("./ EyeSimulation / lena128x128. bmp ",AMPLITUDE ); CudaWaveField source = temp ; // source.

20 20 2 CudaWaveFieldLib 2.4 ExactAsmProp() # include <cwfl.h> using namespace std ; using namespace wfl ; using namespace cwfl ; int main ( void ) { // /////////// /////////// StartCWFL (); // FieldParam :: SetDefault (128,128,2e -6); // double R = 1e -3; // WaveField saver ; // W F // /////// ///////// WaveField temp ; // W F // b m p temp. LoadBmp ("./ EyeSimulation / lena128x128. bmp ",AMPLITUDE ); CudaWaveField source = temp ; // source. Embed (3); //8 6 4 // //////// ////////// CudaWaveField lens ; // C W F lens. Embed (3); // 64 lens. SetGaussian (R,100); // lens. SetQuadraticPhase (10e -3); // 10 m m lens >> saver ; // CWF -> W F saver. SaveAsWf ("./ EyeSimulation / lens.wf"); // W F // ////////// -> //////// source. AsmProp (20e -3); // (20 mm) // ////////// ////////////////// lens *= source ; // source. Dispose (); // s o u r c e // ///////// -> /////// lens. AsmProp (20e -3); // (20 mm) // ///////// //////////////// CudaWaveField & screen = lens ; // screen. Extract (3); // 64 ( ) screen. Normalize (); // // ///////// /////////// screen >> saver ; // G P U C P U saver. SaveAsBmp ("./ EyeSimulation / image. bmp ",AMPLITUDE ); // } return 0;

21 FieldParam::SetDefault() CudaWaveField WaveField

23 WaveField void StartCWFL(iconst char* fname = NULL,int msglevel = 1,bool deviceproperty = true,int devicenumber = 0) CudaWaveField fname msglevel deviceproperty GPU devicenumber (GPU) ID Note CudaWaveFieldLib CUDA int GetCUDADeviceCount(void) GPU PC GPU(CUDA ) Note

24 24 3 PC CUDA CWFL cudadeviceprop GetCUDADeviceInformation(int deviceid) cudadeviceprop deviceid ( ComputeCapability ) cudadeviceprop CUDA4.0 void CudaMallocHost(wfl::Complex* plpointer,long nx,long ny) nx * ny * sizeof(complex) CUDA CPU plpointer CudaWaveField Init() Example // ( C P U ) float2 * host_memory ; // G P U CudaMallocHost ( static_cast < wfl :: Complex * >( host_memory ),1024,1024); // C W F CudaWaveField cwf (1024,1024); // C W F cwf. Init ( host_memory ); // C W F

25 3.2 CudaWaveField CudaWaveField CudaWaveField(void) CudaWaveField(long nxy) CudaWaveField(long nx, long ny) CudaWaveField(long nx, long ny, double pxy) CudaWaveField(long nx, long ny, double px, double py) CudaWaveField(long nx,long ny,double px,double py, double wavelength) CudaWaveField(const wfl::wavefield& wf) CudaWaveField(const wfl::fieldparam& param) 1 2 Note x y nxy 3 x nx y nx x y p [ :m] 4 x y px py 5 wavelength [ :m] 6 wfl::wavefield CPU GPU 7 wfl::fieldparam () N x = N y = 256 P x = P y = 1µm λ = 632.8nm wfl::fieldparam::setdefault*** Example // CudaWaveField cwf1 ; //512 x CudaWaveField cwf2 (512); // 1024 x x y 1 0 u m CudaWaveField cwf3 (1024, 1024, 10e -6);

26 26 3 // 2048 x x 5 u m y 1 0 u m CudaWaveField cwf3 (2048, 2048, 5e -6, 10e -6); // W a v e F i e l d C W F wfl :: WaveField wf1 (1024,1024); CudaWaveField cwf4 ( wf1 ); // wfl :: F i e l d P a r a m C W F wfl :: FieldParam param (4096,4096); CudaWaveField cwf5 ( param ); CudaWaveField(const CudaWaveField& temp) ~CudaWaveField(void) GPU void Init(float2* d = NULL) CudaWaveField GPU d Dispose(void) GPU (Init() NULL ) void Clear(void) 0

27 3.2 CudaWaveField 27 void ThreadsOptimization(void) void ThreadsOptimization(long nx,long ny) CudaWaveField 1 CudaWaveField nx,ny 2 nx,ny void CopyParam(const wfl::fieldparam& param) void CopyParam(const CudaWaveField& param) void CopyParam(const wfl::wavefield& param) nx,ny void CopyParamAll(const wfl::fieldparam& param) void CopyParamAll(const CudaWaveField& param) void CopyParamAll(const wfl::wavefield& param) nx,ny void CopyWinAreaData(CudaWaveField& dest) dest Window dest Window

28 CudaWaveField& Embed(int nxy) CudaWaveField& Embed(int nx,int ny) 1 2 nxy 2 2 nx 2 ny 0 CudaWaveField& Extract(int nxy) CudaWaveField& Extract(int nx,int ny) / nxy / nx 2 1 / ny void AdjustNumPixel(void) CudaWaveField& SetRect(double wx,double wy,float amp = 1.0) wx wy amp CudaWaveField& SetGaussian(double r,double nn = 2.0,double amp = 1.0) r nn amp CudaWaveField& SetPlaneWave(wfl::Vector v,wfl::phase phs = 0.0) v phs

29 3.2 CudaWaveField 29 CudaWaveField& SetSphericalWave(wfl::Point sphwcenter, double amp = 1.0) sphwcenter amp (CudaWaveField::GetOrigin() ) CudaWaveField& MultiplyPlaneWave(Vector v,wfl::phase phs = 0.0) CudaWaveField& MultiplyPlaneWave(double CosA,double CosB, wfl::phase phs = 0.0) 1 v phs 2 x CosA y CosB phs CudaWaveField& AddSphericalWave(wfl::Point p) CudaWaveField& AddSphericalWave(wfl::Point p, wfl::phase phs,double a) p 2 phs a p CudaWaveField& SetRandomPhase(void) CudaWaveField& SetQuadraticPhase(double f) f 2

30 float2 GetDataPointer(void) const float2 CudaWaveField GPU CPU size t GetDataSize(void) const size t size t [byte] void SetParam(wfl::FieldParam& param) wfl::fieldparam void SetBlocks(int x,int y,int z) dim3 GetBlocks(void) dim3 CUDA dim3 void SetThreads(int x,int y,int z) dim3 GetThreads(void) dim3 CUDA dim3

31 3.2 CudaWaveField FFT void FakeFft(void) void Fft(int s, bool beforeshift = true,bool aftershift = true) void Fft(int s,cudafftplan cuplan,bool beforeshift = true,bool aftershift = true) 1 s(-1 1 ) beforeshift,aftershift ( ) 2 CudaFFTPlan FFT CudaFFTPlan void SwitchFs (void) CudaWaveField& RawScaledFft(double s,double t) ScaledFFT CudaWaveField& Rotate(const CudaWaveField& source, wfl::rmatrix& crmat, wfl::sfrequency c) source wfl::rmatrix crmat wfl::sfrequency c NULL c

32 32 3 CudaWaveField& RotateFs(wfl::Rmatrix& crmat,cudawavefield& cpfb,double2 offset) void AsmProp(double d) d[m] void AsmPropFs(double d) d[m] void RawAsmProp(double d) d[m] void RawAsmPropFs(double d) d[m] CudaWaveField& ExactAsmProp(double d) d[m] 4 CudaWaveField& ShiftedAsmProp(const CudaWaveField& source,int precision = 1) source source

33 3.2 CudaWaveField 33 CudaWaveField& FourierProp(double f) f. CudaWaveField& BackFourierProp(double f) f CudaWaveField& ShiftedFresnelProp(const CudaWaveField& source) CudaWaveField& ShiftedFresnelProp(CudaWaveField& source, CudaShiftedFresnelPropDesc& csfpd) 1 source 2 CudaShiftedFresnelPropDesc CudaWaveField& ShiftedFresnelPropAdd(CudaWaveField& source, CudaShiftedFresnelPropDesc& csfpd) source. CudaWaveField& ShiftedFresnelPropEx(CudaWaveField& source) source CudaWaveField& ShiftedFresnelPropAddEx(CudaWaveField& source) source.

34 void PaintTriangle(const wfl::pointarray& p, float2 amp,bool memset) p amp memset PaintPolygonShape(const wfl::pointarray& p,float2 amp) p amp double MaxReduction(void) double double2 MaxMinRedution(void) double2 double2 CUDA double.x double.y void Normalize(double amp = 1.0) amp void AddFrom(CudaWaveField& target) target

35 3.2 CudaWaveField CudaWaveField& operator=(const CudaWaveField& cwf) CudaWaveField& operator+(const CudaWaveField& cwf) CudaWaveField& operator+=(const CudaWaveField& cwf) CudaWaveField& operator-(const CudaWaveField& cwf) CudaWaveField& operator-=(const CudaWaveField& cwf) CudaWaveField& operator (const CudaWaveField& cwf) CudaWaveField& operator =(const CudaWaveField& cwf) CudaWaveField& operator/(const CudaWaveField& cwf) CudaWaveField& operator/=(const CudaWaveField& cwf) CWF CudaWaveField& operator = (double val) CudaWaveField& opeartor = (const wfl::complexdouble& val) CudaWaveField& operator/= (double val) CudaWaveField& operator/= (const wfl::complexdouble& val) CudaWaveField& operator+= (const wfl::complexdouble& val) CudaWaveField& operator-=(const wfl::complexdouble& val) double Real Imaginary val wfl::complexdouble void operator>> (wfl::wavefield& target) CudaWaveField& operator<< (const wfl::wavefield& wf) 1 CudaWaveField WaveField 2 WaveField CudaWaveField CudaWaveField CudaWaveField HDD CudaWaveField. Example CudaWaveField cwf1 ;

36 36 3 WaveField temp ; // c w f 1 cwf1 >> temp ; temp. SaveAsWf (" cwf1.wf"); //(ex. ) temp. LoadBmp (" lena. bmp ",AMPLITUDE ); cwf1 << temp ; Note WaveField CPU CudaWaveField GPU PCI-E PCI-E ( )

Wave Field Library (WFL) Polygon Source Library (PSL)

Wave Field Library (WFL) Polygon Source Library (PSL) 2019 6 5 2 WaveField Rel 3.5.1 PolygonSource Rel 1.9.0 Rel 2.13 WaveField ( WFL) PolygonSource ( PSL) / WFL/PSL WFL/PSL http://www.laser.ee.kansai-u.ac.jp/wavefieldtools