diff options
author | Emile <HanEmile@users.noreply.github.com> | 2018-01-01 20:50:15 +0100 |
---|---|---|
committer | GitHub <noreply@github.com> | 2018-01-01 20:50:15 +0100 |
commit | 19bae8d7db56d4b87ca997b19e1ac1daa68b8555 (patch) | |
tree | 4df16102d71a607e6b31871938238a0166b64717 /src/c/lookup.c | |
parent | 7f409708c3ca3e0af213efa4d6f8390b4f46f78e (diff) | |
parent | eed019682e3d54334fbcb3b9a612086cfb9d4295 (diff) |
Merge pull request #3 from HanEmile/merge_read_lookup
merged the two files langfassung.c and read.c
Diffstat (limited to 'src/c/lookup.c')
-rw-r--r-- | src/c/lookup.c | 74 |
1 files changed, 0 insertions, 74 deletions
diff --git a/src/c/lookup.c b/src/c/lookup.c deleted file mode 100644 index 01b5508..0000000 --- a/src/c/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); - } -} |