about summary refs log tree commit diff
path: root/src/c/lookup.c
blob: 01b5508fb3d6b27cef3b11484803aa887d7147bd (plain)
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
#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);
  }
}