Created
February 28, 2025 12:40
-
-
Save UmarZein/6dd8956e21460653ba90f49a598dead7 to your computer and use it in GitHub Desktop.
proj
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| #include <math.h> | |
| #include <stdio.h> | |
| #include <stdlib.h> | |
| #include <time.h> | |
| int bernoulli_sample(double p) { | |
| return ((double) rand() / RAND_MAX) < p ? 1 : 0; | |
| } | |
| const double SECOND_SCALE=1000000000; | |
| int main(int argc, char *argv[]) { | |
| srand(time(NULL)); | |
| if (argc != 3) { // Check if a filename was provided | |
| fprintf(stderr, "Usage: %s <filename> <corruption rate>\n", argv[0]); | |
| return 1; | |
| } | |
| double p = atof(argv[2]); | |
| FILE *file; | |
| fopen_s(&file, argv[1], "r"); // Open the file specified in the command-line argument | |
| if (!file) { | |
| perror("Error opening file"); | |
| return 1; | |
| } | |
| long long timestamp; | |
| double value; | |
| long long xa=0.0; | |
| long long xb=0.0; | |
| double ya=0.0; | |
| double yb=0.0; | |
| double total_error = 0.0; | |
| unsigned int n = 0; | |
| // mainloop is here. | |
| while (fscanf_s(file, "%lld %lf", ×tamp, &value) == 2) { | |
| if (bernoulli_sample(p)){ | |
| // corrupted | |
| double dx = xb-xa; | |
| double dy = yb-ya; | |
| double grad=dy/dx; | |
| double pred = ya+grad*((double)(timestamp-xb)); | |
| double se = powf(value-pred, 2.0); | |
| total_error += se; | |
| n++; | |
| xa=xb; | |
| xb=timestamp; | |
| ya=yb; | |
| yb=pred; | |
| } else { | |
| // regular | |
| xa=xb; | |
| ya=yb; | |
| xb=timestamp; | |
| yb=value; | |
| } | |
| } | |
| fclose(file); | |
| double final_error = sqrtf(total_error/((double)n)); | |
| printf("RMSE: %lf", final_error); | |
| return 0; | |
| } | |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| #include <math.h> | |
| #include <stdio.h> | |
| #include <stdlib.h> | |
| #include <time.h> | |
| int bernoulli_sample(double p) { | |
| return ((double) rand() / RAND_MAX) < p ? 1 : 0; | |
| } | |
| const double SECOND_SCALE=1000000000; | |
| #define SIZE 1000 | |
| int main(int argc, char *argv[]) { | |
| srand(time(NULL)); | |
| if (argc != 3) { // Check if a filename was provided | |
| fprintf(stderr, "Usage: %s <filename> <corruption rate>\n", argv[0]); | |
| return 1; | |
| } | |
| double p = atof(argv[2]); | |
| FILE *file; | |
| fopen_s(&file, argv[1], "r"); // Open the file specified in the command-line argument | |
| if (!file) { | |
| perror("Error opening file"); | |
| return 1; | |
| } | |
| long long timestamp; | |
| double value; | |
| int isPredicting=0; | |
| long long xL=0.0; | |
| long long xR=0.0; | |
| double yL=0.0; | |
| double yR=0.0; | |
| long long *unknowns = malloc(SIZE * sizeof(long long)); | |
| unsigned int nUnknowns=0; | |
| double *targets = malloc(SIZE * sizeof(double)); | |
| double totalSquaredError = 0.0; | |
| unsigned int n = 0; | |
| // mainloop is here. | |
| while (fscanf_s(file, "%lld %lf", ×tamp, &value) == 2) { | |
| if (bernoulli_sample(p)){ | |
| // corrupted | |
| isPredicting=1; | |
| unknowns[nUnknowns]=timestamp; | |
| targets[nUnknowns]=value; | |
| nUnknowns++; | |
| n++; | |
| } else { | |
| // regular | |
| if (isPredicting){ | |
| xR=timestamp; | |
| yR=value; | |
| long long dx=xR-xL; | |
| double dy=yR-yL; | |
| double grad=dy/dx; | |
| for (int i=0;i<nUnknowns;i++){ | |
| long long x = unknowns[i]; | |
| double y = targets[i]; | |
| double pred = grad*(x-xL)+yL; | |
| totalSquaredError += powf(pred-y, 2.0); | |
| } | |
| nUnknowns=0; | |
| } | |
| xL=timestamp; | |
| yL=value; | |
| isPredicting=0; | |
| } | |
| } | |
| fclose(file); | |
| double final_error = sqrtf(totalSquaredError/((double)n)); | |
| printf("RMSE: %lf", final_error); | |
| return 0; | |
| } | |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| #include <math.h> | |
| #include <stdio.h> | |
| #include <stdlib.h> | |
| #include <time.h> | |
| int bernoulli_sample(double p) { | |
| return ((double) rand() / RAND_MAX) < p ? 1 : 0; | |
| } | |
| const double SECOND_SCALE=1000000000; | |
| int main(int argc, char *argv[]) { | |
| srand(time(NULL)); | |
| if (argc != 4) { // Check if a filename was provided | |
| fprintf(stderr, "Usage: %s <filename> <corruption rate> <momentum>\n", argv[0]); | |
| return 1; | |
| } | |
| double p = atof(argv[2]); | |
| double momentum = atof(argv[3]); | |
| FILE *file; | |
| fopen_s(&file, argv[1], "r"); // Open the file specified in the command-line argument | |
| if (!file) { | |
| perror("Error opening file"); | |
| return 1; | |
| } | |
| long long timestamp; | |
| double value; | |
| long long xa=0.0; | |
| long long xb=0.0; | |
| double ya=0.0; | |
| double yb=0.0; | |
| double total_error = 0.0; | |
| unsigned int n = 0; | |
| // mainloop is here. | |
| while (fscanf_s(file, "%lld %lf", ×tamp, &value) == 2) { | |
| if (bernoulli_sample(p)){ | |
| // corrupted | |
| double dx = xb-xa; | |
| double dy = yb-ya; | |
| double grad=dy/dx; | |
| double pred = ya+grad*((double)(timestamp-xb)); | |
| double se = powf(value-pred, 2.0); | |
| total_error += se; | |
| n++; | |
| xa=xb; | |
| xb=timestamp; | |
| ya=yb; | |
| yb=pred; | |
| } else { | |
| // regular | |
| xa=xb; | |
| ya=yb; | |
| xb=timestamp; | |
| yb=value; | |
| } | |
| } | |
| fclose(file); | |
| double final_error = sqrtf(total_error/((double)n)); | |
| printf("RMSE: %lf", final_error); | |
| return 0; | |
| } | |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment