-
Notifications
You must be signed in to change notification settings - Fork 2
Expand file tree
/
Copy pathhost_window_t.cpp
More file actions
78 lines (71 loc) · 2.09 KB
/
host_window_t.cpp
File metadata and controls
78 lines (71 loc) · 2.09 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
#include <iostream>
#include <vector>
#include "host_window.h"
/**
* @brief Functor returning a gaussian
* \f[
f(x,y) = Ae^{-(\frac{(x-x_0)^2}{2\sigma_x^2} + \frac{(y-y_0)^2}{2\sigma_y^2}}
\f]
*/
struct Gaussian
{
/**
* @brief Functor returning a gaussian
*
* @param x0 x-center-coordinate
* @param y0 y-center-coordinate
* @param sigma_x x - variance
* @param sigma_y y - variance
* @param amp Amplitude
*/
Gaussian( double x0, double y0, double sigma_x, double sigma_y, double amp)
: x00(x0), y00(y0), sigma_x(sigma_x), sigma_y(sigma_y), amplitude(amp){}
/**
* @brief Return the value of the gaussian
*
* \f[
f(x,y) = Ae^{-(\frac{(x-x_0)^2}{2\sigma_x^2} + \frac{(y-y_0)^2}{2\sigma_y^2}}
\f]
* @param x x - coordinate
* @param y y - coordinate
*
* @return gaussian
*/
double operator()(double x, double y)
{
return amplitude*
exp( -((x-x00)*(x-x00)/2./sigma_x/sigma_x +
(y-y00)*(y-y00)/2./sigma_y/sigma_y) );
}
private:
double x00, y00, sigma_x, sigma_y, amplitude;
};
const unsigned Nx = 7000, Ny = 4000;
const double lx = 2., ly = 2.;
const double hx = lx/(double)Nx, hy = ly/(double)Ny;
int main()
{
//Create Window and set window title
GLFWwindow* w = draw::glfwInitAndCreateWindow( 400, 800, "Hello world!");
draw::RenderHostData render( 2, 1);
// generate a vector on the grid to visualize
Gaussian g( 1.2, 0.3, .1, .1, 1);
std::vector<float> visual(Nx*Ny);
for(unsigned i=0; i<Ny; i++)
for( unsigned j=0; j<Nx; j++)
visual[i*Nx+j] = -g( (double)j*hx, (double)i*hy);
//create a colormap
draw::ColorMapRedBlueExt colors( 1.);
draw::ColorMapRedBlueExtMinMax colorsminmax( -1.,2.);
//set scale
colors.scale() = 1.;
while (!glfwWindowShouldClose( w ))
{
render.renderQuad( visual, Nx, Ny, colors);
render.renderQuad( visual, Nx, Ny, colorsminmax);
glfwSwapBuffers(w);
glfwWaitEvents();
}
glfwTerminate();
return 0;
}