Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Segmentation Fault from DDECAL #343

Open
a-benati opened this issue Jun 11, 2024 · 0 comments
Open

Segmentation Fault from DDECAL #343

a-benati opened this issue Jun 11, 2024 · 0 comments

Comments

@a-benati
Copy link

Hi all,

I get a segmentation fault after launching DP3 6.0.0 on my data, for which I want solutions in 4 different directions. The DP3 command is the following (it is launched from within loafer_facetselfcal, therefore some parameters are set by default and not by me):

DP3 numthreads=64 msin=1685906777_sdp_l0-Abell3667-corr_avg.ms.copy msout=. ddecal.mode=sc
alarphase msin.weightcolumn=WEIGHT_SPECTRUM steps=[ddecal] ddecal.type=ddecal msout.storagemanager=dy
sco msout.storagemanager.weightbitrate=16 ddecal.maxiter=100 ddecal.propagatesolutions=True ddecal.so
lint=3 ddecal.solutions_per_direction='[1,1,1,1]' ddecal.nchan=1 ddecal.h5parm=scalarphase0_selfcalcy
le000_1685906777_sdp_l0-Abell3667-corr_avg.ms.copy.h5 ddecal.modeldatacolumns='[MODEL_DATA_DD0,MODEL_
DATA_DD1,MODEL_DATA_DD2,MODEL_DATA_DD3]' ddecal.solveralgorithm=directioniterative msin.datacolumn=DA
TA ddecal.uvlambdamin=350.0 ddecal.smoothnessconstraint=100000000.0 ddecal.smoothnessreffrequency=0.0
 ddecal.smoothnessspectralexponent=-1.0 ddecal.smoothnessrefdistance=0.0 ddecal.tolerance=0.0001

And this is the output:

 MSReader   
  input MS:       /data/abell_3667/ddcal_pipeline/1685906777_sdp_l0-Abell3667-corr_avg.ms.copy
  band            0
  startchan:      0  (0)
  nchan:          512  (0)
  ncorrelations:  4
  nbaselines:     1891
  first time:     2023/06/04/19:52:30
  last time:      2023/06/05/05:12:02
  ntimes:         1402
  time interval:  23.963
  DATA column:    DATA
  WEIGHT column:  WEIGHT_SPECTRUM
  FLAG column:    FLAG
  autoweight:     false
DDECal ddecal.
  mode (constraints):  scalarphase
  algorithm:           directioniterative
  H5Parm:              scalarphase0_selfcalcyle000_1685906777_sdp_l0-Abell3667-corr_avg.ms.copy.h5
  write sol to buffer: false
  solution interval:   3
  nchan:               1
  direction count:     4
  directions:          [[MODEL_DATA_DD0],[MODEL_DATA_DD1],[MODEL_DATA_DD2],[MODEL_DATA_DD3]]
  sols per direction:  [1,1,1,1]
  tolerance:           0.0001
  max iter:            100
  flag unconverged:    false
     diverged only:    false
  propagate solutions: true
       converged only: false
  detect stalling:     true
  step size:           0.2
  smoothnessconstraint:1e+08
  approximate fitter:  false
  only predict:        false
  subtract model:      false
  keep model:          false
Model steps for direction MODEL_DATA_DD0
MsColumnReader ddecal.
  column:      MODEL_DATA_DD0


Model steps for direction MODEL_DATA_DD1
MsColumnReader ddecal.
  column:      MODEL_DATA_DD1


Model steps for direction MODEL_DATA_DD2
MsColumnReader ddecal.
  column:      MODEL_DATA_DD2


Model steps for direction MODEL_DATA_DD3
MsColumnReader ddecal.
  column:      MODEL_DATA_DD3


UVWFlagger ddecal.
  uvm:            []
  um:             []
  vm:             []
  wm:             []
  uvlambda:       [-1e+15,350]
  ulambda:        []
  vlambda:        []
  wlambda:        []
  phasecenter:    []
MSUpdater msout.
  MS:             /data/abell_3667/ddcal_pipeline/1685906777_sdp_l0-Abell3667-corr_avg.ms.copy
  datacolumn:     DATA
  flagcolumn:     FLAG
  weightcolumn    WEIGHT_SPECTRUM
  Compressed:     yes
  Data bitrate:   10
  Weight bitrate: 16
  Dysco mode:     AF TruncatedGaussian(2.5)

  flush:          0

*** WARNING: the following parset keywords were not used ***
             maybe they are misspelled
    ddecal.smoothnessspectralexponent
Processing 1402 time slots ...

0%....10....20....30....40....50....60....70....80....90....100%
Segmentation fault

I attach here the H5Parm.
scalarphase0_h5.tgz

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant