diff options
Diffstat (limited to 'src/lookup.c')
-rw-r--r-- | src/lookup.c | 74 |
1 files changed, 0 insertions, 74 deletions
diff --git a/src/lookup.c b/src/lookup.c deleted file mode 100644 index 01b5508..0000000 --- a/src/lookup.c +++ /dev/null @@ -1,74 +0,0 @@ -#include <stdio.h> -#include <math.h> -#include <stdlib.h> - -// Define some variables -int sigma = 200; -float f_0 = 0.1; -float R_s = 1e4; - -// Define some constants -float pi = M_PI; -float e = M_E; -float G = 4.302e-3; - -// Define the directory where the data shall be stored - -// Prototype -float rho(float r); -float phi(float x); - -int main(int argc, char *argv[] ){ - if( argc == 2 ) { - printf("Generating %s Values...\n", argv[1]); - - FILE * fp = fopen("testFile.csv", "w+"); - if(fp == NULL){ - printf("Error opening the file!\n"); - exit(1); - } - - int x = atoi(argv[1]); - - for(int i = 0; i < x; i++){ - fprintf(fp, "%d, %f\n", i, phi(i)); - } - - // for(int i = 0; i < 1e7; i = i + 1e4){ - // printf("%s, %f\n", i, phi(i)); - // } - - fclose(fp); - - } - else if( argc > 2 ) { - printf("Too many arguments supplied.\n"); - return 0; - } - else { - printf("One argument expected.\n"); - return 0; - } - - return 0; -} - -// Define rho function -float rho(float r){ - float a = (1) / ( sqrt(2 * pi) * sigma); - float b = exp( - (phi(r) / pow(sigma, 2) )); - return a; -} - -// Define phi function -float phi(float x){ - if(x == 0) { - float a = -4 * pi * f_0 * G * pow(R_s, 2); - return(a); - } - else { - float a = - (4 * pi * G * f_0 * pow(R_s, 3) / x); - float b = log(1 + (x / R_s) ); - return(a * b); - } -} |