threads/sem_init/s-c1.c

465 	struct row *Table = NULL;
550 	Table = calloc( N, sizeof( struct row ) );
552 	if ( Table == NULL )
566 		Table[ N ].X = ( long ) cur->nsem;
567 		Table[ N ].LnX = log( ( double ) cur->nsem );
571 			Table[ N ]._x = Table[ N ].X - Xavg ;
572 			Table[ N ]._lnx = Table[ N ].LnX - LnXavg;
573 			Table[ N ].Y_o = cur->_data_open;
574 			Table[ N ]._y_o = Table[ N ].Y_o - Yavg_o ;
575 			Table[ N ].LnY_o = log( ( double ) cur->_data_open );
576 			Table[ N ]._lny_o = Table[ N ].LnY_o - LnYavg_o;
577 			Table[ N ].Y_c = cur->_data_close;
578 			Table[ N ]._y_c = Table[ N ].Y_c - Yavg_c ;
579 			Table[ N ].LnY_c = log( ( double ) cur->_data_close );
580 			Table[ N ]._lny_c = Table[ N ].LnY_c - LnYavg_c;
604 		r1_o += ( ( double ) Table[ r ]._y_o / array_max ) * ( double ) Table[ r ]._y_o;
606 		_q_o[ 0 ] += Table[ r ]._y_o * Table[ r ]._x;
607 		_d_o[ 0 ] += Table[ r ]._x * Table[ r ]._x;
609 		_q_o[ 1 ] += Table[ r ]._lny_o * Table[ r ]._lnx;
610 		_d_o[ 1 ] += Table[ r ]._lnx * Table[ r ]._lnx;
612 		_q_o[ 2 ] += Table[ r ]._lny_o * Table[ r ]._x;
613 		_d_o[ 2 ] += Table[ r ]._x * Table[ r ]._x;
616 		r1_c += ( ( double ) Table[ r ]._y_c / array_max ) * ( double ) Table[ r ]._y_c;
618 		_q_c[ 0 ] += Table[ r ]._y_c * Table[ r ]._x;
619 		_d_c[ 0 ] += Table[ r ]._x * Table[ r ]._x;
621 		_q_c[ 1 ] += Table[ r ]._lny_c * Table[ r ]._lnx;
622 		_d_c[ 1 ] += Table[ r ]._lnx * Table[ r ]._lnx;
624 		_q_c[ 2 ] += Table[ r ]._lny_c * Table[ r ]._x;
625 		_d_c[ 2 ] += Table[ r ]._x * Table[ r ]._x;
642 		t = ( Table[ r ]._y_o - ( ( _q_o[ 0 ] * Table[ r ]._x ) / _d_o[ 0 ] ) );
645 		t = ( Table[ r ]._y_c - ( ( _q_c[ 0 ] * Table[ r ]._x ) / _d_c[ 0 ] ) );
652 		t = ( Table[ r ].Y_o
654 		          * powl( Table[ r ].X, ( _q_o[ 1 ] / _d_o[ 1 ] ) )
658 		t = ( Table[ r ].Y_c
660 		          * powl( Table[ r ].X, ( _q_c[ 1 ] / _d_c[ 1 ] ) )
669 		t = ( Table[ r ].Y_o
670 		      - expl( ( _q_o[ 2 ] / _d_o[ 2 ] ) * Table[ r ]._x + LnYavg_o ) );
673 		t = ( Table[ r ].Y_c
674 		      - expl( ( _q_c[ 2 ] / _d_c[ 2 ] ) * Table[ r ]._x + LnYavg_c ) );
775 	free( Table );