Add source for version 7.1.0.012

SourceCode/DElight/SourceCode/Radiosity.cpp

@@ -19,24 +19,24 @@
  * under contract with Lawrence Berkeley National Laboratory.
  **************************************************************/
 
-// This work was supported by the Assistant Secretary for Energy Efficiency 
-// and Renewable Energy, Office of Building Technologies, 
-// Building Systems and Materials Division of the 
+// This work was supported by the Assistant Secretary for Energy Efficiency
+// and Renewable Energy, Office of Building Technologies,
+// Building Systems and Materials Division of the
 // U.S. Department of Energy under Contract No. DE-AC03-76SF00098.
 
 /*
-NOTICE: The Government is granted for itself and others acting on its behalf 
-a paid-up, nonexclusive, irrevocable worldwide license in this data to reproduce, 
-prepare derivative works, and perform publicly and display publicly. 
+NOTICE: The Government is granted for itself and others acting on its behalf
+a paid-up, nonexclusive, irrevocable worldwide license in this data to reproduce,
+prepare derivative works, and perform publicly and display publicly.
 Beginning five (5) years after (date permission to assert copyright was obtained),
-subject to two possible five year renewals, the Government is granted for itself 
+subject to two possible five year renewals, the Government is granted for itself
 and others acting on its behalf a paid-up, nonexclusive, irrevocable worldwide
-license in this data to reproduce, prepare derivative works, distribute copies to 
-the public, perform publicly and display publicly, and to permit others to do so. 
+license in this data to reproduce, prepare derivative works, distribute copies to
+the public, perform publicly and display publicly, and to permit others to do so.
 NEITHER THE UNITED STATES NOR THE UNITED STATES DEPARTMENT OF ENERGY, NOR ANY OF
-THEIR EMPLOYEES, MAKES ANY WARRANTY, EXPRESS OR IMPLIED, OR ASSUMES ANY LEGAL 
-LIABILITY OR RESPONSIBILITY FOR THE ACCURACY, COMPLETENESS, OR USEFULNESS OF ANY 
-INFORMATION, APPARATUS, PRODUCT, OR PROCESS DISCLOSED, OR REPRESENTS THAT ITS USE 
+THEIR EMPLOYEES, MAKES ANY WARRANTY, EXPRESS OR IMPLIED, OR ASSUMES ANY LEGAL
+LIABILITY OR RESPONSIBILITY FOR THE ACCURACY, COMPLETENESS, OR USEFULNESS OF ANY
+INFORMATION, APPARATUS, PRODUCT, OR PROCESS DISCLOSED, OR REPRESENTS THAT ITS USE
 WOULD NOT INFRINGE PRIVATELY OWNED RIGHTS.
 */
 #pragma warning(disable:4786)
@@ -170,7 +170,7 @@ int slite_interreflect(
 		if ((iRefptIllumRetVal = refpt_total_illum(bldg_ptr,sun_ptr,iz,pofdmpfile)) < 0) {
             // If errors were detected then return now, else register warnings and continue processing
             if (iRefptIllumRetVal != -10) {
-				*pofdmpfile << "ERROR: DElight Bad return from refpt_total_illum()\n"; 
+				*pofdmpfile << "ERROR: DElight Bad return from refpt_total_illum()\n";
 				return(-1);
             }
             else {
@@ -258,7 +258,7 @@ int surf_interreflect(
 				delf_overcast[inode] += fij * bldg_ptr->zone[iz]->surf[jsurf]->skyolum[jnode];
 
 // rjh debug
-//*pofdmpfile << "Configuration Factor: fij = " << fij << " isurf = " << isurf << " inode = " << inode << " jsurf = " << jsurf << " jnode = " << jnode << "\n"; 
+//*pofdmpfile << "Configuration Factor: fij = " << fij << " isurf = " << isurf << " inode = " << inode << " jsurf = " << jsurf << " jnode = " << jnode << "\n";
 
 				/* for each Sun Position Altitude */
 				for (iphs=0; iphs<sun_ptr->nphs; iphs++) {
@@ -387,12 +387,15 @@ int wndo_interreflect(
 	/* improve values for total node luminance for each node on current window */
 	for (inode=0; inode<bldg_ptr->zone[iz]->surf[is]->nnodes; inode++) {
 		/* for overcast sky condition */
+        if ( _isnan(delf_overcast[inode]) )delf_overcast[inode]=0.;
 		bldg_ptr->zone[iz]->surf[is]->wndo[iw]->skyolum[inode] = bldg_ptr->zone[iz]->surf[is]->wndo[iw]->direct_skyolum[inode] + frac * delf_overcast[inode];
 		/* for each Sun Position Altitude */
 		for (iphs=0; iphs<sun_ptr->nphs; iphs++) {
 			/* for each Sun Position Azimuth */
 			for (iths=0; iths<sun_ptr->nths; iths++) {
 				/* for each clear sky sun position */
+				if ( _isnan(delf_skyclear[inode][iphs][iths]) )delf_skyclear[inode][iphs][iths]=0.;
+				if ( _isnan(delf_sunclear[inode][iphs][iths]) )delf_sunclear[inode][iphs][iths]=0.;
 				bldg_ptr->zone[iz]->surf[is]->wndo[iw]->skyclum[inode][iphs][iths] = bldg_ptr->zone[iz]->surf[is]->wndo[iw]->direct_skyclum[inode][iphs][iths] + frac * delf_skyclear[inode][iphs][iths];
 				bldg_ptr->zone[iz]->surf[is]->wndo[iw]->sunclum[inode][iphs][iths] = bldg_ptr->zone[iz]->surf[is]->wndo[iw]->direct_sunclum[inode][iphs][iths] + frac * delf_sunclear[inode][iphs][iths];
 			}
@@ -442,7 +445,7 @@ int refpt_total_illum(
 
 	/* calc direction cosine values */
 	/* Note: see calc_dircos() for complete logic for arbitrary surfaces */
-	// RJH - 8/24/03 - Only Z-axis direction cosines are needed for ref pts 
+	// RJH - 8/24/03 - Only Z-axis direction cosines are needed for ref pts
 	// and they are always 0,0,1 since ref pts are assumed to be on horizontal surface facing upward.
 	for (int idircos=0; idircos<NDC; idircos++)
 	{