crypto/weakcrypto/test/tbigint/tbasicmathsfb.cpp
author hgs
Thu, 24 Jun 2010 15:39:07 +0530
changeset 72 de46a57f75fb
permissions -rw-r--r--
201023_02
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
72
de46a57f75fb 201023_02
hgs
parents:
diff changeset
     1
/*
de46a57f75fb 201023_02
hgs
parents:
diff changeset
     2
* Copyright (c) 2002-2009 Nokia Corporation and/or its subsidiary(-ies).
de46a57f75fb 201023_02
hgs
parents:
diff changeset
     3
* All rights reserved.
de46a57f75fb 201023_02
hgs
parents:
diff changeset
     4
* This component and the accompanying materials are made available
de46a57f75fb 201023_02
hgs
parents:
diff changeset
     5
* under the terms of the License "Eclipse Public License v1.0"
de46a57f75fb 201023_02
hgs
parents:
diff changeset
     6
* which accompanies this distribution, and is available
de46a57f75fb 201023_02
hgs
parents:
diff changeset
     7
* at the URL "http://www.eclipse.org/legal/epl-v10.html".
de46a57f75fb 201023_02
hgs
parents:
diff changeset
     8
*
de46a57f75fb 201023_02
hgs
parents:
diff changeset
     9
* Initial Contributors:
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    10
* Nokia Corporation - initial contribution.
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    11
*
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    12
* Contributors:
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    13
*
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    14
* Description: 
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    15
*
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    16
*/
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    17
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    18
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    19
#include "tbasicmathsfb.h"
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    20
#include "t_input.h"
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    21
#include "t_output.h"
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    22
#include <bigint.h>
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    23
#include <random.h>
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    24
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    25
CTestAction* CBasicMathsFB::NewL(RFs& aFs, CConsoleBase& aConsole, 
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    26
	Output& aOut, const TTestActionSpec& aTestActionSpec)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    27
	{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    28
	CTestAction* self = CBasicMathsFB::NewLC(aFs, aConsole,
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    29
		aOut, aTestActionSpec);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    30
	CleanupStack::Pop();
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    31
	return self;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    32
	}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    33
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    34
CTestAction* CBasicMathsFB::NewLC(RFs& aFs, CConsoleBase& aConsole, 
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    35
	Output& aOut, const TTestActionSpec& aTestActionSpec)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    36
	{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    37
	CBasicMathsFB* self = new(ELeave) CBasicMathsFB(aFs, aConsole, aOut);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    38
	CleanupStack::PushL(self);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    39
	self->ConstructL(aTestActionSpec);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    40
	return self;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    41
	}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    42
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    43
CBasicMathsFB::~CBasicMathsFB()
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    44
	{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    45
	delete iBody;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    46
	}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    47
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    48
CBasicMathsFB::CBasicMathsFB(RFs& aFs, CConsoleBase& aConsole, Output& aOut)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    49
	: CTestAction(aConsole, aOut), iFs(aFs)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    50
	{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    51
	}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    52
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    53
void CBasicMathsFB::ConstructL(const TTestActionSpec& aTestActionSpec)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    54
	{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    55
	CTestAction::ConstructL(aTestActionSpec);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    56
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    57
	iBody = HBufC8::NewL(aTestActionSpec.iActionBody.Length());
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    58
	iBody->Des().Copy(aTestActionSpec.iActionBody);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    59
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    60
	//HBufC8* length = Input::ParseElementHexL(*iBody, _L8("<bits>"));
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    61
	TUint bits = Input::ParseIntElement(*iBody, _L8("<bits>"), _L8("</bits>"));
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    62
	// the final /7 gives the number of times we have to increment by 7 to get
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    63
	// to that number of bytes and hence bits.
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    64
	iIterations = ((bits+7)/8)/7 + 1;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    65
	}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    66
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    67
void CBasicMathsFB::DoPerformPrerequisite(TRequestStatus& aStatus)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    68
	{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    69
	TRequestStatus* status = &aStatus;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    70
	User::RequestComplete(status, KErrNone);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    71
	iActionState = CTestAction::EAction;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    72
	}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    73
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    74
void CBasicMathsFB::DoPerformPostrequisite(TRequestStatus& aStatus)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    75
	{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    76
	TRequestStatus* status = &aStatus;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    77
	iFinished = ETrue;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    78
	User::RequestComplete(status, KErrNone);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    79
	}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    80
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    81
void CBasicMathsFB::DoReportAction(void)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    82
	{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    83
	}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    84
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    85
void CBasicMathsFB::DoCheckResult(TInt)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    86
	{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    87
	}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    88
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    89
void CBasicMathsFB::PerformAction(TRequestStatus& aStatus)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    90
	{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    91
	__UHEAP_MARK;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    92
	TRequestStatus* status = &aStatus;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    93
	iResult = ETrue;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    94
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    95
	//min max values for NewRandomLC call
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    96
	RInteger min = RInteger::NewL(10);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    97
	CleanupStack::PushL(min);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    98
	RInteger max = RInteger::NewL(100);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
    99
	CleanupStack::PushL(max);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   100
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   101
	//Generate iIterations*7 byte random sequences we are using 7 as it's a generator
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   102
	//mod 8.  Thus we'll cycle through every value (0-7) every 8 iterations.
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   103
	//This gives us a better feeling that certain byte lengths (and thus bit
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   104
	//lengths as the byte is chosen randomly) don't have errors.
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   105
	for(TUint i=1; i<iIterations; i++)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   106
		{ 
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   107
		HBufC8* buf = HBufC8::NewMaxLC(i*7);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   108
		TPtr8 ptr = buf->Des();
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   109
		TRandom::RandomL(ptr);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   110
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   111
		//This is this iteration's random number
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   112
		RInteger initial = RInteger::NewL(ptr);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   113
		CleanupStack::PushL(initial);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   114
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   115
		//get a number x | 10 < x < 100
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   116
		RInteger crange = RInteger::NewRandomL(min, max);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   117
		CleanupStack::PushL(crange);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   118
		TUint range = crange.ConvertToLongL();
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   119
		CleanupStack::PopAndDestroy(); //crange
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   120
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   121
		AddSub(initial, range);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   122
		MulDiv(initial, range);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   123
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   124
		//GCD
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   125
		CleanupStack::PopAndDestroy(); //initial
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   126
		CleanupStack::PopAndDestroy();//buf
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   127
		iConsole.Printf(_L("."));
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   128
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   129
	
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   130
	//Test a single iteration where the initial random number is less than a
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   131
	//word so the division and modulo routines that take words rather than
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   132
	//TIntegers can run.
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   133
	//do
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   134
		{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   135
		//This is this iteration's random number
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   136
		RInteger initial = RInteger::NewRandomL(31);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   137
		CleanupStack::PushL(initial);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   138
		//get a number x | 10 < x < 100
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   139
		RInteger crange = RInteger::NewRandomL(min, max);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   140
		CleanupStack::PushL(crange);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   141
		TUint range = crange.ConvertToLongL();
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   142
		CleanupStack::PopAndDestroy(&crange); //crange
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   143
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   144
		AddSub(initial, range);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   145
		MulDiv(initial, range);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   146
		CleanupStack::PopAndDestroy(&initial); //initial
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   147
		iConsole.Printf(_L("."));
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   148
		} //while (0);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   149
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   150
	CleanupStack::PopAndDestroy();//max
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   151
	CleanupStack::PopAndDestroy(); //min
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   152
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   153
	User::RequestComplete(status, KErrNone);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   154
	iActionState = CTestAction::EPostrequisite;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   155
	__UHEAP_MARK;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   156
	}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   157
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   158
void CBasicMathsFB::AddSub(const TInteger& aInitial, TUint aRange)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   159
	{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   160
	__UHEAP_MARK;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   161
	//This is the copy we are going to do stuff to
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   162
	RInteger a = RInteger::NewL(aInitial);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   163
	CleanupStack::PushL(a);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   164
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   165
	// compute a*aRange using doubling
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   166
	TUint j=1;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   167
	for(; j<aRange; j++)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   168
		{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   169
		a += aInitial;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   170
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   171
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   172
	//b = a*aRange;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   173
	RInteger b = RInteger::NewL(a);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   174
	CleanupStack::PushL(b);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   175
	//compute (a*aRange)/aRange using subtraction
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   176
	for(j=1; j<aRange; j++)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   177
		{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   178
		b -= aInitial;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   179
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   180
	// b should be the same as the initial value
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   181
	if( b != aInitial )
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   182
		{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   183
		iResult = EFalse;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   184
		iOut.writeString(_L("AddSub Failure:"));
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   185
		iOut.writeNewLine();
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   186
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   187
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   188
	RInteger c = RInteger::NewL(aInitial);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   189
	CleanupStack::PushL(c);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   190
	// compute a*aRange using normal multiplication
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   191
	c *= aRange;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   192
	
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   193
	// c and a should now be the same
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   194
	if( c != a )
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   195
		{ 
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   196
		iResult = EFalse;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   197
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   198
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   199
	RInteger d = RInteger::NewL(a);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   200
	CleanupStack::PushL(d);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   201
	//compute (a*aRange)/aRange using normal division
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   202
	d /= aRange;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   203
	if( d != aInitial )
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   204
		{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   205
		iResult = EFalse;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   206
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   207
	RInteger e = RInteger::NewL(a);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   208
	CleanupStack::PushL(e);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   209
	e %= aRange;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   210
	// (a*aRange)%aRange == 0
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   211
	if( e != 0 )
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   212
		{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   213
		iResult = EFalse;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   214
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   215
	CleanupStack::PopAndDestroy(5); //e,d,c,b,a
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   216
	__UHEAP_MARKEND;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   217
	}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   218
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   219
void CBasicMathsFB::MulDiv(const TInteger& aInitial, TUint aRange)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   220
	{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   221
	__UHEAP_MARK;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   222
	//This is the copy we are going to do stuff to
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   223
	RInteger a = RInteger::NewL(aInitial);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   224
	CleanupStack::PushL(a);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   225
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   226
	//compute a = aInitial^aRange using repeated multiplication
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   227
	TUint j=1;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   228
	for(; j<aRange; j++)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   229
		{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   230
		a *= aInitial;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   231
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   232
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   233
	//b = a
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   234
	RInteger b = RInteger::NewL(a);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   235
	CleanupStack::PushL(b);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   236
	//try to find aInitial by repeatedly dividing b by aInitial aRange times
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   237
	for(j=1; j<aRange; j++)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   238
		{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   239
		TRAPD(res, b /= aInitial);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   240
		//the first time through aInitial is 0 so this is expected
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   241
		if(res == KErrDivideByZero && aInitial.IsZero())
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   242
			{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   243
			break;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   244
			}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   245
		else if(res == KErrDivideByZero && aInitial.NotZero())
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   246
			{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   247
			iResult = EFalse;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   248
			}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   249
		else if(res != KErrNone)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   250
			{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   251
			User::Leave(res);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   252
			}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   253
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   254
	// b should be the same as the initial value
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   255
	if( b != aInitial )
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   256
		{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   257
		iResult = EFalse;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   258
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   259
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   260
	//tests division by something smaller than a word
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   261
	if(aInitial.WordCount() <= 1)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   262
		{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   263
		RInteger dividend = RInteger::NewL(a);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   264
		CleanupStack::PushL(dividend);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   265
		for(j=1; j<aRange; j++)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   266
			{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   267
			RInteger quotient;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   268
			//try to find aInitial by repeatedly dividing dividend by aInitial aRange times 
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   269
			TRAPD(res, quotient = dividend.DividedByL(aInitial.ConvertToLongL()));
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   270
			//the first time through aInitial is 0 so this is expected
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   271
			if(res == KErrDivideByZero && aInitial.IsZero())
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   272
				{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   273
				break;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   274
				}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   275
			else if(res == KErrDivideByZero && aInitial.NotZero())
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   276
				{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   277
				iResult = EFalse;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   278
				}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   279
			else if(res != KErrNone)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   280
				{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   281
				User::Leave(res);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   282
				}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   283
			dividend.Set(quotient);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   284
			}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   285
		if( dividend != aInitial )
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   286
			{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   287
			iResult = EFalse;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   288
			}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   289
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   290
		TUint remainder=1;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   291
		TRAPD(res, remainder = a.ModuloL(aInitial.ConvertToLongL()));
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   292
		//the first time through aInitial is 0
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   293
		if(res != KErrDivideByZero && res != KErrNone)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   294
			{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   295
			User::Leave(res);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   296
			}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   297
		else if(res == KErrDivideByZero && aInitial.NotZero())
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   298
			{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   299
			iResult = EFalse;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   300
			}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   301
		//else we have an expected divide by zero, ignore it.
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   302
		if(remainder != 0)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   303
			{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   304
			iResult = EFalse;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   305
			}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   306
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   307
		CleanupStack::PopAndDestroy(&dividend);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   308
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   309
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   310
	RInteger c = RInteger::NewL(aRange);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   311
	CleanupStack::PushL(c);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   312
	RInteger d = aInitial.ExponentiateL(c);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   313
	CleanupStack::PushL(d);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   314
	// c and a should now be the same
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   315
	if( d != a )
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   316
		{ 
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   317
		iResult = EFalse;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   318
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   319
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   320
	RInteger e = RInteger::NewL(a);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   321
	CleanupStack::PushL(e);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   322
	TRAPD(res, e %= aInitial);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   323
	//the first time through aInitial is 0
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   324
	if(res != KErrDivideByZero && res != KErrNone)
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   325
		{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   326
		User::Leave(res);
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   327
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   328
	else if(res == KErrDivideByZero && aInitial.NotZero())
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   329
		{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   330
		iResult = EFalse;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   331
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   332
	//else we have an expected divide by zero, ignore it.
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   333
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   334
	// (aInitial^aRange)%aInitial == 0
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   335
	if( e != 0 )
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   336
		{
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   337
		iResult = EFalse;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   338
		}
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   339
	CleanupStack::PopAndDestroy(5);//e,d,c,b,a
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   340
	__UHEAP_MARKEND;
de46a57f75fb 201023_02
hgs
parents:
diff changeset
   341
	}