From 1e55dd14365f9d2c4988bac9a8c7f7267dd40e1d Mon Sep 17 00:00:00 2001
From: "Mulas, Giacomo" <gmulas@oa-cagliari.inaf.it>
Date: Wed, 20 Mar 2024 15:24:22 +0100
Subject: [PATCH] correct first 2 valgrind errors (one to double check by
Giovanni)
---
src/libnptm/clu_subs.cpp | 4 ++--
src/libnptm/sph_subs.cpp | 2 +-
2 files changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/libnptm/clu_subs.cpp b/src/libnptm/clu_subs.cpp
index 72cac991..744eae86 100644
--- a/src/libnptm/clu_subs.cpp
+++ b/src/libnptm/clu_subs.cpp
@@ -834,8 +834,8 @@ void hjv(
const int rfj_size = (lit > lmt) ? lit : lmt;
const int rfn_size = c4->litpo;
double *rfj, *rfn;
- rfj = new double[rfj_size]();
- rfn = new double[rfn_size]();
+ rfj = new double[rfj_size+1]();
+ rfn = new double[rfn_size+1]();
jer = 0;
int ivhb = 0;
for (int nf40 = 1; nf40 <= nsphmo; nf40++) { // GPU portable?
diff --git a/src/libnptm/sph_subs.cpp b/src/libnptm/sph_subs.cpp
index 72f8cb10..f3a7d54d 100644
--- a/src/libnptm/sph_subs.cpp
+++ b/src/libnptm/sph_subs.cpp
@@ -668,7 +668,7 @@ void rbf(int n, double x, int &nm, double sj[]) {
double abs_sb = (sb < 0.0) ? -sb : sb;
if (abs_sa > abs_sb) cs = sa / f;
else cs = sb / f0;
- for (int k = 0; k <= nm; k++) {
+ for (int k = 0; k <= nm ; k++) {
sj[k] = cs * sj[k];
}
}
--
GitLab