static int
write_8(FILE *out, uint8_t val)
{
return write_data(out, &val, sizeof(val));
}
/* Adapted from libavformat/spdifenc.c:
* It seems Dolby TrueHD frames have to be encapsulated in MAT frames before
* they can be encapsulated in IEC 61937.
* Here we encapsulate 24 TrueHD frames in a single MAT frame, padding them
* to achieve constant rate.
* The actual format of a MAT frame is unknown, but the below seems to work.
* However, it seems it is not actually necessary for the 24 TrueHD frames to
* be in an exact alignment with the MAT frame
*/
static int write_buffer_truehd( filter_t *p_filter, block_t *p_in_buf )
{
#define TRUEHD_FRAME_OFFSET 2560
filter_sys_t *p_sys = p_filter->p_sys;
if( !p_sys->p_out_buf
&& write_init( p_filter, p_in_buf, 61440, 61440 / 16 ) )
return SPDIF_ERROR;
int i_padding = 0;
if( p_sys->truehd.i_frame_count == 0 )
{
static const char p_mat_start_code[20] = {
0x07, 0x9E, 0x00, 0x03, 0x84, 0x01, 0x01, 0x01, 0x80, 0x00,
0x56, 0xA5, 0x3B, 0xF4, 0x81, 0x83, 0x49, 0x80, 0x77, 0xE0
};
write_data( p_filter, p_mat_start_code, 20, true );
/* We need to include the S/PDIF header in the first MAT frame */
i_padding = TRUEHD_FRAME_OFFSET - p_in_buf->i_buffer - 20
- SPDIF_HEADER_SIZE;
}
else if( p_sys->truehd.i_frame_count == 11 )
{
/* The middle mat code need to be at the ((2560 * 12) - 4) offset */
i_padding = TRUEHD_FRAME_OFFSET - p_in_buf->i_buffer - 4;
}
else if( p_sys->truehd.i_frame_count == 12 )
{
static const char p_mat_middle_code[12] = {
0xC3, 0xC1, 0x42, 0x49, 0x3B, 0xFA,
0x82, 0x83, 0x49, 0x80, 0x77, 0xE0
};
write_data( p_filter, p_mat_middle_code, 12, true );
i_padding = TRUEHD_FRAME_OFFSET - p_in_buf->i_buffer - ( 12 - 4 );
}
else if( p_sys->truehd.i_frame_count == 23 )
{
static const char p_mat_end_code[16] = {
0xC3, 0xC2, 0xC0, 0xC4, 0x00, 0x00, 0x00, 0x00,
0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x97, 0x11
};
/* The end mat code need to be at the ((2560 * 24) - 24) offset */
i_padding = TRUEHD_FRAME_OFFSET - p_in_buf->i_buffer - 24;
if( i_padding < 0 || p_in_buf->i_buffer + i_padding >
p_sys->p_out_buf->i_buffer - p_sys->i_out_offset )
return SPDIF_ERROR;
write_buffer( p_filter, p_in_buf );
write_padding( p_filter, i_padding );
write_data( p_filter, p_mat_end_code, 16, true );
write_finalize( p_filter, IEC61937_TRUEHD, 1 /* in bytes */ );
p_sys->truehd.i_frame_count = 0;
return SPDIF_SUCCESS;
}
else
i_padding = TRUEHD_FRAME_OFFSET - p_in_buf->i_buffer;
if( i_padding < 0 || p_in_buf->i_buffer + i_padding >
p_sys->p_out_buf->i_buffer - p_sys->i_out_offset )
return SPDIF_ERROR;
write_buffer( p_filter, p_in_buf );
write_padding( p_filter, i_padding );
p_sys->truehd.i_frame_count++;
return SPDIF_MORE_DATA;
}
/* struct smbd_smb2_read_state destructor. Send the SMB2_READ data. */
static int smb2_sendfile_send_data(struct smbd_smb2_read_state *state)
{
struct lock_struct lock;
uint32_t in_length = state->in_length;
uint64_t in_offset = state->in_offset;
files_struct *fsp = state->fsp;
const DATA_BLOB *hdr = state->smb2req->queue_entry.sendfile_header;
NTSTATUS *pstatus = state->smb2req->queue_entry.sendfile_status;
struct smbXsrv_connection *xconn = state->smb2req->xconn;
ssize_t nread;
ssize_t ret;
int saved_errno;
nread = SMB_VFS_SENDFILE(xconn->transport.sock,
fsp,
hdr,
in_offset,
in_length);
DEBUG(10,("smb2_sendfile_send_data: SMB_VFS_SENDFILE returned %d on file %s\n",
(int)nread,
fsp_str_dbg(fsp) ));
if (nread == -1) {
saved_errno = errno;
/*
* Returning ENOSYS means no data at all was sent.
Do this as a normal read. */
if (errno == ENOSYS) {
goto normal_read;
}
if (errno == EINTR) {
/*
* Special hack for broken Linux with no working sendfile. If we
* return EINTR we sent the header but not the rest of the data.
* Fake this up by doing read/write calls.
*/
set_use_sendfile(SNUM(fsp->conn), false);
nread = fake_sendfile(xconn, fsp, in_offset, in_length);
if (nread == -1) {
saved_errno = errno;
DEBUG(0,("smb2_sendfile_send_data: fake_sendfile "
"failed for file %s (%s) for client %s. "
"Terminating\n",
fsp_str_dbg(fsp), strerror(saved_errno),
smbXsrv_connection_dbg(xconn)));
*pstatus = map_nt_error_from_unix_common(saved_errno);
return 0;
}
goto out;
}
DEBUG(0,("smb2_sendfile_send_data: sendfile failed for file "
"%s (%s) for client %s. Terminating\n",
fsp_str_dbg(fsp), strerror(saved_errno),
smbXsrv_connection_dbg(xconn)));
*pstatus = map_nt_error_from_unix_common(saved_errno);
return 0;
} else if (nread == 0) {
/*
* Some sendfile implementations return 0 to indicate
* that there was a short read, but nothing was
* actually written to the socket. In this case,
* fallback to the normal read path so the header gets
* the correct byte count.
*/
DEBUG(3, ("send_file_readX: sendfile sent zero bytes "
"falling back to the normal read: %s\n",
fsp_str_dbg(fsp)));
goto normal_read;
}
/*
* We got a short read
*/
goto out;
normal_read:
/* Send out the header. */
ret = write_data(xconn->transport.sock,
(const char *)hdr->data, hdr->length);
if (ret != hdr->length) {
saved_errno = errno;
DEBUG(0,("smb2_sendfile_send_data: write_data failed for file "
"%s (%s) for client %s. Terminating\n",
fsp_str_dbg(fsp), strerror(saved_errno),
smbXsrv_connection_dbg(xconn)));
*pstatus = map_nt_error_from_unix_common(saved_errno);
return 0;
}
nread = fake_sendfile(xconn, fsp, in_offset, in_length);
if (nread == -1) {
saved_errno = errno;
DEBUG(0,("smb2_sendfile_send_data: fake_sendfile "
"failed for file %s (%s) for client %s. "
"Terminating\n",
fsp_str_dbg(fsp), strerror(saved_errno),
smbXsrv_connection_dbg(xconn)));
*pstatus = map_nt_error_from_unix_common(saved_errno);
return 0;
}
out:
if (nread < in_length) {
ret = sendfile_short_send(xconn, fsp, nread,
hdr->length, in_length);
if (ret == -1) {
saved_errno = errno;
DEBUG(0,("%s: sendfile_short_send "
"failed for file %s (%s) for client %s. "
"Terminating\n",
__func__,
fsp_str_dbg(fsp), strerror(saved_errno),
smbXsrv_connection_dbg(xconn)));
*pstatus = map_nt_error_from_unix_common(saved_errno);
return 0;
}
}
init_strict_lock_struct(fsp,
fsp->op->global->open_persistent_id,
in_offset,
in_length,
READ_LOCK,
&lock);
SMB_VFS_STRICT_UNLOCK(fsp->conn, fsp, &lock);
*pstatus = NT_STATUS_OK;
return 0;
}
MYBOOL __WINAPI write_lpex(lprec *lp, void *userhandle, write_modeldata_func write_modeldata)
{
int i, j, b,
nrows = lp->rows,
ncols = lp->columns,
nchars, maxlen = LP_MAXLINELEN;
MYBOOL ok;
REAL a;
char *ptr;
if(lp->matA->is_roworder) {
report(lp, IMPORTANT, "LP_writefile: Cannot write to LP file while in row entry mode.\n");
return(FALSE);
}
if(!mat_validate(lp->matA)) {
report(lp, IMPORTANT, "LP_writefile: Could not validate the data matrix.\n");
return(FALSE);
}
/* Write name of model */
ptr = get_lp_name(lp);
if(ptr != NULL)
if(*ptr)
write_lpcomment(userhandle, write_modeldata, ptr, FALSE);
else
ptr = NULL;
/* Write the objective function */
write_lpcomment(userhandle, write_modeldata, "Objective function", (MYBOOL) (ptr != NULL));
if(is_maxim(lp))
write_data(userhandle, write_modeldata, "max: ");
else
write_data(userhandle, write_modeldata, "min: ");
write_lprow(lp, 0, userhandle, write_modeldata, maxlen);
a = get_rh(lp, 0);
if(a != 0)
write_data(userhandle, write_modeldata, " %+.12g", a);
write_data(userhandle, write_modeldata, ";\n");
/* Write constraints */
if(nrows > 0)
write_lpcomment(userhandle, write_modeldata, "Constraints", TRUE);
for(j = 1; j <= nrows; j++) {
if(((lp->names_used) && (lp->row_name[j] != NULL)) || (write_lprow(lp, j, userhandle, NULL, maxlen) == 1))
ptr = get_row_name(lp, j);
else
ptr = NULL;
if((ptr != NULL) && (*ptr))
write_data(userhandle, write_modeldata, "%s: ", ptr);
#ifndef SingleBoundedRowInLP
/* Write the ranged part of the constraint, if specified */
if ((lp->orig_upbo[j]) && (lp->orig_upbo[j] < lp->infinite)) {
if(my_chsign(is_chsign(lp, j), lp->orig_rhs[j]) == -lp->infinite)
write_data(userhandle, write_modeldata, "-Inf %s ", (is_chsign(lp, j)) ? ">=" : "<=");
else if(my_chsign(is_chsign(lp, j), lp->orig_rhs[j]) == lp->infinite)
write_data(userhandle, write_modeldata, "+Inf %s ", (is_chsign(lp, j)) ? ">=" : "<=");
else
write_data(userhandle, write_modeldata, "%+.12g %s ",
(lp->orig_upbo[j]-lp->orig_rhs[j]) * (is_chsign(lp, j) ? 1.0 : -1.0) / (lp->scaling_used ? lp->scalars[j] : 1.0),
(is_chsign(lp, j)) ? ">=" : "<=");
}
#endif
if((!write_lprow(lp, j, userhandle, write_modeldata, maxlen)) && (ncols >= 1))
write_data(userhandle, write_modeldata, "0 %s", get_col_name(lp, 1));
if(lp->orig_upbo[j] == 0)
write_data(userhandle, write_modeldata, " =");
else if(is_chsign(lp, j))
write_data(userhandle, write_modeldata, " >=");
else
write_data(userhandle, write_modeldata, " <=");
if(fabs(get_rh(lp, j) + lp->infinite) < 1)
write_data(userhandle, write_modeldata, " -Inf;\n");
else if(fabs(get_rh(lp, j) - lp->infinite) < 1)
write_data(userhandle, write_modeldata, " +Inf;\n");
else
write_data(userhandle, write_modeldata, " %.12g;\n", get_rh(lp, j));
#ifdef SingleBoundedRowInLP
/* Write the ranged part of the constraint, if specified */
if ((lp->orig_upbo[j]) && (lp->orig_upbo[j] < lp->infinite)) {
if(((lp->names_used) && (lp->row_name[j] != NULL)) || (write_lprow(lp, j, userhandle, NULL, maxlen) == 1))
ptr = get_row_name(lp, j);
else
ptr = NULL;
if((ptr != NULL) && (*ptr))
write_data(userhandle, write_modeldata, "%s: ", ptr);
if((!write_lprow(lp, j, userhandle, write_modeldata, maxlen)) && (get_Ncolumns(lp) >= 1))
write_data(userhandle, write_modeldata, "0 %s", get_col_name(lp, 1));
write_data(userhandle, write_modeldata, " %s %g;\n",
(is_chsign(lp, j)) ? "<=" : ">=",
(lp->orig_upbo[j]-lp->orig_rhs[j]) * (is_chsign(lp, j) ? 1.0 : -1.0) / (lp->scaling_used ? lp->scalars[j] : 1.0));
}
#endif
}
/* Write bounds on variables */
ok = FALSE;
for(i = nrows + 1; i <= lp->sum; i++)
if(!is_splitvar(lp, i - nrows)) {
if(lp->orig_lowbo[i] == lp->orig_upbo[i]) {
if(!ok) {
write_lpcomment(userhandle, write_modeldata, "Variable bounds", TRUE);
ok = TRUE;
}
write_data(userhandle, write_modeldata, "%s = %.12g;\n", get_col_name(lp, i - nrows), get_upbo(lp, i - nrows));
}
else {
#ifndef SingleBoundedRowInLP
if((lp->orig_lowbo[i] != 0) && (lp->orig_upbo[i] < lp->infinite)) {
if(!ok) {
write_lpcomment(userhandle, write_modeldata, "Variable bounds", TRUE);
ok = TRUE;
}
if(lp->orig_lowbo[i] == -lp->infinite)
write_data(userhandle, write_modeldata, "-Inf");
else
write_data(userhandle, write_modeldata, "%.12g", get_lowbo(lp, i - nrows));
write_data(userhandle, write_modeldata, " <= %s <= ", get_col_name(lp, i - nrows));
if(lp->orig_lowbo[i] == lp->infinite)
write_data(userhandle, write_modeldata, "+Inf");
else
write_data(userhandle, write_modeldata, "%.12g", get_upbo(lp, i - nrows));
write_data(userhandle, write_modeldata, ";\n");
}
else
#endif
{
if(lp->orig_lowbo[i] != 0) {
if(!ok) {
write_lpcomment(userhandle, write_modeldata, "Variable bounds", TRUE);
ok = TRUE;
}
if(lp->orig_lowbo[i] == -lp->infinite)
write_data(userhandle, write_modeldata, "%s >= -Inf;\n", get_col_name(lp, i - nrows));
else if(lp->orig_lowbo[i] == lp->infinite)
write_data(userhandle, write_modeldata, "%s >= +Inf;\n", get_col_name(lp, i - nrows));
else
write_data(userhandle, write_modeldata, "%s >= %.12g;\n",
get_col_name(lp, i - nrows), get_lowbo(lp, i - nrows));
}
if(lp->orig_upbo[i] != lp->infinite) {
if(!ok) {
write_lpcomment(userhandle, write_modeldata, "Variable bounds", TRUE);
ok = TRUE;
}
write_data(userhandle, write_modeldata, "%s <= %.12g;\n",
get_col_name(lp, i - nrows), get_upbo(lp, i - nrows));
}
}
}
}
/* Write optional integer section */
if(lp->int_vars > 0) {
write_lpcomment(userhandle, write_modeldata, "Integer definitions", TRUE);
i = 1;
while((i <= ncols) && !is_int(lp, i))
i++;
if(i <= ncols) {
nchars = write_data(userhandle, write_modeldata, "int %s", get_col_name(lp, i));
i++;
for(; i <= ncols; i++)
if((!is_splitvar(lp, i)) && (is_int(lp, i))) {
if((maxlen!= 0) && (nchars > maxlen)) {
write_data(userhandle, write_modeldata, "%s", "\n");
nchars = 0;
}
write_data(userhandle, write_modeldata, ",%s", get_col_name(lp, i));
}
write_data(userhandle, write_modeldata, ";\n");
}
}
/* Write optional SEC section */
if(lp->sc_vars > 0) {
write_lpcomment(userhandle, write_modeldata, "Semi-continuous variables", TRUE);
i = 1;
while((i <= ncols) && !is_semicont(lp, i))
i++;
if(i <= ncols) {
nchars = write_data(userhandle, write_modeldata, "sec %s", get_col_name(lp, i));
i++;
for(; i <= ncols; i++)
if((!is_splitvar(lp, i)) && (is_semicont(lp, i))) {
if((maxlen != 0) && (nchars > maxlen)) {
write_data(userhandle, write_modeldata, "%s", "\n");
nchars = 0;
}
nchars += write_data(userhandle, write_modeldata, ",%s", get_col_name(lp, i));
}
write_data(userhandle, write_modeldata, ";\n");
}
}
/* Write optional SOS section */
if(SOS_count(lp) > 0) {
SOSgroup *SOS = lp->SOS;
write_lpcomment(userhandle, write_modeldata, "SOS definitions", TRUE);
write_data(userhandle, write_modeldata, "SOS\n");
for(b = 0, i = 0; i < SOS->sos_count; b = SOS->sos_list[i]->priority, i++) {
nchars = write_data(userhandle, write_modeldata, "%s: ",
(SOS->sos_list[i]->name == NULL) ||
(*SOS->sos_list[i]->name==0) ? "SOS" : SOS->sos_list[i]->name); /* formatnumber12((double) lp->sos_list[i]->priority) */
for(a = 0.0, j = 1; j <= SOS->sos_list[i]->size; a = SOS->sos_list[i]->weights[j], j++) {
if((maxlen != 0) && (nchars > maxlen)) {
write_data(userhandle, write_modeldata, "%s", "\n");
nchars = 0;
}
if(SOS->sos_list[i]->weights[j] == ++a)
nchars += write_data(userhandle, write_modeldata, "%s%s",
(j > 1) ? "," : "",
get_col_name(lp, SOS->sos_list[i]->members[j]));
else
nchars += write_data(userhandle, write_modeldata, "%s%s:%.12g",
(j > 1) ? "," : "",
get_col_name(lp, SOS->sos_list[i]->members[j]),
SOS->sos_list[i]->weights[j]);
}
if(SOS->sos_list[i]->priority == ++b)
nchars += write_data(userhandle, write_modeldata, " <= %d;\n", SOS->sos_list[i]->type);
else
nchars += write_data(userhandle, write_modeldata, " <= %d:%d;\n", SOS->sos_list[i]->type, SOS->sos_list[i]->priority);
}
}
ok = TRUE;
return(ok);
}
static void command_list(void)
{
write_command(ILI9341_PWCTRB);
write_data(0x00);
write_data(0XC1);
write_data(0X30);
write_command(ILI9341_TIMCTRA);
write_data(0x85);
write_data(0x00);
write_data(0x78);
write_command(ILI9341_PWCTRSEQ);
write_data(0x39);
write_data(0x2C);
write_data(0x00);
write_data(0x34);
write_data(0x02);
write_command(ILI9341_PUMP);
write_data(0x20);
write_command(ILI9341_TIMCTRB);
write_data(0x00);
write_data(0x00);
write_command(ILI9341_PWCTR1);
write_data(0x23);
write_command(ILI9341_PWCTR2);
write_data(0x10);
write_command(ILI9341_VMCTR1);
write_data(0x3e);
write_data(0x28);
write_command(ILI9341_VMCTR2);
write_data(0x86);
write_command(ILI9341_MADCTL);
write_data(0x48);
write_command(ILI9341_PIXFMT);
write_data(0x55);
write_command(ILI9341_FRMCTR1);
write_data(0x00);
write_data(0x18);
write_command(ILI9341_DFUNCTR);
write_data(0x08);
write_data(0x82);
write_data(0x27);
write_command(ILI9341_ENGMCTR);
write_data(0x00);
write_command(ILI9341_GAMMASET);
write_data(0x01);
write_command(ILI9341_GMCTRP1);
write_data(0x0F);
write_data(0x31);
write_data(0x2B);
write_data(0x0C);
write_data(0x0E);
write_data(0x08);
write_data(0x4E);
write_data(0xF1);
write_data(0x37);
write_data(0x07);
write_data(0x10);
write_data(0x03);
write_data(0x0E);
write_data(0x09);
write_data(0x00);
write_command(ILI9341_GMCTRN1);
write_data(0x00);
write_data(0x0E);
write_data(0x14);
write_data(0x03);
write_data(0x11);
write_data(0x07);
write_data(0x31);
write_data(0xC1);
write_data(0x48);
write_data(0x08);
write_data(0x0F);
write_data(0x0C);
write_data(0x31);
write_data(0x36);
write_data(0x0F);
write_command(ILI9341_SLPOUT);
nrf_delay_ms(120);
write_command(ILI9341_DISPON);
}
void
do_work(unsigned int start)
{
unsigned int stack_array[STACK_ARRAY_ELEMS];
unsigned int i = 0;
unsigned int array1_start = start;
unsigned int array2_start = start + (ELEMS / (ELEMS_PER_PAGE)) + 10;
printf("Checking uninitialized array1\n");
/* check the uninitialized array1 before initialization */
for (i=0; i<ELEMS; i++) {
if (array1[i] != 0) {
printf("FAILED in file %s at line %d: array1[%d] = %u != %d\n", __FILE__, __LINE__, i, array1[i], 0);
exit(1);
}
}
printf("Checking uninitialized array2\n");
/* check the uninitialized array2 before initialization */
for (i=0; i<ELEMS; i++) {
if (array2[i] != 0) {
printf("FAILED in file %s at line %d: array2[%d] = %u != %d\n", __FILE__, __LINE__, i, array2[i], 0);
exit(1);
}
}
for (i=0; i<STACK_ARRAY_ELEMS; i++) {
stack_array[i] = i * 1000;
}
for (i=0; i<2; i++) {
call_all();
write_data(array1, array1_start);
call_all();
printf("Checking initialized array1\n");
read_data(array1, array1_start, "array1");
}
/* check the uninitialized array2 again before initialization */
printf("Checking initialized array2 again\n");
for (i=0; i<ELEMS; i++) {
if (array2[i] != 0) {
printf("FAILED in file %s at line %d: array2[%d] = %u != %d\n", __FILE__, __LINE__, i, array2[i], 0);
exit(1);
}
}
printf("Checking initialized stack_array\n");
for (i=0; i<STACK_ARRAY_ELEMS; i++) {
if (stack_array[i] != i * 1000) {
printf("FAILED in file %s at line %d: stack_array[%d] = %u != %d\n", __FILE__, __LINE__, i, stack_array[i], i);
exit(1);
}
}
printf("Checking initialized init\n");
/* check the initialized array */
for (i=0; i<INIT_ARRAY_ELEMS; i++) {
if (init[i] != i) {
printf("FAILED in file %s at line %d: init[%d] = %u != %d\n", __FILE__, __LINE__, i, init[i], i);
exit(1);
}
}
for (i=0; i<2; i++) {
call_all();
write_data(array2, array2_start);
call_all();
printf("Checking initialized array2\n");
read_data(array2, array2_start, "array2");
}
printf("Checking initialized stack_array\n");
for (i=0; i<STACK_ARRAY_ELEMS; i++) {
if (stack_array[i] != i * 1000) {
printf("FAILED in file %s at line %d: stack_array[%d] = %u != %d\n", __FILE__, __LINE__, i, stack_array[i], i);
exit(1);
}
}
/* check the initialized array */
printf("Checking initialized init\n");
for (i=0; i<INIT_ARRAY_ELEMS; i++) {
if (init[i] != i) {
printf("FAILED in file %s at line %d: init[%d] = %u != %d\n", __FILE__, __LINE__, i, init[i], i);
exit(1);
}
}
printf("Checking initialized array1 for the last time\n");
read_data(array1, array1_start, "array1");
printf("Checking initialized array2 for the last time\n");
read_data(array2, array2_start, "array2");
printf("SUCCEEDED\n");
}
/* This implements the state machine defined in the IPMI manual, see
that for details on how this works. Divide that flowchart into
sections delimited by "Wait for IBF" and this will become clear. */
enum kcs_result kcs_event(struct kcs_data *kcs, long time)
{
unsigned char status;
unsigned char state;
status = read_status(kcs);
#ifdef DEBUG_STATE
printk(" State = %d, %x\n", kcs->state, status);
#endif
/* All states wait for ibf, so just do it here. */
if (!check_ibf(kcs, status, time))
return KCS_CALL_WITH_DELAY;
/* Just about everything looks at the KCS state, so grab that, too. */
state = GET_STATUS_STATE(status);
switch (kcs->state) {
case KCS_IDLE:
/* If there's and interrupt source, turn it off. */
clear_obf(kcs, status);
if (GET_STATUS_ATN(status))
return KCS_ATTN;
else
return KCS_SM_IDLE;
case KCS_START_OP:
if (state != KCS_IDLE) {
start_error_recovery(kcs,
"State machine not idle at start");
break;
}
clear_obf(kcs, status);
write_cmd(kcs, KCS_WRITE_START);
kcs->state = KCS_WAIT_WRITE_START;
break;
case KCS_WAIT_WRITE_START:
if (state != KCS_WRITE_STATE) {
start_error_recovery(
kcs,
"Not in write state at write start");
break;
}
read_data(kcs);
if (kcs->write_count == 1) {
write_cmd(kcs, KCS_WRITE_END);
kcs->state = KCS_WAIT_WRITE_END;
} else {
write_next_byte(kcs);
kcs->state = KCS_WAIT_WRITE;
}
break;
case KCS_WAIT_WRITE:
if (state != KCS_WRITE_STATE) {
start_error_recovery(kcs,
"Not in write state for write");
break;
}
clear_obf(kcs, status);
if (kcs->write_count == 1) {
write_cmd(kcs, KCS_WRITE_END);
kcs->state = KCS_WAIT_WRITE_END;
} else {
write_next_byte(kcs);
}
break;
case KCS_WAIT_WRITE_END:
if (state != KCS_WRITE_STATE) {
start_error_recovery(kcs,
"Not in write state for write end");
break;
}
clear_obf(kcs, status);
write_next_byte(kcs);
kcs->state = KCS_WAIT_READ;
break;
case KCS_WAIT_READ:
if ((state != KCS_READ_STATE) && (state != KCS_IDLE_STATE)) {
start_error_recovery(
kcs,
"Not in read or idle in read state");
break;
}
if (state == KCS_READ_STATE) {
if (! check_obf(kcs, status, time))
return KCS_CALL_WITH_DELAY;
read_next_byte(kcs);
} else {
/* We don't implement this exactly like the state
machine in the spec. Some broken hardware
does not write the final dummy byte to the
read register. Thus obf will never go high
here. We just go straight to idle, and we
handle clearing out obf in idle state if it
happens to come in. */
clear_obf(kcs, status);
kcs->orig_write_count = 0;
kcs->state = KCS_IDLE;
return KCS_TRANSACTION_COMPLETE;
}
break;
case KCS_ERROR0:
clear_obf(kcs, status);
write_cmd(kcs, KCS_GET_STATUS_ABORT);
kcs->state = KCS_ERROR1;
break;
case KCS_ERROR1:
clear_obf(kcs, status);
write_data(kcs, 0);
kcs->state = KCS_ERROR2;
break;
case KCS_ERROR2:
if (state != KCS_READ_STATE) {
start_error_recovery(kcs,
"Not in read state for error2");
break;
}
if (! check_obf(kcs, status, time))
return KCS_CALL_WITH_DELAY;
clear_obf(kcs, status);
write_data(kcs, KCS_READ_BYTE);
kcs->state = KCS_ERROR3;
break;
case KCS_ERROR3:
if (state != KCS_IDLE_STATE) {
start_error_recovery(kcs,
"Not in idle state for error3");
break;
}
if (! check_obf(kcs, status, time))
return KCS_CALL_WITH_DELAY;
clear_obf(kcs, status);
if (kcs->orig_write_count) {
restart_kcs_transaction(kcs);
} else {
kcs->state = KCS_IDLE;
return KCS_TRANSACTION_COMPLETE;
}
break;
case KCS_HOSED:
break;
}
if (kcs->state == KCS_HOSED) {
init_kcs_data(kcs, kcs->port, kcs->addr);
return KCS_SM_HOSED;
}
return KCS_CALL_WITHOUT_DELAY;
}
/*
* move the file pointer to the point where the next key can be written and return
* that offset
*/
int
write_key_init(int fd, UINT32 pub_data_size, UINT32 blob_size, UINT32 vendor_data_size)
{
UINT32 num_keys;
BYTE version;
int rc, offset;
/* seek to the PS version */
rc = lseek(fd, TSSPS_VERSION_OFFSET, SEEK_SET);
if (rc == ((off_t) - 1)) {
LogError("lseek: %s", strerror(errno));
return -1;
}
/* go to NUM_KEYS */
rc = lseek(fd, TSSPS_NUM_KEYS_OFFSET, SEEK_SET);
if (rc == ((off_t) - 1)) {
LogError("lseek: %s", strerror(errno));
return -1;
}
/* read the number of keys */
rc = read(fd, &num_keys, sizeof(UINT32));
num_keys = LE_32(num_keys);
if (rc == -1) {
LogError("read of %zd bytes: %s", sizeof(UINT32), strerror(errno));
return -1;
} else if (rc == 0) {
/* This is the first key being written */
num_keys = 1;
version = 1;
/* seek to the PS version */
rc = lseek(fd, TSSPS_VERSION_OFFSET, SEEK_SET);
if (rc == ((off_t) - 1)) {
LogError("lseek: %s", strerror(errno));
return -1;
}
/* write out the version info byte */
if ((rc = write_data(fd, &version, sizeof(BYTE)))) {
LogError("%s", __FUNCTION__);
return rc;
}
rc = lseek(fd, TSSPS_NUM_KEYS_OFFSET, SEEK_SET);
if (rc == ((off_t) - 1)) {
LogError("lseek: %s", strerror(errno));
return -1;
}
num_keys = LE_32(num_keys);
if ((rc = write_data(fd, &num_keys, sizeof(UINT32)))) {
LogError("%s", __FUNCTION__);
return rc;
}
/* return the offset */
return (TSSPS_NUM_KEYS_OFFSET + sizeof(UINT32));
}
/* if there is a hole in the file we can write to, find it */
offset = find_write_offset(pub_data_size, blob_size, vendor_data_size);
if (offset != -1) {
/* we found a hole, seek to it and don't increment the # of keys on disk */
rc = lseek(fd, offset, SEEK_SET);
} else {
/* we didn't find a hole, increment the number of keys on disk and seek
* to the end of the file
*/
num_keys++;
/* go to the beginning */
rc = lseek(fd, TSSPS_NUM_KEYS_OFFSET, SEEK_SET);
if (rc == ((off_t) - 1)) {
LogError("lseek: %s", strerror(errno));
return -1;
}
num_keys = LE_32(num_keys);
if ((rc = write_data(fd, &num_keys, sizeof(UINT32)))) {
LogError("%s", __FUNCTION__);
return rc;
}
rc = lseek(fd, 0, SEEK_END);
}
if (rc == ((off_t) - 1)) {
LogError("lseek: %s", strerror(errno));
return -1;
}
/* lseek returns the number of bytes of offset from the beginning of the file */
return rc;
}
void lcd_initial()
{
int i;
// hard reset
for (i=0; i<120; i++)
SPI_LCD_REST_LOW();
for (i=0; i<120; i++)
SPI_LCD_REST_HIGH();
//-------------------Software Reset-------------------------------//
write_command(0x11);//Sleep exit
//ST7735R Frame Rate
write_command(0xB1);
write_data(0x01); write_data(0x2C); write_data(0x2D);
write_command(0xB2);
write_data(0x01); write_data(0x2C); write_data(0x2D);
write_command(0xB3);
write_data(0x01); write_data(0x2C); write_data(0x2D);
write_data(0x01); write_data(0x2C); write_data(0x2D);
write_command(0xB4); //Column inversion
write_data(0x07);
//ST7735R Power Sequence
write_command(0xC0);
write_data(0xA2); write_data(0x02); write_data(0x84);
write_command(0xC1); write_data(0xC5);
write_command(0xC2);
write_data(0x0A); write_data(0x00);
write_command(0xC3);
write_data(0x8A); write_data(0x2A);
write_command(0xC4);
write_data(0x8A); write_data(0xEE);
write_command(0xC5); //VCOM
write_data(0x0E);
write_command(0x36); //MX, MY, RGB mode
#if ROTATION == 90
write_data(0xa0 | SBIT_RGB);
#elif ROTATION == 180
write_data(0xc0 | SBIT_RGB);
#elif ROTATION == 270
write_data(0x60 | SBIT_RGB);
#else
write_data(0x00 | SBIT_RGB);
#endif
#if ROTATION == 90 || ROTATION == 270
write_command(0x2a);
write_data(0x00);write_data(0x00);
write_data(0x00);write_data(0x9f);
write_command(0x2b);
write_data(0x00);write_data(0x00);
write_data(0x00);write_data(0x7f);
#else
write_command(0x2a);
write_data(0x00);write_data(0x00);
write_data(0x00);write_data(0x7f);
write_command(0x2b);
write_data(0x00);write_data(0x00);
write_data(0x00);write_data(0x9f);
#endif
//ST7735R Gamma Sequence
write_command(0xe0);
write_data(0x0f); write_data(0x1a);
write_data(0x0f); write_data(0x18);
write_data(0x2f); write_data(0x28);
write_data(0x20); write_data(0x22);
write_data(0x1f); write_data(0x1b);
write_data(0x23); write_data(0x37);
write_data(0x00); write_data(0x07);
write_data(0x02); write_data(0x10);
write_command(0xe1);
write_data(0x0f); write_data(0x1b);
write_data(0x0f); write_data(0x17);
write_data(0x33); write_data(0x2c);
write_data(0x29); write_data(0x2e);
write_data(0x30); write_data(0x30);
write_data(0x39); write_data(0x3f);
write_data(0x00); write_data(0x07);
write_data(0x03); write_data(0x10);
write_command(0xF0); //Enable test command
write_data(0x01);
write_command(0xF6); //Disable ram power save mode
write_data(0x00);
write_command(0x3A); //65k mode
write_data(0x05);
write_command(0x29);//Display on
}
void write_item(uint8_t *buf)
{
write_data(buf+1, buf[0]);
}
void isapnp_write_byte(unsigned char idx, unsigned char val)
{
write_address(idx);
write_data(val);
}
void RequestMonitoredPages(void)
{
DIR *dir;
struct dirent* ent;
/* Open the monitor subdirectory. */
if(chdir("monitor"))
{PrintMessage(Warning,"Cannot change to directory 'monitor'; [%!s] no files monitored.");return;}
dir=opendir(".");
if(!dir)
{PrintMessage(Warning,"Cannot open directory 'monitor'; [%!s] no files monitored.");ChangeBackToSpoolDir();return;}
ent=readdir(dir);
if(!ent)
{PrintMessage(Warning,"Cannot read directory 'monitor'; [%!s] no files monitored.");closedir(dir);ChangeBackToSpoolDir();return;}
/* Scan through all of the files. */
do
{
struct stat buf;
if(ent->d_name[0]=='.' && (ent->d_name[1]==0 || (ent->d_name[1]=='.' && ent->d_name[2]==0)))
continue; /* skip . & .. */
if(stat(ent->d_name,&buf))
{PrintMessage(Inform,"Cannot stat file 'monitor/%s'; [%!s] race condition?",ent->d_name);return;}
else if(S_ISREG(buf.st_mode) && *ent->d_name=='O')
{
URL *Url=FileNameToURL(ent->d_name);
int last,next;
ChangeBackToSpoolDir();
if(!Url)
continue;
MonitorTimes(Url,&last,&next);
PrintMessage(Debug,"Monitoring '%s' last=%dh next=%dh => %s",Url->name,last,next,next?"No":"Yes");
chdir("monitor");
if(next==0)
{
int ifd=open(ent->d_name,O_RDONLY|O_BINARY);
if(ifd==-1)
PrintMessage(Warning,"Cannot open monitored file 'monitor/%s' to read; [%!s].",ent->d_name);
else
{
int ofd;
init_io(ifd);
ChangeBackToSpoolDir();
ofd=OpenNewOutgoingSpoolFile();
if(ofd==-1)
PrintMessage(Warning,"Cannot open outgoing spool file for monitored URL '%s'; [%!s].",Url->name);
else
{
char *contents=(char*)malloc(buf.st_size+1);
init_io(ofd);
read_data(ifd,contents,buf.st_size);
if(write_data(ofd,contents,buf.st_size)==-1)
PrintMessage(Warning,"Cannot write to outgoing file; disk full?");
finish_io(ofd);
CloseNewOutgoingSpoolFile(ofd,Url);
free(contents);
}
chdir("monitor");
finish_io(ifd);
close(ifd);
if(utime(URLToFileName(Url,'M',0),NULL))
PrintMessage(Warning,"Cannot change timestamp of monitored file 'monitor/%s'; [%!s].",URLToFileName(Url,'M',0));
}
}
FreeURL(Url);
}
}
while((ent=readdir(dir)));
ChangeBackToSpoolDir();
closedir(dir);
}
int main(int argc, char *argv[]){
DWORD status = 0;
DWORD errorstatus = 0;
string error_str;
scu_mil myscu;
/* int test = 0;
BYTE blub = 0x38;
cin>>hex>>test;
printf("test 0x%x \n", (unsigned char)test);
cout<<"HEX :"<<hex<<test<<endl;
cout<<"DEC :"<<dec<<test<<endl;
cout<<"BYTE HEX :"<<hex<<char(test)<<endl;
cout<<"BYTE DEC :"<<dec<<BYTE(test)<<endl;
//blub = test;
cout<<"HEX blub: "<<hex<<blub<<endl;*/
if(argc == 1){
cout<<"No target device defined !"<<endl;
cout<<"example : scumil tcp/scuxl0089.acc"<<endl;
cout<<"end"<<endl;
return 0;
}
int choice = 0;
bool run = true;
// generator starten
srand((unsigned)time(NULL));
// status = myscu.scu_milbusopen("tcp/scuxl0089.acc", errorstatus);
status = myscu.scu_milbusopen(argv[1], errorstatus);
error_str = myscu.scu_milerror(status);
cout <<"open scu : "<<error_str<<endl;
if (status != status_ok){
return 0;
}
// -------------------------
do{
choice = print_menue();
switch (choice)
{
case 1:
scan_milbus_loop(myscu);
break;
case 2:
write_data(myscu);
break;
case 3:
write_cmd(myscu);
break;
case 4:
read_data(myscu);
break;
case 5:
write_ifk(myscu);
break;
case 6:
read_ifk(myscu);
break;
case 7:
cout<<endl<<"IFKs online:"<<endl;
cout<<"----------------"<<endl;
echo_test(myscu,choose_ifk(myscu));
break;
default : run = false;
break;
}
}while (run);
// -------------------------
status = myscu.scu_milbusclose(errorstatus);
error_str = myscu.scu_milerror(status);
cout <<"close scu : "<<error_str<<endl;
return 0;
}
int main(int argc, char *argv[])
{
/* Any variable that begins with 't' means topocentric */
/* Any variable that begins with 'b' means barycentric */
FILE **outfiles = NULL;
float **outdata;
double dtmp, *dms, avgdm = 0.0, dsdt = 0, maxdm;
double *dispdt, tlotoa = 0.0, blotoa = 0.0, BW_ddelay = 0.0;
double max = -9.9E30, min = 9.9E30, var = 0.0, avg = 0.0;
double *btoa = NULL, *ttoa = NULL, avgvoverc = 0.0;
char obs[3], ephem[10], rastring[50], decstring[50];
int totnumtowrite, *idispdt, **offsets;
int ii, jj, numadded = 0, numremoved = 0, padding = 0, good_inputs = 1;
int numbarypts = 0, numread = 0, numtowrite = 0, totwrote = 0, datawrote = 0;
int padwrote = 0, padtowrite = 0, statnum = 0;
int numdiffbins = 0, *diffbins = NULL, *diffbinptr = NULL, good_padvals = 0;
double local_lodm;
char *datafilenm, *outpath, *outfilenm, *hostname;
struct spectra_info s;
infodata idata;
mask obsmask;
MPI_Init(&argc, &argv);
MPI_Comm_size(MPI_COMM_WORLD, &numprocs);
MPI_Comm_rank(MPI_COMM_WORLD, &myid);
set_using_MPI();
{
FILE *hostfile;
char tmpname[100];
int retval;
hostfile = chkfopen("/etc/hostname", "r");
retval = fscanf(hostfile, "%s\n", tmpname);
if (retval==0) {
printf("Warning: error reading /etc/hostname on proc %d\n", myid);
}
hostname = (char *) calloc(strlen(tmpname) + 1, 1);
memcpy(hostname, tmpname, strlen(tmpname));
fclose(hostfile);
}
/* Call usage() if we have no command line arguments */
if (argc == 1) {
if (myid == 0) {
Program = argv[0];
usage();
}
MPI_Finalize();
exit(1);
}
make_maskbase_struct();
make_spectra_info_struct();
/* Parse the command line using the excellent program Clig */
cmd = parseCmdline(argc, argv);
spectra_info_set_defaults(&s);
// If we are zeroDMing, make sure that clipping is off.
if (cmd->zerodmP) cmd->noclipP = 1;
s.clip_sigma = cmd->clip;
if (cmd->noclipP) {
cmd->clip = 0.0;
s.clip_sigma = 0.0;
}
if (cmd->ifsP) {
// 0 = default or summed, 1-4 are possible also
s.use_poln = cmd->ifs + 1;
}
if (!cmd->numoutP)
cmd->numout = INT_MAX;
#ifdef DEBUG
showOptionValues();
#endif
if (myid == 0) { /* Master node only */
printf("\n\n");
printf(" Parallel Pulsar Subband De-dispersion Routine\n");
printf(" by Scott M. Ransom\n\n");
s.filenames = cmd->argv;
s.num_files = cmd->argc;
s.clip_sigma = cmd->clip;
// -1 causes the data to determine if we use weights, scales, &
// offsets for PSRFITS or flip the band for any data type where
// we can figure that out with the data
s.apply_flipband = (cmd->invertP) ? 1 : -1;
s.apply_weight = (cmd->noweightsP) ? 0 : -1;
s.apply_scale = (cmd->noscalesP) ? 0 : -1;
s.apply_offset = (cmd->nooffsetsP) ? 0 : -1;
s.remove_zerodm = (cmd->zerodmP) ? 1 : 0;
if (RAWDATA) {
if (cmd->filterbankP) s.datatype = SIGPROCFB;
else if (cmd->psrfitsP) s.datatype = PSRFITS;
else if (cmd->pkmbP) s.datatype = SCAMP;
else if (cmd->bcpmP) s.datatype = BPP;
else if (cmd->wappP) s.datatype = WAPP;
else if (cmd->spigotP) s.datatype = SPIGOT;
} else { // Attempt to auto-identify the data
identify_psrdatatype(&s, 1);
if (s.datatype==SIGPROCFB) cmd->filterbankP = 1;
else if (s.datatype==PSRFITS) cmd->psrfitsP = 1;
else if (s.datatype==SCAMP) cmd->pkmbP = 1;
else if (s.datatype==BPP) cmd->bcpmP = 1;
else if (s.datatype==WAPP) cmd->wappP = 1;
else if (s.datatype==SPIGOT) cmd->spigotP = 1;
else if (s.datatype==SUBBAND) insubs = 1;
else {
printf("\nError: Unable to identify input data files. Please specify type.\n\n");
good_inputs = 0;
}
}
// So far we can only handle PSRFITS, filterbank, and subbands
if (s.datatype!=PSRFITS &&
s.datatype!=SIGPROCFB &&
s.datatype!=SUBBAND) good_inputs = 0;
// For subbanded data
if (!RAWDATA) s.files = (FILE **)malloc(sizeof(FILE *) * s.num_files);
if (good_inputs && (RAWDATA || insubs)) {
char description[40];
psrdatatype_description(description, s.datatype);
if (s.num_files > 1)
printf("Reading %s data from %d files:\n", description, s.num_files);
else
printf("Reading %s data from 1 file:\n", description);
for (ii = 0; ii < s.num_files; ii++) {
printf(" '%s'\n", cmd->argv[ii]);
if (insubs) s.files[ii] = chkfopen(s.filenames[ii], "rb");
}
printf("\n");
if (RAWDATA) {
read_rawdata_files(&s);
print_spectra_info_summary(&s);
spectra_info_to_inf(&s, &idata);
} else { // insubs
char *root, *suffix;
cmd->nsub = s.num_files;
s.N = chkfilelen(s.files[0], sizeof(short));
s.start_subint = gen_ivect(1);
s.num_subint = gen_ivect(1);
s.start_MJD = (long double *)malloc(sizeof(long double));
s.start_spec = (long long *)malloc(sizeof(long long));
s.num_spec = (long long *)malloc(sizeof(long long));
s.num_pad = (long long *)malloc(sizeof(long long));
s.start_spec[0] = 0L;
s.start_subint[0] = 0;
s.num_spec[0] = s.N;
s.num_subint[0] = s.N / SUBSBLOCKLEN;
s.num_pad[0] = 0L;
s.padvals = gen_fvect(s.num_files);
for (ii = 0 ; ii < ii ; ii++)
s.padvals[ii] = 0.0;
if (split_root_suffix(s.filenames[0], &root, &suffix) == 0) {
printf("\nError: The input filename (%s) must have a suffix!\n\n", s.filenames[0]);
exit(1);
}
if (strncmp(suffix, "sub", 3) == 0) {
char *tmpname;
tmpname = calloc(strlen(root) + 10, 1);
sprintf(tmpname, "%s.sub", root);
readinf(&idata, tmpname);
free(tmpname);
strncpy(s.telescope, idata.telescope, 40);
strncpy(s.backend, idata.instrument, 40);
strncpy(s.observer, idata.observer, 40);
strncpy(s.source, idata.object, 40);
s.ra2000 = hms2rad(idata.ra_h, idata.ra_m,
idata.ra_s) * RADTODEG;
s.dec2000 = dms2rad(idata.dec_d, idata.dec_m,
idata.dec_s) * RADTODEG;
ra_dec_to_string(s.ra_str,
idata.ra_h, idata.ra_m, idata.ra_s);
ra_dec_to_string(s.dec_str,
idata.dec_d, idata.dec_m, idata.dec_s);
s.num_channels = idata.num_chan;
s.start_MJD[0] = idata.mjd_i + idata.mjd_f;
s.dt = idata.dt;
s.T = s.N * s.dt;
s.lo_freq = idata.freq;
s.df = idata.chan_wid;
s.hi_freq = s.lo_freq + (s.num_channels - 1.0) * s.df;
s.BW = s.num_channels * s.df;
s.fctr = s.lo_freq - 0.5 * s.df + 0.5 * s.BW;
s.beam_FWHM = idata.fov / 3600.0;
s.spectra_per_subint = SUBSBLOCKLEN;
print_spectra_info_summary(&s);
} else {
printf("\nThe input files (%s) must be subbands! (i.e. *.sub##)\n\n",
cmd->argv[0]);
MPI_Finalize();
exit(1);
}
free(root);
free(suffix);
}
}
}
// If we don't have good input data, exit
MPI_Bcast(&good_inputs, 1, MPI_INT, 0, MPI_COMM_WORLD);
if (!good_inputs) {
MPI_Finalize();
exit(1);
}
MPI_Bcast(&insubs, 1, MPI_INT, 0, MPI_COMM_WORLD);
if (insubs)
cmd->nsub = cmd->argc;
/* Determine the output file names and open them */
local_numdms = cmd->numdms / (numprocs - 1);
dms = gen_dvect(local_numdms);
if (cmd->numdms % (numprocs - 1)) {
if (myid == 0)
printf
("\nThe number of DMs must be divisible by (the number of processors - 1).\n\n");
MPI_Finalize();
exit(1);
}
local_lodm = cmd->lodm + (myid - 1) * local_numdms * cmd->dmstep;
split_path_file(cmd->outfile, &outpath, &outfilenm);
datafilenm = (char *) calloc(strlen(outfilenm) + 20, 1);
if (myid > 0) {
if (chdir(outpath) == -1) {
printf("\nProcess %d on %s cannot chdir() to '%s'. Exiting.\n\n",
myid, hostname, outpath);
MPI_Finalize();
exit(1);
}
outfiles = (FILE **) malloc(local_numdms * sizeof(FILE *));
for (ii = 0; ii < local_numdms; ii++) {
dms[ii] = local_lodm + ii * cmd->dmstep;
avgdm += dms[ii];
sprintf(datafilenm, "%s_DM%.2f.dat", outfilenm, dms[ii]);
outfiles[ii] = chkfopen(datafilenm, "wb");
}
avgdm /= local_numdms;
}
// Broadcast the raw data information
broadcast_spectra_info(&s, myid);
if (myid > 0) {
spectra_info_to_inf(&s, &idata);
if (s.datatype==SIGPROCFB) cmd->filterbankP = 1;
else if (s.datatype==PSRFITS) cmd->psrfitsP = 1;
else if (s.datatype==SCAMP) cmd->pkmbP = 1;
else if (s.datatype==BPP) cmd->bcpmP = 1;
else if (s.datatype==WAPP) cmd->wappP = 1;
else if (s.datatype==SPIGOT) cmd->spigotP = 1;
else if (s.datatype==SUBBAND) insubs = 1;
}
s.filenames = cmd->argv;
/* Read an input mask if wanted */
if (myid > 0) {
int numpad = s.num_channels;
if (insubs)
numpad = s.num_files;
s.padvals = gen_fvect(numpad);
for (ii = 0 ; ii < numpad ; ii++)
s.padvals[ii] = 0.0;
}
if (cmd->maskfileP) {
if (myid == 0) {
read_mask(cmd->maskfile, &obsmask);
printf("Read mask information from '%s'\n\n", cmd->maskfile);
good_padvals = determine_padvals(cmd->maskfile, &obsmask, s.padvals);
}
broadcast_mask(&obsmask, myid);
MPI_Bcast(&good_padvals, 1, MPI_INT, 0, MPI_COMM_WORLD);
MPI_Bcast(s.padvals, obsmask.numchan, MPI_FLOAT, 0, MPI_COMM_WORLD);
} else {
obsmask.numchan = obsmask.numint = 0;
MPI_Bcast(&good_padvals, 1, MPI_INT, 0, MPI_COMM_WORLD);
}
// The number of topo to bary time points to generate with TEMPO
numbarypts = (int) (s.T * 1.1 / TDT + 5.5) + 1;
// Identify the TEMPO observatory code
{
char *outscope = (char *) calloc(40, sizeof(char));
telescope_to_tempocode(idata.telescope, outscope, obs);
free(outscope);
}
// Broadcast or calculate a few extra important values
if (insubs) avgdm = idata.dm;
idata.dm = avgdm;
dsdt = cmd->downsamp * idata.dt;
maxdm = cmd->lodm + cmd->numdms * cmd->dmstep;
BW_ddelay = delay_from_dm(maxdm, idata.freq) -
delay_from_dm(maxdm, idata.freq + (idata.num_chan-1) * idata.chan_wid);
blocksperread = ((int) (BW_ddelay / idata.dt) / s.spectra_per_subint + 1);
worklen = s.spectra_per_subint * blocksperread;
if (cmd->nsub > s.num_channels) {
printf
("Warning: The number of requested subbands (%d) is larger than the number of channels (%d).\n",
cmd->nsub, s.num_channels);
printf(" Re-setting the number of subbands to %d.\n\n", s.num_channels);
cmd->nsub = s.num_channels;
}
if (s.spectra_per_subint % cmd->downsamp) {
if (myid == 0) {
printf
("\nError: The downsample factor (%d) must be a factor of the\n",
cmd->downsamp);
printf(" blocklength (%d). Exiting.\n\n", s.spectra_per_subint);
}
MPI_Finalize();
exit(1);
}
tlotoa = idata.mjd_i + idata.mjd_f; /* Topocentric epoch */
if (cmd->numoutP)
totnumtowrite = cmd->numout;
else
totnumtowrite = (int) idata.N / cmd->downsamp;
if (cmd->nobaryP) { /* Main loop if we are not barycentering... */
/* Dispersion delays (in bins). The high freq gets no delay */
/* All other delays are positive fractions of bin length (dt) */
dispdt = subband_search_delays(s.num_channels, cmd->nsub, avgdm,
idata.freq, idata.chan_wid, 0.0);
idispdt = gen_ivect(s.num_channels);
for (ii = 0; ii < s.num_channels; ii++)
idispdt[ii] = NEAREST_INT(dispdt[ii] / idata.dt);
vect_free(dispdt);
/* The subband dispersion delays (see note above) */
offsets = gen_imatrix(local_numdms, cmd->nsub);
for (ii = 0; ii < local_numdms; ii++) {
double *subdispdt;
subdispdt = subband_delays(s.num_channels, cmd->nsub, dms[ii],
idata.freq, idata.chan_wid, 0.0);
dtmp = subdispdt[cmd->nsub - 1];
for (jj = 0; jj < cmd->nsub; jj++)
offsets[ii][jj] = NEAREST_INT((subdispdt[jj] - dtmp) / dsdt);
vect_free(subdispdt);
}
/* Allocate our data array and start getting data */
if (myid == 0) {
printf("De-dispersing using %d subbands.\n", cmd->nsub);
if (cmd->downsamp > 1)
printf("Downsampling by a factor of %d (new dt = %.10g)\n",
cmd->downsamp, dsdt);
printf("\n");
}
/* Print the nodes and the DMs they are handling */
print_dms(hostname, myid, numprocs, local_numdms, dms);
outdata = gen_fmatrix(local_numdms, worklen / cmd->downsamp);
numread = get_data(outdata, blocksperread, &s,
&obsmask, idispdt, offsets, &padding);
while (numread == worklen) {
numread /= cmd->downsamp;
if (myid == 0)
print_percent_complete(totwrote, totnumtowrite);
/* Write the latest chunk of data, but don't */
/* write more than cmd->numout points. */
numtowrite = numread;
if (cmd->numoutP && (totwrote + numtowrite) > cmd->numout)
numtowrite = cmd->numout - totwrote;
if (myid > 0) {
write_data(outfiles, local_numdms, outdata, 0, numtowrite);
/* Update the statistics */
if (!padding) {
for (ii = 0; ii < numtowrite; ii++)
update_stats(statnum + ii, outdata[0][ii], &min, &max, &avg, &var);
statnum += numtowrite;
}
}
totwrote += numtowrite;
/* Stop if we have written out all the data we need to */
if (cmd->numoutP && (totwrote == cmd->numout))
break;
numread = get_data(outdata, blocksperread, &s,
&obsmask, idispdt, offsets, &padding);
}
datawrote = totwrote;
} else { /* Main loop if we are barycentering... */
/* What ephemeris will we use? (Default is DE200) */
if (cmd->de405P)
strcpy(ephem, "DE405");
else
strcpy(ephem, "DE200");
/* Define the RA and DEC of the observation */
ra_dec_to_string(rastring, idata.ra_h, idata.ra_m, idata.ra_s);
ra_dec_to_string(decstring, idata.dec_d, idata.dec_m, idata.dec_s);
/* Allocate some arrays */
btoa = gen_dvect(numbarypts);
ttoa = gen_dvect(numbarypts);
for (ii = 0; ii < numbarypts; ii++)
ttoa[ii] = tlotoa + TDT * ii / SECPERDAY;
/* Call TEMPO for the barycentering */
if (myid == 0) {
double maxvoverc = -1.0, minvoverc = 1.0, *voverc = NULL;
printf("\nGenerating barycentric corrections...\n");
voverc = gen_dvect(numbarypts);
barycenter(ttoa, btoa, voverc, numbarypts, rastring, decstring, obs, ephem);
for (ii = 0; ii < numbarypts; ii++) {
if (voverc[ii] > maxvoverc)
maxvoverc = voverc[ii];
if (voverc[ii] < minvoverc)
minvoverc = voverc[ii];
avgvoverc += voverc[ii];
}
avgvoverc /= numbarypts;
vect_free(voverc);
printf(" Average topocentric velocity (c) = %.7g\n", avgvoverc);
printf(" Maximum topocentric velocity (c) = %.7g\n", maxvoverc);
printf(" Minimum topocentric velocity (c) = %.7g\n\n", minvoverc);
printf("De-dispersing using %d subbands.\n", cmd->nsub);
if (cmd->downsamp > 1) {
printf(" Downsample = %d\n", cmd->downsamp);
printf(" New sample dt = %.10g\n", dsdt);
}
printf("\n");
}
/* Print the nodes and the DMs they are handling */
print_dms(hostname, myid, numprocs, local_numdms, dms);
MPI_Bcast(btoa, numbarypts, MPI_DOUBLE, 0, MPI_COMM_WORLD);
MPI_Bcast(&avgvoverc, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD);
blotoa = btoa[0];
/* Dispersion delays (in bins). The high freq gets no delay */
/* All other delays are positive fractions of bin length (dt) */
dispdt = subband_search_delays(s.num_channels, cmd->nsub, avgdm,
idata.freq, idata.chan_wid, avgvoverc);
idispdt = gen_ivect(s.num_channels);
for (ii = 0; ii < s.num_channels; ii++)
idispdt[ii] = NEAREST_INT(dispdt[ii] / idata.dt);
vect_free(dispdt);
/* The subband dispersion delays (see note above) */
offsets = gen_imatrix(local_numdms, cmd->nsub);
for (ii = 0; ii < local_numdms; ii++) {
double *subdispdt;
subdispdt = subband_delays(s.num_channels, cmd->nsub, dms[ii],
idata.freq, idata.chan_wid, avgvoverc);
dtmp = subdispdt[cmd->nsub - 1];
for (jj = 0; jj < cmd->nsub; jj++)
offsets[ii][jj] = NEAREST_INT((subdispdt[jj] - dtmp) / dsdt);
vect_free(subdispdt);
}
/* Convert the bary TOAs to differences from the topo TOAs in */
/* units of bin length (dt) rounded to the nearest integer. */
dtmp = (btoa[0] - ttoa[0]);
for (ii = 0; ii < numbarypts; ii++)
btoa[ii] = ((btoa[ii] - ttoa[ii]) - dtmp) * SECPERDAY / dsdt;
/* Find the points where we need to add or remove bins */
{
int oldbin = 0, currentbin;
double lobin, hibin, calcpt;
numdiffbins = abs(NEAREST_INT(btoa[numbarypts - 1])) + 1;
diffbins = gen_ivect(numdiffbins);
diffbinptr = diffbins;
for (ii = 1; ii < numbarypts; ii++) {
currentbin = NEAREST_INT(btoa[ii]);
if (currentbin != oldbin) {
if (currentbin > 0) {
calcpt = oldbin + 0.5;
lobin = (ii - 1) * TDT / dsdt;
hibin = ii * TDT / dsdt;
} else {
calcpt = oldbin - 0.5;
lobin = -((ii - 1) * TDT / dsdt);
hibin = -(ii * TDT / dsdt);
}
while (fabs(calcpt) < fabs(btoa[ii])) {
/* Negative bin number means remove that bin */
/* Positive bin number means add a bin there */
*diffbinptr =
NEAREST_INT(LININTERP
(calcpt, btoa[ii - 1], btoa[ii], lobin, hibin));
diffbinptr++;
calcpt = (currentbin > 0) ? calcpt + 1.0 : calcpt - 1.0;
}
oldbin = currentbin;
}
}
*diffbinptr = cmd->numout; /* Used as a marker */
}
diffbinptr = diffbins;
/* Now perform the barycentering */
outdata = gen_fmatrix(local_numdms, worklen / cmd->downsamp);
numread = get_data(outdata, blocksperread, &s,
&obsmask, idispdt, offsets, &padding);
while (numread == worklen) { /* Loop to read and write the data */
int numwritten = 0;
double block_avg, block_var;
numread /= cmd->downsamp;
/* Determine the approximate local average */
avg_var(outdata[0], numread, &block_avg, &block_var);
if (myid == 0)
print_percent_complete(totwrote, totnumtowrite);
/* Simply write the data if we don't have to add or */
/* remove any bins from this batch. */
/* OR write the amount of data up to cmd->numout or */
/* the next bin that will be added or removed. */
numtowrite = abs(*diffbinptr) - datawrote;
if (cmd->numoutP && (totwrote + numtowrite) > cmd->numout)
numtowrite = cmd->numout - totwrote;
if (numtowrite > numread)
numtowrite = numread;
if (myid > 0) {
write_data(outfiles, local_numdms, outdata, 0, numtowrite);
/* Update the statistics */
if (!padding) {
for (ii = 0; ii < numtowrite; ii++)
update_stats(statnum + ii, outdata[0][ii], &min, &max, &avg, &var);
statnum += numtowrite;
}
}
datawrote += numtowrite;
totwrote += numtowrite;
numwritten += numtowrite;
if ((datawrote == abs(*diffbinptr)) &&
(numwritten != numread) &&
(totwrote < cmd->numout)) { /* Add/remove a bin */
int skip, nextdiffbin;
skip = numtowrite;
/* Write the rest of the data after adding/removing a bin */
do {
if (*diffbinptr > 0) {
/* Add a bin */
if (myid > 0)
write_padding(outfiles, local_numdms, block_avg, 1);
numadded++;
totwrote++;
} else {
/* Remove a bin */
numremoved++;
datawrote++;
numwritten++;
skip++;
}
diffbinptr++;
/* Write the part after the diffbin */
numtowrite = numread - numwritten;
if (cmd->numoutP && (totwrote + numtowrite) > cmd->numout)
numtowrite = cmd->numout - totwrote;
nextdiffbin = abs(*diffbinptr) - datawrote;
if (numtowrite > nextdiffbin)
numtowrite = nextdiffbin;
if (myid > 0) {
write_data(outfiles, local_numdms, outdata, skip, numtowrite);
/* Update the statistics and counters */
if (!padding) {
for (ii = 0; ii < numtowrite; ii++)
update_stats(statnum + ii,
outdata[0][skip + ii], &min, &max, &avg, &var);
statnum += numtowrite;
}
}
numwritten += numtowrite;
datawrote += numtowrite;
totwrote += numtowrite;
skip += numtowrite;
/* Stop if we have written out all the data we need to */
if (cmd->numoutP && (totwrote == cmd->numout))
break;
} while (numwritten < numread);
}
/* Stop if we have written out all the data we need to */
if (cmd->numoutP && (totwrote == cmd->numout))
break;
numread = get_data(outdata, blocksperread, &s,
&obsmask, idispdt, offsets, &padding);
}
}
if (myid > 0) {
/* Calculate the amount of padding we need */
if (cmd->numoutP && (cmd->numout > totwrote))
padwrote = padtowrite = cmd->numout - totwrote;
/* Write the new info file for the output data */
idata.dt = dsdt;
update_infodata(&idata, totwrote, padtowrite, diffbins,
numdiffbins, cmd->downsamp);
for (ii = 0; ii < local_numdms; ii++) {
idata.dm = dms[ii];
if (!cmd->nobaryP) {
double baryepoch, barydispdt, baryhifreq;
baryhifreq = idata.freq + (s.num_channels - 1) * idata.chan_wid;
barydispdt = delay_from_dm(dms[ii], doppler(baryhifreq, avgvoverc));
baryepoch = blotoa - (barydispdt / SECPERDAY);
idata.bary = 1;
idata.mjd_i = (int) floor(baryepoch);
idata.mjd_f = baryepoch - idata.mjd_i;
}
sprintf(idata.name, "%s_DM%.2f", outfilenm, dms[ii]);
writeinf(&idata);
}
/* Set the padded points equal to the average data point */
if (idata.numonoff >= 1) {
int index, startpad, endpad;
for (ii = 0; ii < local_numdms; ii++) {
fclose(outfiles[ii]);
sprintf(datafilenm, "%s_DM%.2f.dat", outfilenm, dms[ii]);
outfiles[ii] = chkfopen(datafilenm, "rb+");
}
for (ii = 0; ii < idata.numonoff; ii++) {
index = 2 * ii;
startpad = idata.onoff[index + 1];
if (ii == idata.numonoff - 1)
endpad = idata.N - 1;
else
endpad = idata.onoff[index + 2];
for (jj = 0; jj < local_numdms; jj++)
chkfseek(outfiles[jj], (startpad + 1) * sizeof(float), SEEK_SET);
padtowrite = endpad - startpad;
write_padding(outfiles, local_numdms, avg, padtowrite);
}
}
}
/* Print simple stats and results */
var /= (datawrote - 1);
if (myid == 0)
print_percent_complete(1, 1);
if (myid == 1) {
printf("\n\nDone.\n\nSimple statistics of the output data:\n");
printf(" Data points written: %d\n", totwrote);
if (padwrote)
printf(" Padding points written: %d\n", padwrote);
if (!cmd->nobaryP) {
if (numadded)
printf(" Bins added for barycentering: %d\n", numadded);
if (numremoved)
printf(" Bins removed for barycentering: %d\n", numremoved);
}
printf(" Maximum value of data: %.2f\n", max);
printf(" Minimum value of data: %.2f\n", min);
printf(" Data average value: %.2f\n", avg);
printf(" Data standard deviation: %.2f\n", sqrt(var));
printf("\n");
}
/* Close the files and cleanup */
if (cmd->maskfileP)
free_mask(obsmask);
if (myid > 0) {
for (ii = 0; ii < local_numdms; ii++)
fclose(outfiles[ii]);
free(outfiles);
}
vect_free(outdata[0]);
vect_free(outdata);
vect_free(dms);
free(hostname);
vect_free(idispdt);
vect_free(offsets[0]);
vect_free(offsets);
free(datafilenm);
free(outfilenm);
free(outpath);
if (!cmd->nobaryP) {
vect_free(btoa);
vect_free(ttoa);
vect_free(diffbins);
}
MPI_Finalize();
return (0);
}
ssize_t
mifare_desfire_write_record_ex (MifareTag tag, uint8_t file_no, off_t offset, size_t length, void *data, int cs)
{
return write_data (tag, 0x3B, file_no, offset, length, data, cs);
}
int
fd_read_body (int fd, FILE *out, wgint toread, wgint startpos,
wgint *qtyread, wgint *qtywritten, double *elapsed, int flags)
{
int ret = 0;
#undef max
#define max(a,b) ((a) > (b) ? (a) : (b))
int dlbufsize = max (BUFSIZ, 8 * 1024);
char *dlbuf = xmalloc (dlbufsize);
struct ptimer *timer = NULL;
double last_successful_read_tm = 0;
/* The progress gauge, set according to the user preferences. */
void *progress = NULL;
/* Non-zero if the progress gauge is interactive, i.e. if it can
continually update the display. When true, smaller timeout
values are used so that the gauge can update the display when
data arrives slowly. */
bool progress_interactive = false;
bool exact = !!(flags & rb_read_exactly);
/* Used only by HTTP/HTTPS chunked transfer encoding. */
bool chunked = flags & rb_chunked_transfer_encoding;
wgint skip = 0;
/* How much data we've read/written. */
wgint sum_read = 0;
wgint sum_written = 0;
wgint remaining_chunk_size = 0;
if (flags & rb_skip_startpos)
skip = startpos;
if (opt.verbose)
{
/* If we're skipping STARTPOS bytes, pass 0 as the INITIAL
argument to progress_create because the indicator doesn't
(yet) know about "skipping" data. */
wgint start = skip ? 0 : startpos;
progress = progress_create (start, start + toread);
progress_interactive = progress_interactive_p (progress);
}
if (opt.limit_rate)
limit_bandwidth_reset ();
/* A timer is needed for tracking progress, for throttling, and for
tracking elapsed time. If either of these are requested, start
the timer. */
if (progress || opt.limit_rate || elapsed)
{
timer = ptimer_new ();
last_successful_read_tm = 0;
}
/* Use a smaller buffer for low requested bandwidths. For example,
with --limit-rate=2k, it doesn't make sense to slurp in 16K of
data and then sleep for 8s. With buffer size equal to the limit,
we never have to sleep for more than one second. */
if (opt.limit_rate && opt.limit_rate < dlbufsize)
dlbufsize = opt.limit_rate;
/* Read from FD while there is data to read. Normally toread==0
means that it is unknown how much data is to arrive. However, if
EXACT is set, then toread==0 means what it says: that no data
should be read. */
while (!exact || (sum_read < toread))
{
int rdsize;
double tmout = opt.read_timeout;
if (chunked)
{
if (remaining_chunk_size == 0)
{
char *line = fd_read_line (fd);
char *endl;
if (line == NULL)
{
ret = -1;
break;
}
remaining_chunk_size = strtol (line, &endl, 16);
if (remaining_chunk_size == 0)
{
ret = 0;
if (fd_read_line (fd) == NULL)
ret = -1;
break;
}
}
rdsize = MIN (remaining_chunk_size, dlbufsize);
}
else
rdsize = exact ? MIN (toread - sum_read, dlbufsize) : dlbufsize;
if (progress_interactive)
{
/* For interactive progress gauges, always specify a ~1s
timeout, so that the gauge can be updated regularly even
when the data arrives very slowly or stalls. */
tmout = 0.95;
if (opt.read_timeout)
{
double waittm;
waittm = ptimer_read (timer) - last_successful_read_tm;
if (waittm + tmout > opt.read_timeout)
{
/* Don't let total idle time exceed read timeout. */
tmout = opt.read_timeout - waittm;
if (tmout < 0)
{
/* We've already exceeded the timeout. */
ret = -1, errno = ETIMEDOUT;
break;
}
}
}
}
ret = fd_read (fd, dlbuf, rdsize, tmout);
if (progress_interactive && ret < 0 && errno == ETIMEDOUT)
ret = 0; /* interactive timeout, handled above */
else if (ret <= 0)
break; /* EOF or read error */
if (progress || opt.limit_rate || elapsed)
{
ptimer_measure (timer);
if (ret > 0)
last_successful_read_tm = ptimer_read (timer);
}
if (ret > 0)
{
sum_read += ret;
if (!write_data (out, dlbuf, ret, &skip, &sum_written))
{
ret = -2;
goto out;
}
if (chunked)
{
remaining_chunk_size -= ret;
if (remaining_chunk_size == 0)
if (fd_read_line (fd) == NULL)
{
ret = -1;
break;
}
}
}
if (opt.limit_rate)
limit_bandwidth (ret, timer);
if (progress)
progress_update (progress, ret, ptimer_read (timer));
#ifdef WINDOWS
if (toread > 0 && !opt.quiet)
ws_percenttitle (100.0 *
(startpos + sum_read) / (startpos + toread));
#endif
}
if (ret < -1)
ret = -1;
out:
if (progress)
progress_finish (progress, ptimer_read (timer));
if (elapsed)
*elapsed = ptimer_read (timer);
if (timer)
ptimer_destroy (timer);
if (qtyread)
*qtyread += sum_read;
if (qtywritten)
*qtywritten += sum_written;
free (dlbuf);
return ret;
}
/* relay traffic from source to destination socket
*
*/
static int
relay_traffic( int ssockfd, int dsockfd, struct server_ctx* ctx,
int dfilefd, const struct in_addr* mifaddr )
{
volatile sig_atomic_t quit = 0;
int rc = 0;
ssize_t nmsgs = -1;
ssize_t nrcv = 0, nsent = 0, nwr = 0,
lrcv = 0, lsent = 0;
char* data = NULL;
size_t data_len = g_uopt.rbuf_len;
struct rdata_opt ropt;
time_t pause_time = 0, rfr_tm = time(NULL);
sigset_t ubset;
const int ALLOW_PAUSES = get_flagval( "UDPXY_ALLOW_PAUSES", 0 );
const ssize_t MAX_PAUSE_MSEC =
get_sizeval( "UDPXY_PAUSE_MSEC", 1000);
/* permissible variation in data-packet size */
static const ssize_t t_delta = 0x20;
struct dstream_ctx ds;
static const int SET_PID = 1;
struct tps_data tps;
assert( ctx && mifaddr && MAX_PAUSE_MSEC > 0 );
(void) sigemptyset (&ubset);
sigaddset (&ubset, SIGINT);
sigaddset (&ubset, SIGQUIT);
sigaddset (&ubset, SIGTERM);
/* restore the ability to receive *quit* signals */
rc = sigprocmask (SIG_UNBLOCK, &ubset, NULL);
if (0 != rc) {
mperror (g_flog, errno, "%s: sigprocmask", __func__);
return -1;
}
/* NOPs to eliminate warnings in lean version */
(void)&lrcv; (void)&lsent; (void)&t_delta;
check_fragments( NULL, 0, 0, 0, 0, g_flog );
/* INIT
*/
rc = calc_buf_settings( &nmsgs, NULL );
if (0 != rc) return -1;
TRACE( (void)tmfprintf( g_flog, "Data buffer will hold up to "
"[%d] messages\n", nmsgs ) );
rc = init_dstream_ctx( &ds, ctx->rq.cmd, g_uopt.srcfile, nmsgs );
if( 0 != rc ) return -1;
(void) set_nice( g_uopt.nice_incr, g_flog );
do {
if( NULL == g_uopt.srcfile ) {
rc = set_timeouts( ssockfd, dsockfd,
ctx->rcv_tmout, 0,
ctx->snd_tmout, 0 );
if( 0 != rc ) break;
}
if( dsockfd > 0 ) {
rc = sync_dsockbuf_len( ssockfd, dsockfd );
if( 0 != rc ) break;
rc = send_http_response( dsockfd, 200, "OK" );
if( 0 != rc ) break;
/* timeshift: to detect PAUSE make destination
* socket non-blocking, otherwise make it blocking
* (since it might have been set unblocking earlier)
*/
rc = set_nblock( dsockfd, (ALLOW_PAUSES ? 1 : 0) );
if( 0 != rc ) break;
}
data = malloc(data_len);
if( NULL == data ) {
mperror( g_flog, errno, "%s: malloc", __func__ );
break;
}
if( g_uopt.cl_tpstat )
tpstat_init( &tps, SET_PID );
} while(0);
TRACE( (void)tmfprintf( g_flog, "Relaying traffic from socket[%d] "
"to socket[%d], buffer size=[%d], Rmsgs=[%d], pauses=[%d]\n",
ssockfd, dsockfd, data_len, g_uopt.rbuf_msgs, ALLOW_PAUSES) );
/* RELAY LOOP
*/
ropt.max_frgs = g_uopt.rbuf_msgs;
ropt.buf_tmout = g_uopt.dhold_tmout;
pause_time = 0;
while( (0 == rc) && !(quit = must_quit()) ) {
if( g_uopt.mcast_refresh > 0 ) {
check_mcast_refresh( ssockfd, &rfr_tm, mifaddr );
}
nrcv = read_data( &ds, ssockfd, data, data_len, &ropt );
if( -1 == nrcv ) break;
TRACE( check_fragments( "received new", data_len,
lrcv, nrcv, t_delta, g_flog ) );
lrcv = nrcv;
if( dsockfd && (nrcv > 0) ) {
nsent = write_data( &ds, data, nrcv, dsockfd );
if( -1 == nsent ) break;
if ( nsent < 0 ) {
if ( !ALLOW_PAUSES ) break;
if ( 0 != pause_detect( nsent, MAX_PAUSE_MSEC, &pause_time ) )
break;
}
TRACE( check_fragments("sent", nrcv,
lsent, nsent, t_delta, g_flog) );
lsent = nsent;
}
if( (dfilefd > 0) && (nrcv > 0) ) {
nwr = write_data( &ds, data, nrcv, dfilefd );
if( -1 == nwr )
break;
TRACE( check_fragments( "wrote to file",
nrcv, lsent, nwr, t_delta, g_flog ) );
lsent = nwr;
}
if( ds.flags & F_SCATTERED ) reset_pkt_registry( &ds );
if( uf_TRUE == g_uopt.cl_tpstat )
tpstat_update( ctx, &tps, nsent );
} /* end of RELAY LOOP */
/* CLEANUP
*/
TRACE( (void)tmfprintf( g_flog, "Exited relay loop: received=[%ld], "
"sent=[%ld], quit=[%ld]\n", (long)nrcv, (long)nsent, (long)quit ) );
free_dstream_ctx( &ds );
if( NULL != data ) free( data );
if( 0 != (quit = must_quit()) ) {
TRACE( (void)tmfprintf( g_flog, "Child process=[%d] must quit\n",
getpid()) );
}
return rc;
}
//********************************************************
int read_satpower()
{
int i, j, k, b, n;
int acqnum;
char acqstr[12];
char timemark[3];
char dayofweek[3];
char strmonth[3];
int date;
int hour;
int minute;
int second;
int year;
float fftbin1[515], fftbin2[515], fftbin3[515], fftbin4[515];
char tmp1[8], tmp2[10], tmp3[10], tmp4[10], tmp5[10], tmp6[10], tmp7[10], tmp8[10], tmp9[10];
double dyear;
int month=0, mlen=0, maccum=0;
int centerbin;
int fftcenterbinsofvisiblesats[20];
float temp1,temp2,temp3,temp4;
int visiblesatnum;
struct
{
char satname[10];
int satnum;
float azimuth;
float elevation;
float range;
float velocity;
float meananom;
long orbitnum;
} satinfo[10];
if ((satpowerfile=fopen(satpowerfilename,"r"))==NULL)
{
printf("Cannot open sat power file. \n");
exit(1);
}
while (feof(satpowerfile)==0)
{
fscanf(satpowerfile, "%s%d",acqstr, &acqnum);
fscanf(satpowerfile, "%s%s%s%d%d:%d:%d%d", timemark, dayofweek, strmonth, &date, &hour, &minute, &second, &year);
if (strcmp(strmonth, "Jan")==0) { month=1; mlen=31; maccum=0;}
else if (strcmp(strmonth, "Feb")==0){month=2; mlen=28; maccum=31;}
else if (strcmp(strmonth, "Mar")==0){month=3; mlen=31; maccum=59;}
else if (strcmp(strmonth, "Apr")==0){month=4; mlen=30; maccum=90;}
else if (strcmp(strmonth, "May")==0){month=5; mlen=31; maccum=120;}
else if (strcmp(strmonth, "Jun")==0){month=6; mlen=30; maccum=151;}
else if (strcmp(strmonth, "Jul")==0){month=7; mlen=31; maccum=181;}
else if (strcmp(strmonth, "Aug")==0){month=8; mlen=31; maccum=212;}
else if (strcmp(strmonth, "Sep")==0){month=9; mlen=30; maccum=243;}
else if (strcmp(strmonth, "Oct")==0){month=10;mlen=31; maccum=273;}
else if (strcmp(strmonth, "Nov")==0){month=11;mlen=30; maccum=304;}
else if (strcmp(strmonth, "Dec")==0){month=12;mlen=31; maccum=334;}
else
{
printf("No option selected.");
}
dyear=year+ ((double)maccum + (double)date + ((double)hour+(double)minute/60+(double)second/3600)/24)/365;
i=1;
fscanf(satpowerfile, "%s", satinfo[i].satname);
while (strcmp(satinfo[i].satname, "Bin")!=0)
{
//printf("reading sat num from string: '%s'\n",satinfo[i].satname);
char stuff[256];
int david;
sprintf(stuff,"%s",satinfo[i].satname);
printf("stuff: %s\n",stuff);
sscanf(stuff, "FM%d",&david);
satinfo[i].satnum = david;
printf("david: %d\tsatinfo[i].satnum = %d\n",david,satinfo[i].satnum);
//sscanf(stuff, "FM%d",&satinfo[i].satnum);
//printf("satnum: %d\n",&satinfo[i].satnum);
fscanf(satpowerfile, "%f%f%f%f%f%ld", &satinfo[i].azimuth, &satinfo[i].elevation, &satinfo[i].range, &satinfo[i].velocity, &satinfo[i].meananom, &satinfo[i].orbitnum);
i++;
fscanf(satpowerfile, "%s", satinfo[i].satname);
}
// printf("Reading junk ...\n");
fscanf(satpowerfile,"%s\n %s %s %s %s %s %s %s %s\n", tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7, tmp8, tmp9);
aq.tot1=0; aq.tot2=0; aq.tot3=0; aq.tot4=0;
for (b=1; b<513; b++)
{
fscanf(satpowerfile, "%f %f %f %f", &fftbin1[b], &fftbin2[b], &fftbin3[b], &fftbin4[b]);
aq.tot1 += fftbin1[b]*fftbin1[b];
aq.tot2 += fftbin2[b]*fftbin2[b];
aq.tot3 += fftbin3[b]*fftbin3[b];
aq.tot4 += fftbin4[b]*fftbin4[b];
}
// prepare data for storage
// printf("Preparing data storage ...\n");
aq.dyear=dyear;
aq.month = month;
aq.day = date;
aq.hour = hour;
aq.min = minute;
aq.sec = second;
for (j=1; j<(totsatpassid+1); j++){ // for all sats
for (k=1; k<i; k++){ // for all passes
// check if the current time is between the pass start and end time
if ((dyear>satpasses[j].wstart) && (dyear< satpasses[j].wstop) && (satpasses[j].satid==satinfo[k].satnum)){
aq.satfmnum=satinfo[k].satnum;
aq.satorbitnum=satinfo[k].orbitnum;
aq.meananom=satinfo[k].meananom;
sprintf(aq.station, "GB");
aq.azimuth=satinfo[k].azimuth;
aq.elevation=satinfo[k].elevation;
aq.range=satinfo[k].range;
aq.satchannel=satpasses[j].chan;
centerbin = (int)((float) aq.satchannel*2.5*512/1000);
aq.fftbincenter=centerbin;
aq.bin1[1]=fftbin1[centerbin-3];
aq.bin1[2]=fftbin1[centerbin-2];
aq.bin1[3]=fftbin1[centerbin-1];
aq.bin1[4]=fftbin1[centerbin];
aq.bin1[5]=fftbin1[centerbin+1];
aq.bin1[6]=fftbin1[centerbin+2];
aq.bin1[7]=fftbin1[centerbin+3];
aq.bin2[1]=fftbin2[centerbin-3];
aq.bin2[2]=fftbin2[centerbin-2];
aq.bin2[3]=fftbin2[centerbin-1];
aq.bin2[4]=fftbin2[centerbin];
aq.bin2[5]=fftbin2[centerbin+1];
aq.bin2[6]=fftbin2[centerbin+2];
aq.bin2[7]=fftbin2[centerbin+3];
aq.bin3[1]=fftbin3[centerbin-3];
aq.bin3[2]=fftbin3[centerbin-2];
aq.bin3[3]=fftbin3[centerbin-1];
aq.bin3[4]=fftbin3[centerbin];
aq.bin3[5]=fftbin3[centerbin+1];
aq.bin3[6]=fftbin3[centerbin+2];
aq.bin3[7]=fftbin3[centerbin+3];
aq.bin4[1]=fftbin4[centerbin-3];
aq.bin4[2]=fftbin4[centerbin-2];
aq.bin4[3]=fftbin4[centerbin-1];
aq.bin4[4]=fftbin4[centerbin];
aq.bin4[5]=fftbin4[centerbin+1];
aq.bin4[6]=fftbin4[centerbin+2];
aq.bin4[7]=fftbin4[centerbin+3];
write_data();
channel1=0;
channel2=0;
channel3=0;
channel4=0;
for (n=1; n<8; n++) channel1=channel1+aq.bin1[n];
for (n=1; n<8; n++) channel2=channel2+aq.bin2[n];
for (n=1; n<8; n++) channel3=channel3+aq.bin3[n];
for (n=1; n<8; n++) channel4=channel4+aq.bin4[n];
azIndex=(int) rintf(10*aq.azimuth);
elIndex=(int) rintf(10*aq.elevation);
//printf("%d %d\n", azIndex, elIndex);
mapAzimuth[azIndex].ratio1[elIndex]=(channel1*channel1/(50*2))/(channel2*channel2/(50*2));
mapAzimuth[azIndex].ratio2[elIndex]=(channel3*channel3/(50*2))/(channel4*channel4/(50*2));
// printf("Processing satellite %d \n", satinfo[k].satnum);
}
}
}
}
fclose(satpowerfile);
return 0;
}
static inline void write_next_byte(struct kcs_data *kcs)
{
write_data(kcs, kcs->write_data[kcs->write_pos]);
(kcs->write_pos)++;
(kcs->write_count)--;
}
const
size_t assemble(
uint8_t *const bytecode,
const size_t bytecode_sz,
const char *const assembly,
const size_t asm_sz)
{
char path[PATH_MAX];
snprintf(path, sizeof(path), "/tmp/rappel-input.XXXXXX");
const int t = mkstemp(path);
if (t == -1) {
perror("mkstemp");
exit(EXIT_FAILURE);
}
REQUIRE (unlink(path) == 0);
dprintf(t, BITSTR);
dprintf(t, "section .text vstart=%ld\n", options.start);
write_data(t, (uint8_t *)assembly, asm_sz);
int fildes[2];
REQUIRE (pipe(fildes) == 0);
const pid_t asm_pid = fork();
if (asm_pid < 0) {
perror("fork");
exit(EXIT_FAILURE);
} else if (asm_pid == 0) {
_child(fildes, t);
// Not reached
abort();
}
verbose_printf("nasm is pid %d\n", asm_pid);
REQUIRE (close(fildes[1]) == 0);
mem_assign(bytecode, bytecode_sz, TRAP, TRAP_SZ);
size_t sz = read_data(fildes[0], bytecode, bytecode_sz);
if (sz >= bytecode_sz) {
fprintf(stderr, "Too much bytecode to handle, exiting...\n");
exit(EXIT_FAILURE);
}
int status;
REQUIRE (waitpid(asm_pid, &status, 0) != -1);
if (WIFSIGNALED(status))
fprintf(stderr, "nasm exited with signal %d.\n", WTERMSIG(status));
else if (WIFEXITED(status) && WEXITSTATUS(status))
fprintf(stderr, "nasm exited %d.\n", WEXITSTATUS(status));
REQUIRE (close(t) == 0);
return sz;
}
void *Importobs(void *argument)
{
int operatorID;
char line[MAX_LINE_LEN];
int streamID;
int tsID;
int gridID, zaxisID, taxisID, vlistID;
int i, j;
int nvars = MAX_VARS;
int vdate = 0, vtime = 0;
int vdate0 = 0, vtime0 = 0;
int gridsize, xsize, ysize;
double *xvals = NULL, *yvals = NULL;
double *data[MAX_VARS];
FILE *fp;
char dummy[32], station[32], datetime[32];
float lat, lon, height1, pressure, height2, value;
double latmin = 90, latmax = -90, lonmin = 360, lonmax = -360;
int code;
int index;
double dx, dy;
char *pstation;
cdoInitialize(argument);
cdoOperatorAdd("import_obs", 0, 0, "grid description file or name");
operatorID = cdoOperatorID();
operatorInputArg(cdoOperatorEnter(operatorID));
gridID = cdoDefineGrid(operatorArgv()[0]);
if ( gridInqType(gridID) != GRID_LONLAT )
cdoAbort("Unsupported grid type: %s", gridNamePtr(gridInqType(gridID)));
gridsize = gridInqSize(gridID);
xsize = gridInqXsize(gridID);
ysize = gridInqYsize(gridID);
// printf("gridsize=%d, xsize=%d, ysize=%d\n", gridsize, xsize, ysize);
xvals = (double*) malloc(gridsize*sizeof(double));
yvals = (double*) malloc(gridsize*sizeof(double));
gridInqXvals(gridID, xvals);
gridInqYvals(gridID, yvals);
/* Convert lat/lon units if required */
{
char units[CDI_MAX_NAME];
gridInqXunits(gridID, units);
grid_to_degree(units, gridsize, xvals, "grid center lon");
gridInqYunits(gridID, units);
grid_to_degree(units, gridsize, yvals, "grid center lat");
}
fp = fopen(cdoStreamName(0)->args, "r");
if ( fp == NULL ) { perror(cdoStreamName(0)->args); exit(EXIT_FAILURE); }
vdate = getDate(cdoStreamName(0)->args);
if ( vdate <= 999999 ) vdate = vdate*100 + 1;
streamID = streamOpenWrite(cdoStreamName(1), cdoFiletype());
zaxisID = zaxisCreate(ZAXIS_SURFACE, 1);
taxisID = taxisCreate(TAXIS_ABSOLUTE);
vlistID = vlistCreate();
vlistDefTaxis(vlistID, taxisID);
{
for ( i = 0; i < nvars; ++i ) data[i] = (double*) malloc(gridsize*sizeof(double));
init_vars(vlistID, gridID, zaxisID, nvars);
streamDefVlist(streamID, vlistID);
vdate0 = 0;
vtime0 = 0;
//ntime = 0;
tsID = 0;
while ( readline(fp, line, MAX_LINE_LEN) )
{
sscanf(line, "%s %s %s %g %g %g %d %g %g %g",
dummy, station, datetime, &lat, &lon, &height1, &code, &pressure, &height2, &value);
sscanf(datetime, "%d_%d", &vdate, &vtime);
if ( vdate != vdate0 || vtime != vtime0 )
{
if ( tsID > 0 ) write_data(streamID, vlistID, nvars, data);
vdate0 = vdate;
vtime0 = vtime;
/*
printf("%s %d %d %g %g %g %d %g %g %g\n",
station, vdate, vtime, lat, lon, height1, code, pressure, height2, value);
*/
taxisDefVdate(taxisID, vdate);
taxisDefVtime(taxisID, vtime);
streamDefTimestep(streamID, tsID);
init_data(vlistID, nvars, data);
tsID++;
}
if ( lon < lonmin ) lonmin = lon;
if ( lon > lonmax ) lonmax = lon;
if ( lat < latmin ) latmin = lat;
if ( lat > latmax ) latmax = lat;
dy = yvals[1] - yvals[0];
for ( j = 0; j < ysize; ++j )
if ( lat >= (yvals[j]-dy/2) && lat < (yvals[j]+dy/2) ) break;
dx = xvals[1] - xvals[0];
if ( lon < (xvals[0] - dx/2) && lon < 0 ) lon+=360;
for ( i = 0; i < xsize; ++i )
if ( lon >= (xvals[i]-dx/2) && lon < (xvals[i]+dx/2) ) break;
index = -1;
if ( code == 11 ) index = 0;
if ( code == 17 ) index = 1;
if ( code == 33 ) index = 2;
if ( code == 34 ) index = 3;
//printf("%d %d %d %g %g %g %g\n", i, j, index, dx, dy, lon, lat);
if ( i < xsize && j < ysize && index >= 0 )
{
pstation = station;
while (isalpha(*pstation)) pstation++;
// printf("station %s %d\n", pstation, atoi(pstation));
data[index][j*xsize+i] = value;
data[ 4][j*xsize+i] = height1;
data[ 5][j*xsize+i] = pressure;
// data[ 6][j*xsize+i] = atoi(pstation);
}
/*
printf("%s %d %d %g %g %g %d %g %g %g\n",
station, vdate, vtime, lat, lon, height1, code, pressure, height2, value);
*/
}
write_data(streamID, vlistID, nvars, data);
for ( i = 0; i < nvars; ++i ) free(data[i]);
}
printf("lonmin=%g, lonmax=%g, latmin=%g, latmax=%g\n", lonmin, lonmax, latmin, latmax);
processDefVarNum(vlistNvars(vlistID), streamID);
streamClose(streamID);
fclose(fp);
vlistDestroy(vlistID);
gridDestroy(gridID);
zaxisDestroy(zaxisID);
taxisDestroy(taxisID);
free(xvals);
free(yvals);
cdoFinish();
return (0);
}
int main(int argc,char *argv[])
{
if(argc<2)
{
fprintf(stderr,"usage:%s port\n",argv[0]);
exit(1);
}
// 创建一个socket套接字
sockfd=socket(AF_INET,SOCK_STREAM,0);
// AF_INET表示使用IPv4地址族(AF_INET6表示使用IPv6地址族)
// SOCK_STREAM是数据流（tcp协议）;若传递报文则是udp协议
// 0由tcp决定
if(sockfd<0)
{
fprintf(stderr,"socket:%s\n",strerror(errno));
exit(1);
}
// 设置ip和端口号
struct sockaddr_in addr;
memset(&addr,0,sizeof(addr)); // 目的：把ip地址中的后8个字节置为0
addr.sin_family=AF_INET;
addr.sin_port=htons(atoi(argv[1]));// argv[1]是传入的端口地址,atio将其转为整形，htons将这个整形转为网络字节序
addr.sin_addr.s_addr=INADDR_ANY; // INADDR_ANY是本地的一个宏，代表本地回环地址的网络字节序
//绑定ip和端口号
int len=sizeof(addr);
if(bind(sockfd,(struct sockaddr*)&addr,len)<0)
{
fprintf(stderr,"bind:%s\n",strerror(errno));
exit(1);
}
// 开始监听
if(listen(sockfd,10)<0) // 10 代表客户端的个数，假定有10个
{
fprintf(stderr,"listen:%s\n",strerror(errno));
exit(1);
}
while(1) // 循环获取监听对象
{
struct sockaddr_in client_addr;
int c_len=sizeof(client_addr);
int fd=accept(sockfd,(struct sockaddr*)&client_addr,&len); // 返回一个文件描述符
if(fd<0)
{
fprintf(stderr,"accept:%s\n",strerror(errno));
continue;
}
print_client_information(client_addr); // 将客户端的连接信息打印出来
write_data(fd);//向fd中写入要传入的信息
close(fd);
}
return 0;
}
STATIC void write_lpcomment(void *userhandle, write_modeldata_func write_modeldata, char *string, MYBOOL newlinebefore)
{
write_data(userhandle, write_modeldata, "%s/* %s */\n", (newlinebefore) ? "\n" : "", string);
}
int main (void)
{
IPTR args[ARG_NUM] = { (IPTR) NULL, (IPTR) NULL, (IPTR) NULL, FALSE, FALSE};
struct RDArgs *rda;
ULONG error = 0;
locale = OpenLocale(NULL);
if (!locale)
{
PutStr("Could not open locale!\n");
return -1;
}
rda = ReadArgs(TEMPLATE, args, NULL);
if (rda)
{
BPTR lock_in;
lock_in = Lock((STRPTR)args[ARG_FROM], ACCESS_READ);
if (lock_in)
{
BPTR file_out = Open((STRPTR)args[ARG_TO], MODE_NEWFILE);
if (NULL != file_out)
{
struct FileInfoBlock fib;
UBYTE * data = NULL;
BOOL success = Examine(lock_in, &fib);
/*
** Read the input file into memory
*/
if (fib.fib_Size && DOSTRUE == success)
data = AllocVec(fib.fib_Size, MEMF_ANY);
if (data)
{
ULONG read = Read(lock_in, data, fib.fib_Size);
if (-1 != read)
{
struct sorted_data * sd;
sd = sort(data,
fib.fib_Size,
(STRPTR)args[ARG_COLSTART],
(BOOL)args[ARG_CASE],
(BOOL)args[ARG_NUMERIC]);
error = write_data(sd, file_out);
}
FreeVec(data);
}/* if (data) */
Close(file_out);
} /* if (file_out) */
UnLock(lock_in);
} /* if (lock_in) */
FreeArgs(rda);
}
else
error=RETURN_FAIL;
if (error)
PrintFault(IoErr(), "Sort");
return error;
}
static void write_buffer( filter_t *p_filter, block_t *p_in_buf )
{
write_data( p_filter, p_in_buf->p_buffer, p_in_buf->i_buffer,
is_big_endian( p_filter, p_in_buf ) );
p_filter->p_sys->p_out_buf->i_length += p_in_buf->i_length;
}
void RelaxedVariables::write(std::ostream& s) const
{ write_data(s, allContinuousVars, all_continuous_variable_labels()); }
/* return value:
* 0 success
* -1 fopen failed
* -2 curl initialization failed
* -2 faild to parse the URL (InternetCrackUrl) (windows)
* -2 scheme is neither http nor https
* -3 https responce status is not HTTP_STATUS_OK (windows)
* -4 HttpQueryInfo failed (windows)
*/
int download_simple (char* uri,char* path,int verbose) {
#ifndef HAVE_WINDOWS_H
CURL *curl;
CURLcode res=!CURLE_OK;
char* path_partial=cat(path,".partial",NULL);
curl = curl_easy_init();
if(curl) {
FILE *bodyfile;
char* current=get_opt("ros.proxy",1);
if(current) {
/*<[protocol://][user:password@]proxyhost[:port]>*/
char *reserve=current,*protocol=NULL,*userpwd=NULL,*port=NULL,*uri=NULL;
int pos=position_char("/",current);
if(pos>0 && current[pos-1]==':' && current[pos+1]=='/')
protocol=current,current[pos-1]='\0',current=current+pos+2;
pos=position_char("@",current);
if(pos!=-1)
userpwd=current,current[pos]='\0',current=current+pos+1;
pos=position_char(":",current);
if(pos!=-1)
current[pos]='\0',port=current+pos+1,uri=current;
curl_easy_setopt(curl, CURLOPT_PROXY, uri);
if(port)
curl_easy_setopt(curl, CURLOPT_PROXYPORT,atoi(port));
if(userpwd)
curl_easy_setopt(curl, CURLOPT_PROXYUSERPWD, userpwd);
s(reserve);
}
count=0,content_length=0;
curl_easy_setopt(curl, CURLOPT_URL, uri);
curl_easy_setopt(curl, CURLOPT_FOLLOWLOCATION, 1L);
curl_easy_setopt(curl, CURLOPT_WRITEFUNCTION, write_data);
curl_easy_setopt(curl, CURLOPT_HEADERFUNCTION, header_callback);
curl_easy_setopt(curl, CURLOPT_HEADERDATA, &verbose);
bodyfile = fopen(path_partial,"wb");
if (bodyfile == NULL) {
curl_easy_cleanup(curl);
s(path_partial);
return -1;
}
curl_easy_setopt(curl,CURLOPT_WRITEDATA,bodyfile);
res = curl_easy_perform(curl);
if(res != CURLE_OK && verbose) {
fprintf(stderr, "curl_easy_perform() failed: %s\n",
curl_easy_strerror(res));
}
curl_easy_cleanup(curl);
fclose(bodyfile);
}
if(res != CURLE_OK)
return -2;
#else
URL_COMPONENTS u;
TCHAR szHostName[4096];
TCHAR szUrlPath[4096];
FILE *bodyfile;
char* path_partial=cat(path,".partial",NULL);
bodyfile = fopen(path_partial,"wb");
if (bodyfile == NULL) {
s(path_partial);
return -1;
}
u.dwStructSize = sizeof( u );
u.dwSchemeLength = 1;
u.dwHostNameLength = 4096;
u.dwUserNameLength = 1;
u.dwPasswordLength = 1;
u.dwUrlPathLength = 4096;
u.dwExtraInfoLength = 1;
u.lpszScheme = NULL;
u.lpszHostName = szHostName;
u.lpszUserName = NULL;
u.lpszPassword = NULL;
u.lpszUrlPath = szUrlPath;
u.lpszExtraInfo = NULL;
if(!InternetCrackUrl(uri,(DWORD)strlen(uri),0,&u)) {
fclose(bodyfile);
return -2;
}
HINTERNET hSession = InternetOpen("WinInet",INTERNET_OPEN_TYPE_PRECONFIG,NULL,NULL,0);
HINTERNET hConnection = InternetConnect(hSession,szHostName,u.nPort,NULL,NULL,INTERNET_SERVICE_HTTP,0,0);
DWORD dwFlags = INTERNET_FLAG_RELOAD | INTERNET_FLAG_DONT_CACHE;
if(INTERNET_SCHEME_HTTP == u.nScheme) {
}else if( INTERNET_SCHEME_HTTPS == u.nScheme ) {
dwFlags = dwFlags | INTERNET_FLAG_SECURE| INTERNET_FLAG_IGNORE_CERT_DATE_INVALID| INTERNET_FLAG_IGNORE_CERT_CN_INVALID;
}else {
fclose(bodyfile);
return -2;
}
HINTERNET hRequest = HttpOpenRequest(hConnection,"GET",szUrlPath,NULL,NULL,NULL,dwFlags,0);
HttpSendRequest(hRequest,NULL,0,NULL,0);
DWORD dwStatusCode,dwContentLen;
DWORD dwLength = sizeof(DWORD);
if(HttpQueryInfo(hRequest,HTTP_QUERY_CONTENT_LENGTH | HTTP_QUERY_FLAG_NUMBER,&dwContentLen,&dwLength,0))
content_length=dwContentLen;
if(!HttpQueryInfo(hRequest,HTTP_QUERY_STATUS_CODE|HTTP_QUERY_FLAG_NUMBER,&dwStatusCode,&dwLength,0)) {
fclose(bodyfile);
return -4;
}
if(HTTP_STATUS_OK != dwStatusCode) {
fclose(bodyfile);
return -3;
}
char pData[10000];
DWORD dwBytesRead = 1;
count=0;
while (dwBytesRead) {
InternetReadFile(hRequest, pData, 99, &dwBytesRead);
pData[dwBytesRead] = 0;
write_data(pData,dwBytesRead,1,bodyfile);
}
fclose(bodyfile);
#endif
fprintf(stderr, "\n");
int ret=rename_file(path_partial,path);
s(path_partial);
return ret?0:-3;
}
/****************************************************************************
Desc: This function handle the details of extracting the parameters
needed to interpret the request and then generating the response
HTML page
****************************************************************************/
RCODE F_FFilePage::display(
FLMUINT uiNumParams,
const char ** ppszParams)
{
RCODE rc = FERR_OK;
#define GENERIC_SIZE_B 20
char szFrom[ GENERIC_SIZE_B];
char szBucket[ 4];
FLMUINT uiBucket;
FFILE localFFile;
FFILE * pFile;
FLMBOOL bRefresh;
void * pvAddress;
char szAddress[GENERIC_SIZE_B];
char szLink[GENERIC_SIZE_B];
FLMBOOL bFlmLocked = FALSE;
DATASTRUCT DataStruct;
FLMBYTE * pszTemp = NULL;
FLMBYTE * pszTemp1 = NULL;
if( RC_BAD( rc = f_alloc( 150, &pszTemp)))
{
printErrorPage( rc, TRUE, (char *)"Failed to allocate temporary buffer");
goto Exit;
}
if( RC_BAD( rc = f_alloc( 150, &pszTemp1)))
{
printErrorPage( rc, TRUE, (char *)"Failed to allocate temporary buffer");
goto Exit;
}
// Initialize a few variables first...
szFrom[0] = '\0';
szBucket[0] = '\0';
pFile = NULL;
// Get the "From" parameter. We use this to determine everything else.
if (RC_BAD( rc = ExtractParameter( uiNumParams,
ppszParams,
"From",
sizeof( szFrom),
szFrom)))
{
goto Exit;
}
f_mutexLock( gv_FlmSysData.hShareMutex);
bFlmLocked = TRUE;
if (!f_stricmp( szFrom, "FileHashTbl"))
{
// Get the hash bucket index
if (RC_BAD( rc = ExtractParameter( uiNumParams,
ppszParams,
"Bucket",
sizeof( szBucket),
szBucket)))
{
goto Exit;
}
uiBucket = f_atoud( szBucket);
pFile = (FFILE *)gv_FlmSysData.pFileHashTbl[uiBucket].pFirstInBucket;
}
else if ( (f_stricmp( szFrom, "SCacheBlock") == 0) ||
(f_stricmp( szFrom, "RCache") == 0) ||
(f_stricmp( szFrom, "FDB") == 0))
{
// Get the FFile address and the Hash Bucket
if (RC_BAD( rc = ExtractParameter( uiNumParams,
ppszParams,
"Bucket",
sizeof( szBucket),
szBucket)))
{
goto Exit;
}
uiBucket = f_atoud( szBucket);
if (RC_BAD( rc = ExtractParameter( uiNumParams,
ppszParams,
"Address",
sizeof( szAddress),
szAddress)))
{
goto Exit;
}
pvAddress = (void *)f_atoud( szAddress);
pFile = (FFILE *)gv_FlmSysData.pFileHashTbl[uiBucket].pFirstInBucket;
while (pFile && (void *)pFile != pvAddress)
{
pFile = pFile->pNext;
}
}
else if (f_stricmp( szFrom, "FlmSysData") == 0)
{
// Get the Link and the FFile address
if (RC_BAD( rc = ExtractParameter( uiNumParams,
ppszParams,
"Link",
sizeof( szLink),
szLink)))
{
goto Exit;
}
if (RC_BAD( rc = ExtractParameter( uiNumParams,
ppszParams,
"Address",
sizeof( szAddress),
szAddress)))
{
goto Exit;
}
pvAddress = (void *)f_atoud( szAddress);
if (f_stricmp( szLink, "pMrnuFile") == 0)
{
pFile = gv_FlmSysData.pMrnuFile;
// Now let's make sure we are looking at the right FFile...
while (pFile && (void *)pFile != pvAddress)
{
pFile = pFile->pNextNUFile;
}
}
else if (f_stricmp( szLink, "pLrnuFile") == 0)
{
pFile = gv_FlmSysData.pLrnuFile;
// Now let's make sure we are looking at the right FFile...
while (pFile && (void *)pFile != pvAddress)
{
pFile = pFile->pPrevNUFile;
}
}
}
else if (f_stricmp( szFrom, "FFile") == 0)
{
// We need to get the Link, Bucket & Address
if (RC_BAD(rc = ExtractParameter( uiNumParams,
ppszParams,
"Link",
sizeof( szLink),
szLink)))
{
goto Exit;
}
if (RC_BAD(rc = ExtractParameter( uiNumParams,
ppszParams,
"Address",
sizeof( szAddress),
szAddress)))
{
goto Exit;
}
pvAddress = (void *)f_atoud( szAddress);
if (RC_BAD(rc = ExtractParameter( uiNumParams,
ppszParams,
"Bucket",
sizeof( szBucket),
szBucket)))
{
goto Exit;
}
uiBucket = f_atoud( szBucket);
// First, let's get a reference to an FFile from the specified bucket
if (gv_FlmSysData.pFileHashTbl[uiBucket].pFirstInBucket)
{
pFile = (FFILE *)gv_FlmSysData.pFileHashTbl[uiBucket].pFirstInBucket;
}
// Now let's make sure we are looking at the right FFile...
while (pFile && (void *)pFile != pvAddress)
{
pFile = pFile->pNext;
}
// Now what link are we supposed to follow?
if (f_stricmp( szLink, "pNext") == 0)
{
pFile = pFile->pNext;
}
else if (f_stricmp( szLink, "pPrev") == 0)
{
pFile = pFile->pPrev;
}
else if (f_stricmp( szLink, "pNextNUFile") == 0)
{
pFile = pFile->pNextNUFile;
}
else if (f_stricmp( szLink, "pPrevNUFile") == 0)
{
pFile = pFile->pPrevNUFile;
}
}
// Gather additional data if present. Initialize the structure before
// using it.
f_memset( &DataStruct, 0, sizeof(DataStruct));
if (pFile)
{
f_memcpy( &localFFile, pFile, sizeof(localFFile));
if (pFile->pSCacheList)
{
DataStruct.SCacheBlkAddress = pFile->pSCacheList->uiBlkAddress;
DataStruct.SCacheLowTransID = scaGetLowTransID( pFile->pSCacheList),
DataStruct.SCacheHighTransID = pFile->pSCacheList->uiHighTransID;
}
if (pFile->pPendingWriteList)
{
DataStruct.PendingWriteBlkAddress = pFile->pPendingWriteList->uiBlkAddress;
DataStruct.PendingWriteLowTransID = scaGetLowTransID( pFile->pPendingWriteList),
DataStruct.PendingWriteHighTransID = pFile->pPendingWriteList->uiHighTransID;
}
if (pFile->pLastDirtyBlk)
{
DataStruct.LastDirtyBlkAddress = pFile->pLastDirtyBlk->uiBlkAddress;
DataStruct.LastDirtyLowTransID = scaGetLowTransID( pFile->pLastDirtyBlk),
DataStruct.LastDirtyHighTransID = pFile->pLastDirtyBlk->uiHighTransID;
}
if (pFile->pFirstRecord)
{
DataStruct.FirstRecordContainer = pFile->pFirstRecord->uiContainer;
DataStruct.FirstRecordDrn = pFile->pFirstRecord->uiDrn;
DataStruct.FirstRecordLowTransId = pFile->pFirstRecord->uiLowTransId;
}
if (pFile->pLastRecord)
{
DataStruct.LastRecordContainer = pFile->pLastRecord->uiContainer;
DataStruct.LastRecordDrn = pFile->pLastRecord->uiDrn;
DataStruct.LastRecordLowTransId = pFile->pLastRecord->uiLowTransId;
}
}
f_mutexUnlock( gv_FlmSysData.hShareMutex);
bFlmLocked = FALSE;
stdHdr();
fnPrintf( m_pHRequest, HTML_DOCTYPE);
fnPrintf( m_pHRequest, "<html>\n");
// Determine if we are being requested to refresh this page or not.
if ((bRefresh = DetectParameter( uiNumParams,
ppszParams,
"Refresh")) == TRUE)
{
// Send back the page with a refresh command in the header
f_sprintf( (char *)pszTemp, "%s/FFile?Refresh&From=%s&Bucket=%s",
m_pszURLString,
szFrom, szBucket);
fnPrintf( m_pHRequest,
"<HEAD>"
"<META http-equiv=\"refresh\" content=\"5; url=%s\">"
"<TITLE>FFile Structure</TITLE>\n", pszTemp);
}
else
{
fnPrintf( m_pHRequest, "<HEAD><TITLE>FFile Structure</TITLE>\n");
}
printStyle();
fnPrintf( m_pHRequest, "</HEAD>\n");
fnPrintf( m_pHRequest, "<body>\n");
// If we are not to refresh this page, then don't include the
// refresh meta command
if (!bRefresh)
{
f_sprintf( (char *)pszTemp,
"<A HREF=%s/FFile?Refresh&From=%s&Bucket=%s>Start Auto-refresh (5 sec.)</A>",
m_pszURLString, szFrom, szBucket);
}
else
{
f_sprintf( (char *)pszTemp,
"<A HREF=%s/FFile?From=%s&Bucket=%s>Stop Auto-refresh</A>",
m_pszURLString, szFrom, szBucket);
}
// Prepare the refresh link.
f_sprintf( (char *)pszTemp1,
"<A HREF=%s/FFile?From=%s&Bucket=%s>Refresh</A>",
m_pszURLString, szFrom, szBucket);
// Show the table headings and the refresh option.
if (pFile)
{
// Write out the table headings
printTableStart( "FFile Structure", 4, 100);
printTableRowStart();
printColumnHeading( "", JUSTIFY_LEFT, FLM_IMON_COLOR_PUTTY_1, 4, 1, FALSE);
fnPrintf( m_pHRequest, "%s, ", pszTemp1);
fnPrintf( m_pHRequest, "%s\n", pszTemp);
printColumnHeadingClose();
printTableRowEnd();
// Write out the table headings.
printTableRowStart();
printColumnHeading( "Byte Offset (hex)");
printColumnHeading( "Field Name");
printColumnHeading( "Field Type");
printColumnHeading( "Value");
printTableRowEnd();
write_data( (pFile ? &localFFile: NULL), (void *)pFile, &DataStruct);
}
else
{
// Write out an error page...
fnPrintf( m_pHRequest,
"<P>Unable to find the FFile structure that you requested."
" This is probably because the state of the cache changed between "
"the time that you displayed the previous page and the time that you "
"clicked on the link that brought you here.\n"
"<P>Click on your browser's \"Back\" button, then click \"Reload\" "
"and then try the link again.\n");
}
fnPrintf( m_pHRequest, "</body></html>\n");
fnEmit();
Exit:
if (bFlmLocked)
{
f_mutexUnlock( gv_FlmSysData.hShareMutex);
bFlmLocked = FALSE;
}
if (pszTemp)
{
f_free( &pszTemp);
}
if (pszTemp1)
{
f_free( &pszTemp1);
}
return( rc);
}
static void
io_write (SoupSocket *sock, SoupMessage *msg)
{
SoupMessagePrivate *priv = SOUP_MESSAGE_GET_PRIVATE (msg);
SoupMessageIOData *io = priv->io_data;
write_more:
switch (io->write_state) {
case SOUP_MESSAGE_IO_STATE_NOT_STARTED:
return;
case SOUP_MESSAGE_IO_STATE_HEADERS:
if (!io->write_buf->len) {
io->get_headers_cb (msg, io->write_buf,
&io->write_encoding,
io->user_data);
if (!io->write_buf->len) {
soup_message_io_pause (msg);
return;
}
}
if (!write_data (msg, io->write_buf->str,
io->write_buf->len, FALSE))
return;
g_string_truncate (io->write_buf, 0);
if (io->write_encoding == SOUP_ENCODING_CONTENT_LENGTH) {
SoupMessageHeaders *hdrs =
(io->mode == SOUP_MESSAGE_IO_CLIENT) ?
msg->request_headers : msg->response_headers;
io->write_length = soup_message_headers_get_content_length (hdrs);
}
if (io->mode == SOUP_MESSAGE_IO_SERVER &&
SOUP_STATUS_IS_INFORMATIONAL (msg->status_code)) {
if (msg->status_code == SOUP_STATUS_CONTINUE) {
/* Stop and wait for the body now */
io->write_state =
SOUP_MESSAGE_IO_STATE_BLOCKING;
io->read_state = io_body_state (io->read_encoding);
} else {
/* We just wrote a 1xx response
* header, so stay in STATE_HEADERS.
* (The caller will pause us from the
* wrote_informational callback if he
* is not ready to send the final
* response.)
*/
}
} else if (io->mode == SOUP_MESSAGE_IO_CLIENT &&
soup_message_headers_get_expectations (msg->request_headers) & SOUP_EXPECTATION_CONTINUE) {
/* Need to wait for the Continue response */
io->write_state = SOUP_MESSAGE_IO_STATE_BLOCKING;
io->read_state = SOUP_MESSAGE_IO_STATE_HEADERS;
} else {
io->write_state = io_body_state (io->write_encoding);
/* If the client was waiting for a Continue
* but we sent something else, then they're
* now done writing.
*/
if (io->mode == SOUP_MESSAGE_IO_SERVER &&
io->read_state == SOUP_MESSAGE_IO_STATE_BLOCKING)
io->read_state = SOUP_MESSAGE_IO_STATE_FINISHING;
}
SOUP_MESSAGE_IO_PREPARE_FOR_CALLBACK;
if (SOUP_STATUS_IS_INFORMATIONAL (msg->status_code)) {
soup_message_wrote_informational (msg);
soup_message_cleanup_response (msg);
} else
soup_message_wrote_headers (msg);
SOUP_MESSAGE_IO_RETURN_IF_CANCELLED_OR_PAUSED;
break;
case SOUP_MESSAGE_IO_STATE_BLOCKING:
io_read (sock, msg);
/* If io_read reached a point where we could write
* again, it would have recursively called io_write.
* So (a) we don't need to try to keep writing, and
* (b) we can't anyway, because msg may have been
* destroyed.
*/
return;
case SOUP_MESSAGE_IO_STATE_BODY:
if (!io->write_length && io->write_encoding != SOUP_ENCODING_EOF) {
wrote_body:
io->write_state = SOUP_MESSAGE_IO_STATE_FINISHING;
SOUP_MESSAGE_IO_PREPARE_FOR_CALLBACK;
soup_message_wrote_body (msg);
SOUP_MESSAGE_IO_RETURN_IF_CANCELLED_OR_PAUSED;
break;
}
if (!io->write_chunk) {
io->write_chunk = soup_message_body_get_chunk (io->write_body, io->write_body_offset);
if (!io->write_chunk) {
soup_message_io_pause (msg);
return;
}
if (io->write_chunk->length > io->write_length &&
io->write_encoding != SOUP_ENCODING_EOF) {
/* App is trying to write more than it
* claimed it would; we have to truncate.
*/
SoupBuffer *truncated =
soup_buffer_new_subbuffer (io->write_chunk,
0, io->write_length);
soup_buffer_free (io->write_chunk);
io->write_chunk = truncated;
} else if (io->write_encoding == SOUP_ENCODING_EOF &&
!io->write_chunk->length)
goto wrote_body;
}
if (!write_data (msg, io->write_chunk->data,
io->write_chunk->length, TRUE))
return;
soup_message_body_wrote_chunk (io->write_body, io->write_chunk);
io->write_body_offset += io->write_chunk->length;
soup_buffer_free (io->write_chunk);
io->write_chunk = NULL;
SOUP_MESSAGE_IO_PREPARE_FOR_CALLBACK;
soup_message_wrote_chunk (msg);
SOUP_MESSAGE_IO_RETURN_IF_CANCELLED_OR_PAUSED;
break;
case SOUP_MESSAGE_IO_STATE_CHUNK_SIZE:
if (!io->write_chunk) {
io->write_chunk = soup_message_body_get_chunk (io->write_body, io->write_body_offset);
if (!io->write_chunk) {
soup_message_io_pause (msg);
return;
}
g_string_append_printf (io->write_buf, "%lx\r\n",
(unsigned long) io->write_chunk->length);
io->write_body_offset += io->write_chunk->length;
}
if (!write_data (msg, io->write_buf->str,
io->write_buf->len, FALSE))
return;
g_string_truncate (io->write_buf, 0);
if (io->write_chunk->length == 0) {
/* The last chunk has no CHUNK_END... */
io->write_state = SOUP_MESSAGE_IO_STATE_TRAILERS;
break;
}
io->write_state = SOUP_MESSAGE_IO_STATE_CHUNK;
/* fall through */
case SOUP_MESSAGE_IO_STATE_CHUNK:
if (!write_data (msg, io->write_chunk->data,
io->write_chunk->length, TRUE))
return;
soup_message_body_wrote_chunk (io->write_body, io->write_chunk);
soup_buffer_free (io->write_chunk);
io->write_chunk = NULL;
io->write_state = SOUP_MESSAGE_IO_STATE_CHUNK_END;
SOUP_MESSAGE_IO_PREPARE_FOR_CALLBACK;
soup_message_wrote_chunk (msg);
SOUP_MESSAGE_IO_RETURN_IF_CANCELLED_OR_PAUSED;
/* fall through */
case SOUP_MESSAGE_IO_STATE_CHUNK_END:
if (!write_data (msg, SOUP_MESSAGE_IO_EOL,
SOUP_MESSAGE_IO_EOL_LEN, FALSE))
return;
io->write_state = SOUP_MESSAGE_IO_STATE_CHUNK_SIZE;
break;
case SOUP_MESSAGE_IO_STATE_TRAILERS:
if (!write_data (msg, SOUP_MESSAGE_IO_EOL,
SOUP_MESSAGE_IO_EOL_LEN, FALSE))
return;
io->write_state = SOUP_MESSAGE_IO_STATE_FINISHING;
SOUP_MESSAGE_IO_PREPARE_FOR_CALLBACK;
soup_message_wrote_body (msg);
SOUP_MESSAGE_IO_RETURN_IF_CANCELLED_OR_PAUSED;
/* fall through */
case SOUP_MESSAGE_IO_STATE_FINISHING:
if (io->write_tag) {
g_signal_handler_disconnect (io->sock, io->write_tag);
io->write_tag = 0;
}
io->write_state = SOUP_MESSAGE_IO_STATE_DONE;
if (io->mode == SOUP_MESSAGE_IO_CLIENT) {
io->read_state = SOUP_MESSAGE_IO_STATE_HEADERS;
io_read (sock, msg);
} else
soup_message_io_finished (msg);
return;
case SOUP_MESSAGE_IO_STATE_DONE:
default:
g_return_if_reached ();
}
goto write_more;
}
static int
write_16(FILE *out, uint16_t val)
{
val = htons(val);
return write_data(out, &val, sizeof(val));
}
