void Prepreprocess::magic_quote(const std::string &line)
{
int end=line.length();
bool use_suffix=true;
state=START;
for (int pos=0; pos<end; ++pos) {
if (line[pos]=='$') {
if (pos == end-1) {
use_suffix=false;
break;
} else {
if (id0(line[pos+1])) {
++pos;
int begin=pos;
while (pos < end && id1(line[pos])) ++pos;
emit_expression(line.substr(begin,(pos-begin)));
--pos;
continue;
}
if (line[pos+1] == '(') {
++pos;
int begin=pos;
int parens=0;
while (pos < end) {
if (line[pos]=='(') ++parens;
if (line[pos]==')') {--parens; if (parens==0) { ++pos; break; } }
++pos;
}
emit_expression(line.substr(begin,(pos-begin)));
--pos;
continue;
}
if (line[pos+1]=='$') {
++pos;
}
}
}
emit_char(line[pos]);
}
emit_end(use_suffix);
}
int CallableWrapT<arr_xt<T> >::__getarray(Executor& ewsl,int pm)
{
if(pm==0)
{
ewsl.push(this);
}
else if(pm==1)
{
size_t sz=value.size();
Variant* _bp=ewsl.ci1.nsp;
idx_1t id0;
if(id0.update(_bp[0],0,sz)!=0)
{
ewsl.kerror("invalid array index");
}
if(id0.size==1)
{
_bp[0].reset(value[id0(0)]);
}
else
{
arr_type tmp;
tmp.resize(id0.size);
for(size_t k0=0;k0<id0.size;k0++)
{
tmp[k0]=value[id0(k0)];
}
_bp[0].reset(tmp);
}
}
else if(pm<=6)
{
sub_xt<arr_type> sub(value);
arr_xt_dims xtd=value.size_ptr();
Variant* _bp=ewsl.ci1.nsp-pm+1;
int flag=0;
for(int i=0;i<pm;i++)
{
flag|=sub.idx[i].update(_bp[i],0,xtd[i]);
xtd[i]=sub.idx[i].size;
}
for(int i=pm;i<6;i++)
{
sub.idx[i].size=xtd[i];
}
if(flag!=0)
{
ewsl.kerror("invalid array index");
}
if(xtd.size()==1)
{
_bp[0].reset(sub(0,0,0,0,0,0));
}
else
{
arr_type tmp;
tmp.resize(xtd);
for(size_t k5=0;k5<sub.idx[5].size;k5++)
for(size_t k4=0;k4<sub.idx[4].size;k4++)
for(size_t k3=0;k3<sub.idx[3].size;k3++)
for(size_t k2=0;k2<sub.idx[2].size;k2++)
for(size_t k1=0;k1<sub.idx[1].size;k1++)
for(size_t k0=0;k0<sub.idx[0].size;k0++)
{
tmp(k0,k1,k2,k3,k4,k5)=sub(k0,k1,k2,k3,k4,k5);
}
_bp[0].reset(tmp);
}
}
else
{
ewsl.kerror("too many dims");
}
return 1;
}
// this is d0 transformed
pardisoMatrix DDGMatrices::dual_d1( meshMetaInfo & aMesh )
{
pardisoMatrix d_0= d0(aMesh);
d_0*=(-1);//^1
return (id0(aMesh) % d_0);
}
int main(int argc, char **argv)
{
bool fail = false;
//
// now create a domain and a modelbuilder
// and build the model
//
FEM_ObjectBroker theBroker;
Domain *theDomain = new Domain();
FileDatastore *theDatabase
= new FileDatastore("/tmp/database/test1",*theDomain,theBroker);
FileDatastore &theDb = *theDatabase;
opserr << "TESTING IDs: \n";
ID id1(2);
id1(0) = 1;
id1(1) = 1;
ID id0(2);
id0(0) = 0;
id0(1) = 0;
ID id2(2);
id2(0) = 2;
id2(1) = 2;
ID id3(2);
id3(0) = 3;
id3(1) = 3;
ID id4(2);
id4(0) = 4;
id4(1) = 4;
ID id5(2);
id5(0) = 5;
id5(1) = 5;
ID id6(2);
id6(0) = 6;
id6(1) = 6;
theDb.sendID(1,1,id1);
theDb.sendID(1,2,id2);
theDb.sendID(2,1,id3);
theDb.sendID(2,2,id4);
opserr << "RESULTS\n";
ID recvID(2);
theDb.recvID(1,1,recvID);
opserr << "1: " << recvID;
theDb.recvID(1,2,recvID);
opserr << "2: " << recvID;
theDb.recvID(2,1,recvID);
opserr << "3: " << recvID;
theDb.recvID(2,2,recvID);
opserr << "4: " << recvID;
theDb.sendID(1,1,id1);
theDb.sendID(3,1,id3);
theDb.sendID(2,1,id2);
theDb.sendID(0,1,id0);
theDb.sendID(1,2,id1);
theDb.sendID(2,2,id2);
theDb.sendID(3,1,id3);
theDb.sendID(5,1,id5);
theDb.sendID(4,1,id4);
theDb.sendID(1,1,id1);
theDb.recvID(3,1,id5);
opserr << "3: " << id5;
theDb.recvID(1,1,id5);
opserr << "1: " << id5;
theDb.recvID(2,1,id5);
opserr << "2: " << id5;
theDb.recvID(1,2,id5);
opserr << "1: " << id5;
theDb.recvID(2,2,id5);
opserr << "3: " << id5;
theDb.recvID(3,1,id5);
opserr << "3: " << id5;
theDb.recvID(4,1,id5);
opserr << "4: " << id5;
theDb.recvID(5,1,id5);
opserr << "5: " << id5;
opserr << "FAILURE: " << theDb.recvID(6,1,id5) << " returned\n";
opserr << "FAILURE " << theDb.recvID(6,1,id5) << " returned\n";
theDb.recvID(0,1,id5);
opserr << "0: " << id5;
theDb.recvID(5,1,id5);
opserr << "5: " << id5;
ID id64(4);
id64(0) = 6;
id64(1) = 6;
id64(2) = 6;
id64(3) = 6;
theDb.sendID(6,1,id64);
theDb.recvID(6,1,id64);
opserr << id64;
opserr << "TESTING MATRICES: \n";
Matrix mat1(2,2);
mat1(0,0) = 1.1;
mat1(0,1) = 11.1;
mat1(1,0) = 111.1;
mat1(1,1) = 1111.1;
Matrix mat2(2,2);
mat2(0,0) = 2.2;
mat2(0,1) = 22.2;
mat2(1,0) = 222.2;
mat2(1,1) = 2222.2;
theDb.sendMatrix(2,1,mat2);
theDb.sendMatrix(1,1,mat1);
theDb.sendMatrix(3,2,mat2);
theDb.sendMatrix(3,1,mat1);
Matrix mat3(2,2);
theDb.recvMatrix(1,1,mat3);
opserr << mat1 << mat3 << endln;
theDb.recvMatrix(2,1,mat3);
opserr << mat2 << mat3 << endln;
theDb.recvMatrix(3,2,mat3);
opserr << mat2 << mat3 << endln;
theDb.recvMatrix(3,1,mat3);
opserr << mat1 << mat3 << endln;
// theDb.sendMatrix(2,1,mat1);
theDb.recvMatrix(2,1,mat3);
opserr << mat2 << mat3;
opserr << "TESTING VECTORS: \n";
Vector vect1(2);
vect1(0) = 1.1;
vect1(1) = 2.22;
Vector vect2(2);
vect2(0) = 3;
vect2(1) = 4.2;
Vector vect3(2);
vect3(0) = 5;
vect3(1) = 6;
Vector vect4(2);
vect4(0) = 7;
vect4(1) = 8.8e12;
theDb.sendVector(1,1,vect1);
theDb.sendVector(1,2,vect2);
theDb.sendVector(2,1,vect3);
theDb.sendVector(2,2,vect4);
opserr << "RESULTS\n";
Vector vect5(2);
theDb.recvVector(1,1,vect5);
opserr << vect1 << vect5 << endln;
theDb.recvVector(1,2,vect5);
opserr << vect2 << vect5 << endln;
theDb.recvVector(2,1,vect5);
opserr << vect3 << vect5 << endln;
theDb.recvVector(2,2,vect5);
opserr << vect4 << vect5 << endln;
theDb.sendVector(2,2,vect1);
theDb.sendVector(2,1,vect2);
theDb.sendVector(1,2,vect3);
theDb.sendVector(1,1,vect4);
theDb.recvVector(1,1,vect5);
opserr << vect4 << vect5 << endln;
theDb.recvVector(1,2,vect5);
opserr << vect3 << vect5 << endln;
theDb.recvVector(2,1,vect5);
opserr << vect2 << vect5 << endln;
theDb.recvVector(2,2,vect5);
opserr << vect1 << vect5 << endln;
theDb.sendVector(4,4,vect5);
theDb.recvVector(4,4,vect5);
opserr << vect5 << vect5 << endln;
theDb.recvVector(5,5,vect5);
opserr << "FAIL\n";
theDatabase->commitState(0);
/* */
/*
theDb.sendID(2,2,id1);
theDb.sendID(2,1,id2);
theDb.sendID(1,2,id3);
theDb.sendID(1,1,id4);
theDb.recvID(1,1,id5);
opserr << id5;
theDb.recvID(1,2,id5);
opserr << id5;
theDb.recvID(2,1,id5);
opserr << id5;
theDb.recvID(2,2,id5);
opserr << id5;
theDb.sendID(4,4,id5);
theDb.recvID(4,4,id5);
opserr << id5;
theDb.recvID(5,5,id5);
opserr << id5;
*/
/**************************
**************************/
/*
*/
// theModelBuilder->buildFE_Model();
// theDb.commitState(0);
// theDb.restoreState(0);
// theDb.restoreElements(0);
/*
theDb.restoreNode(1,0);
theDb.restoreNode(2,0);
theDb.restoreNode(3,0);
theDb.restoreNode(4,0);
theDb.restoreElement(0,0);
theDb.restoreElement(1,0);
theDb.restoreElement(2,0);
*/
// opserr << *theDomain;
delete theDatabase;
exit(0);
}
