/*****************************************************************************
**
**    consistency.c                   NQ                       Werner Nickel
**                                         Werner.Nickel@math.rwth-aachen.de
*/


#include "nq.h"

static
void	printEv( ev )
expvec	ev;

{	long	i;

	for( i = 1; i <= NrPcGens+NrCenGens; i++ )
	    printf( "%2d ", ev[i] );
}

do_cba( c, b, a )
gen	c, b, a;

{	int	i;
	expvec	ev1, ev2;

	if( Wt(c)+Wt(b)+Wt(abs(a)) > Class+1 ) return;

	/* the left hand side first : c (b a) = c a b^a */
	ev1 = ExpVecWord( Generators[c] );
	Collect( ev1, Generators[a], 1 );
	Collect( ev1, Conjugate[b][a], 1 );

	/* then the right hand side : (c b) a */
	ev2 = ExpVecWord( Generators[c] );
	Collect( ev2, Generators[b], 1 );
	Collect( ev2, Generators[a], 1 );

	for( i = 1; i <= NrPcGens+NrCenGens; i++ ) ev1[i] -= ev2[i];
	free( ev2 );

	if( Debug ) {
	    printf( "cba( %2d %2d %2d ) : ",c,b,a);
	    printEv( ev1 );
	    printf("\n");
	}
	addRow( ev1 );
}

do_cbb( c, b )
gen	c, b;

{	expvec	ev;

	if( Wt(c)+Wt(b) > Class+1 ) return;

	/* (c b) b^-1 */
	ev = ExpVecWord( Generators[c] );
	Collect( ev, Generators[ b], 1 );
	Collect( ev, Generators[-b], 1 );
	ev[ c ] -= 1;

	if( Debug ) {
	    printf( "cbb( %2d %2d    ) : ", c, b );
	    printEv( ev );
	    printf("\n");
	}
	addRow( ev );

        if( EarlyStop ) return;

	/* (c b^-1) b */
	ev = ExpVecWord( Generators[c] );
	Collect( ev, Generators[-b], 1 );
	Collect( ev, Generators[ b], 1 );
	ev[ c ] -= 1;

	if( Debug ) {
	    printf( "cbb( %2d %2d    ) : ", c, -b );
	    printEv( ev );
	    printf("\n");
	}
	addRow( ev );
}

do_ccb( c, b )
gen	c, b;

{	expvec	ev;

	if( Wt(c)+Wt(abs(b)) > Class+1 ) return;

	/* c^-1 (c b) = c^-1 b c^b */
	ev = ExpVecWord( Generators[-c] );
	Collect( ev, Generators[b], 1 );
	Collect( ev, Conjugate[c][b], 1 );
	ev[abs(b)] -= sgn(b);

	if( Debug ) {
	    printf( "ccb( %2d %2d    ) : ", c, b );
	    printEv( ev );
	    printf("\n");
	}
	addRow( ev );
}

do_cbn( c, b )
gen	c, b;

{	int	i;
	expvec	ev1, ev2;

	if( Wt(c)+Wt(b) > Class+1 ) return;

	/* c (b^n) */
	ev1 = ExpVecWord( Generators[c] );
	Collect( ev1, Power[b], 1 );

	/* (c b) b^(n-1) */
	ev2 = ExpVecWord( Generators[c] );
	Collect( ev2, Generators[b], Exponent[b] );

	for( i = 1; i <= NrPcGens+NrCenGens; i++ ) ev1[i] -= ev2[i];

	free( ev2 );
	if( Debug ) {
	    printf( "cbn( %2d %2d    ) : ", c, b );
	    printEv( ev1 );
	    printf("\n");
	}
	addRow( ev1 );
}

do_cnb( c, b )
gen	c, b;

