libs/dtgf/detrend_C99.c

// $Header: $

#include <stdio.h>
#include <stdlib.h>
#include <assert.h>
#include <string.h>
#include <complex.h>
#include <math.h>
#include "ctypes.h"
#include "detrend_C99.h"
//#include "xmem.h"

//#define DEBUGFILE  //uncomment this to write debugging files

#ifdef FLOAT
  #define TYPE FLOAT
  #include "detrend_code_C99.h" 
  #undef TYPE
#endif

#ifdef DOUBLE
  #define TYPE DOUBLE
  #include "detrend_code_C99.h" 
  #undef TYPE
#endif

#ifdef COMPLEXFLOAT
static int count2=0;
void cdetrend_discontig( int n, _Complex float *data, int *isgood, 
                         int degree, int length, int skip, 
                         int m, int *sect, int detrend_first)
{
  int i,first,last;
  // Seperately detrend m sections of data separated by 
  // discontinuities.
  // The first section begins on data[0], 
  // the m-1 first sections end on data[sect_last[i]], i=0,1,...,m-2,
  // and the last section ends on data[n-1].

  if (m==0)
    cdetrend(n,data,isgood,degree,length,skip,detrend_first);
  else
  {
    for (i=0; i<m; i++)
    {
      first = sect[2*i];
      last = sect[2*i+1];
#ifdef DEBUGOUT
      printf("Detrending [%d:%d]\n",first,last);
#endif
      cdetrend(last-first+1,&data[first],&isgood[first],degree,length,skip,detrend_first);
    }
  }
}

void cdetrend( int n, _Complex float *data, int *isgood, int degree, 
               int length, int skip, int detrend_first)
{
  int i;

  float *real, *imag;
  real = (float *) malloc(n*sizeof(float));
  for (i=0;i<n;i++)
    real[i] = crealf(data[i]);
  sdetrend( n, real, isgood, degree, length, skip, detrend_first);

  imag = (float *) malloc(n*sizeof(float));
  for (i=0;i<n;i++)
    imag[i] = cimagf(data[i]);
  sdetrend( n, imag, isgood, degree, length, skip, detrend_first);
  for (i=0;i<n;i++)
    data[i] = real[i] + _Complex_I*imag[i];

#ifdef DEBUGFILE
 {
   FILE *fh;
   char name[1024];
  printf("writing debug2 file #%d\n",count2);
  sprintf(name,"detrend_debug2_%d.out",count2++);
  fh = fopen(name,"w");
  assert(fh);
  for (i=0; i<n; i++)
    fprintf(fh,"%e %e\n",crealf(data[i]),cimagf(data[i]));
  fclose(fh);
 }
#endif

  free(real);
  free(imag);
}
#endif


#ifdef COMPLEXDOUBLE
void zdetrend_discontig( int n, _Complex double *data, int *isgood, 
                         int degree, int length, int skip, 
                         int m, int *sect, int detrend_first)
{
  int i,first,last;
  // Seperately detrend m sections of data separated by 
  // discontinuities.
  // The first section begins on data[0], 
  // the m-1 first sections end on data[sect_last[i]], i=0,1,...,m-2,
  // and the last section ends on data[n-1].

  if (m==0)
    zdetrend(n,data,isgood,degree,length,skip,detrend_first);
  else
  {
    first = 0;
    for (i=0; i<m; i++)
    {
      first = sect[2*i];
      last = sect[2*i+1];
      zdetrend(last-first+1,&data[first],&isgood[first],degree,length,skip,detrend_first);
    }
  }
}


