Corrected calculation of the gradient of the WendlandQuinticC2Kernel. See also Line 591. (#347)

holgergoe · web-flow · commit e773d3e26655 · 2025-02-04T08:39:32.000+01:00
diff --git a/SPlisHSPlasH/SPHKernels.h b/SPlisHSPlasH/SPHKernels.h
@@ -315,7 +315,7 @@ namespace SPH
 			if (q <= 1.0)
 			{
 				const Vector3r gradq = r * (static_cast<Real>(1.0) / (rl*m_radius));
-				res = m_l*gradq*pow(static_cast<Real>(1.0) - q, 3);
+				res = m_l*q*pow(static_cast<Real>(1.0) - q, static_cast<Real>(3.0))*gradq;
 			}
 			else
 				res.setZero();

Original file line number	Diff line number	Diff line change
`@@ -315,7 +315,7 @@ namespace SPH`
`315`	`315`	`if (q <= 1.0)`
`316`	`316`	`{`
`317`	`317`	`const Vector3r gradq = r * (static_cast<Real>(1.0) / (rl*m_radius));`
`318`		`- res = m_lgradqpow(static_cast<Real>(1.0) - q, 3);`
	`318`	`+ res = m_lqpow(static_cast<Real>(1.0) - q, static_cast<Real>(3.0))*gradq;`
`319`	`319`	`}`
`320`	`320`	`else`
`321`	`321`	`res.setZero();`