shoot_boundary.c
View Options

			/* $Id: shoot_boundary.c 4875 2012-08-30 20:04:30Z john $ */
			/* (c) John Ashburner (2011) */

			#include "mex.h"
			#include "shoot_boundary.h"

			/* Neumann boundary condition */
			static mwSignedIndex neumann_boundary(mwSignedIndex i, mwSize m)
			{
			mwSignedIndex m2 = m*2;
			i = (i<0) ? m2-((-i-1)%m2)-1 : (i%m2);
			return((m<=i)? m2-i-1: i);
			}

			static mwSignedIndex circulant_boundary(mwSignedIndex i, mwSize m)
			{
			return((i>=0) ? i%((signed)m) : (((signed)m)+i%((signed)m))%(signed)m);
			}

			mwSignedIndex (*bound)() = circulant_boundary;
			static int bound_type = BOUND_CIRCULANT;

			void set_bound(int t)
			{
			bound_type = t;
			if (t==BOUND_CIRCULANT)
			bound = circulant_boundary;
			else if (t==BOUND_NEUMANN)
			bound = neumann_boundary;
			else
			mexErrMsgTxt("Undefined boundary condition.");
			}

			int get_bound()
			{
			return(bound_type);
			}

Slepi/Code/Utilities/spm12/src/shoot_boundary.c8a186cd82567master