From e5fcaffec1692618aa8a541a25a74776304cb6c6 Mon Sep 17 00:00:00 2001
From: Francois Gygi <fgygi@ucdavis.edu>
Date: Wed, 8 Aug 2018 14:42:21 -0700
Subject: [PATCH] Cleanup headers, add assert on basis size (inactive)

---
 src/ExternalPotential.C | 13 +++----------
 1 file changed, 3 insertions(+), 10 deletions(-)

diff --git a/src/ExternalPotential.C b/src/ExternalPotential.C
index 8f15a16..ac6ca35 100644
--- a/src/ExternalPotential.C
+++ b/src/ExternalPotential.C
@@ -16,14 +16,9 @@
 //
 ////////////////////////////////////////////////////////////////////////////////
 
-//#include <unistd.h>
-
 #include <iostream>
 #include <fstream>
 #include <cassert>
-#include <algorithm>
-#include <numeric>
-#include <functional>
 using namespace std;
 
 #include "Basis.h"
@@ -130,16 +125,14 @@ void ExternalPotential::update(const ChargeDensity& cd)
   basis.resize(cell,cell,ecut_);
 
   FourierTransform *vft = cd.vft();
-  // Following assertions are found to cause compatibility issues
-  // in Qbox-WEST coupling calculations
-  // assert(basis.np(0)<=vft->np0());
-  // assert(basis.np(1)<=vft->np1());
-  // assert(basis.np(2)<=vft->np2());
 
   FourierTransform ft1(basis,n_[0],n_[1],n_[2]);
   vext_read_loc.resize(ft1.np012loc());
   vector<complex<double> > vext_g(basis.localsize());
 
+  // check that the basis fits in the vft grid
+  //assert(basis.fits_in_grid(vft->np0(),vft->np1(),vft->np2()));
+
   FourierTransform ft2(basis,vft->np0(),vft->np1(),vft->np2());
   vext_r_.resize(ft2.np012loc());
 
--
libgit2 0.26.0