/*
* Checkup routine
*/
int md2_self_test( int verbose )
{
int i;
unsigned char md2sum[16];
for( i = 0; i < 7; i++ )
{
if( verbose != 0 )
printf( " MD2 test #%d: ", i + 1 );
md2( (unsigned char *) md2_test_str[i],
strlen( md2_test_str[i] ), md2sum );
if( memcmp( md2sum, md2_test_sum[i], 16 ) != 0 )
{
if( verbose != 0 )
printf( "failed\n" );
return( 1 );
}
if( verbose != 0 )
printf( "passed\n" );
}
if( verbose != 0 )
printf( "\n" );
return( 0 );
}
/*
* Wrapper for x509 hashes.
*/
static void x509_hash( const unsigned char *in, size_t len, int alg,
unsigned char *out )
{
switch( alg )
{
#if defined(POLARSSL_MD2_C)
case SIG_RSA_MD2 : md2( in, len, out ); break;
#endif
#if defined(POLARSSL_MD4_C)
case SIG_RSA_MD4 : md4( in, len, out ); break;
#endif
#if defined(POLARSSL_MD5_C)
case SIG_RSA_MD5 : md5( in, len, out ); break;
#endif
#if defined(POLARSSL_SHA1_C)
case SIG_RSA_SHA1 : sha1( in, len, out ); break;
#endif
#if defined(POLARSSL_SHA2_C)
case SIG_RSA_SHA224 : sha2( in, len, out, 1 ); break;
case SIG_RSA_SHA256 : sha2( in, len, out, 0 ); break;
#endif
#if defined(POLARSSL_SHA4_C)
case SIG_RSA_SHA384 : sha4( in, len, out, 1 ); break;
case SIG_RSA_SHA512 : sha4( in, len, out, 0 ); break;
#endif
default:
memset( out, '\xFF', 64 );
break;
}
}
/*
* MD2 HMAC context setup
*/
void md2_hmac_starts( md2_context *ctx, unsigned char *key, int keylen )
{
int i;
unsigned char sum[16];
if( keylen > 64 )
{
md2( key, keylen, sum );
keylen = 16;
key = sum;
}
memset( ctx->ipad, 0x36, 64 );
memset( ctx->opad, 0x5C, 64 );
for( i = 0; i < keylen; i++ )
{
ctx->ipad[i] = (unsigned char)( ctx->ipad[i] ^ key[i] );
ctx->opad[i] = (unsigned char)( ctx->opad[i] ^ key[i] );
}
md2_starts( ctx );
md2_update( ctx, ctx->ipad, 64 );
memset( sum, 0, sizeof( sum ) );
}
/*
* MD2 HMAC context setup
*/
void md2_hmac_starts( md2_context *ctx, const unsigned char *key,
size_t keylen )
{
size_t i;
unsigned char sum[16];
if( keylen > 16 )
{
md2( key, keylen, sum );
keylen = 16;
key = sum;
}
memset( ctx->ipad, 0x36, 16 );
memset( ctx->opad, 0x5C, 16 );
for( i = 0; i < keylen; i++ )
{
ctx->ipad[i] = (unsigned char)( ctx->ipad[i] ^ key[i] );
ctx->opad[i] = (unsigned char)( ctx->opad[i] ^ key[i] );
}
md2_starts( ctx );
md2_update( ctx, ctx->ipad, 16 );
polarssl_zeroize( sum, sizeof( sum ) );
}
void *
carmen_mdalloc(int ndim, int width, ...)
{
va_list argp;
int *dims, i;
char ***tip;
va_start(argp, width);
/* allocate storage for variable args (dimensions) */
dims = malloc(ndim*sizeof(int));
/* check_alloc checked */
if(dims == NULL)
return NULL;
/* initialize dimensions array for subsequent calls */
for(i=0; i<ndim; i++)
dims[i] = va_arg(argp,int);
w_units = width; /* global used by md2 and md3 */
/* allocate required pointer and array element storage */
tip = (char ***)md2(dims[0], ndim, &dims[1]);
if(ndim>1 && tip)
md3(tip, dims[0], ndim-1, &dims[1]); /* init pointers */
free(dims);
return tip;
}
/* ****************************************************************************
*
* check_json -
*/
TEST(AppendContextElementRequest, check_json)
{
AppendContextElementRequest acer;
std::string out;
ContextAttribute ca("caName", "caType", "121");
Metadata md("mdName", "mdType", "122");
const char* outfile1 = "ngsi10.appendContextElementResponse.predetectedError.valid.json";
const char* outfile2 = "ngsi10.appendContextElementResponse.missingAttributeName.valid.json";
const char* outfile3 = "ngsi10.appendContextElementResponse.missingMetadataName.valid.json";
ConnectionInfo ci;
utInit();
acer.attributeDomainName.set("ADN");
acer.contextAttributeVector.push_back(&ca);
acer.domainMetadataVector.push_back(&md);
// 1. ok
ci.outMimeType = JSON;
out = acer.check(&ci, AppendContextElement, "", "", 0);
EXPECT_STREQ("OK", out.c_str());
// 2. Predetected error
EXPECT_EQ("OK", testDataFromFile(expectedBuf, sizeof(expectedBuf), outfile1)) << "Error getting test data from '" << outfile1 << "'";
out = acer.check(&ci, AppendContextElement, "", "Error is predetected", 0);
EXPECT_STREQ(expectedBuf, out.c_str());
// 3. bad ContextAttribute
ContextAttribute ca2("", "caType", "121");
acer.contextAttributeVector.push_back(&ca2);
out = acer.check(&ci, AppendContextElement, "", "", 0);
EXPECT_EQ("OK", testDataFromFile(expectedBuf, sizeof(expectedBuf), outfile2)) << "Error getting test data from '" << outfile2 << "'";
EXPECT_STREQ(expectedBuf, out.c_str());
ca2.name = "ca2Name";
// 4. Bad domainMetadata
Metadata md2("", "mdType", "122");
acer.domainMetadataVector.push_back(&md2);
out = acer.check(&ci, AppendContextElement, "", "", 0);
EXPECT_EQ("OK", testDataFromFile(expectedBuf, sizeof(expectedBuf), outfile3)) << "Error getting test data from '" << outfile3 << "'";
EXPECT_STREQ(expectedBuf, out.c_str());
// 5. Bad attributeDomainName
// FIXME P3: AttributeDomainName::check always returns "OK"
utExit();
}
void testMediaDescriptor_TypeAndDirectionality()
{
const char *sdp =
"v=0\r\n"
"o=mhandley 2890844526 2890842807 IN IP4 126.16.64.4\r\n"
"s=SDP Seminar\r\n"
"i=A Seminar on the session description protocol\r\n"
"u=http://www.cs.ucl.ac.uk/staff/M.Handley/sdp.03.ps\r\n"
"[email protected] (Mark Handley)\r\n"
"c=IN IP4 224.2.17.12/127\r\n"
"t=2873397496 2873404696\r\n"
"m=audio 49170 RTP/AVP 0\r\n" // media description 0
"c=IN IP4 224.2.17.12/127\r\n"
"a=recvonly\r\n"
"m=video 51372 RTP/AVP 31\r\n" // media description 1
"a=inactive\r\n"
"m=application 32416 udp wb\r\n" // media description 2
"a=sendonly\r\n"
"m=audio 55554 RTP/AVP 0\r\n" // media description 3
"c=IN IP4 224.2.17.12/127\r\n"
"a=sendrecv\r\n"
"m=audio 55560 RTP/AVP 0\r\n" // media description 4
;
SdpBody body(sdp);
MediaDescriptor md0( body, 0, CALLER );
MediaDescriptor md1( body, 1, CALLER );
MediaDescriptor md2( body, 2, CALLER );
MediaDescriptor md3( body, 3, CALLER );
MediaDescriptor md4( body, 4, CALLER );
CPPUNIT_ASSERT( md0.getType() == "audio" );
CPPUNIT_ASSERT( md0.getDirectionality() == RECV_ONLY );
CPPUNIT_ASSERT( md0.getDirectionalityOverride() == NOT_A_DIRECTION );
md0.setDirectionalityOverride( INACTIVE );
CPPUNIT_ASSERT( md0.getDirectionalityOverride() == INACTIVE );
CPPUNIT_ASSERT( md1.getType() == "video" );
CPPUNIT_ASSERT( md1.getDirectionality() == INACTIVE );
CPPUNIT_ASSERT( md2.getType() == "application" );
CPPUNIT_ASSERT( md2.getDirectionality() == SEND_ONLY );
CPPUNIT_ASSERT( md3.getType() == "audio" );
CPPUNIT_ASSERT( md3.getDirectionality() == SEND_RECV );
CPPUNIT_ASSERT( md4.getType() == "audio" );
CPPUNIT_ASSERT( md4.getDirectionality() == SEND_RECV );
}
virtual void SetUp()
{
n = 10;
vmf::initialize();
TEST_SCHEMA_NAME = "TEST_SCHEMA_NAME";
TEST_PROPERTY_NAME1 = "TEST_PROPERTY_NAME1";
TEST_PROPERTY_NAME2 = "TEST_PROPERTY_NAME2";
TEST_FIELD_NAME = "TEST_FIELD_NAME";
TEST_FIELD_NAME2 = "TEST_FIELD_NAME2";
schema = std::shared_ptr<vmf::MetadataSchema>(new vmf::MetadataSchema(TEST_SCHEMA_NAME));
VMF_METADATA_BEGIN(TEST_PROPERTY_NAME1)
VMF_FIELD_INT(TEST_FIELD_NAME)
VMF_FIELD_INT(TEST_FIELD_NAME2)
VMF_METADATA_END(schema);
VMF_METADATA_BEGIN(TEST_PROPERTY_NAME2)
VMF_FIELD_INT(TEST_FIELD_NAME)
VMF_METADATA_END(schema);
copyFile(TEST_FILE_SRC, TEST_FILE);
if (!stream.open(TEST_FILE, vmf::MetadataStream::ReadWrite))
std::cout << "Cann't open stream!" << std::endl;
stream.addSchema(schema);
for (int i = 0; i < n; i++)
{
std::shared_ptr<vmf::Metadata> md(new vmf::Metadata(schema->findMetadataDesc(TEST_PROPERTY_NAME1)));
md->setFieldValue(TEST_FIELD_NAME, i);
md->setFieldValue(TEST_FIELD_NAME2, 2*i);
stream.add(md);
std::shared_ptr<vmf::Metadata> md2(new vmf::Metadata(schema->findMetadataDesc(TEST_PROPERTY_NAME2)));
md2->setFieldValue(TEST_FIELD_NAME, n + i);
stream.add(md2);
}
stream.save();
stream.close();
stream.clear();
}
TEST(MetaDataTest, write)
{
PDMlib::MetaData md("MetaDataTest.txt");
//単位系情報を追加
PDMlib::UnitElem p = {"Pressure", "Pa", 2.0, 0.3, "true"};
PDMlib::UnitElem v = {"velocity", "m/s", 0.9, 30, "false"};
md.AddUnit(p);
md.AddUnit(v);
double bbox[6]={0.1, 0.2, 0.3, 0.4, 0.5, 0.6};
md.SetBoundingBox(bbox);
//コンテナを追加
PDMlib::ContainerInfo PV = {"ParticleVerocity", "", "zip", PDMlib::FLOAT, "vel", 3, PDMlib::NIJK};
PDMlib::ContainerInfo T = {"temperature", "", "zip", PDMlib::FLOAT, "temp", 1};
PDMlib::ContainerInfo id = {"ParticleID", "", "RLE", PDMlib::INT64, "id", 1};
md.AddContainer(PV);
md.AddContainer(T);
md.AddContainer(id);
EXPECT_EQ(0, md.Write());
// md.WriteTimeSlice();
PDMlib::MetaData md2("MetaDataTest.txt");
EXPECT_EQ(0, md2.Read());
double bbox2[6];
md2.GetBoundingBox(bbox2);
//mdとmd2が保持する値が等しいことを確認する
EXPECT_TRUE(md.Compare(md2));
//bounding boxはgoogle testのマクロで検査
for (int i =0; i<6; i++)
{
EXPECT_DOUBLE_EQ(bbox[i], bbox2[i]);
}
}
int main()
{
uint8_t hash[16];
char hexhash[33];
uint8_t *test_plain[] = {
(uint8_t *)"",
(uint8_t *)"a",
(uint8_t *)"abc",
(uint8_t *)"message digest",
(uint8_t *)"abcdefghijklmnopqrstuvwxyz",
(uint8_t *)"ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwxyz0123456789",
(uint8_t *)"12345678901234567890123456789012345678901234567890123456789012345678901234567890",
NULL};
char *test_hash[] = {
"8350e5a3e24c153df2275c9f80692773",
"32ec01ec4a6dac72c0ab96fb34c0b5d1",
"da853b0d3f88d99b30283a69e6ded6bb",
"ab4f496bfb2a530b219ff33031fe06b0",
"4e8ddff3650292ab5a4108c3aa47940b",
"da33def2a42df13975352846c30338cd",
"d5976f79d83d3a0dc9806c3c66f3efd8",
NULL};
for(int i = 0; test_plain[i] && test_hash[i]; i++)
{
md2(test_plain[i], strlen((char *)test_plain[i]), hash);
hex_encode(hexhash, hash, 16);
printf("%s - \"%s\"", hexhash, test_plain[i]);
if(strcmp(hexhash, test_hash[i]) == 0)
printf("\n");
else
printf("\n ERROR! Expected %s\n\n", test_hash[i]);
}
}
// *****************************************************************************
// Main
int main(int argc, char* const argv[])
{
try {
if (argc != 2) {
std::cout << "Usage: " << argv[0] << " file\n";
return 1;
}
std::string file(argv[1]);
std::cout <<"----- One IFD0 tag\n";
Exiv2::ExifData ed1;
Exiv2::Exifdatum md1(Exiv2::ExifKey("Exif.Image.Model"));
md1.setValue("Test 1");
ed1.add(md1);
write(file, ed1);
print(file);
std::cout <<"\n----- One Exif tag\n";
Exiv2::ExifData ed2;
Exiv2::Exifdatum md2(Exiv2::ExifKey("Exif.Photo.DateTimeOriginal"));
md2.setValue("Test 2");
ed2.add(md2);
write(file, ed2);
print(file);
// Todo: One Makernote tag for each Makernote
std::cout <<"\n----- One IOP tag\n";
Exiv2::ExifData ed3;
Exiv2::Exifdatum md3(Exiv2::ExifKey("Exif.Iop.InteroperabilityVersion"));
md3.setValue("Test 3");
ed3.add(md3);
write(file, ed3);
print(file);
std::cout <<"\n----- One GPS tag\n";
Exiv2::ExifData ed4;
Exiv2::Exifdatum md4(Exiv2::ExifKey("Exif.GPSInfo.GPSVersionID"));
md4.setValue("Test 4");
ed4.add(md4);
write(file, ed4);
print(file);
// Todo: Fix this
std::cout <<"\n----- One IFD1 tag\n";
Exiv2::ExifData ed5;
Exiv2::Exifdatum md5(Exiv2::ExifKey("Exif.Thumbnail.Artist"));
md5.setValue("Test 5");
ed5.add(md5);
Exiv2::Exifdatum md6(Exiv2::ExifKey("Exif.Image.Model"));
md6.setValue("Test 5 (Fix me!)");
ed5.add(md6);
write(file, ed5);
print(file);
return 0;
}
catch (Exiv2::Error& e) {
std::cout << "Caught Exiv2 exception '" << e << "'\n";
return -1;
}
}
void testMediaDescriptor_GetEndpointData()
{
const char *sdpOffer =
"v=0\r\n"
"o=mhandley 2890844526 2890842807 IN IP4 126.16.64.4\r\n"
"s=SDP Seminar\r\n"
"i=A Seminar on the session description protocol\r\n"
"u=http://www.cs.ucl.ac.uk/staff/M.Handley/sdp.03.ps\r\n"
"[email protected] (Mark Handley)\r\n"
"c=IN IP4 10.10.10.1\r\n"
"t=2873397496 2873404696\r\n"
"m=audio 10000 RTP/AVP 0\r\n" // media description 0
"c=IN IP4 10.10.10.2\r\n"
"a=recvonly\r\n"
"m=video 10002 RTP/AVP 31\r\n" // media description 1
"a=inactive\r\n"
"m=application 10004 udp wb\r\n" // media description 2
"a=sendonly\r\n"
"m=audio 10006 RTP/AVP 0\r\n" // media description 3
"c=IN IP4 10.10.10.3\r\n"
"a=sendrecv\r\n"
"m=audio 10008 RTP/AVP 0\r\n" // media description 4
;
SdpBody offerBody(sdpOffer);
const char *sdpAnswer =
"v=0\r\n"
"o=mhandley 2890844526 2890842807 IN IP4 126.16.64.4\r\n"
"s=SDP Seminar\r\n"
"i=A Seminar on the session description protocol\r\n"
"u=http://www.cs.ucl.ac.uk/staff/M.Handley/sdp.03.ps\r\n"
"[email protected] (Mark Handley)\r\n"
"c=IN IP4 20.10.10.1\r\n"
"t=2873397496 2873404696\r\n"
"m=audio 20000 RTP/AVP 0\r\n" // media description 0
"c=IN IP4 20.10.10.2\r\n"
"a=recvonly\r\n"
"m=video 20002 RTP/AVP 31\r\n" // media description 1
"a=inactive\r\n"
"m=application 20004 udp wb\r\n" // media description 2
"a=sendonly\r\n"
"m=audio 20006 RTP/AVP 0\r\n" // media description 3
"c=IN IP4 20.10.10.3\r\n"
"a=sendrecv\r\n"
"m=audio 20008 RTP/AVP 0\r\n" // media description 4
;
SdpBody answerBody(sdpAnswer);
MediaDescriptor md0( offerBody, 0, CALLER );
MediaDescriptor md1( offerBody, 1, CALLER );
MediaDescriptor md2( offerBody, 2, CALLER );
MediaDescriptor md3( offerBody, 3, CALLER );
MediaDescriptor md4( offerBody, 4, CALLER );
md0.setEndpointData( answerBody, 0, CALLEE );
md1.setEndpointData( answerBody, 1, CALLEE );
md2.setEndpointData( answerBody, 2, CALLEE );
md3.setEndpointData( answerBody, 3, CALLEE );
md4.setEndpointData( answerBody, 4, CALLEE );
MediaEndpoint mediaEndpoint;
// Media Descriptor 0
mediaEndpoint = md0.getEndpoint( CALLER );
CPPUNIT_ASSERT( mediaEndpoint.getAddress() == "10.10.10.2" );
CPPUNIT_ASSERT( mediaEndpoint.getRtpPort() == 10000 );
CPPUNIT_ASSERT( mediaEndpoint.getRtcpPort() == 10001 );
mediaEndpoint = md0.getEndpoint( CALLEE );
CPPUNIT_ASSERT( mediaEndpoint.getAddress() == "20.10.10.2" );
CPPUNIT_ASSERT( mediaEndpoint.getRtpPort() == 20000 );
CPPUNIT_ASSERT( mediaEndpoint.getRtcpPort() == 20001 );
// Media Descriptor 1
mediaEndpoint = md1.getEndpoint( CALLER );
CPPUNIT_ASSERT( mediaEndpoint.getAddress() == "10.10.10.1" );
CPPUNIT_ASSERT( mediaEndpoint.getRtpPort() == 10002 );
CPPUNIT_ASSERT( mediaEndpoint.getRtcpPort() == 10003 );
mediaEndpoint = md1.getEndpoint( CALLEE );
CPPUNIT_ASSERT( mediaEndpoint.getAddress() == "20.10.10.1" );
CPPUNIT_ASSERT( mediaEndpoint.getRtpPort() == 20002 );
CPPUNIT_ASSERT( mediaEndpoint.getRtcpPort() == 20003 );
// Media Descriptor 2
mediaEndpoint = md2.getEndpoint( CALLER );
CPPUNIT_ASSERT( mediaEndpoint.getAddress() == "10.10.10.1" );
CPPUNIT_ASSERT( mediaEndpoint.getRtpPort() == 10004 );
CPPUNIT_ASSERT( mediaEndpoint.getRtcpPort() == 10005 );
mediaEndpoint = md2.getEndpoint( CALLEE );
CPPUNIT_ASSERT( mediaEndpoint.getAddress() == "20.10.10.1" );
CPPUNIT_ASSERT( mediaEndpoint.getRtpPort() == 20004 );
CPPUNIT_ASSERT( mediaEndpoint.getRtcpPort() == 20005 );
// Media Descriptor 3
mediaEndpoint = md3.getEndpoint( CALLER );
CPPUNIT_ASSERT( mediaEndpoint.getAddress() == "10.10.10.3" );
CPPUNIT_ASSERT( mediaEndpoint.getRtpPort() == 10006 );
CPPUNIT_ASSERT( mediaEndpoint.getRtcpPort() == 10007 );
mediaEndpoint = md3.getEndpoint( CALLEE );
CPPUNIT_ASSERT( mediaEndpoint.getAddress() == "20.10.10.3" );
CPPUNIT_ASSERT( mediaEndpoint.getRtpPort() == 20006 );
CPPUNIT_ASSERT( mediaEndpoint.getRtcpPort() == 20007 );
// Media Descriptor 4
mediaEndpoint = md4.getEndpoint( CALLER );
CPPUNIT_ASSERT( mediaEndpoint.getAddress() == "10.10.10.1" );
CPPUNIT_ASSERT( mediaEndpoint.getRtpPort() == 10008 );
CPPUNIT_ASSERT( mediaEndpoint.getRtcpPort() == 10009 );
mediaEndpoint = md4.getEndpoint( CALLEE );
CPPUNIT_ASSERT( mediaEndpoint.getAddress() == "20.10.10.1" );
CPPUNIT_ASSERT( mediaEndpoint.getRtpPort() == 20008 );
CPPUNIT_ASSERT( mediaEndpoint.getRtcpPort() == 20009 );
}
int main(void) {
gslpp::complex zi = gslpp::complex::i();
std::vector<double> sd = {10., 5., 1.};
std::vector<gslpp::complex> sc = {30. + zi, 2. + 3. * zi, 1. + 4. * zi};
gslpp::matrix<double> md1(2, 2);
md1(0, 0) = 10.;
md1(1, 0) = 20.;
md1(0, 1) = 5.;
md1(1, 1) = 1.;
gslpp::matrix<double> md2(2, 2);
md2(0, 0) = -3.;
md2(1, 0) = 30.;
md2(0, 1) = -5.;
md2(1, 1) = 4.;
gslpp::matrix<gslpp::complex> mc1(2, 2);
mc1.assign(0, 0, 9. + 2 * zi);
mc1.assign(1, 0, 19. - 4 * zi);
mc1.assign(0, 1, 4. - 6 * zi);
mc1.assign(1, 1, 3. + 2 * zi);
gslpp::matrix<gslpp::complex> mc2(2, 2);
mc2.assign(0, 0, -8. + 3 * zi);
mc2.assign(1, 0, 11. - 5 * zi);
mc2.assign(0, 1, 2. + 5 * zi);
mc2.assign(1, 1, -3. + 4 * zi);
gslpp::vector<double> vd1(2);
gslpp::vector<double> vd2(2);
vd1(0) = 14.;
vd1(1) = -5;
vd2(0) = -9.;
vd2(1) = 3;
gslpp::vector<gslpp::complex> vc1(2);
gslpp::vector<gslpp::complex> vc2(2);
vc1.assign(0, 4. - 2. * zi);
vc1.assign(1, 6. - 8. * zi);
vc2.assign(0, 5. + 3. * zi);
vc2.assign(1, 4. - 12. * zi);
Expanded<double> esd(sd);
Expanded<gslpp::complex> esc(sc);
std::vector<gslpp::matrix<double> > mdv = {md1, md2};
Expanded<gslpp::matrix<double> > emd(mdv);
std::vector<gslpp::matrix<gslpp::complex> > mcv = {mc1, mc2};
Expanded<gslpp::matrix<gslpp::complex> > emc(mcv);
std::vector<gslpp::vector<double> > vdv = {vd1, vd2};
Expanded<gslpp::vector<double> > evd(vdv);
std::vector<gslpp::vector<gslpp::complex> > vcv = {vc1, vc2};
Expanded<gslpp::vector<gslpp::complex> > evc(vcv);
// Print Input
std::cout << std::endl;
std::cout << "esd " << esd << std::endl;
std::cout << "esc " << esc << std::endl;
std::cout << "evd " << evd << std::endl;
std::cout << "evc " << evc << std::endl;
std::cout << "emd " << emd << std::endl;
std::cout << "emc " << emc << std::endl;
std::cout << "-------------" << std::endl;
std::cout << "-sd " << -esd << std::endl;
std::cout << "-sc " << -esc << std::endl;
std::cout << "-vd " << -evd << std::endl;
std::cout << "-vc " << -evc << std::endl;
std::cout << "-md " << -emd << std::endl;
std::cout << "-mc " << -emc << std::endl;
std::cout << "-------------" << std::endl;
// Expanded * Expanded
std::cout << "sd*sd " << esd * esd << std::endl;
std::cout << "sd*sc " << esd * esc << std::endl;
std::cout << "sd*vd " << esd * evd << std::endl;
std::cout << "sd*vc " << esd * evc << std::endl;
std::cout << "sd*md " << esd * emd << std::endl;
std::cout << "sd*mc " << esd * emc << std::endl;
std::cout << "sc*sd " << esc * esd << std::endl;
std::cout << "sc*sc " << esc * esc << std::endl;
std::cout << "sc*vd " << esc * evd << std::endl;
std::cout << "sc*vc " << esc * evc << std::endl;
std::cout << "sc*md " << esc * emd << std::endl;
std::cout << "sc*mc " << esc * emc << std::endl;
std::cout << "vd*sd " << evd * esd << std::endl;
std::cout << "vd*sc " << evd * esc << std::endl;
std::cout << "vd*vd " << evd * evd << std::endl;
std::cout << "vd*vc " << evd * evc << std::endl;
std::cout << "vd*md " << evd * emd << std::endl;
std::cout << "vd*mc " << evd * emc << std::endl;
std::cout << "vc*sd " << evc * esd << std::endl;
std::cout << "vc*sc " << evc * esc << std::endl;
std::cout << "vc*vd " << evc * evd << std::endl;
std::cout << "vc*vc " << evc * evc << std::endl;
std::cout << "vc*md " << evc * emd << std::endl;
std::cout << "vc*mc " << evc * emc << std::endl;
std::cout << "md*sd " << emd * esd << std::endl;
std::cout << "md*sc " << emd * esc << std::endl;
std::cout << "md*vd " << emd * evd << std::endl;
std::cout << "md*vc " << emd * evc << std::endl;
std::cout << "md*md " << emd * emd << std::endl;
std::cout << "md*mc " << emd * emc << std::endl;
std::cout << "mc*sd " << emc * esd << std::endl;
std::cout << "mc*sc " << emc * esc << std::endl;
std::cout << "mc*vd " << emc * evd << std::endl;
std::cout << "mc*vc " << emc * evc << std::endl;
std::cout << "mc*md " << emc * emd << std::endl;
std::cout << "mc*mc " << emc * emc << std::endl;
std::cout << "-------------" << std::endl;
// Expanded + Expanded
std::cout << "sd + sd " << esd + esd << std::endl;
std::cout << "sd + sc " << esd + esc << std::endl;
std::cout << "sc + sd " << esc + esd << std::endl;
std::cout << "sc + sc " << esc + esc << std::endl;
std::cout << "vd + vd " << evd + evd << std::endl;
std::cout << "vd + vc " << evd + evc << std::endl;
std::cout << "vc + vd " << evc + evd << std::endl;
std::cout << "vc + vc " << evc + evc << std::endl;
std::cout << "md + md " << emd + emd << std::endl;
std::cout << "md + mc " << emd + emc << std::endl;
std::cout << "mc + md " << emc + emd << std::endl;
std::cout << "mc + mc " << emc + emc << std::endl;
//
std::cout << "-------------" << std::endl;
// Expanded - Expanded
std::cout << "sd - sd " << esd - esd << std::endl;
std::cout << "sd - sc " << esd - esc << std::endl;
std::cout << "sc - sd " << esc - esd << std::endl;
std::cout << "sc - sc " << esc - esc << std::endl;
std::cout << "vd - vd " << evd - evd << std::endl;
std::cout << "vd - vc " << evd - evc << std::endl;
std::cout << "vc - vd " << evc - evd << std::endl;
std::cout << "vc - vc " << evc - evc << std::endl;
std::cout << "md - md " << emd - emd << std::endl;
std::cout << "md - mc " << emd - emc << std::endl;
std::cout << "mc - md " << emc - emd << std::endl;
std::cout << "mc - mc " << emc - emc << std::endl;
std::cout << "-------------" << std::endl;
// Expanded * UnExpanded
std::cout << "esd*5 = " << esd * 5. << std::endl;
std::cout << "esd*(4-2I) = " << esd * (4. - 2. * zi) << std::endl;
std::cout << "esd*vd1 = " << esd * vd1 << std::endl;
std::cout << "esd*vc1 = " << esd * vc1 << std::endl;
std::cout << "esd*md1 = " << esd * md1 << std::endl;
std::cout << "esd*mc1 = " << esd * mc1 << std::endl;
//
std::cout << "esc*5 = " << esc * 5. << std::endl;
std::cout << "esc*(4-2I) = " << esc * (4. - 2. * zi) << std::endl;
std::cout << "esc*vd1 = " << esc * vd1 << std::endl;
std::cout << "esc*vc1 = " << esc * vc1 << std::endl;
std::cout << "esc*md1 = " << esc * md1 << std::endl;
std::cout << "esc*mc1 = " << esc * mc1 << std::endl;
//
std::cout << "evd*5 = " << evd * 5. << std::endl;
std::cout << "evd*(4-2I) = " << evd * (4. - 2. * zi) << std::endl;
std::cout << "evd*vd1 = " << evd * vd1 << std::endl;
std::cout << "evd*vc1 = " << evd * vc1 << std::endl;
std::cout << "evd*md1 = " << evd * md1 << std::endl;
std::cout << "evd*mc1 = " << evd * mc1 << std::endl;
//
std::cout << "evc*5 = " << evc * 5. << std::endl;
std::cout << "evc*(4-2I) = " << evc * (4. - 2. * zi) << std::endl;
std::cout << "evc*vd1 = " << evc * vd1 << std::endl;
std::cout << "evc*vc1 = " << evc * vc1 << std::endl;
std::cout << "evc*md1 = " << evc * md1 << std::endl;
std::cout << "evc*mc1 = " << evc * mc1 << std::endl;
//
std::cout << "emd*5 = " << emd * 5. << std::endl;
std::cout << "emd*(4-2I) = " << emd * (4. - 2. * zi) << std::endl;
std::cout << "emd*vd1 = " << emd * vd1 << std::endl;
std::cout << "emd*vc1 = " << emd * vc1 << std::endl;
std::cout << "emd*md1 = " << emd * md1 << std::endl;
std::cout << "emd*mc1 = " << emd * mc1 << std::endl;
//
std::cout << "emc*5 = " << emc * 5. << std::endl;
std::cout << "emc*(4-2I) = " << emc * (4. - 2. * zi) << std::endl;
std::cout << "emc*vd1 = " << emc * vd1 << std::endl;
std::cout << "emc*vc1 = " << emc * vc1 << std::endl;
std::cout << "emc*md1 = " << emc * md1 << std::endl;
std::cout << "emc*mc1 = " << emc * mc1 << std::endl;
std::cout << "-------------" << std::endl;
// // Expanded / UnExpanded-Scalar
std::cout << "esd/5 = " << esd / 5. << std::endl;
std::cout << "esd/(4-2I) = " << esd / (4. - 2. * zi) << std::endl;
std::cout << "esc/5 = " << esc / 5. << std::endl;
std::cout << "esc/(4-2I) = " << esc / (4. - 2. * zi) << std::endl;
std::cout << "evd/5 = " << evd / 5. << std::endl;
std::cout << "evd/(4-2I) = " << evd / (4. - 2. * zi) << std::endl;
std::cout << "evc/5 = " << evc / 5. << std::endl;
std::cout << "evc/(4-2I) = " << evc / (4. - 2. * zi) << std::endl;
std::cout << "emd/5 = " << emd / 5. << std::endl;
std::cout << "emd/(4-2I) = " << emd / (4. - 2. * zi) << std::endl;
std::cout << "emc/5 = " << emc / 5. << std::endl;
std::cout << "emc/(4-2I) = " << emc / (4. - 2. * zi) << std::endl;
std::cout << "-------------" << std::endl;
// UnExpanded * Expanded easy-check (must be 0)
std::cout << "zero1 = " << 5. * esd - esd * 5. << std::endl;
std::cout << "zero1 = " << 5. * esc - esc * 5. << std::endl;
std::cout << "zero1 = " << 5. * evd - evd * 5. << std::endl;
std::cout << "zero1 = " << 5. * evc - evc * 5. << std::endl;
std::cout << "zero1 = " << 5. * emd - emd * 5. << std::endl;
std::cout << "zero1 =" << 5. * emc - emc * 5. << std::endl;
//
std::cout << "zero2 = " << (4 - 2 * zi) * esd - esd * (4 - 2 * zi) << std::endl;
std::cout << "zero2 = " << (4 - 2 * zi) * esc - esc * (4 - 2 * zi) << std::endl;
std::cout << "zero2 = " << (4 - 2 * zi) * evd - evd * (4 - 2 * zi) << std::endl;
std::cout << "zero2 = " << (4 - 2 * zi) * evc - evc * (4 - 2 * zi) << std::endl;
std::cout << "zero2 = " << (4 - 2 * zi) * emd - emd * (4 - 2 * zi) << std::endl;
std::cout << "zero2 = " << (4 - 2 * zi) * emc - emc * (4 - 2 * zi) << std::endl;
std::cout << "zero3 = " << vd1 * esd - esd * vd1 << std::endl;
std::cout << "zero3 = " << vd1 * esc - esc * vd1 << std::endl;
std::cout << "zero3 = " << vd1 * evd - evd * vd1 << std::endl;
std::cout << "zero3 = " << vd1 * evc - evc * vd1 << std::endl;
std::cout << "zero3 = " << vd1 * emd - emd.transpose() * vd1 << std::endl;
std::cout << "zero3 = " << vd1 * emc - emc.transpose() * vd1 << std::endl;
//
std::cout << "zero4 = " << vc1 * esd - esd * vc1 << std::endl;
std::cout << "zero4 = " << vc1 * esc - esc * vc1 << std::endl;
std::cout << "zero4 = " << vc1 * evd - evd * vc1 << std::endl;
std::cout << "zero4 = " << vc1 * evc - evc * vc1 << std::endl;
std::cout << "zero4 = " << vc1 * emd - emd.transpose() * vc1 << std::endl;
std::cout << "zero4 = " << vc1 * emc - emc.transpose() * vc1 << std::endl;
//
std::cout << "zero5 = " << md1 * esd - esd * md1 << std::endl;
std::cout << "zero5 = " << md1 * esc - esc * md1 << std::endl;
std::cout << "zero5 = " << md1 * evd - evd * md1.transpose() << std::endl;
std::cout << "zero5 = " << md1 * evc - evc * md1.transpose() << std::endl;
std::cout << "zero5 = " << md1 * emd - (emd.transpose() * md1.transpose()).transpose() << std::endl;
std::cout << "zero5 = " << md1 * emc - (emc.transpose() * md1.transpose()).transpose() << std::endl;
//
std::cout << "zero6 = " << mc1 * esd - esd * mc1 << std::endl;
std::cout << "zero6 = " << mc1 * esc - esc * mc1 << std::endl;
std::cout << "zero6 = " << mc1 * evd - evd * mc1.transpose() << std::endl;
std::cout << "zero6 = " << mc1 * evc - evc * mc1.transpose() << std::endl;
std::cout << "zero6 = " << mc1 * emd - (emd.transpose() * mc1.transpose()).transpose() << std::endl;
std::cout << "zero6 = " << mc1 * emc - (emc.transpose() * mc1.transpose()).transpose() << std::endl;
std::cout << "true = " << (esd == esd) << std::endl;
std::cout << "true = " << (esc == esc) << std::endl;
std::cout << "true = " << (evd == evd) << std::endl;
std::cout << "true = " << (evc == evc) << std::endl;
std::cout << "true = " << (emd == emd) << std::endl;
std::cout << "true = " << (emc == emc) << std::endl;
std::cout << "false = " << (esd == esd * esd) << std::endl;
std::cout << "false = " << (esc == esc * esc) << std::endl;
std::cout << "false = " << (evd == 5 * evd) << std::endl;
std::cout << "false = " << (evc == 7 * evc) << std::endl;
std::cout << "false = " << (emd == emd * emd) << std::endl;
std::cout << "false = " << (emc == emc * emc) << std::endl;
Expanded<double> esd1(sd, 2);
std::vector<double> sdx = {10., 5., 1.,-7};
Expanded<double> esd2(sdx, 1);
std::cout << "esd1 * esd2 = " << esd1*esd2 << std::endl;
std::cout << "esd1 + esd2 = " << esd1+esd2 << std::endl;
std::cout << "esd1 - esd2 = " << esd1-esd2 << std::endl;
std::vector<double> q1={1.,0.};
std::vector<double> q2={5.,7.};
std::vector<std::vector<double> > v1 = {q1,q2};
std::vector<double> w1={3.,4.};
std::vector<double> w2={8.};
std::vector<std::vector<double> > v2 = {w1,w2};
Expanded<double> ev1(v1, 1, 0);
Expanded<double> ev2(v2);
std::cout << "ev1 = " << ev1 << std::endl;
std::cout << "ev2 = " << ev2 << std::endl;
std::cout << "ev1 * ev2 = " << ev1*ev2 << std::endl;
std::cout << "ev1 + ev2 = " << ev1+ev2 << std::endl;
std::cout << "1/ev1 = " << ev1.inverse() << std::endl;
std::cout << "1/ev2 = " << ev2.inverse() << std::endl;
complex ii(0, 1);
std::vector<complex> r1={1.+2.*ii, 3.+4.*ii};
std::vector<complex> r2={5. + 6.*ii, 7.+8.*ii, 9.+10.*ii};
std::vector<std::vector<complex> > vv1 = {r1,r2};
Expanded<complex> evv1(vv1, 1, 1);
std::cout << "evv1 = " << evv1 << std::endl;
std::cout << "1/evv1 = " << evv1.inverse() << std::endl;
std::cout << "ev1/evv1 = " << ev1/evv1 << std::endl;
std::cout << "evv1/ev1 = " << evv1/ev1 << std::endl;
Expanded<complex> tmp(evv1.truncate(std::vector<int>(2,2),2));
std::cout << "evv1 trunc at (std::vector<int>(2,2),2)" << tmp << std::endl;
std::cout << "evv1 series = " << evv1.Series(std::vector<int>(2,2),2) << std::endl;
Expanded<double> pp;
return (0);
}
