int TestConnProbe::outputState(double timed){
//Grab weights of probe and test for the value of .625/1.5, or .4166666
HyPerConn* conn = getTargetHyPerConn();
int numPreExt = conn->preSynapticLayer()->getNumExtended();
int syw = conn->yPatchStride(); // stride in patch
for(int kPre = 0; kPre < numPreExt; kPre++){
PVPatch * weights = conn->getWeights(kPre, 0);
int nk = conn->fPatchSize() * weights->nx;
pvwdata_t * data = conn->get_wData(0,kPre);
int ny = weights->ny;
for (int y = 0; y < ny; y++) {
pvwdata_t * dataYStart = data + y * syw;
for(int k = 0; k < nk; k++){
if(fabs(timed - 0) < (parent->getDeltaTime()/2)){
if(fabs(dataYStart[k] - 1) > .01){
std::cout << "dataYStart[k]: " << dataYStart[k] << "\n";
}
assert(fabs(dataYStart[k] - 1) <= .01);
}
else if(fabs(timed - 1) < (parent->getDeltaTime()/2)){
if(fabs(dataYStart[k] - 1.375) > .01){
std::cout << "dataYStart[k]: " << dataYStart[k] << "\n";
}
assert(fabs(dataYStart[k] - 1.375) <= .01);
}
}
}
}
return PV_SUCCESS;
}
int MomentumTestConnProbe::outputState(double timed){
//Grab weights of probe and test for the value of .625/1.5, or .4166666
HyPerConn* conn = getTargetHyPerConn();
int numPreExt = conn->preSynapticLayer()->getNumExtended();
int syw = conn->yPatchStride(); // stride in patch
for(int kPre = 0; kPre < numPreExt; kPre++){
PVPatch * weights = conn->getWeights(kPre, 0);
int nk = conn->fPatchSize() * weights->nx;
pvwdata_t * data = conn->get_wData(0,kPre);
int ny = weights->ny;
pvdata_t wCorrect;
for (int y = 0; y < ny; y++) {
pvwdata_t * dataYStart = data + y * syw;
for(int k = 0; k < nk; k++){
pvdata_t wObserved = dataYStart[k];
if(timed < 3){
wCorrect = 0;
}
else{
wCorrect = .376471;
for(int i = 0; i < (timed-3); i++){
wCorrect += .376471 * exp(-(2*(i+1)));
}
}
assert(fabs(wObserved - wCorrect) <= 1e-4);
}
}
}
return PV_SUCCESS;
}