Added missing return values

openmmapi/src/CMAPTorsionForceImpl.cpp

@@ -55,6 +55,7 @@ void CMAPTorsionForceImpl::initialize(ContextImpl& context) {
 double CMAPTorsionForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcCMAPTorsionForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 vector<string> CMAPTorsionForceImpl::getKernelNames() {

openmmapi/src/CustomAngleForceImpl.cpp

@@ -91,6 +91,7 @@ void CustomAngleForceImpl::initialize(ContextImpl& context) {
 double CustomAngleForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcCustomAngleForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 vector<string> CustomAngleForceImpl::getKernelNames() {

openmmapi/src/CustomBondForceImpl.cpp

@@ -85,6 +85,7 @@ void CustomBondForceImpl::initialize(ContextImpl& context) {
 double CustomBondForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcCustomBondForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 vector<string> CustomBondForceImpl::getKernelNames() {

openmmapi/src/CustomExternalForceImpl.cpp

@@ -79,6 +79,7 @@ void CustomExternalForceImpl::initialize(ContextImpl& context) {
 double CustomExternalForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcCustomExternalForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 vector<string> CustomExternalForceImpl::getKernelNames() {

openmmapi/src/CustomGBForceImpl.cpp

@@ -108,6 +108,7 @@ void CustomGBForceImpl::initialize(ContextImpl& context) {
 double CustomGBForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcCustomGBForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 vector<string> CustomGBForceImpl::getKernelNames() {

openmmapi/src/CustomHbondForceImpl.cpp

@@ -181,6 +181,7 @@ void CustomHbondForceImpl::initialize(ContextImpl& context) {
 double CustomHbondForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcCustomHbondForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 vector<string> CustomHbondForceImpl::getKernelNames() {

openmmapi/src/CustomNonbondedForceImpl.cpp

@@ -108,6 +108,7 @@ void CustomNonbondedForceImpl::initialize(ContextImpl& context) {
 double CustomNonbondedForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcCustomNonbondedForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 vector<string> CustomNonbondedForceImpl::getKernelNames() {

openmmapi/src/CustomTorsionForceImpl.cpp

@@ -97,6 +97,7 @@ void CustomTorsionForceImpl::initialize(ContextImpl& context) {
 double CustomTorsionForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcCustomTorsionForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 vector<string> CustomTorsionForceImpl::getKernelNames() {

openmmapi/src/GBSAOBCForceImpl.cpp

@@ -58,6 +58,7 @@ void GBSAOBCForceImpl::initialize(ContextImpl& context) {
 double GBSAOBCForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcGBSAOBCForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 std::vector<std::string> GBSAOBCForceImpl::getKernelNames() {

openmmapi/src/GBVIForceImpl.cpp

@@ -242,6 +242,7 @@ void GBVIForceImpl::findScaledRadii( int numberOfParticles, const std::vector<st
 double GBVIForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcGBVIForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 std::vector<std::string> GBVIForceImpl::getKernelNames() {

openmmapi/src/HarmonicAngleForceImpl.cpp

@@ -52,6 +52,7 @@ void HarmonicAngleForceImpl::initialize(ContextImpl& context) {
 double HarmonicAngleForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcHarmonicAngleForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 std::vector<std::string> HarmonicAngleForceImpl::getKernelNames() {

openmmapi/src/HarmonicBondForceImpl.cpp

@@ -52,6 +52,7 @@ void HarmonicBondForceImpl::initialize(ContextImpl& context) {
 double HarmonicBondForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcHarmonicBondForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 std::vector<std::string> HarmonicBondForceImpl::getKernelNames() {

openmmapi/src/PeriodicTorsionForceImpl.cpp

@@ -52,6 +52,7 @@ void PeriodicTorsionForceImpl::initialize(ContextImpl& context) {
 double PeriodicTorsionForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcPeriodicTorsionForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 std::vector<std::string> PeriodicTorsionForceImpl::getKernelNames() {

openmmapi/src/RBTorsionForceImpl.cpp

@@ -52,6 +52,7 @@ void RBTorsionForceImpl::initialize(ContextImpl& context) {
 double RBTorsionForceImpl::calcForcesAndEnergy(ContextImpl& context, bool includeForces, bool includeEnergy, int groups) {
     if ((groups&(1<<owner.getForceGroup())) != 0)
         return dynamic_cast<CalcRBTorsionForceKernel&>(kernel.getImpl()).execute(context, includeForces, includeEnergy);
+    return 0.0;
 }
 
 std::vector<std::string> RBTorsionForceImpl::getKernelNames() {