implement_new_mate…

solver_common.cc
View Options

			/**
			* file solver_common.cc
			*
			* @author Till Junge <till.junge@epfl.ch>
			*
			* @date 15 May 2018
			*
			* @brief implementation for solver utilities
			*
			* Copyright © 2018 Till Junge
			*
			* µSpectre is free software; you can redistribute it and/or
			* modify it under the terms of the GNU General Public License as
			* published by the Free Software Foundation, either version 3, or (at
			* your option) any later version.
			*
			* µSpectre is distributed in the hope that it will be useful, but
			* WITHOUT ANY WARRANTY; without even the implied warranty of
			* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
			* General Public License for more details.
			*
			* You should have received a copy of the GNU General Public License
			* along with GNU Emacs; see the file COPYING. If not, write to the
			* Free Software Foundation, Inc., 59 Temple Place - Suite 330,
			* Boston, MA 02111-1307, USA.
			*/

			#include "solver/solver_common.hh"
			namespace muSpectre {

			/* ---------------------------------------------------------------------- */
			bool check_symmetry(const Eigen::Ref<const Eigen::ArrayXXd>& eps,
			Real rel_tol){
			return (rel_tol >= (eps-eps.transpose()).matrix().norm()/eps.matrix().norm() \|\|
			rel_tol >= eps.matrix().norm());
			}

			} // muSpectre