Commit 2a0a5f14 authored by Giovanni La Mura's avatar Giovanni La Mura
Browse files

Fix harmonic indexing in wmamp()

parent b4393d99

src/include/sph_subs.h

+1 −1
Original line number Original line Diff line number Diff line
@@ -1239,7 +1239,7 @@ void wmamp( 
) {

) {

	std::complex<double> *ylm = new std::complex<double>[1682];

	std::complex<double> *ylm = new std::complex<double>[1682];

	int nlmp = lm * (lm + 2) + 2;

	int nlmp = lm * (lm + 2) + 2;

	ylm[nlmp] = std::complex<double>(0.0, 0.0);

	ylm[nlmp - 1] = std::complex<double>(0.0, 0.0);

	if (idot != 0) {

	if (idot != 0) {

		if (idot != 1) {

		if (idot != 1) {

			for (int n40 = 0; n40 < nsph; n40++) {

			for (int n40 = 0; n40 < nsph; n40++) {