Esempio n. 1
0
void splint2_clip(double *x1a,double *x2a,double *ya,double *y2a,int32 m,int32 n,double x2,double *x1,int32 nx1,double *y)
{
double *y1temp;
double *yy1temp;
double *y2temp;
double *yy2temp;
int i,j,k;



y2temp = new double[m];
yy2temp = new double[m];


for (i=0; i<m; i++)
	{
	splint_0_clip(x2a,ya+i*n,y2a+i*n,n,x2,yy2temp + i);
	}
spline_0(x1a,yy2temp,m,1.0e30,1.0e30,y2temp);
for (k=0; k<nx1; k++)
	{
	splint_0_clip(x1a,yy2temp,y2temp,m,x1[k],y+k);
	}
	

delete y2temp;
delete yy2temp;
} 
Esempio n. 2
0
void spline2(double *x1a,double *x2a,double *ya,int32 m,int32 n, double *y2a)
{
int i,j;

for (i=0; i<m; i++)
	{
	spline_0(x2a,ya+i*n,n,1.0e30,1.0e30,y2a+i*n);
	}
}
Esempio n. 3
0
//save to input black table when window is closed
McoStatus	PrinterType::setblack(McoHandle blacktableH)
{
	short 	i;
	double *black;
	double *y2_hands;
	double ypn;
	
	y2_hands = new double[num_hands];
	
	black = (double*)McoLockHandle(blacktableH);
	
	ypn=(y_hands[num_hands-1]-y_hands[num_hands-2])/
		(x_hands[num_hands-1]-x_hands[num_hands-2]);
	
	spline_0(x_hands,y_hands,num_hands,1E30,1E30,y2_hands);
//L to Black table	
	for( i = 0; i <= 100; i++)
		splint_0(x_hands,y_hands,y2_hands,num_hands,i, black+i);

//fix on 8/23
//clipping the _yd_points which on the left of letfmost hand
//to the value of the leftmost hand
//on the right of rightmost hand to the value of the rightmost hand
	for( i = 0; i <= 100; i++){
		if( i <= x_hands[0] ) //left of the leftmost hand
			black[i] = y_hands[0];
		if( i >= x_hands[num_hands-1] ) //left of the leftmost hand
			black[i] = y_hands[num_hands-1];
	}

//clipping the table and normaliz to 1
	for( i = 0; i <= 100; i++){
		if(black[i] < 0)
			black[i] = 0;
		else if(black[i] > 100)
			black[i] = 1.0;
		else
			black[i] /= 100.0;
	}

				
	McoUnlockHandle(blacktableH);
	
	delete y2_hands;
	
	return MCO_SUCCESS;
}	
Esempio n. 4
0
// use splines to linearize tables
// the polaroid version (data is in density format)
McoStatus LinearData::setUp(double *data,double *vals,int array_mult,int nc,int *s)
{

int i,j;
double t,t2;
McoStatus state = MCO_SUCCESS;
double z,Ts,Tt;

CleanUp();

numc = nc;

sizes = new int32[numc];
curves = new double*[numc];
values = new double*[numc];
y2 = new double*[numc];
y3 = new double*[numc];
if ((sizes == 0L) || (curves == 0L) || (values == 0L) || (y2 == 0L) || (y3 == 0L)) return MCO_MEM_ALLOC_ERROR;

for (i=0; i<numc; i++) 
	{
	sizes[i] = s[i];
	curves[i] = new double[sizes[i]];
	values[i] = new double[sizes[i]];
	y2[i] = new double[sizes[i]];
	y3[i] = new double[sizes[i]];
	if ((sizes[i] == 0L) || (curves[i] == 0L) || (values[i] == 0L) || (y2[i] == 0L) || (y3[i] == 0L)) return MCO_MEM_ALLOC_ERROR;
	for (j=0; j<sizes[i]; j++)
		{
		values[i][j] = vals[i*array_mult+j];
		}
	}

// build the curves
for (i=0; i<numc; i++)
	{
	for (j=0; j<sizes[i]; j++)
		{
		curves[i][j] = data[i*array_mult+j];
		}
	// remve any reversals from the data
	state = cleanReversals(curves[i],values[i],&sizes[i]);
	if (state != MCO_SUCCESS) goto bail;
	// use %dot

	z = curves[i][0];
	Ts = pow(10,(-curves[i][sizes[i]-1]+z));
	for (j=0; j<sizes[i]; j++)
		{
		 Tt = pow(10,(-curves[i][j]+z));
		 curves[i][j] = 100*(1-Tt)/(1-Ts);
		}	 		
			
	// set up the splines
	spline_0(values[i],curves[i],sizes[i],1.0e30,1.0e30,y2[i]);	
	
	// set up the splines
	spline_0(curves[i],values[i],sizes[i],1.0e30,1.0e30,y3[i]);	
	}
	
	
bail:
return state;
}
Esempio n. 5
0
// use splines to linearize tables
McoStatus LinearData::setUp(double *data,double *vals, double *paper,int nc,int s)
{

int i,j;
double t,t2;
McoStatus state = MCO_SUCCESS;
int		LAB = 0;
double  z,Ts,Tt;

CleanUp();

numc = nc;

sizes = new int32[numc];
curves = new double*[numc];
values = new double*[numc];
y2 = new double*[numc];
y3 = new double*[numc];
if ((sizes == 0L) || (curves == 0L) || (values == 0L) || (y2 == 0L) || (y3 == 0L)) return MCO_MEM_ALLOC_ERROR;

for (i=0; i<numc; i++) 
	{
	sizes[i] = s;
	curves[i] = new double[s];
	values[i] = new double[s];
	y2[i] = new double[s];
	y3[i] = new double[s];
	if ((sizes[i] == 0L) || (curves[i] == 0L) || (values[i] == 0L) || (y2[i] == 0L) || (y3[i] == 0L)) return MCO_MEM_ALLOC_ERROR;
	for (j=0; j<s; j++)
		{
		values[i][j] = vals[j];
		}
	}

// determine if data is in density format
// must have paper density subtracted
for (i=0; i<numc; i++)
	{
	for (j=0; j<sizes[i]-1; j++)
		{
		if (fabs(data[i*(s-1)*3+j*3+1]) > 0.01) LAB = 1;
		if (fabs(data[i*(s-1)*3+j*3+2]) > 0.01) LAB = 1;
		}
	}

// build the curves
for (i=0; i<numc; i++)
	{
	curves[i][0] = 0;	// 0%
	for (j=0; j<sizes[i]-1; j++)
		{
		if (LAB)
			{
			t = data[i*(s-1)*3+j*3]-paper[0];	
			t2 = t*t;
			t = data[i*(s-1)*3+j*3+1]-paper[1];	
			t2 += t*t;
			t = data[i*(s-1)*3+j*3+2]-paper[2];	
			t2 += t*t;
			curves[i][j+1] = sqrt(t2);
			}
		else curves[i][j+1] = data[i*(s-1)*3+j*3];
		}
	// remve any reversals from the data
	state = cleanReversals(curves[i],values[i],&sizes[i]);
	if (state != MCO_SUCCESS) goto bail;
	// scale from 0 to 100 if lab
	if (LAB) for (j=0; j<sizes[i]; j++) curves[i][j] = 100*curves[i][j]/curves[i][sizes[i]-1];	
	else 
		{
		z = curves[i][0];
		Ts = pow(10,(-curves[i][sizes[i]-1]+z));
		for (j=0; j<sizes[i]; j++)
			{
			 Tt = pow(10,(-curves[i][j]+z));
			 curves[i][j] = 100*(1-Tt)/(1-Ts);
			}	 		
		}
	
	// set up the splines
	spline_0(values[i],curves[i],sizes[i],1.0e30,1.0e30,y2[i]);	
	
	// set up the splines
	spline_0(curves[i],values[i],sizes[i],1.0e30,1.0e30,y3[i]);	
	}
	
	
bail:
return state;
}