From 5cf415700b117e5b79a17694dd4e90001e02aaa0 Mon Sep 17 00:00:00 2001
From: Adrian Hurtado <adrian.hurtado@stonybrook.edu>
Date: Wed, 20 Sep 2023 12:48:38 -0400
Subject: [PATCH] reduce the box size and number of points

---
 src/apps/molresponse/Plot_VTK.cc           | 6 +++---
 src/apps/molresponse/response_parameters.h | 2 +-
 2 files changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/apps/molresponse/Plot_VTK.cc b/src/apps/molresponse/Plot_VTK.cc
index 14faa7b50f2..99589e5916c 100644
--- a/src/apps/molresponse/Plot_VTK.cc
+++ b/src/apps/molresponse/Plot_VTK.cc
@@ -57,8 +57,8 @@ namespace madness {
                                        const response_matrix &responseMatrix) {
         // Stuff needed to plot
         //
-        Vector<double, 3> box_lo{-L / 4, -L / 4, -L / 4};
-        Vector<double, 3> box_hi{L / 4, L / 4, L / 4};
+        Vector<double, 3> box_lo{-L / 10, -L / 10, -L / 10};
+        Vector<double, 3> box_hi{L / 10, L / 10, L / 10};
 
 
         std::string vtk_dir = "vtk_plots";
@@ -159,7 +159,7 @@ namespace madness {
         //
         auto L_box = L / 10;
         Vector<double, 3> box_lo{-L_box, -L_box, -L_box};
-        Vector<double, 3> box_hi{L_box , L_box, L_box};
+        Vector<double, 3> box_hi{L_box, L_box, L_box};
 
 
         std::string vtk_dir = "vtk_plots";
diff --git a/src/apps/molresponse/response_parameters.h b/src/apps/molresponse/response_parameters.h
index f3cccc677e4..01ba38c0eb8 100644
--- a/src/apps/molresponse/response_parameters.h
+++ b/src/apps/molresponse/response_parameters.h
@@ -38,7 +38,7 @@ namespace madness {
             initialize<std::vector<double>>("plot_cell", std::vector<double>(),
                                             "lo-hi plot cell (default is all space)");
             initialize<double>("plot_l", -1.0, "Controls the plotting box size");
-            initialize<size_t>("plot_pts", 81, "Controls number of points in plots");
+            initialize<size_t>("plot_pts", 51, "Controls number of points in plots");
             initialize<bool>("plot_all_orbitals", false, "Turn on 2D plotting of response orbitals ");
             initialize<size_t>("maxiter", 25, "maximum number of iterations");
             initialize<double>("dconv", 1.e-4, "recommended values: 1.e-4 < dconv < 1.e-8");