Changeset 15741 for issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.cpp

Timestamp:

08/07/13 11:25:05 (12 years ago)

Author:

Mathieu Morlighem

Message:

CHG: keep the exponent as an integer may execute faster (compiler dependent)

File:

• issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.cpp

@@ -1452 +1452 @@
                                           *(5.+3.*egaus[j]-zgaus[k]-3.*egaus[j]*zgaus[k]);
                                         zeta      =sqrt(2./3.)*(1.+zgaus[k]);
-                                        (*pwgt)[ipt]=xwgt[i]*ewgt[j]*zwgt[k]
-                                          *(SQRT2/16.
-                                                                  *(1.-egaus[j])*pow(1.-zgaus[k],2.));
+                                        (*pwgt)[ipt]=xwgt[i]*ewgt[j]*zwgt[k]*(SQRT2/16.*(1.-egaus[j])*pow(1.-zgaus[k],2));
 
                                         (*pl1 )[ipt]=(1.-xi-eta/SQRT3-zeta/sqrt(6.))/2.;