manifold_ellipsoid.cpp
030f98e98a2alammps-icms

manifold_ellipsoid.cpp
View Options

			#include "manifold_ellipsoid.h"

			using namespace LAMMPS_NS;

			using namespace user_manifold;

			manifold_ellipsoid::manifold_ellipsoid( LAMMPS lmp, int narg, char *argv ) : manifold(lmp)
			{}


			double manifold_ellipsoid::g( const double *x )
			{
			const double ai = 1.0 / params[0];
			const double bi = 1.0 / params[1];
			const double ci = 1.0 / params[2];
			return x[0]x[0]aiai + x[1]x[1]bibi + x[2]x[2]ci*ci - 1.0;
			}

			void manifold_ellipsoid::n( const double x, double n )
			{
			const double ai = 1.0 / params[0];
			const double bi = 1.0 / params[1];
			const double ci = 1.0 / params[2];
			n[0] = 2x[0]ai*ai;
			n[1] = 2x[1]bi*bi;
			n[2] = 2x[2]ci*ci;
			}