aboutsummaryrefslogtreecommitdiff
path: root/init.c
diff options
context:
space:
mode:
authorAnton Khirnov <anton@khirnov.net>2023-01-11 16:19:41 +0100
committerAnton Khirnov <anton@khirnov.net>2023-01-11 16:19:41 +0100
commit5b39260bc65720620cf21438fe49d31ea01fa7da (patch)
tree6bdb313c6ff879cc1d2438cd79bd8de493c2029d /init.c
parentbc9a2906fc30e2d73d6b4c23783df5e15c90278c (diff)
init: drop unused cruft
Diffstat (limited to 'init.c')
-rw-r--r--init.c37
1 files changed, 0 insertions, 37 deletions
diff --git a/init.c b/init.c
index 82f738b..b6b27b5 100644
--- a/init.c
+++ b/init.c
@@ -37,15 +37,6 @@
#include "td_constraints.h"
#include "teukolsky_data.h"
-double tdi_scalarproduct_metric_fma3(size_t len1, size_t len2, double *mat,
- double *vec1, double *vec2);
-double tdi_scalarproduct_metric_avx(size_t len1, size_t len2, double *mat,
- double *vec1, double *vec2);
-double tdi_scalarproduct_metric_sse3(size_t len1, size_t len2, double *mat,
- double *vec1, double *vec2);
-double tdi_scalarproduct_metric_c(size_t len1, size_t len2, double *mat,
- double *vec1, double *vec2);
-
static const enum BasisFamily basis_sets[NB_EQUATIONS][2] = {
{ BASIS_FAMILY_SB_EVEN, BASIS_FAMILY_COS_EVEN },
{ BASIS_FAMILY_SB_EVEN, BASIS_FAMILY_COS_EVEN },
@@ -63,8 +54,6 @@ static int teukolsky_init_check_options(TDContext *td)
TDPriv *s = td->priv;
int ret;
- s->cpu_flags = tdi_init_cpu_flags();
-
if (!td->nb_threads) {
td->nb_threads = tdi_cpu_count();
if (!td->nb_threads)
@@ -78,32 +67,6 @@ static int teukolsky_init_check_options(TDContext *td)
s->logger.log = log_callback;
s->logger.opaque = td;
- //s->scalarproduct_metric = tdi_scalarproduct_metric_c;
- //if (EXTERNAL_SSE3(s->cpu_flags))
- // s->scalarproduct_metric = tdi_scalarproduct_metric_sse3;
- //if (EXTERNAL_AVX(s->cpu_flags))
- // s->scalarproduct_metric = tdi_scalarproduct_metric_avx;
- //if (EXTERNAL_FMA3(s->cpu_flags))
- // s->scalarproduct_metric = tdi_scalarproduct_metric_fma3;
-// for (int i = 0; i < ctx->nb_equations; i++)
-// for (int j = 0; j < 2; j++) {
-// s->ps_ctx->basis[i][j] = ctx->basis[i][j];
-// s->ps_ctx->solve_order[i][j] = basis_order[i][j];
-// max_order = MAX(max_order, basis_order[i][j]);
-// }
-//
- //for (int i = 0; i < max_order; i++) {
- // tdi_log(&s->logger, 2, "%d ", i);
- // for (int j = 0; j < ctx->nb_equations; j++)
- // for (int k = 0; k < 2; k++) {
- // if (i < s->ps_ctx->solve_order[j][k])
- // tdi_log(&s->logger, 2, "%8.8g\t", s->ps_ctx->colloc_grid[j][k][i]);
- // else
- // tdi_log(&s->logger, 2, " ");
- // }
- // tdi_log(&s->logger, 2, "\n");
- //}
-
s->basis_order[0][0] = td->nb_coeffs[0];
s->basis_order[0][1] = td->nb_coeffs[1];
s->basis_order[1][0] = td->nb_coeffs[0];