diff options
| author | JP Appel <jeanpierre.appel01@gmail.com> | 2024-04-20 03:01:37 -0400 |
|---|---|---|
| committer | JP Appel <jeanpierre.appel01@gmail.com> | 2024-04-20 03:02:29 -0400 |
| commit | 627f4641269d273f3bbf76e679381924115aaf16 (patch) | |
| tree | 2220e884d6b51b84f0102554822f8d2a92230718 /src/util.c | |
| parent | 5cba5bc63de403c779fd4fa762216280d4396444 (diff) | |
moved lattice to complex
Diffstat (limited to 'src/util.c')
| -rw-r--r-- | src/util.c | 22 |
1 files changed, 22 insertions, 0 deletions
diff --git a/src/util.c b/src/util.c new file mode 100644 index 0000000..b76da0b --- /dev/null +++ b/src/util.c @@ -0,0 +1,22 @@ +#include "util.h" + +/* + * Converts a grid point into an complex number + */ +double complex lattice_to_complex(const size_t index, const vec2 resolution) { + const double x_min = -2.0; + const double x_max = 2.0; + const double y_min = -2.0; + const double y_max = 2.0; + + const double x_step = (x_max - x_min) / (double)resolution.x; + const double y_step = (y_max - y_min) / (double)resolution.y; + + const size_t x_index = index % resolution.x; + const size_t y_index = index / resolution.x; + + const double x = x_min + x_index * x_step; + const double y = y_min + y_index * y_step; + + return x + y * I; +} |
