ClangBug/main.cpp at master · modkin/ClangBug · GitHub

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
#include <iostream>
#include <cmath>
#include <vector>

inline size_t num_microvertices_per_edge(size_t level)
{
  return static_cast<size_t>(std::pow(2, level) + 1);
}


inline size_t num_microvertices_per_face(size_t level)
{
  return static_cast<size_t>((std::pow(2, level)+1) * (std::pow(2, level-1) + 1));
}

inline void prolongate(size_t level)
{
  size_t rowsize_coarse = num_microvertices_per_edge(level-2);
  size_t rowsize_fine = num_microvertices_per_edge(level-2+1);

//  double* face_data_f = (double *) malloc(num_microvertices_per_face(level-2+1)*sizeof(double));
//  for(int i = 0; i < num_microvertices_per_face(level-2+1)*sizeof(double); ++i){
//    face_data_f[i] = 0.1;
//  }
//  double* face_data_c = (double *) malloc(num_microvertices_per_face(level-2)*sizeof(double));
//  for(int i = 0; i < num_microvertices_per_face(level-2)*sizeof(double); ++i){
//    face_data_c[i] = 5;
//  }

  std::vector<double> face_data_f(num_microvertices_per_face(level-2+1),0);
  std::vector<double> face_data_c(num_microvertices_per_face(level-2),1);


  size_t mr_c = 1;
  size_t mr_f = rowsize_fine + 2;

  size_t i_rowsize_coarse = rowsize_coarse;

  for (size_t i_coarse = 0; i_coarse < rowsize_coarse-2; ++i_coarse)
  {
    for (size_t j_coarse = 0; j_coarse < i_rowsize_coarse-3; ++j_coarse)
    {
      face_data_f[mr_f] = 0.5 * (face_data_c[mr_c] + face_data_c[mr_c + i_rowsize_coarse]);
      face_data_f[mr_f-1] = 0.5 * (face_data_c[mr_c] + face_data_c[mr_c + i_rowsize_coarse - 1]);
      //size_t OMG = mr_f + rowsize_fine - 1 - 1;
      face_data_f[mr_f + rowsize_fine - 1 - 1] = 0.5 * (face_data_c[mr_c + i_rowsize_coarse] + face_data_c[mr_c + i_rowsize_coarse - 1]);

      face_data_f[mr_f + rowsize_fine - 1] = face_data_c[mr_c + i_rowsize_coarse];

      mr_c += 1;
      mr_f += 2;
    }

    face_data_f[mr_f] = 0.5 * (face_data_c[mr_c] + face_data_c[mr_c + i_rowsize_coarse]);
    face_data_f[mr_f-1] = 0.5 * (face_data_c[mr_c] + face_data_c[mr_c + i_rowsize_coarse - 1]);
    face_data_f[mr_f + rowsize_fine - 1 - 1] = 0.5 * (face_data_c[mr_c + i_rowsize_coarse] + face_data_c[mr_c + i_rowsize_coarse - 1]);

    mr_c += 3;
    mr_f += rowsize_fine - 1 + 3;

    rowsize_fine -= 2;
    i_rowsize_coarse -= 1;
  }
  double sum = 0;
  for(size_t i = 0; i < face_data_f.size(); ++i){
    sum += face_data_f[i];
  }
  std::cout << sum/face_data_f.size() << std::endl;
}


int main()
{
  prolongate(8);
  return 0;
}