Commits
Ville Suoranta authored 853e1e55261 Merge
283 283 | sdbs.aggregateScan(), |
284 284 | sdbs.aggregateTime(), |
285 285 | sdbs.aggregateSpw(), |
286 286 | sdbs.freqs(), |
287 287 | sdbs.aggregateFld()); |
288 288 | Dsol.sizeSolveParCurrSpw(sdbs.nChannels()); |
289 289 | |
290 290 | if (DJONES_TEST_VERBOSE) |
291 291 | Dsol.state(); |
292 292 | |
293 - | Dsol.guessPar(sdbs); |
293 + | // call guessPar using generic signature, tests CAS-13179 |
294 + | //Dsol.guessPar(sdbs); |
295 + | SolvableVisCal* svc_=&Dsol; |
296 + | svc_->guessPar(sdbs,False); |
294 297 | |
295 - | Dsol.reportSolvedQU(); |
296 - | |
297 298 | if (DJONES_TEST_VERBOSE) |
298 299 | cout << "srcPolPar() = " << real(Dsol.srcPolPar()) << endl; |
299 300 | ASSERT_NEAR(0.04,real(Dsol.srcPolPar()(0)),1e-3); |
300 301 | ASSERT_NEAR(0.03,real(Dsol.srcPolPar()(1)),1e-3); |
301 302 | |
302 303 | /* // explore time-dep QU signal... |
303 304 | for (Int i=0;i<sdbs.nSDB();++i) { |
304 305 | const Complex& RL(sdbs(i).visCubeModel()(1,0,0)); |
305 306 | cout << i << " " << RL << " " << abs(RL) << " " << arg(RL)*180.0/C::pi << endl; |
306 307 | } |
307 308 | cout << endl << endl; |
308 309 | */ |
309 310 | |
310 311 | // Now solve for the Ds |
312 + | // This calls LLS version since Dsol is DfllsJones |
311 313 | Dsol.selfSolveOne(sdbs); |
312 314 | |
313 315 | // Apply refant (ad hoc) |
314 316 | Dref(Dsol.solveAllCPar(),Dsol.refant(),d(0,0,Dsol.refant())); |
315 317 | |
316 318 | |
317 319 | Cube<Float> soldiff=amplitude(Dsol.solveAllCPar()-d); |
318 320 | Cube<Float> reldiff=soldiff/amplitude(d); |
319 321 | |
320 322 | /* |