diff options
author | Gabe Black <gabeblack@google.com> | 2018-05-24 01:37:55 -0700 |
---|---|---|
committer | Gabe Black <gabeblack@google.com> | 2018-08-08 10:09:54 +0000 |
commit | 16fa8d7cc8c92f5ab879e4cf9c6c0bbb3567860f (patch) | |
tree | 7b6faaacb4574a555e561534aa4a8508c0624c32 /src/systemc/tests/systemc/misc/v1.0/blv | |
parent | 7235d3b5211d0ba8f528d930a4c1e7ad62eec51a (diff) | |
download | gem5-16fa8d7cc8c92f5ab879e4cf9c6c0bbb3567860f.tar.xz |
systemc: Import tests from the Accellera systemc distribution.
Change-Id: Iad76b398949a55d768a34d027a2d8e3739953da6
Reviewed-on: https://gem5-review.googlesource.com/10845
Reviewed-by: Giacomo Travaglini <giacomo.travaglini@arm.com>
Maintainer: Gabe Black <gabeblack@google.com>
Diffstat (limited to 'src/systemc/tests/systemc/misc/v1.0/blv')
-rw-r--r-- | src/systemc/tests/systemc/misc/v1.0/blv/blv.cpp | 148 | ||||
-rw-r--r-- | src/systemc/tests/systemc/misc/v1.0/blv/golden/blv.log | 5 | ||||
-rw-r--r-- | src/systemc/tests/systemc/misc/v1.0/blv/isaac.h | 272 |
3 files changed, 425 insertions, 0 deletions
diff --git a/src/systemc/tests/systemc/misc/v1.0/blv/blv.cpp b/src/systemc/tests/systemc/misc/v1.0/blv/blv.cpp new file mode 100644 index 000000000..e2135fbf9 --- /dev/null +++ b/src/systemc/tests/systemc/misc/v1.0/blv/blv.cpp @@ -0,0 +1,148 @@ +/***************************************************************************** + + Licensed to Accellera Systems Initiative Inc. (Accellera) under one or + more contributor license agreements. See the NOTICE file distributed + with this work for additional information regarding copyright ownership. + Accellera licenses this file to you under the Apache License, Version 2.0 + (the "License"); you may not use this file except in compliance with the + License. You may obtain a copy of the License at + + http://www.apache.org/licenses/LICENSE-2.0 + + Unless required by applicable law or agreed to in writing, software + distributed under the License is distributed on an "AS IS" BASIS, + WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or + implied. See the License for the specific language governing + permissions and limitations under the License. + + *****************************************************************************/ + +/***************************************************************************** + + blv.cpp -- + + Original Author: Martin Janssen, Synopsys, Inc., 2002-02-15 + + *****************************************************************************/ + +/***************************************************************************** + + MODIFICATION LOG - modifiers, enter your name, affiliation, date and + changes you are making here. + + Name, Affiliation, Date: + Description of Modification: + + *****************************************************************************/ + +//---------------------------------------------------------- +// test the basic functionality of sc_bv<> and sc_lv<> +//---------------------------------------------------------- +#include "systemc.h" +#include <time.h> +#include "isaac.h" + +QTIsaac<8> rng; + +template<class T> +void compare(int W) +{ // functionality verification + sc_bv_base x(W); + T st; + // initialize + for(int i=0; i<W; i++) + { + bool la = (rng.rand()&1) == 0; + x[i] = la; + st[i] = la; + } + if(st.to_string()!=x.to_string()) + cout<<"\nERROR: x= "<<x<<" st= "<<st<<"\n"<< + "st="<<st<<"\nx="<<x<< endl; + if( (int) st.xor_reduce()!= x.xor_reduce()) + cout<<"\nERROR: st.xor_reduce="<<(short)st.xor_reduce()<<"; x.xor_reduce="<< + (short)x.xor_reduce()<<"\n"<< + "st="<<st<<"\nx="<<x<< endl; + if( (int) x.or_reduce()!=st.or_reduce()) + cout<<"\nERROR: st.or_reduce="<<(short)st.or_reduce()<< + "; x.or_reduce="<<(short)x.or_reduce()<<"\n"<< + "st="<<st<<"\nx="<<x<< endl; + if( (int) x.and_reduce()!=st.and_reduce()) + cout<<"\nERROR: st.and_reduce="<<(short)st.and_reduce()<< + "; x.and_reduce="<<(short)x.and_reduce()<<"\n"<< + "st="<<st<<"\nx="<<x<< endl; + + if((st,st).to_string()!=(x,x).to_string()) + cout<<"\nERROR: st,st="<<(st,st)<< + "; x,x="<<(x,x)<<"\n"<< + "st="<<st<<"\nx="<<x<< endl; + int first = (int) ( (double) W * ((double)rng.rand() / (double)0x7fffffff)); + int second = (int) ( (double) W * ((double)rng.rand() / (double)0x7fffffff)); + if(st.range(first,second).to_string()!=x.range(first,second).to_string()) + cout<<"st.range("<<first<<","<<second<<")="<<st.range(first,second)<< + "; x.range("<<first<<","<<second<<")="<<x.range(first,second)<<"\n"<< + "st="<<st<<"\nx="<<x<< endl; + + sc_bv_base bv(2*W); + sc_bv_base xv(2*W); + bv = (st,st); + xv = (x,x); + if(bv.to_string()!=xv.to_string()) + cout<<"\nERROR: bv(st,st)="<<bv<<"; xv(x,x)="<<xv<<"\n"<< + "st="<<st<<"\nx="<<x<< endl; + int Len=0; + if(first>second) + Len = first-second; + else + Len = second-first; + sc_bv_base br(Len+1); + sc_bv_base xr(Len+1); + br = st.range(first,second); + xr = x.range(first,second); + if(br.to_string()!=xr.to_string()) + { + cout<<"\nERROR: br("<<first<<","<<second<<")!= xr("<<first<<","<<second + <<")" << endl; + br = st.range(first,second); + xr = x.range(first,second); + cout<<"br="<<br.to_string()<<" xr="<<xr.to_string()<< + "st.range="<<st.range(first,second)<<" x.range="<<x.range(first,second)<< + "st="<<st<<"\nx="<<x<< endl; + } + + // verify assignments + long ra = rng.rand(); + x = ra; + st = ra; + if(st.to_string()!=x.to_string()) + cout<<"\nERROR (assignment): x= "<<x<<" st= "<<st<<" original long="<<ra + << endl; +} + +int sc_main(int, char**) +{ + const int N = 2000; + int Seed = rng.rand(); + cout<<"\nverifying sc_bv<"<<N<<">" << endl; + try{ + for(int i=0; i<1000; i++) + compare<sc_bv<N> >(N); + } + catch(...) + { + cout<<"Test failed due to exception in sc_bv. Seed = "<<Seed<< endl; + throw; + } + try{ + cout<<"\nverifying sc_lv<"<<N<<">" << endl; + for(int i=0; i<1000; i++) + compare<sc_lv<N> >(N); + } + catch(...) + { + cout<<"Test failed due to exception in sc_lv. Seed = "<<Seed<< endl; + throw; + } + + return 0; +} diff --git a/src/systemc/tests/systemc/misc/v1.0/blv/golden/blv.log b/src/systemc/tests/systemc/misc/v1.0/blv/golden/blv.log new file mode 100644 index 000000000..12908a228 --- /dev/null +++ b/src/systemc/tests/systemc/misc/v1.0/blv/golden/blv.log @@ -0,0 +1,5 @@ +SystemC Simulation + +verifying sc_bv<2000> + +verifying sc_lv<2000> diff --git a/src/systemc/tests/systemc/misc/v1.0/blv/isaac.h b/src/systemc/tests/systemc/misc/v1.0/blv/isaac.h new file mode 100644 index 000000000..9625be2ac --- /dev/null +++ b/src/systemc/tests/systemc/misc/v1.0/blv/isaac.h @@ -0,0 +1,272 @@ +#ifndef __ISAAC_HPP +#define __ISAAC_HPP + + +/* + + C++ TEMPLATE VERSION OF Robert J. Jenkins Jr.'s + ISAAC Random Number Generator. + + Ported from vanilla C to to template C++ class + by Quinn Tyler Jackson on 16-23 July 1998. + + quinn@qtj.net + + The function for the expected period of this + random number generator, according to Jenkins is: + + f(a,b) = 2**((a+b*(3+2^^a)-1) + + (where a is ALPHA and b is bitwidth) + + So, for a bitwidth of 32 and an ALPHA of 8, + the expected period of ISAAC is: + + 2^^(8+32*(3+2^^8)-1) = 2^^8295 + + Jackson has been able to run implementations + with an ALPHA as high as 16, or + + 2^^2097263 + +*/ + + +typedef unsigned int UINT32; +const UINT32 GOLDEN_RATIO = UINT32(0x9e3779b9); + + +template <UINT32 ALPHA = (8)> +class QTIsaac +{ + public: + + typedef unsigned char byte; + + struct randctx + { + randctx(void) + { + randrsl = new UINT32[N]; + randmem = new UINT32[N]; + } + + ~randctx(void) + { + delete [] randrsl; + delete [] randmem; + } + + UINT32 randcnt; + UINT32* randrsl; + UINT32* randmem; + UINT32 randa; + UINT32 randb; + UINT32 randc; + }; + + QTIsaac(UINT32 a = 0, UINT32 b = 0, UINT32 c = 0); + virtual ~QTIsaac(void); + + UINT32 rand(void); + virtual void randinit(randctx* ctx, bool bUseSeed); + virtual void srand( + UINT32 a = 0, UINT32 b = 0, UINT32 c = 0, UINT32* s = NULL); + + enum {N = (1<<ALPHA)}; + + protected: + + virtual void isaac(randctx* ctx); + + UINT32 ind(UINT32* mm, UINT32 x); + void rngstep( + UINT32 mix, UINT32& a, UINT32& b, UINT32*& mm, UINT32*& m, + UINT32*& m2, UINT32*& r, UINT32& x, UINT32& y); + virtual void shuffle( + UINT32& a, UINT32& b, UINT32& c, UINT32& d, UINT32& e, UINT32& f, + UINT32& g, UINT32& h); + + private: + randctx m_rc; +}; + + +template<UINT32 ALPHA> +QTIsaac<ALPHA>::QTIsaac(UINT32 a, UINT32 b, UINT32 c) : m_rc() +{ + srand(a, b, c); +} + + +template<UINT32 ALPHA> +QTIsaac<ALPHA>::~QTIsaac(void) +{ + // DO NOTHING +} + + +template<UINT32 ALPHA> +void QTIsaac<ALPHA>::srand(UINT32 a, UINT32 b, UINT32 c, UINT32* s) +{ + for(int i = 0; i < N; i++) + { + m_rc.randrsl[i] = s != NULL ? s[i] : 0; + } + + m_rc.randa = a; + m_rc.randb = b; + m_rc.randc = c; + + randinit(&m_rc, true); +} + + +template<UINT32 ALPHA> +inline UINT32 QTIsaac<ALPHA>::rand(void) +{ + return 0x7fffffff & (!m_rc.randcnt-- ? + (isaac(&m_rc), m_rc.randcnt=(N-1), m_rc.randrsl[m_rc.randcnt]) : + m_rc.randrsl[m_rc.randcnt]); +} + + +template<UINT32 ALPHA> +inline void QTIsaac<ALPHA>::randinit(randctx* ctx, bool bUseSeed) +{ + UINT32 a,b,c,d,e,f,g,h; + int i; + + a = b = c = d = e = f = g = h = GOLDEN_RATIO; + + UINT32* m = (ctx->randmem); + UINT32* r = (ctx->randrsl); + + if(!bUseSeed) + { + ctx->randa = 0; + ctx->randb = 0; + ctx->randc = 0; + } + + // scramble it + for(i=0; i < 4; ++i) + { + shuffle(a,b,c,d,e,f,g,h); + } + + if(bUseSeed) + { + // initialize using the contents of r[] as the seed + + for(i=0; i < N; i+=8) + { + a+=r[i ]; b+=r[i+1]; c+=r[i+2]; d+=r[i+3]; + e+=r[i+4]; f+=r[i+5]; g+=r[i+6]; h+=r[i+7]; + + shuffle(a,b,c,d,e,f,g,h); + + m[i ]=a; m[i+1]=b; m[i+2]=c; m[i+3]=d; + m[i+4]=e; m[i+5]=f; m[i+6]=g; m[i+7]=h; + } + + //do a second pass to make all of the seed affect all of m + + for(i=0; i < N; i += 8) + { + a+=m[i ]; b+=m[i+1]; c+=m[i+2]; d+=m[i+3]; + e+=m[i+4]; f+=m[i+5]; g+=m[i+6]; h+=m[i+7]; + + shuffle(a,b,c,d,e,f,g,h); + + m[i ]=a; m[i+1]=b; m[i+2]=c; m[i+3]=d; + m[i+4]=e; m[i+5]=f; m[i+6]=g; m[i+7]=h; + } + } + else + { + // fill in mm[] with messy stuff + + shuffle(a,b,c,d,e,f,g,h); + + m[i ]=a; m[i+1]=b; m[i+2]=c; m[i+3]=d; + m[i+4]=e; m[i+5]=f; m[i+6]=g; m[i+7]=h; + + } + + isaac(ctx); // fill in the first set of results + ctx->randcnt = N; // prepare to use the first set of results +} + + +template<UINT32 ALPHA> +inline UINT32 QTIsaac<ALPHA>::ind(UINT32* mm, UINT32 x) +{ + return (*(UINT32*)((byte*)(mm) + ((x) & ((N-1)<<2)))); +} + + +template<UINT32 ALPHA> +inline void QTIsaac<ALPHA>::rngstep(UINT32 mix, UINT32& a, UINT32& b, UINT32*& mm, UINT32*& m, UINT32*& m2, UINT32*& r, UINT32& x, UINT32& y) +{ + x = *m; + a = (a^(mix)) + *(m2++); + *(m++) = y = ind(mm,x) + a + b; + *(r++) = b = ind(mm,y>>ALPHA) + x; +} + + +template<UINT32 ALPHA> +inline void QTIsaac<ALPHA>::shuffle(UINT32& a, UINT32& b, UINT32& c, UINT32& d, UINT32& e, UINT32& f, UINT32& g, UINT32& h) +{ + a^=b<<11; d+=a; b+=c; + b^=c>>2; e+=b; c+=d; + c^=d<<8; f+=c; d+=e; + d^=e>>16; g+=d; e+=f; + e^=f<<10; h+=e; f+=g; + f^=g>>4; a+=f; g+=h; + g^=h<<8; b+=g; h+=a; + h^=a>>9; c+=h; a+=b; +} + + +template<UINT32 ALPHA> +inline void QTIsaac<ALPHA>::isaac(randctx* ctx) +{ + UINT32 x,y; + + UINT32* mm = ctx->randmem; + UINT32* r = ctx->randrsl; + + UINT32 a = (ctx->randa); + UINT32 b = (ctx->randb + (++ctx->randc)); + + UINT32* m = mm; + UINT32* m2 = (m+(N/2)); + UINT32* mend = m2; + + for(; m<mend; ) + { + rngstep((a<<13), a, b, mm, m, m2, r, x, y); + rngstep((a>>6) , a, b, mm, m, m2, r, x, y); + rngstep((a<<2) , a, b, mm, m, m2, r, x, y); + rngstep((a>>16), a, b, mm, m, m2, r, x, y); + } + + m2 = mm; + + for(; m2<mend; ) + { + rngstep((a<<13), a, b, mm, m, m2, r, x, y); + rngstep((a>>6) , a, b, mm, m, m2, r, x, y); + rngstep((a<<2) , a, b, mm, m, m2, r, x, y); + rngstep((a>>16), a, b, mm, m, m2, r, x, y); + } + + ctx->randb = b; + ctx->randa = a; +} + + +#endif // __ISAAC_HPP + |