From 1145cf6569db95d307f1e7a3d2d52ec080452584 Mon Sep 17 00:00:00 2001
From: Ardavan Oskooi <ardavan.oskooi@gmail.com>
Date: Tue, 22 Sep 2020 18:19:20 -0700
Subject: [PATCH] remove unused function output_mode_fields (#1367)

---
 doc/docs/Mode_Decomposition.md |  5 +----
 src/dft.cpp                    | 14 --------------
 src/meep.hpp                   |  1 -
 3 files changed, 1 insertion(+), 19 deletions(-)

diff --git a/doc/docs/Mode_Decomposition.md b/doc/docs/Mode_Decomposition.md
index 694812308..1b2f791eb 100644
--- a/doc/docs/Mode_Decomposition.md
+++ b/doc/docs/Mode_Decomposition.md
@@ -153,14 +153,11 @@ These functions are implemented in [src/mpb.cpp](https://github.com/NanoComp/mee
 ````
   void output_dft(dft_flux flux, const char *HDF5FileName);
 
-  void output_mode_fields(void *mode_data, dft_flux flux, const char *HDF5FileName);
 ````
 
 `output_dft` exports the components of the frequency-domain fields stored in `flux` to an HDF5 file with the given filename In general, `flux` will store data for fields at multiple frequencies.
 
-`output_mode_fields` is similar, but instead exports the components of the eigenmode described by `mode_data` which should be the return value of a call to `get_eigenmode`.
-
-These functions are implemented in [src/dft.cpp](https://github.com/NanoComp/meep/blob/master/src/dft.cpp#L982-L1021).
+This function is implemented in [src/dft.cpp](https://github.com/NanoComp/meep/blob/master/src/dft.cpp#L1184-L1196).
 
 ### Computing Overlap Integrals
 ````
diff --git a/src/dft.cpp b/src/dft.cpp
index 7d6f7314b..c9db81fac 100644
--- a/src/dft.cpp
+++ b/src/dft.cpp
@@ -1181,20 +1181,6 @@ void fields::output_dft(dft_fields fdft, const char *HDF5FileName) {
   output_dft_components(chunklists, 1, fdft.where, HDF5FileName);
 }
 
-/***************************************************************/
-/* does the same thing as output_dft(flux ...), but using      */
-/* eigenmode fields instead of dft_flux fields.                */
-/***************************************************************/
-void fields::output_mode_fields(void *mode_data, dft_flux flux, const char *HDF5FileName) {
-  h5file *file = open_h5file(HDF5FileName, h5file::WRITE);
-  delete file;
-
-  dft_chunk *chunklists[2];
-  chunklists[0] = flux.E;
-  chunklists[1] = flux.H;
-  FOR_E_AND_H(c) { process_dft_component(chunklists, 2, 0, c, 0, 0, 0, 0, 0, mode_data, 0, c); }
-}
-
 /***************************************************************/
 /***************************************************************/
 /***************************************************************/
diff --git a/src/meep.hpp b/src/meep.hpp
index f374eccb9..6ecc2b6c3 100644
--- a/src/meep.hpp
+++ b/src/meep.hpp
@@ -1897,7 +1897,6 @@ class fields {
   void output_dft(dft_force force, const char *HDF5FileName);
   void output_dft(dft_near2far n2f, const char *HDF5FileName);
   void output_dft(dft_fields fdft, const char *HDF5FileName);
-  void output_mode_fields(void *mode_data, dft_flux flux, const char *HDF5FileName);
 
   // get array of DFT field values
   std::complex<double> *get_dft_array(dft_flux flux, component c, int num_freq, int *rank,