Update JonswapSpectrum.cpp

Should be frequency squared in the numerator of the gamma power
mhscott · Nov 16, 2024 · f979699 · f979699
1 parent 1474bc4
commit f979699
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/SRC/reliability/domain/spectrum/JonswapSpectrum.cpp b/SRC/reliability/domain/spectrum/JonswapSpectrum.cpp
@@ -83,7 +83,7 @@ JonswapSpectrum::getAmplitude(double frequency)
 		sigma = 0.09;
 	}
 
-	double power = exp(-(frequency-wp)/(2.0*sigma*sigma*wp*wp));
+	double power = exp(-(frequency-wp)*(frequency-wp)/(2.0*sigma*sigma*wp*wp));
 	double GAMMA = pow(gamma,power);
 
 	return GAMMA*alpha*pow(frequency,-5.0)*exp(-5.0/4.0*pow((frequency/wp),-4.0));