Commits
479 479 | if jwidth == -1: |
480 480 | jwidth = 1.55 |
481 481 | gwidth_arcsec = self.__parse_width(gwidth, cell_arcsec) |
482 482 | jwidth_arcsec = self.__parse_width(jwidth, cell_arcsec) |
483 483 | mygjinc = self.trunc(self.gjinc(axis, gwidth=gwidth_arcsec, |
484 484 | jwidth=jwidth_arcsec, |
485 485 | useCasaJinc=True, normalize=False)) |
486 486 | return mygjinc |
487 487 | |
488 488 | def gjinc(self, x, gwidth, jwidth, useCasaJinc=False, normalize=False): |
489 - | """ |
490 - | Migrated from AnalysisUtils (revision 1.2204, 2015/02/18) |
491 - | """ |
489 + | """Migrated from AnalysisUtils (revision 1.2204, 2015/02/18)""" |
492 490 | if (useCasaJinc): |
493 491 | result = self.grdjinc1(x, jwidth, normalize) * self.gjincGauss(x, gwidth) |
494 492 | else: |
495 493 | result = self.jinc(x, jwidth) * self.gjincGauss(x, gwidth) |
496 494 | return result |
497 495 | |
498 496 | def grdjinc1(self, val, c, normalize=True): |
499 - | """ |
500 - | Migrated from AnalysisUtils (revision 1.2204, 2015/02/18) |
501 - | """ |
497 + | """Migrated from AnalysisUtils (revision 1.2204, 2015/02/18)""" |
502 498 | # Casa's function |
503 499 | # // Calculate J_1(x) using approximate formula |
504 500 | xs = np.pi * val / c |
505 501 | result = [] |
506 502 | for x in xs: |
507 503 | x = abs(x) # I added this to make it symmetric |
508 504 | ax = abs(x) |
509 505 | if (ax < 8.0): |
510 506 | y = x * x |
511 507 | ans1 = x * (72362614232.0 + y * |