ibgdsplitl.c

有限元学习研究用源代码(老外的),供科研人员参考

/* last edit: Ilja Schmelzer -------------- 17-OCT-1994 19:20:32.33	*/
/************************************************************************/
/*                                                                      */
/*  <<< I B G >>> - Intersection - Based Grid generation package 	*/
/*                                                                      */
/*  Version 1.1 by Ilja Schmelzer   [e-mail]       */
/*                                                                      */
/*  to be distributed under IBG license conditions (see "readme.ibg")	*/
/*                                                                      */
/************************************************************************/
/* <<< IBGDSplitLine >>> - Intersection-Based Geometry Definition
				Splitting a Line
*/
#define const
#include "ibglib.h"
#include "ibgd.h"
#include "ibgi.h"
#include "ibgd0.h"
#include "ibgdefault.h"

typedef struct{
	ibGeometryRec	g;
	ibGeometry	gold;
	ibgSegment		lold,lnew,node;
	ibgFloat	(*func)(ibgPtObject fThis, ibgPoint *x);
       	ibgPtObject	fThis;
	ibgFloat	epsilon;
}ibgSplitLineRec,*ibgSplitLine;

static int RegionOfPoint(ibGeometry geom, ibgPoint *nnew, ibgPoint *nold)
{ibgSplitLine splt= (ibgSplitLine)geom;
 return ibgiRegionOfPoint(splt->gold,nnew,nold);
}

static int FaceWithEdge(ibGeometry g, ibgPoint *nint,
					ibgPoint *n1, ibgPoint *n2)
{ibgSplitLine splt = (ibgSplitLine) g;
 return ibgiFaceWithEdge(splt->gold,nint,n1,n2);
}

static int LineWithTriangle(ibGeometry g, ibgPoint *nint, ibgPoint *nface,
		    	ibgPoint *n1, ibgPoint *n2, ibgPoint *n3)
{ibgSplitLine splt = (ibgSplitLine ) g;
 int rc = ibgiLineWithTriangle(splt->gold,nint,nface,n1,n2,n3);
 if(ibgpSegment(*nint) == splt->lold){
 	if(splt->func(splt->fThis,nint)>0){
		ibgpSegment(*nint) = splt->lnew;
	}
 }
 return rc;
}

static int Free(ibGeometry g)
{ibgSplitLine splt = (ibgSplitLine ) g;
 ibgdFree(splt->gold);
 return ibgSuccess;
}

static ibGeometryClassRec	ibgSplitLineClass;

ibGeometry ibgdSplitLine(ibGeometry 	gold,
			ibgSegment		lnew,
			ibgSegment		lold,
			ibgFloat	(*func)(ibgPtObject fThis, ibgPoint *n),
			ibgPtObject	fThis,
			ibgFloat	Dfunc)
{ibgSplitLine splt=(ibgSplitLine) malloc(sizeof(ibgSplitLineRec));
 ibgdInitialize((ibGeometry)splt,&ibgSplitLineClass);
 ibgtCopy(&(splt->g.top),&(gold->top));
 splt->gold=gold;
 splt->lold=lold; splt->lnew=lnew; 
 splt->node=ibgDefaultNodeNr;
 splt->func=func; splt->fThis=fThis;
 splt->epsilon = 0.5*splt->g.Delta*Dfunc;
 ibgSplitLineClass.RegionOfPoint	= RegionOfPoint;
 ibgSplitLineClass.FaceWithEdge = FaceWithEdge;
 ibgSplitLineClass.LineWithTriangle = LineWithTriangle;
 ibgSplitLineClass.Free = Free;
 return (ibGeometry) splt;
}