bool Multivector::AssignInnerProduct( const Multivector& multivectorA, const Multivector& multivectorB )
{
if( multivectorA.CountProductTypes( Term::GEOMETRIC_PRODUCT ) > 0 )
if( !const_cast< Multivector* >( &multivectorA )->CollectTerms( Term::OUTER_PRODUCT ) )
return false;
if( multivectorB.CountProductTypes( Term::GEOMETRIC_PRODUCT ) > 0 )
if( !const_cast< Multivector* >( &multivectorB )->CollectTerms( Term::OUTER_PRODUCT ) )
return false;
Multivector* multivectorResult = this;
Multivector multivectorStorage;
if( this == &multivectorA || this == &multivectorB )
multivectorResult = &multivectorStorage;
for( const SumOfTerms::Node* nodeA = multivectorA.sumOfTerms.Head(); nodeA; nodeA = nodeA->Next() )
{
const Term* termA = nodeA->data;
for( const SumOfTerms::Node* nodeB = multivectorB.sumOfTerms.Head(); nodeB; nodeB = nodeB->Next() )
{
const Term* termB = nodeB->data;
Multivector innerProductMultivector;
if( !termA->InnerProductMultiply( termB, innerProductMultivector ) )
return false;
Scalar scalar;
if( !scalar.AssignProduct( *termA->coeficient, *termB->coeficient ) )
return false;
if( !innerProductMultivector.AssignScalarProduct( scalar, innerProductMultivector ) )
return false;
multivectorResult->sumOfTerms.Absorb( &innerProductMultivector.sumOfTerms );
}
}
if( multivectorResult == &multivectorStorage )
return Assign( multivectorStorage );
return true;
}
