Lines Matching refs:n_uv
425 …tion (struct GMTAPI_CTRL *API, double **X, double *u, double *v, uint64_t n_uv, unsigned int mode,… in gpsgridder_do_gps_normalization() argument
437 for (i = 0; i < n_uv; i++) { /* Find mean u and v-values */ in gpsgridder_do_gps_normalization()
444 coeff[GPSGRIDDER_MEAN_U] /= n_uv; /* Average u value to remove/restore */ in gpsgridder_do_gps_normalization()
445 coeff[GPSGRIDDER_MEAN_V] /= n_uv; /* Average v value to remove/restore */ in gpsgridder_do_gps_normalization()
450 coeff[GPSGRIDDER_MEAN_X] /= n_uv; /* Mean x */ in gpsgridder_do_gps_normalization()
451 coeff[GPSGRIDDER_MEAN_Y] /= n_uv; /* Mean y */ in gpsgridder_do_gps_normalization()
452 for (i = 0; i < n_uv; i++) { in gpsgridder_do_gps_normalization()
479 for (i = 0; i < n_uv; i++) { /* Also find min/max or residuals in the process */ in gpsgridder_do_gps_normalization()
503 for (i = 0; i < n_uv; i++) { /* Normalize 0-1 */ in gpsgridder_do_gps_normalization()
601 uint64_t col, row, n_read, p, k, i, j, seg, n_uv, n_params, n_ok = 0, ij; in GMT_gpsgridder() local
683 n_uv = n_read = 0; in GMT_gpsgridder()
688 for (p = 0; p < n_uv; p++) gmt_M_free (GMT, X[p]); in GMT_gpsgridder()
711 X[n_uv] = gmt_M_memory (GMT, NULL, n_cols, double); /* Allocate space for this constraint */ in GMT_gpsgridder()
712 X[n_uv][GMT_X] = in[GMT_X]; /* Save x,y */ in GMT_gpsgridder()
713 X[n_uv][GMT_Y] = in[GMT_Y]; in GMT_gpsgridder()
716 for (i = 0; !skip && i < n_uv; i++) { in GMT_gpsgridder()
717 r = gpsgridder_get_gps_radius (GMT, X[i], X[n_uv]); in GMT_gpsgridder()
736 u[n_uv] = in[GPSGRIDDER_U]; v[n_uv] = in[GPSGRIDDER_V]; /* Save current u,v data pair */ in GMT_gpsgridder()
738 X[n_uv][GPSGRIDDER_WU] = in[4]; /* First just copy over what we got */ in GMT_gpsgridder()
739 X[n_uv][GPSGRIDDER_WV] = in[5]; in GMT_gpsgridder()
741 err_sum_u += X[n_uv][GPSGRIDDER_WU]*X[n_uv][GPSGRIDDER_WU]; /* Update u variance */ in GMT_gpsgridder()
742 err_sum_v += X[n_uv][GPSGRIDDER_WV]*X[n_uv][GPSGRIDDER_WV]; /* Update v variance */ in GMT_gpsgridder()
743 …err_sum += X[n_uv][GPSGRIDDER_WU]*X[n_uv][GPSGRIDDER_WU] + X[n_uv][GPSGRIDDER_WV]*X[n_uv][GPSGRIDD… in GMT_gpsgridder()
744 X[n_uv][GPSGRIDDER_WU] = 1.0 / X[n_uv][GPSGRIDDER_WU]; /* We will square these weights later */ in GMT_gpsgridder()
745 X[n_uv][GPSGRIDDER_WV] = 1.0 / X[n_uv][GPSGRIDDER_WV]; /* We will square these weights later */ in GMT_gpsgridder()
748 var_sum += u[n_uv] * u[n_uv] + v[n_uv] * v[n_uv]; in GMT_gpsgridder()
749 n_uv++; /* Added a new data constraint */ in GMT_gpsgridder()
750 if (n_uv == n_alloc) { /* Get more memory */ in GMT_gpsgridder()
759 for (p = 0; p < n_uv; p++) gmt_M_free (GMT, X[p]); in GMT_gpsgridder()
764 …n_params = 2 * n_uv; /* Dimension of array is twice the size since using both u & v as separate ob… in GMT_gpsgridder()
765 X = gmt_M_memory (GMT, X, n_uv, double *); /* Realloc to exact size */ in GMT_gpsgridder()
767 v = gmt_M_memory (GMT, v, n_uv, double); in GMT_gpsgridder()
768 GMT_Report (API, GMT_MSG_INFORMATION, "Found %" PRIu64 " unique data constraints\n", n_uv); in GMT_gpsgridder()
772 err_sum_u = sqrt (err_sum_u / n_uv); in GMT_gpsgridder()
773 err_sum_v = sqrt (err_sum_v / n_uv); in GMT_gpsgridder()
774 err_sum = sqrt (0.5 * err_sum / n_uv); in GMT_gpsgridder()
789 for (p = 0; p < n_uv; p++) gmt_M_free (GMT, X[p]); in GMT_gpsgridder()
798 … pairs, yielding a %" PRIu64 " by %" PRIu64 " set of linear equations\n", n_uv, n_params, n_params… in GMT_gpsgridder()
803 for (p = 0; p < n_uv; p++) gmt_M_free (GMT, X[p]); in GMT_gpsgridder()
839 for (p = 0; p < n_uv; p++) gmt_M_free (GMT, X[p]); in GMT_gpsgridder()
855 for (p = 0; p < n_uv; p++) gmt_M_free (GMT, X[p]); in GMT_gpsgridder()
864 orig_u = gmt_M_memory (GMT, NULL, n_uv, double); in GMT_gpsgridder()
865 orig_v = gmt_M_memory (GMT, NULL, n_uv, double); in GMT_gpsgridder()
866 gmt_M_memcpy (orig_u, u, n_uv, double); in GMT_gpsgridder()
867 gmt_M_memcpy (orig_v, v, n_uv, double); in GMT_gpsgridder()
881 gpsgridder_do_gps_normalization (API, X, u, v, n_uv, normalize, norm); in GMT_gpsgridder()
894 …off = n_uv * n_params; /* Separation in 1-D index between rows evaluating u and v for same column … in GMT_gpsgridder()
895 for (row = 0; row < n_uv; row++) { /* For each data constraint pair (u,v)_row */ in GMT_gpsgridder()
896 for (col = 0; col < n_uv; col++) { /* For each body force pair (f_x,f_y)_col */ in GMT_gpsgridder()
898 Guv_ij = Gu_ij + n_uv; /* Index for Guv term in equation for u */ in GMT_gpsgridder()
908 gmt_M_memcpy (&u[n_uv], v, n_uv, double); /* Place v array at end of u array */ in GMT_gpsgridder()
936 …for (row = 0; row < n_uv; row++) { /* Set S, the diagonal squared weights (=1/sigma^2) matrix if g… in GMT_gpsgridder()
938 …S[row+n_uv] = (Ctrl->W.mode == GPSGRIDDER_GOT_SIG) ? X[row][GPSGRIDDER_WV] * X[row][GPSGRIDDER_WV]… in GMT_gpsgridder()
1012 for (p = 0; p < n_uv; p++) gmt_M_free (GMT, X[p]); in GMT_gpsgridder()
1057 f_y = &obs[n_uv]; /* Halfway down the array we find the start of the f_y body forces */ in GMT_gpsgridder()
1060 for (p = 0; p < n_uv; p++) fprintf (fp, "%g\t%g\n", f_x[p], f_y[p]); in GMT_gpsgridder()
1068 uint64_t e_dim[GMT_DIM_SIZE] = {1, 1, n_uv, 8+2*Ctrl->W.active}; in GMT_gpsgridder()
1087 S->n_rows = n_uv; in GMT_gpsgridder()
1089 for (j = 0; j < n_uv; j++) { /* For each data constraint pair (u,v) */ in GMT_gpsgridder()
1093 here[GPSGRIDDER_V] = predicted[j+n_uv]; in GMT_gpsgridder()
1096 for (p = 0; p < n_uv; p++) { in GMT_gpsgridder()
1147 rms_u = sqrt (rms_u / n_uv); in GMT_gpsgridder()
1148 rms_v = sqrt (rms_v / n_uv); in GMT_gpsgridder()
1154 …arate u Misfit: N = %u\tMean = %g\tStd.dev = %g\tRMS = %g\tChi^2 = %g\n", n_uv, mean_u, std_u, rms… in GMT_gpsgridder()
1155 …arate v Misfit: N = %u\tMean = %g\tStd.dev = %g\tRMS = %g\tChi^2 = %g\n", n_uv, mean_v, std_v, rms… in GMT_gpsgridder()
1159 …MATION, "Separate u Misfit: N = %u\tMean = %g\tStd.dev = %g\tRMS = %g\n", n_uv, mean_u, std_u, rms… in GMT_gpsgridder()
1160 …MATION, "Separate v Misfit: N = %u\tMean = %g\tStd.dev = %g\tRMS = %g\n", n_uv, mean_v, std_v, rms… in GMT_gpsgridder()
1206 for (p = 0; p < n_uv; p++) { in GMT_gpsgridder()
1264 for (j = 0; j < n_uv; j++) { /* For each data constraint pair (u,v) */ in GMT_gpsgridder()
1267 for (p = 0; p < (int64_t)n_uv; p++) { /* Initialize before adding up all body forces */ in GMT_gpsgridder()
1285 rms = sqrt (rms / n_uv); in GMT_gpsgridder()
1286 rms_u = sqrt (rms_u / n_uv); in GMT_gpsgridder()
1287 rms_v = sqrt (rms_v / n_uv); in GMT_gpsgridder()
1308 #pragma omp parallel for private(row,V,col,ij,p,G) shared(Out,yp,xp,n_uv,GMT,X,par,geo,f_x,f_y,norm… in GMT_gpsgridder()
1317 …for (p = 0, V[GPSGRIDDER_U] = V[GPSGRIDDER_V] = 0.0; p < (int64_t)n_uv; p++) { /* Initialize befor… in GMT_gpsgridder()
1381 #pragma omp parallel for private(row,V,col,ij,p,G) shared(Out,yp,xp,n_uv,GMT,X,par,geo,f_x,f_y,norm… in GMT_gpsgridder()
1390 …for (p = 0, V[GPSGRIDDER_U] = V[GPSGRIDDER_V] = 0.0; p < (int64_t)n_uv; p++) { /* Initialize befor… in GMT_gpsgridder()
1418 for (p = 0; p < n_uv; p++) gmt_M_free (GMT, X[p]); in GMT_gpsgridder()