sleef-sys 0.1.2

Rust FFI bindings to the SLEEF Vectorized Math Library 
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
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106

//          Copyright Naoki Shibata 2010 - 2017.
// Distributed under the Boost Software License, Version 1.0.
//    (See accompanying file LICENSE.txt or copy at
//          http://www.boost.org/LICENSE_1_0.txt)

#include <stdio.h>
#include <stdlib.h>
#include <stdint.h>
#include <string.h>

#define CONFIGMAX 4

char *replaceAll(const char *in, const char *pat, const char *replace) {
  const int replaceLen = strlen(replace);
  const int patLen = strlen(pat);

  char *str = malloc(strlen(in)+1);
  strcpy(str, in);

  for(;;) {
    char *p = strstr(str, pat);
    if (p == NULL) return str;

    int replace_pos = p - str;
    int tail_len = strlen(p + patLen);

    char *newstr = malloc(strlen(str) + (replaceLen - patLen) + 1);

    memcpy(newstr, str, replace_pos);
    memcpy(newstr + replace_pos, replace, replaceLen);
    memcpy(newstr + replace_pos + replaceLen, str + replace_pos + patLen, tail_len+1);

    free(str);
    str = newstr;
  }

  return str;
}

#define LEN 1024
char line[LEN+10];

int main(int argc, char **argv) {
  if (argc < 2) {
    fprintf(stderr, "Usage : %s <Base type> <ISA> ...\n", argv[0]);
    exit(-1);
  }

  const char *baseType = argv[1];
  const int isastart = 2;
  const int isamax = argc - isastart;
  const int maxbutwidth = 6;

  for(int config=0;config<CONFIGMAX;config++) {
#if ENABLE_STREAM == 0
    if ((config & 1) != 0) continue;
#endif
    for(int isa=isastart;isa<argc;isa++) {
      char *isaString = argv[isa];
      char configString[100];
      sprintf(configString, "%d", config);
      
      FILE *fpin = fopen("unroll0.org", "r");

      sprintf(line, "unroll_%d_%s.c", config, isaString);
      FILE *fpout = fopen(line, "w");
      fputs("#include \"vectortype.h\"\n\n", fpout);
      fprintf(fpout, "extern %s ctbl_%s[];\n", baseType, baseType);
      fprintf(fpout, "#define ctbl ctbl_%s\n\n", baseType);

      for(;;) {
	if (fgets(line, LEN, fpin) == NULL) break;
	char *s;
	if ((config & 1) == 0) {
	  char *s0 = replaceAll(line, "%ISA%", isaString);
	  s = replaceAll(s0, "%CONFIG%", configString);
	  free(s0);
	} else {
	  char *s0 = replaceAll(line, "%ISA%", isaString);
	  char *s1 = replaceAll(s0, "%CONFIG%", configString);
	  char *s2 = replaceAll(s1, "store(", "stream(");
	  s = replaceAll(s2, "scatter(", "scstream(");
	  free(s0); free(s1); free(s2);
	}

	if ((config & 2) == 0) {
	  char *s0 = replaceAll(s, "#pragma", "//");
	  free(s);
	  s = s0;
	}

	if (config == 0) {
	  char *s0 = replaceAll(s, "#undef EMITREALSUB", "#define EMITREALSUB");
	  free(s);
	  s = s0;
	}
	
	fputs(s, fpout);
	free(s);
      }
    
      fclose(fpin);
      fclose(fpout);
    }
  }
}