Skip to content
Snippets Groups Projects
Commit 74f28a91 authored by meijeren's avatar meijeren
Browse files

BugID:345 

Added 3C343 test, using the entire BBS3 system.
parent 738e4bdd
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
CEP/BB/BBS3/test/BB3C343.input1 0 → 100644
+ 38
− 0
View file @ 74f28a91
# 3C343 test with:
# /data/PSS3/diepen/BBS2/3C343/10008336_p1.MS
# nrsrc=2
# nrstat=14
# solve for StokesI as function of freq
# model=""
# time domain is all times
# solve for entire time domain
nrPrediffers=1 # Number of Prediffers
CTRLparams.nrStrategies=1 # Number of Strategy Controllers
CTRLparams.SC1params.writeParms=T # Does the Controller write the final solution in parmeter table?
CTRLparams.SC1params.strategy="Simple" # Strategy type
CTRLparams.SC1params.maxNrIterations=3 # Maximum number of iterations per time interval
CTRLparams.SC1params.fitCriterion=-1 # Iterate until this absolute value is obtained for the fit (-1 means no fit criterion, iterate 'maxNrIterations' times)
CTRLparams.SC1params.controlParmUpdate=F # Does the controller update the parameters?
CTRLparams.SC1params.solvableParams=["StokesI.*"]
CTRLparams.SC1params.timeInterval=1.0e12 # Time interval
CTRLparams.SC1params.startTime=-1 # Start time (-1 means take beginning of measurement set)
CTRLparams.SC1params.endTime=-1 # End time (-1 means take end of measurement set)
CTRLparams.SC1params.antennas=[0,1,2,3,4,5,6,7,8,9,10,11,12,13]
CTRLparams.SC1params.startFreq=1170703125 # Start frequency
CTRLparams.SC1params.freqLength=8750000 # Size of frequency interval
CTRLparams.SC1params.sources=[] # Solvable source numbers (0 based)
CTRLparams.SC1params.modelType="" # Model type
CTRLparams.SC1params.calcUVW=F # Calculate UVW coordinates?
CTRLparams.SC1params.useAutoCorr=F # Use auto correlations
CTRLparams.SC1params.correlations=[0,3] # Only use XX and YY correlations
CTRLparams.SC1params.useSVD=F # Use Singular Value Decomposition in solve?
CTRLparams.SC1params.MSDBparams.MSName="10008336_p1.MS" # Name of the Measurement Set
CTRLparams.SC1params.MSDBparams.DBHost="localhost" # Parameter database host name
CTRLparams.SC1params.MSDBparams.DBType="aips" # Parameter database type
CTRLparams.SC1params.MSDBparams.DBPwd="" # Parameter database password
CTRLparams.SC1params.MSDBparams.DBMasterPort=13157 # port on which the master db is listening
CTRLparams.SC1params.MSDBparams.meqTableName="10008336" # Meq model table name *
CTRLparams.SC1params.MSDBparams.skyTableName="10008336_gsm" # Sky model table name *
CTRLparams.SC1params.MSDBparams.sourceGroups = ["srcs=[1,2]"]
CEP/BB/BBS3/test/BB3C343.input2 0 → 100644
+ 39
− 0
View file @ 74f28a91
# 3C343 test with:
# /data/PSS3/diepen/BBS2/3C343/10008336_p1.MS
# use middle 56 channels
# nrsrc=2
# nrstat=14
# solve for phase
# model="AP"
# solve for steps in time
nrPrediffers=1 # Number of Prediffers
CTRLparams.nrStrategies=1 # Number of Strategy Controllers
CTRLparams.SC1params.writeParms=T # Does the Controller write the final solution in parmeter table?
CTRLparams.SC1params.strategy="Simple" # Strategy type
CTRLparams.SC1params.maxNrIterations=3 # Maximum number of iterations per time interval
CTRLparams.SC1params.fitCriterion=-1 # Iterate until this absolute value is obtained for the fit (-1 means no fit criterion, iterate 'maxNrIterations' times)
CTRLparams.SC1params.controlParmUpdate=F # Does the controller update the parameters?
CTRLparams.SC1params.solvableParams=["EJ11.phase.*","EJ22.phase.*"]
CTRLparams.SC1params.excludeParams=["EJ11.phase.SR1.SG1","EJ11.phase.SR1","EJ22.phase.SR1.SG1","EJ22.phase.SR1"]
CTRLparams.SC1params.timeInterval=60 # Time interval
CTRLparams.SC1params.startTime=4472025720 # Start time (-1 means take beginning of measurement set)
CTRLparams.SC1params.endTime=4472068820 # End time (-1 means take end of measurement set)
CTRLparams.SC1params.antennas=[0,1,2,3,4,5,6,7,8,9,10,11,12,13]
CTRLparams.SC1params.startFreq=1170703125 # Start frequency
CTRLparams.SC1params.freqLength=8750000 # Size of frequency interval
CTRLparams.SC1params.sources=[] # Solvable source numbers (0 based)
CTRLparams.SC1params.modelType="AP" # Model type
CTRLparams.SC1params.calcUVW=F # Calculate UVW coordinates?
CTRLparams.SC1params.useAutoCorr=F # Use auto correlations
CTRLparams.SC1params.correlations=[0,3] # Only use XX and YY correlations
CTRLparams.SC1params.useSVD=F # Use Singular Value Decomposition in solve?
CTRLparams.SC1params.MSDBparams.MSName="10008336_p1.MS" # Name of the Measurement Set
CTRLparams.SC1params.MSDBparams.DBHost="localhost" # Parameter database host name
CTRLparams.SC1params.MSDBparams.DBType="aips" # Parameter database type
CTRLparams.SC1params.MSDBparams.DBPwd="" # Parameter database password
CTRLparams.SC1params.MSDBparams.DBMasterPort=13157 # port on which the master db is listening
CTRLparams.SC1params.MSDBparams.meqTableName="10008336" # Meq model table name *
CTRLparams.SC1params.MSDBparams.skyTableName="10008336_gsm" # Sky model table name *
CTRLparams.SC1params.MSDBparams.sourceGroups = ["srcs=[1,2]"]
CEP/BB/BBS3/test/run.BB3C343 0 → 100755
+ 127
− 0
View file @ 74f28a91
#!/bin/sh
# this script must be called with nriter and usernm
# This script does the same as the run.3C343 script. However,
# this script uses the whole BBS3 system, instead of just
# the Prediffer and Solver class.
# It uses the input parameters sets in files: BB3C343.input1
# and BB3C343.input2.
if [ "$#" != "2" ]; then
echo "run as: run.3C343 usernm dbtype"
exit 1
fi
usernm=$1
dbtype=$2
# Define functions to create the parm tables.
beginParmdbInputFile () {
rm -f fill.parmdb
add2file "create $*"
add2file "connect $*"
add2file "clear"
}
add2file () {
echo $* >> fill.parmdb
}
applyfile () {
add2file "quit"
./parmdb < fill.parmdb
}
makeMeq () {
diff="1e-6"
if [ "$#" -eq "1" ]; then
diff=$1
fi
add2file "adddef frot values=0, diff=${diff}"
add2file "adddef drot values=0, diff=${diff}"
add2file "adddef dell values=0, diff=${diff}"
add2file "adddef gain.11 values=1, diff=${diff}"
add2file "adddef gain.22 values=1, diff=${diff}"
add2file "adddef EJ11.ampl values=1, diff=${diff}"
add2file "adddef EJ12.ampl values=0, diff=${diff}"
add2file "adddef EJ21.ampl values=0, diff=${diff}"
add2file "adddef EJ22.ampl values=1, diff=${diff}"
add2file "adddef EJ11.phase values=0, diff=${diff}"
add2file "adddef EJ12.phase values=0, diff=${diff}"
add2file "adddef EJ21.phase values=0, diff=${diff}"
add2file "adddef EJ22.phase values=0, diff=${diff}"
}
makeGsm () {
diff="1e-6"
if [ "$#" -eq "1" ]; then
diff=$1
fi
add2file "adddef RA.CP1 values=4.356645791155902, srcnr=0, diff=${diff}, diff_rel=F"
add2file "adddef DEC.CP1 values=1.092208429052697, srcnr=0, diff=${diff}, diff_rel=F"
# add2file "adddef RA.CP1 values=4.356645, srcnr=0, diff=${diff}, diff_rel=F"
# add2file "adddef DEC.CP1 values=1.092208, srcnr=0, diff=${diff}, diff_rel=F"
add2file "adddef StokesI.CP1 values=1.0, srcnr=0, diff=${diff} "
add2file "adddef StokesQ.CP1 values=0, srcnr=0, diff=${diff} "
add2file "adddef StokesU.CP1 values=0, srcnr=0, diff=${diff} "
add2file "adddef StokesV.CP1 values=0, srcnr=0, diff=${diff} "
add2file "adddef RA.CP2 values=4.3396003966265599, srcnr=1, diff=${diff}, diff_rel=F"
add2file "adddef DEC.CP2 values=1.0953677174056471, srcnr=1, diff=${diff}, diff_rel=F"
# add2file "adddef RA.CP2 values=4.339600, srcnr=1, diff=${diff}, diff_rel=F"
# add2file "adddef DEC.CP2 values=1.095367, srcnr=1, diff=${diff}, diff_rel=F"
add2file "adddef StokesI.CP2 values=1.0, srcnr=1, diff=${diff} "
add2file "adddef StokesQ.CP2 values=0, srcnr=1, diff=${diff} "
add2file "adddef StokesU.CP2 values=0, srcnr=1, diff=${diff} "
add2file "adddef StokesV.CP2 values=0, srcnr=1, diff=${diff} "
}
makeGsm2 () {
diff="1e-6"
if [ "$#" -eq "1" ]; then
diff=$1
fi
add2file "adddef RA.CP1 values=4.356645791155902, srcnr=0, diff=${diff}, diff_rel=F"
add2file "adddef DEC.CP1 values=1.092208429052697, srcnr=0, diff=${diff}, diff_rel=F"
# add2file "adddef RA.CP1 values=4.356645, srcnr=0, diff=${diff}, diff_rel=F"
# add2file "adddef DEC.CP1 values=1.092208, srcnr=0, diff=${diff}, diff_rel=F"
add2file "adddef StokesI.CP1 nx=2, values=[1.0,0.], srcnr=0, diff=${diff} "
add2file "adddef StokesQ.CP1 values=0, srcnr=0, diff=${diff} "
add2file "adddef StokesU.CP1 values=0, srcnr=0, diff=${diff} "
add2file "adddef StokesV.CP1 values=0, srcnr=0, diff=${diff} "
add2file "adddef RA.CP2 values=4.3396003966265599, srcnr=1, diff=${diff}, diff_rel=F"
add2file "adddef DEC.CP2 values=1.0953677174056471, srcnr=1, diff=${diff}, diff_rel=F"
# add2file "adddef RA.CP2 values=4.339600, srcnr=1, diff=${diff}, diff_rel=F"
# add2file "adddef DEC.CP2 values=1.095367, srcnr=1, diff=${diff}, diff_rel=F"
add2file "adddef StokesI.CP2 nx=4, values=[1.0,0.,0.,0.], srcnr=1, diff=${diff} "
add2file "adddef StokesQ.CP2 values=0, srcnr=1, diff=${diff} "
add2file "adddef StokesU.CP2 values=0, srcnr=1, diff=${diff} "
add2file "adddef StokesV.CP2 values=0, srcnr=1, diff=${diff} "
}
echo " Generating parm tables ..."
# use functions to initialize parmtables
beginParmdbInputFile "dbtype='$dbtype', user='$usernm', tablename='10008336'"
makeMeq
applyfile
beginParmdbInputFile "dbtype='$dbtype', user='$usernm', tablename='10008336_gsm'"
makeGsm2
applyfile
ln -s /data/PSS3/diepen/BBS2/3C343/10008336_p1.MS
./clearBB
# Solve for StokesI as function of freq.
time ./BBS3 ../../test/BB3C343.input1 $usernm
# Create a new gsm table.
beginParmdbInputFile "dbtype='$dbtype', user='$usernm', tablename='10008336_gsm'"
makeGsm
applyfile
./clearBB
# Solve for StokesI (and save that).
time ./BBS3 ../../test/BB3C343.input1 $usernm
./clearBB
# Solve for phase
time ./BBS3 ../../test/BB3C343.input2 $usernm
CEP/BB/BBS3/test/t3C343.cc
+ 1
− 1
View file @ 74f28a91
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment