Method optcost

constr.cpp:70–121

	}

	virtual double optcost( const double* const p )
	{
		double s1 = 0.0;
		double s2 = 0.0;
		double s3 = 0.0;
		int i;

		for( i = 0; i < 4; i++ )
		{
			s1 += p[ i ];
			s2 += sqr( p[ i ]);
		}

		for( i = 4; i < 13; i++ )
		{
			s3 += p[ i ];
		}

		double cost = 5.0*s1-5.0*s2-s3;

		const int n_con = 9;
		double pn[ n_con ];
		con_notmet = 0;

		pn[ 0 ] = penalty( 2.0*p[0]+2.0*p[1]+p[9]+p[10]-10.0 );
		pn[ 1 ] = penalty( 2.0*p[0]+2.0*p[2]+p[9]+p[11]-10.0 );
		pn[ 2 ] = penalty( 2.0*p[1]+2.0*p[2]+p[10]+p[11]-10.0 );
		pn[ 3 ] = penalty( -8.0*p[0]+p[9] );
		pn[ 4 ] = penalty( -8.0*p[1]+p[10] );
		pn[ 5 ] = penalty( -8.0*p[2]+p[11] );
		pn[ 6 ] = penalty( -2.0*p[3]-p[4]+p[9] );
		pn[ 7 ] = penalty( -2.0*p[5]-p[6]+p[10] );
		pn[ 8 ] = penalty( -2.0*p[7]-p[8]+p[11] );

		real_value = cost;

		if( con_notmet > 0 )
		{
			const double ps = 1.0 + 1.0 / n_con;
			double pns = 0.0;

			for( i = 0; i < n_con; i++ )
			{
				const double v = pn[ i ];
				pns = ( 1.0 + pns ) * ps + ( v + v * v + v * v * v ) * 0.33333;
			}

			cost += 1e10 * pns;
		}

		return( cost );
	}
};

int main()