diff --git a/src/aqdata/aq_base.cc b/src/aqdata/aq_base.cc
index 18d6cd98c..9a689003e 100644
--- a/src/aqdata/aq_base.cc
+++ b/src/aqdata/aq_base.cc
@@ -189,11 +189,5 @@ std::vector<Tensor<1, dim> > AQBase<dim>::get_all_directions ()
   return omega_i;
 }
 
-template <int dim>
-std::vector<double> AQBase<dim>::get_tensor_norms ()
-{
-  return tensor_norms;
-}
-
 template class AQBase<2>;
 template class AQBase<3>;
diff --git a/src/aqdata/aq_base.h b/src/aqdata/aq_base.h
index 40b9c7111..a4f7607c4 100644
--- a/src/aqdata/aq_base.h
+++ b/src/aqdata/aq_base.h
@@ -26,7 +26,6 @@ class AQBase
   unsigned int get_n_dir ();
   unsigned int get_n_total_ho_vars ();
   std::vector<double> get_angular_weights ();
-  std::vector<double> get_tensor_norms ();
   std::vector<Tensor<1, dim> > get_all_directions ();
   std::map<std::pair<unsigned int, unsigned int>, unsigned int>
   get_component_index_map ();
diff --git a/src/aqdata/lsgc.cc b/src/aqdata/lsgc.cc
index 027af6faf..ccfdf80e2 100644
--- a/src/aqdata/lsgc.cc
+++ b/src/aqdata/lsgc.cc
@@ -59,14 +59,6 @@ void LSGC<dim>::produce_angular_quad ()
                ExcMessage("calculated number of angles should be the same as number of angles"));
 
   this->n_total_ho_vars = this->n_dir * this->n_group;
-  // estimate tensor norm to do penalty method for EP
-  if (this->transport_model_name=="ep" &&
-      this->discretization=="dfem")
-    for (unsigned int i=0; i<this->n_dir; ++i)
-    {
-      Tensor<2, dim> tensor_tmp = outer_product(this->omega_i[i], this->omega_i[i]);
-      this->tensor_norms.push_back(tensor_tmp.norm());
-    }
 }
 
 template class LSGC<2>;
diff --git a/src/common/bart_driver.h b/src/common/bart_driver.h
index 5619d0790..8cc19d7c4 100644
--- a/src/common/bart_driver.h
+++ b/src/common/bart_driver.h
@@ -71,28 +71,10 @@ class BartDriver
   void assemble_lo_system ();
   void prepare_correction_aflx ();
   void initialize_ho_preconditioners ();
-  void ho_solve ();
-  void lo_solve ();
   void refine_grid ();
   void output_results () const;
-  void power_iteration ();
-  void initialize_fiss_process ();
-  void update_ho_moments_in_fiss ();
-  void update_fiss_source_keff ();
-  void source_iteration ();
-  void scale_fiss_transfer_matrices ();
-  void renormalize_sflx (std::vector<LA::MPI::Vector*> &target_sflxes);
-  void NDA_PI ();
-  void NDA_SI ();
   void initialize_aq (ParameterHandler &prm);
   
-  double estimate_k (double &fiss_source,
-                     double &fiss_source_prev_gen,
-                     double &k_prev_gen);
-  double estimate_fiss_source (std::vector<Vector<double> > &phis_this_process);
-  double estimate_phi_diff (std::vector<LA::MPI::Vector*> &phis_newer,
-                            std::vector<LA::MPI::Vector*> &phis_older);
-  
   const ParameterHandler &prm;
   
   std_cxx11::shared_ptr<ProblemDefinition> def_ptr;
@@ -108,24 +90,7 @@ class BartDriver
   std::string discretization;
   std::string namebase;
   std::string aq_name;
-  
-  unsigned int get_component_index (unsigned int incident_angle_index, unsigned int g);
-  unsigned int get_component_direction (unsigned int comp_ind);
-  unsigned int get_component_group (unsigned int comp_ind);
-  
-  unsigned int get_reflective_direction_index (unsigned int boundary_id,
-                                               unsigned int incident_angle_index);
-  
-  void radio (std::string str);
-  void radio (std::string str1, std::string str2);
-  void radio (std::string str1, unsigned int num1,
-              std::string str2, unsigned int num2,
-              std::string str3, unsigned int num3);
-  void radio (std::string str, double num);
-  void radio (std::string str, unsigned int num);
-  void radio (std::string str, bool boolean);
-  void radio ();
-  
+
   std::vector<typename DoFHandler<dim>::active_cell_iterator> local_cells;
   std::vector<typename DoFHandler<dim>::active_cell_iterator> ref_bd_cells;
   std::vector<bool> is_cell_at_bd;
@@ -138,9 +103,6 @@ class BartDriver
   std_cxx11::shared_ptr<FEFaceValues<dim> > fvf;
   std_cxx11::shared_ptr<FEFaceValues<dim> > fvf_nei;
   
-  
-  MPI_Comm mpi_communicator;
-  
   parallel::distributed::Triangulation<dim> triangulation;
   
   DoFHandler<dim> dof_handler;
@@ -148,17 +110,9 @@ class BartDriver
   IndexSet local_dofs;
   IndexSet relevant_dofs;
   
-  const double err_k_tol;
-  const double err_phi_tol;
-  const double err_phi_eigen_tol;
-  
   double ssor_omega;
   double keff;
-  double keff_prev_gen;
   double total_angle;
-  double c_penalty;
-  double fission_source;
-  double fission_source_prev_gen;
   
   bool is_eigen_problem;
   bool do_nda;
diff --git a/src/transport/even_parity.cc b/src/transport/even_parity.cc
index f3ba6700b..ecc580f84 100644
--- a/src/transport/even_parity.cc
+++ b/src/transport/even_parity.cc
@@ -11,7 +11,11 @@ TransportBase<dim>(prm, msh_ptr, aqd_ptr, mat_ptr)
 {
   if (this->discretization=="dfem")
   {
-    tensor_norms = aqd_ptr->get_tensor_norms ();
+    for (unsigned int i=0; i<this->n_dir; ++i)
+    {
+      Tensor<2, dim> tensor_tmp = outer_product(this->omega_i[i], this->omega_i[i]);
+      tensor_norms.push_back(tensor_tmp.norm());
+    }
     // note that the first constant on rhs is tunable and should be cautious
     // when tuning it.
     c_penalty = 1.0 * this->p_order * (this->p_order + 1.0);