{	int	i;
	expvec	ev1, ev2;

	if( Wt(c)+Wt(abs(b)) > Class+1 ) return;

	/* (c^n) b */
	ev1 = ExpVecWord( Power[c] );
	Collect( ev1, Generators[b], 1 );

	/* c^(n-1) (c b) = c^(n-1) b c^b */
	ev2 = ExpVecWord( Generators[c] );
	if( Exponent[c] > 2 ) Collect( ev2, Generators[c], Exponent[c]-2 );
	Collect( ev2, Generators[b], 1 );
	Collect( ev2, Conjugate[c][b], 1 );

	for( i = 1; i <= NrPcGens+NrCenGens; i++ ) ev1[i] -= ev2[i];

	free( ev2 );
	if( Debug ) {
	    printf( "cnb( %2d %2d    ) : ", c, b );
	    printEv( ev1 );
	    printf("\n");
	}
	addRow( ev1 );
}

do_cnc( c )
gen	c;

{	int	i;
	expvec	ev1, ev2;

	if( 2*Wt(c) > Class+1 ) return;

	/* c^n c */
	ev1 = ExpVecWord( Power[c] );
	Collect( ev1, Generators[c], 1 );

	/* c c^n */
	ev2 = ExpVecWord( Generators[c] );
	Collect( ev2, Power[c], 1 );

	for( i = 1; i <= NrPcGens+NrCenGens; i++ ) ev1[i] -= ev2[i];

	free( ev2 );
	if( Debug ) {
	    printf( "cnb( %2d       ) : ", c );
	    printEv( ev1 );
	    printf("\n");
	}
	addRow( ev1 );
}

void	Consistency() {

	long	t;
	gen	a, b, c;

	if( Verbose ) t = RunTime();

	/*    c * ( b * a ) = ( c * b ) * a  for all generators c > b > a.
	*/
	for( a = 1; !EarlyStop && a <= Dimension[1]; a++ )
	    for( b = a+1; !EarlyStop && b <= NrPcGens; b++ )
		for( c = b+1; !EarlyStop && c <= NrPcGens; c++ )
		    do_cba( c, b, a );


	/*
	**    c * ( b  * b' ) = ( c * b  ) * b' and
	**    c * ( b' * b  ) = ( c * b' ) * b  for all c > b,
	**                                      Exponent[b] == 0
	*/
	for( b = 1; !EarlyStop && b <= NrPcGens; b++ )
	    if( Exponent[b] == 0 )
		for( c = b+1; !EarlyStop && c <= NrPcGens; c++ )
		    do_cbb( c, b );

	/*
	**    c * ( c' * b ) = ( c * c' ) * b  and
	**    c * ( c' * b ) = ( c * c' ) * b' for all generators c > b,
	**                                     Exponent[c] == 0.
	*/
	for( c = 1; !EarlyStop && c <= NrPcGens; c++ )
	    if( Exponent[c] == 0 ) {
		for( b = 1; !EarlyStop && b <= min(c-1,Dimension[1]); b++ )
		    do_ccb( c, b );
		for( b = 1; !EarlyStop && b < c; b++ )
		    if( Exponent[b] == 0 ) do_ccb( c, -b );
	    }

	/*
	**    c * b^n = ( c * b ) * b^(n-1) for all generators c > b,
	**                                  Exponent[b] == n > 0.
	*/
	for( b = 1; !EarlyStop && b <= NrPcGens; b++ )
	    if( Exponent[b] > 0 )
		for( c = b+1; !EarlyStop && c <= NrPcGens; c++ )
		    do_cbn( c, b );

	/*
	**    c^n * b = c^(n-1) * ( c * b ) for all generators c > b,
	**                                  Exponent[c] == n > 0.
	*/
	for( c = 1; !EarlyStop && c <= NrPcGens; c++ )
	    if( Exponent[c] > 0 ) {
		for( b = 1; !EarlyStop && b <= min(c-1,Dimension[1]); b++ )
		    do_cnb( c, b );
		for( b = 1; !EarlyStop && b < c; b++ )
		    if( Exponent[b] == 0 ) do_cnb( c, -b );
	    }

	/*
	**    c^n * c = c * c^n for all generators c, Exponent[c] == n > 0.
	*/
	for( c = 1; !EarlyStop && c <= NrPcGens; c++ )
	    if( Exponent[c] > 0 )
		do_cnc( c );

	if( Verbose )
	    printf("#    Checked consistency (%d msec).\n",RunTime()-t);
}