void zdetrend( int n, _Complex double *data, int *isgood, int degree, 
               int length, int skip, int detrend_first)
{
  int i;
  double *real, *imag;

  real = (double *) malloc(n*sizeof(double));
  for (i=0;i<n;i++)
    real[i] = creal(data[i]);
  ddetrend( n, real, isgood, degree, length, skip, detrend_first);

  imag = (double *) malloc(n*sizeof(double));
  for (i=0;i<n;i++)
    imag[i] = cimag(data[i]);
  ddetrend( n, imag, isgood, degree, length, skip, detrend_first);
  for (i=0;i<n;i++)
    data[i] = real[i] + _Complex_I*imag[i];

  free(real);
  free(imag);
}
#endif
Revision:	1.1
Committed:	Sun Apr 28 07:46:58 2013 UTC (10 years, 4 months ago) by tplarson
Content type:	text/plain
Branch:	MAIN
CVS Tags:	globalhs_version_5, Ver_8-7, Ver_8-5, globalhs_version_23, globalhs_version_22, globalhs_version_21, globalhs_version_20, Ver_LATEST, globalhs_version_24, Ver_8-3, globalhs_version_8, globalhs_version_9, globalhs_version_0, globalhs_version_1, globalhs_version_2, globalhs_version_3, globalhs_version_4, Ver_9-41, globalhs_version_6, globalhs_version_7, Ver_9-5, Ver_8-8, globalhs_version_19, Ver_8-2, Ver_8-10, Ver_8-1, Ver_8-6, Ver_9-1, Ver_8-4, Ver_9-2, globalhs_version_12, globalhs_version_13, globalhs_version_10, globalhs_version_11, globalhs_version_16, globalhs_version_17, globalhs_version_14, globalhs_version_15, globalhs_version_18, Ver_9-4, Ver_9-3, Ver_8-11, Ver_8-12, Ver_9-0, HEAD
Log Message:	functions needed for detrending and gapfilling
#	Content
1	// $Header: $
2
3	#include <stdio.h>
4	#include <stdlib.h>
5	#include <assert.h>
6	#include <string.h>
7	#include <complex.h>
8	#include <math.h>
9	#include "ctypes.h"
10	#include "detrend_C99.h"
11	//#include "xmem.h"
12
13	//#define DEBUGFILE //uncomment this to write debugging files
14
15	#ifdef FLOAT
16	#define TYPE FLOAT
17	#include "detrend_code_C99.h"
18	#undef TYPE
19	#endif
20
21	#ifdef DOUBLE
22	#define TYPE DOUBLE
23	#include "detrend_code_C99.h"
24	#undef TYPE
25	#endif
26
27	#ifdef COMPLEXFLOAT
28	static int count2=0;
29	void cdetrend_discontig( int n, _Complex float data, int isgood,
30	int degree, int length, int skip,
31	int m, int *sect, int detrend_first)
32	{
33	int i,first,last;
34	// Seperately detrend m sections of data separated by
35	// discontinuities.
36	// The first section begins on data[0],
37	// the m-1 first sections end on data[sect_last[i]], i=0,1,...,m-2,
38	// and the last section ends on data[n-1].
39
40	if (m==0)
41	cdetrend(n,data,isgood,degree,length,skip,detrend_first);
42	else
43	{
44	for (i=0; i<m; i++)
45	{
46	first = sect[2*i];
47	last = sect[2*i+1];
48	#ifdef DEBUGOUT
49	printf("Detrending [%d:%d]\n",first,last);
50	#endif
51	cdetrend(last-first+1,&data[first],&isgood[first],degree,length,skip,detrend_first);
52	}
53	}
54	}
55
56	void cdetrend( int n, _Complex float data, int isgood, int degree,
57	int length, int skip, int detrend_first)
58	{
59	int i;
60
61	float real, imag;
62	real = (float ) malloc(nsizeof(float));
63	for (i=0;i<n;i++)
64	real[i] = crealf(data[i]);
65	sdetrend( n, real, isgood, degree, length, skip, detrend_first);
66
67	imag = (float ) malloc(nsizeof(float));
68	for (i=0;i<n;i++)
69	imag[i] = cimagf(data[i]);
70	sdetrend( n, imag, isgood, degree, length, skip, detrend_first);
71	for (i=0;i<n;i++)
72	data[i] = real[i] + _Complex_I*imag[i];
73
74	#ifdef DEBUGFILE
75	{
76	FILE *fh;
77	char name[1024];
78	printf("writing debug2 file #%d\n",count2);
79	sprintf(name,"detrend_debug2_%d.out",count2++);
80	fh = fopen(name,"w");
81	assert(fh);
82	for (i=0; i<n; i++)
83	fprintf(fh,"%e %e\n",crealf(data[i]),cimagf(data[i]));
84	fclose(fh);
85	}
86	#endif
87
88	free(real);
89	free(imag);
90	}
91	#endif
92
93
94	#ifdef COMPLEXDOUBLE
95	void zdetrend_discontig( int n, _Complex double data, int isgood,
96	int degree, int length, int skip,
97	int m, int *sect, int detrend_first)
98	{
99	int i,first,last;
100	// Seperately detrend m sections of data separated by
101	// discontinuities.
102	// The first section begins on data[0],
103	// the m-1 first sections end on data[sect_last[i]], i=0,1,...,m-2,
104	// and the last section ends on data[n-1].
105
106	if (m==0)
107	zdetrend(n,data,isgood,degree,length,skip,detrend_first);
108	else
109	{
110	first = 0;
111	for (i=0; i<m; i++)
112	{
113	first = sect[2*i];
114	last = sect[2*i+1];
115	zdetrend(last-first+1,&data[first],&isgood[first],degree,length,skip,detrend_first);
116	}
117	}
118	}
119
120
121	void zdetrend( int n, _Complex double data, int isgood, int degree,
122	int length, int skip, int detrend_first)
123	{
124	int i;
125	double real, imag;
126
127	real = (double ) malloc(nsizeof(double));
128	for (i=0;i<n;i++)
129	real[i] = creal(data[i]);
130	ddetrend( n, real, isgood, degree, length, skip, detrend_first);
131
132	imag = (double ) malloc(nsizeof(double));
133	for (i=0;i<n;i++)
134	imag[i] = cimag(data[i]);
135	ddetrend( n, imag, isgood, degree, length, skip, detrend_first);
136	for (i=0;i<n;i++)
137	data[i] = real[i] + _Complex_I*imag[i];
138
139	free(real);
140	free(imag);
141	}
142	#endif