add reference C program

Author: scivision

CMakeLists.txt

@@ -1,12 +1,15 @@
 cmake_minimum_required(VERSION 3.14)
 
-project(DonutFortran LANGUAGES Fortran)
+project(DonutFortran
+LANGUAGES C Fortran)
 
 if(CMAKE_Fortran_COMPILER_ID STREQUAL GNU)
-  add_compile_options(-Wall -Werror=array-bounds -fcheck=bounds
-  "$<$<CONFIG:Debug,RelWithDebInfo>:-Wextra;-Wpedantic>"
-  "$<$<CONFIG:Release>:-fno-backtrace>"
+  add_compile_options(-Wall -Wextra -Wpedantic
+  "$<$<COMPILE_LANGUAGE:Fortran>:-Werror=array-bounds;-fcheck=bounds>"
+  "$<$<AND:$<CONFIG:Release>,$<COMPILE_LANGUAGE:Fortran>>:-fno-backtrace>"
   )
 endif()
 
 add_executable(donut app/main.f90)
+
+add_executable(donut_ref app/main.c)

app/main.c

@@ -0,0 +1,45 @@
+// https://gist.github.com/gcr/1075131
+
+#include <stdio.h>
+#include <string.h>
+#include <math.h>
+
+int k;
+
+int main(){
+    float A=0, B=0, i, j, z[1760];
+    char b[1760];
+    printf("\x1b[2J");
+    for(; ; ) {
+        memset(b,32,1760);
+        memset(z,0,7040);
+        for(j=0; 6.28>j; j+=0.07) {
+            for(i=0; 6.28 >i; i+=0.02) {
+                float sini=sin(i),
+                      cosj=cos(j),
+                      sinA=sin(A),
+                      sinj=sin(j),
+                      cosA=cos(A),
+                      cosj2=cosj+2,
+                      mess=1/(sini*cosj2*sinA+sinj*cosA+5),
+                      cosi=cos(i),
+                      cosB=cos(B),
+                      sinB=sin(B),
+                      t=sini*cosj2*cosA-sinj* sinA;
+                int x=40+30*mess*(cosi*cosj2*cosB-t*sinB),
+                    y= 12+15*mess*(cosi*cosj2*sinB +t*cosB),
+                    o=x+80*y,
+                    N=8*((sinj*sinA-sini*cosj*cosA)*cosB-sini*cosj*sinA-sinj*cosA-cosi *cosj*sinB);
+                if(22>y&&y>0&&x>0&&80>x&&mess>z[o]){
+                    z[o]=mess;
+                    b[o]=".,-~:;=!*#$@"[N>0?N:0];
+                }
+            }
+        }
+        printf("\x1b[d");
+        for(k=0; 1761>k; k++)
+            putchar(k%80?b[k]:10);
+        A+=0.04;
+        B+= 0.02;
+    }
+}