# exafmm

committed d649325

Fix for double precision.

# include/serialfmm.h

`   void setBounds(Bodies &bodies) {`
`     startTimer("Set bounds");                                   // Start timer`
`     vec3Pair bounds = getBounds(bodies.begin(), bodies.end());  // Get Xmin (first) and Xmax (second) of domain`
`-    localXmin = bounds.first;                                   // Set local Xmin`
`-    localXmax = bounds.second;                                  // Set local Xmax`
`-    localCenter = (localXmax + localXmin) / 2;                  // Calculate center of domain`
`+    for (int d=0; d<3; d++) localXmin[d] = bounds.first[d];     // Set local Xmin`
`+    for (int d=0; d<3; d++) localXmax[d] = bounds.second[d];    // Set local Xmax`
`+    for (int d=0; d<3; d++) localCenter[d] = (localXmax[d] + localXmin[d]) / 2;// Calculate center of domain`
`     localRadius = 0;                                            // Initialize localRadius`
`     for (int d=0; d<3; d++) {                                   // Loop over dimensions`
`       localRadius = std::max(localCenter[d] - localXmin[d], localRadius);// Calculate min distance from center`

# include/types.h

` typedef std::pair<vec3,vec3> vec3Pair;                          //!< Pair of vec3`
` `
` // Compile-time parameters`
`-const int P = 3;                                                //!< Order of expansions`
`+const int P = 5;                                                //!< Order of expansions`
` const float EPS2 = .0;                                          //!< Softening parameter (squared)`
` #if COMkernel`
` const int MTERM = P*(P+1)*(P+2)/6-3;                            //!< Number of Cartesian mutlipole terms`
Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for ProjectModifiedEvent.java.
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.