/[escript]/branches/clazy/escriptcore/src/DataLazy.cpp
ViewVC logotype

Diff of /branches/clazy/escriptcore/src/DataLazy.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 6082 by jfenwick, Tue Mar 22 02:57:49 2016 UTC revision 6083 by jfenwick, Tue Mar 22 03:30:18 2016 UTC
# Line 1191  DataLazy::resolveNodeReduction(int tid, Line 1191  DataLazy::resolveNodeReduction(int tid,
1191            for (unsigned int z=0;z<ndpps;++z)            for (unsigned int z=0;z<ndpps;++z)
1192            {            {
1193              FMin op;              FMin op;
1194              *result=DataMaths::reductionOpVector(*leftres, m_left->getShape(), loffset, op, numeric_limits<double>::max());              *result=escript::reductionOpVector(*leftres, m_left->getShape(), loffset, op, numeric_limits<double>::max());
1195              loffset+=psize;              loffset+=psize;
1196              result++;              result++;
1197            }            }
# Line 1202  DataLazy::resolveNodeReduction(int tid, Line 1202  DataLazy::resolveNodeReduction(int tid,
1202            for (unsigned int z=0;z<ndpps;++z)            for (unsigned int z=0;z<ndpps;++z)
1203            {            {
1204            FMax op;            FMax op;
1205            *result=DataMaths::reductionOpVector(*leftres, m_left->getShape(), loffset, op, numeric_limits<double>::max()*-1);            *result=escript::reductionOpVector(*leftres, m_left->getShape(), loffset, op, numeric_limits<double>::max()*-1);
1206            loffset+=psize;            loffset+=psize;
1207            result++;            result++;
1208            }            }
# Line 1240  DataLazy::resolveNodeNP1OUT(int tid, int Line 1240  DataLazy::resolveNodeNP1OUT(int tid, int
1240      case SYM:      case SYM:
1241          for (loop=0;loop<numsteps;++loop)          for (loop=0;loop<numsteps;++loop)
1242          {          {
1243              DataMaths::symmetric(*leftres,m_left->getShape(),subroffset, m_samples, getShape(), offset);              escript::symmetric(*leftres,m_left->getShape(),subroffset, m_samples, getShape(), offset);
1244              subroffset+=step;              subroffset+=step;
1245              offset+=step;              offset+=step;
1246          }          }
# Line 1248  DataLazy::resolveNodeNP1OUT(int tid, int Line 1248  DataLazy::resolveNodeNP1OUT(int tid, int
1248      case NSYM:      case NSYM:
1249          for (loop=0;loop<numsteps;++loop)          for (loop=0;loop<numsteps;++loop)
1250          {          {
1251              DataMaths::nonsymmetric(*leftres,m_left->getShape(),subroffset, m_samples, getShape(), offset);              escript::nonsymmetric(*leftres,m_left->getShape(),subroffset, m_samples, getShape(), offset);
1252              subroffset+=step;              subroffset+=step;
1253              offset+=step;              offset+=step;
1254          }          }
# Line 1287  DataLazy::resolveNodeNP1OUT_P(int tid, i Line 1287  DataLazy::resolveNodeNP1OUT_P(int tid, i
1287      case TRACE:      case TRACE:
1288          for (loop=0;loop<numsteps;++loop)          for (loop=0;loop<numsteps;++loop)
1289          {          {
1290              DataMaths::trace(*leftres,m_left->getShape(),subroffset, m_samples ,getShape(),offset,m_axis_offset);              escript::trace(*leftres,m_left->getShape(),subroffset, m_samples ,getShape(),offset,m_axis_offset);
1291              subroffset+=instep;              subroffset+=instep;
1292              offset+=outstep;              offset+=outstep;
1293          }          }
# Line 1295  DataLazy::resolveNodeNP1OUT_P(int tid, i Line 1295  DataLazy::resolveNodeNP1OUT_P(int tid, i
1295      case TRANS:      case TRANS:
1296          for (loop=0;loop<numsteps;++loop)          for (loop=0;loop<numsteps;++loop)
1297          {          {
1298              DataMaths::transpose(*leftres,m_left->getShape(),subroffset, m_samples, getShape(),offset,m_axis_offset);              escript::transpose(*leftres,m_left->getShape(),subroffset, m_samples, getShape(),offset,m_axis_offset);
1299              subroffset+=instep;              subroffset+=instep;
1300              offset+=outstep;              offset+=outstep;
1301          }          }
# Line 1332  DataLazy::resolveNodeNP1OUT_2P(int tid, Line 1332  DataLazy::resolveNodeNP1OUT_2P(int tid,
1332      case SWAP:      case SWAP:
1333          for (loop=0;loop<numsteps;++loop)          for (loop=0;loop<numsteps;++loop)
1334          {          {
1335              DataMaths::swapaxes(*leftres,m_left->getShape(),subroffset, m_samples, getShape(),offset, m_axis_offset, m_transpose);              escript::swapaxes(*leftres,m_left->getShape(),subroffset, m_samples, getShape(),offset, m_axis_offset, m_transpose);
1336              subroffset+=instep;              subroffset+=instep;
1337              offset+=outstep;              offset+=outstep;
1338          }          }
# Line 1520  LAZYDEBUG(cout << "Right res["<< rroffse Line 1520  LAZYDEBUG(cout << "Right res["<< rroffse
1520    {    {
1521      case ADD:      case ADD:
1522          //PROC_OP(NO_ARG,plus<double>());          //PROC_OP(NO_ARG,plus<double>());
1523        DataMaths::binaryOpVectorLazyHelper<real_t, real_t, real_t>(resultp,        escript::binaryOpVectorLazyHelper<real_t, real_t, real_t>(resultp,
1524               &(*left)[0],               &(*left)[0],
1525               &(*right)[0],               &(*right)[0],
1526               chunksize,               chunksize,
# Line 1536  LAZYDEBUG(cout << "Right res["<< rroffse Line 1536  LAZYDEBUG(cout << "Right res["<< rroffse
1536               escript::ESFunction::PLUSF);                 escript::ESFunction::PLUSF);  
1537          break;          break;
1538      case SUB:      case SUB:
1539        DataMaths::binaryOpVectorLazyHelper<real_t, real_t, real_t>(resultp,        escript::binaryOpVectorLazyHelper<real_t, real_t, real_t>(resultp,
1540               &(*left)[0],               &(*left)[0],
1541               &(*right)[0],               &(*right)[0],
1542               chunksize,               chunksize,
# Line 1554  LAZYDEBUG(cout << "Right res["<< rroffse Line 1554  LAZYDEBUG(cout << "Right res["<< rroffse
1554          break;          break;
1555      case MUL:      case MUL:
1556          //PROC_OP(NO_ARG,multiplies<double>());          //PROC_OP(NO_ARG,multiplies<double>());
1557        DataMaths::binaryOpVectorLazyHelper<real_t, real_t, real_t>(resultp,        escript::binaryOpVectorLazyHelper<real_t, real_t, real_t>(resultp,
1558               &(*left)[0],               &(*left)[0],
1559               &(*right)[0],               &(*right)[0],
1560               chunksize,               chunksize,
# Line 1571  LAZYDEBUG(cout << "Right res["<< rroffse Line 1571  LAZYDEBUG(cout << "Right res["<< rroffse
1571          break;          break;
1572      case DIV:      case DIV:
1573          //PROC_OP(NO_ARG,divides<double>());          //PROC_OP(NO_ARG,divides<double>());
1574        DataMaths::binaryOpVectorLazyHelper<real_t, real_t, real_t>(resultp,        escript::binaryOpVectorLazyHelper<real_t, real_t, real_t>(resultp,
1575               &(*left)[0],               &(*left)[0],
1576               &(*right)[0],               &(*right)[0],
1577               chunksize,               chunksize,
# Line 1588  LAZYDEBUG(cout << "Right res["<< rroffse Line 1588  LAZYDEBUG(cout << "Right res["<< rroffse
1588          break;          break;
1589      case POW:      case POW:
1590         //PROC_OP(double (double,double),::pow);         //PROC_OP(double (double,double),::pow);
1591        DataMaths::binaryOpVectorLazyHelper<real_t, real_t, real_t>(resultp,        escript::binaryOpVectorLazyHelper<real_t, real_t, real_t>(resultp,
1592               &(*left)[0],               &(*left)[0],
1593               &(*right)[0],               &(*right)[0],
1594               chunksize,               chunksize,

Legend:
Removed from v.6082  
changed lines
  Added in v.6083

  ViewVC Help
Powered by ViewVC 1.1.26