Skip to content

Segmentation Fault from DDECAL #343

@a-benati

Description

@a-benati

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

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions