Commits
691 691 | |
692 692 | convFunctions_p[actualConvIndex_p]= new Array<Complex>(); |
693 693 | convWeights_p[actualConvIndex_p]= new Array<Complex>(); |
694 694 | convSupportBlock_p[actualConvIndex_p]=new Vector<Int>(); |
695 695 | Int newConvSize=2*(max(convSupport_p)+2)*convSampling; |
696 696 | Int newRealConvSize=newConvSize* Int(Double(convSamp)/Double(convSampling)); |
697 697 | Int lattSize=convFuncTemp.shape()(0); |
698 698 | (*convSupportBlock_p[actualConvIndex_p])=convSupport_p; |
699 699 | |
700 700 | //cerr << "convSize " << newConvSize << " " << newRealConvSize<< " lattSize " << lattSize << endl; |
701 + | LogIO os(LogOrigin("HetArrConvFunc", "findConvFunction", WHERE)); |
702 + | os << "convolution function support: " << convSupport_p << LogIO::POST; |
701 703 | |
702 704 | if(newConvSize < lattSize){ |
703 705 | IPosition blc(5, (lattSize/2)-(newConvSize/2), |
704 706 | (lattSize/2)-(newConvSize/2),0,0,0); |
705 707 | IPosition trc(5, (lattSize/2)+(newConvSize/2-1), |
706 708 | (lattSize/2)+(newConvSize/2-1), nBeamPols-1, nBeamChans-1,ndishpair-1); |
707 709 | IPosition shp(5, newConvSize, newConvSize, nBeamPols, nBeamChans, ndishpair); |
708 710 | |
709 711 | convFunctions_p[actualConvIndex_p]= new Array<Complex>(IPosition(5, newRealConvSize, newRealConvSize, nBeamPols, nBeamChans, ndishpair )); |
710 712 | convWeights_p[actualConvIndex_p]= new Array<Complex>(IPosition(5, newRealConvSize, newRealConvSize, nBeamPols, nBeamChans, ndishpair )); |