1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
|
/*********************************************************************
MakeProfiles - Create mock astronomical profiles.
MakeProfiles is part of GNU Astronomy Utilities (Gnuastro) package.
Original author:
Mohammad Akhlaghi <mohammad@akhlaghi.org>
Contributing author(s):
Copyright (C) 2015-2025 Free Software Foundation, Inc.
Gnuastro is free software: you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by the
Free Software Foundation, either version 3 of the License, or (at your
option) any later version.
Gnuastro is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
General Public License for more details.
You should have received a copy of the GNU General Public License
along with Gnuastro. If not, see <http://www.gnu.org/licenses/>.
**********************************************************************/
#ifndef MOCKGALS_H
#define MOCKGALS_H
#include <gnuastro/threads.h>
#include "main.h"
struct mkonthread
{
/* General parameters: */
double r; /* Elliptical radius at this point. */
double coord[3]; /* Pixel coord. (from profile center). */
double lower[3]; /* Coordinates of lower pixel position. */
double higher[3]; /* Coordinates of higher pixel position. */
double c[3]; /* Cosine of position angle(s). */
double s[3]; /* Sine of position angle(s). */
double q[2]; /* Axis ratio(s). */
double center[3]; /* Center (in FITS) in oversampled image.*/
double (*profile)(struct mkonthread *); /* Function to use. */
double truncr; /* Truncation radius in pixels. */
double intruncr; /* Inner truncation radius in pixels. */
long width[3]; /* Enclosing box in FITS axes, not C. */
float peakflux; /* Flux at profile peak. */
float brightness; /* The brightness of the profile. */
uint8_t func; /* Radial function code of the profile. */
long *onaxes; /* Sides of the unover-sampled image. */
long fpixel_i[3]; /* fpixel_i before running overlap. */
int correction; /* ==1: correct the pixels afterwards. */
unsigned long rng_seed; /* Seed used to generate this profile. */
gal_data_t *customimg; /* Custom image for this profile. */
/* Random number generator: */
gsl_rng *rng; /* Copy of main random number generator. */
/* Profile specific parameters: */
double sersic_re; /* r/re in Sersic profile. */
double sersic_inv_n; /* Sersic index of Sersic profile. */
double sersic_nb; /* Negative of b(n) constant. */
double moffat_alphasq; /* r divided by alpha in Moffat. */
double moffat_nb; /* Negative beta in the Moffat. */
double gaussian_c; /* Constant value in Gaussian. */
double fixedvalue; /* Value of a point source. */
/* General parameters */
struct mkprofparams *p; /* Pointer to the main.h structure. */
size_t *indexs; /* Indexs to build on this thread. */
pthread_barrier_t *b; /* Pthread barrier pointer. */
struct builtqueue *ibq; /* Internally built queue. */
};
void
mkprof(struct mkprofparams *p);
#endif
|
