Форум программистов, компьютерный форум, киберфорум
Наши страницы

Исправить программу для спрямления (гиперболических) кривых - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Определить прямоугольник в плоскости, в который попадают все указанные точки http://www.cyberforum.ru/cpp/thread1709844.html
Таблица содержит геометрические точки, заданные в полярных координатах (R ,α) . Определить прямоугольник в плоскости Oxy , в который попадают все указанные точки x=R⋅cosα , y=R⋅sinα . Стороны...
C++ Шахматный ИИ Ищу примеры реализации поиска самого выгодного хода для двухмерного массива интов 8 на 8. В моем случае: 0 - пустая клетка, 1 - Король, 2 - Ферзь, 3 - Офицер, 4 - Конь, 5 - Ладья, 6 - Пешка ... http://www.cyberforum.ru/cpp/thread1708996.html
C++ Метод Монте-Карло для систем многих частиц
Есть задача реализовать метод Монте-Карло для систем многих частиц. На данный момент сделан расчет взаимодействия двух частиц, а также граничные условия. Далее надо 1) реализовать схему...
Компилятор g++ выдает ошибку: "Версия этого файла несовместима с используемой версией Windows" C++
Когда запускаю программу с указателями, то виндовс ее сразу закрывает. Почему и как исправить? Все на скришоте
C++ Проверка скорости кода. Обмен опытом http://www.cyberforum.ru/cpp/thread1707297.html
Вводные данные: - C++ стандарта 11 - gcc Работаю над ускорением кода для работы со строками. К примеру, сравниваю между собой скорость работы: -strncpy -strcpy -через ассемблер movs ...
C++ Реализации преобразований с геометрической фигурой (сферой) Товарищи подскажите кто чем сможет. Надо написать программу, позволяющую делать различные преобразования со сферой (перемещение, сжатие, растяжение, изменение цвета, освещение и т. д.); управление... подробнее

Показать сообщение отдельно
Kerim_Geophysic
0 / 0 / 0
Регистрация: 08.01.2016
Сообщений: 137

Исправить программу для спрямления (гиперболических) кривых - C++

13.04.2016, 03:17. Просмотров 225. Ответов 0
Метки (Все метки)

впервые встречаюсь с необходимостью что-то изменить в программе вроде C++
немного знаком с Матлабом
значит есть программа, ее описание есть здесь http://ahay.org/blog/2013/04/08/program-of-the-month-sfnmo/
суть в том, что есть три формулы, которые можно найти по ссылке выше. По умолчанию используется первая формула.
в этой формуле to- значение по вертикальной оси на рисунке (время), х- значение по горизонтали (расстояние).
входным параметрами являются to(время) и V(скорость).
программа написана так, что to(время) и V(скорость) надо задавать в файле, который непонятным для меня образом формируется. Я хочу задавать эти параметры двумя столбцами: первый- to(время), второй: V. То есть каждому времени соответствует скорость.
В приведенном ниже коде этот входной файл именуется как velocity
код программы:
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
#include <math.h>
 
#include <rsf.h>
 
#include "fint1.h"
 
static float *vel, *par, h;
static bool slow, squared;
 
static float nmo_map(float t, int it) {
    float v;
    
    v = vel[it];
    if (!squared) v *= v;
    v = slow ? h*v : h/v;
    t = t*t + v;
    if (t > 0.) t=sqrtf(t);
    return t;
}
 
static float shifted_nmo_map(float t, int it) {
    float v, q, s;
    
    v = vel[it];
    s = par[it];
 
    if (!squared) v *= v;
    v = slow ? h*v : h/v;
    q = t*t + s*v;
    if (q > 0.) {
    t=t*(1.-1./s)+sqrtf(q)/s;
    return t;
    } else {
    return q;
    }
}
 
static float taner_nmo_map(float t, int it) {
    float v, a;
    
    v = vel[it];
    a = par[it];
 
    if (!squared) v *= v;
    v = slow ? h*v : h/v;
    t = t*t + v/(1+a*v);
    if (t > 0.) t=sqrtf(t);
    return t;
}
 
int main (int argc, char* argv[])
{
    fint1 nmo;
    bool half;
    int ix,ih, nt,nx,nw, nh, noff, nmask, CDPtype, mute, *mask;
    float dt, t0, h0, dh, dy, str;
    float *trace, *off;
    mapfunc nmofunc;
    sf_file cmp, nmod, velocity, offset, msk, het;
 
    sf_init (argc,argv);
    cmp      = sf_input("in");
    velocity = sf_input("velocity");
    nmod     = sf_output("out");
 
    if (NULL != sf_getstring("s")) {
        het = sf_input("s");
    nmofunc = shifted_nmo_map;  
    } else if (NULL != sf_getstring("a")) {
        het = sf_input("a");
        nmofunc = taner_nmo_map;
    } else {
    het = NULL;
    nmofunc =  nmo_map;
    } 
 
    if (SF_FLOAT != sf_gettype(cmp)) sf_error("Need float input");
    if (!sf_histint  (cmp,"n1",&nt)) sf_error("No n1= in input");
    if (!sf_histfloat(cmp,"d1",&dt)) sf_error("No d1= in input");
    if (!sf_histfloat(cmp,"o1",&t0)) sf_error("No o1= in input");
    if (!sf_histint  (cmp,"n2",&nh)) sf_error("No n2= in input");
    nx = sf_leftsize(cmp,2);
 
    if (!sf_getbool("half",&half)) half=true;
    /* if y, the second axis is half-offset instead of full offset */
    if (!sf_getfloat("str",&str)) str=0.5;
    /* maximum stretch allowed */
 
    if (!sf_getint("mute",&mute)) mute=12;
    /* mute zone */
 
    CDPtype=1;
    off = sf_floatalloc(nh);
 
    if (NULL != sf_getstring("offset")) {
    offset = sf_input("offset");
    if (SF_FLOAT != sf_gettype(offset)) sf_error("Need float offset");
    noff = sf_filesize(offset);
    if (noff == nh) {
        sf_floatread (off,nh,offset);
    } else if (noff != nh*nx) {
        sf_error("Wrong dimensions in offset");
    }
    } else {
    if (!sf_histfloat(cmp,"d2",&dh)) sf_error("No d2= in input");
    if (!sf_histfloat(cmp,"o2",&h0)) sf_error("No o2= in input");
 
    if (sf_histfloat(cmp,"d3",&dy) && !sf_getint("CDPtype",&CDPtype)) {
        CDPtype=half? 0.5+dh/dy : 0.5+0.5*dh/dy;
        if (CDPtype < 1) {
        CDPtype=1;
        } else if (1 != CDPtype) {
        sf_histint(cmp,"CDPtype",&CDPtype);
            sf_warning("CDPtype=%d",CDPtype);
        }
    }       
 
    for (ih = 0; ih < nh; ih++) {
        off[ih] = h0 + ih*dh; 
    }
 
    noff = nh;
    offset = NULL;
    }
    
    if (NULL != sf_getstring("mask")) {
    msk = sf_input("mask");
 
    if (SF_INT != sf_gettype(msk)) sf_error("Need integer mask");
    nmask = sf_filesize(msk);
    mask = sf_intalloc(nh);
 
    if (nmask == nh) {
        sf_intread (mask,nh,msk);
    } else if (nmask != nh*nx) {
        sf_error("Wrong dimensions in mask");
    }
    } else {
    nmask = nh;
 
    msk = NULL;
    mask = NULL;
    }
 
    if (!sf_getbool("slowness",&slow)) slow=false;
    /* if y, use slowness instead of velocity */
 
    if (!sf_getbool("squared",&squared)) squared=false;
    /* if y, the slowness or velocity is squared */
 
    if (!sf_getfloat ("h0",&h0)) h0=0.;
    /* reference offset */
    if (half) h0 *= 2.;
    if (!sf_getint("extend",&nw)) nw=4;
    /* trace extension */
 
    trace = sf_floatalloc(nt);
    vel   = sf_floatalloc(nt);
    if (NULL != het) {
    par = sf_floatalloc(nt);
    } else {
    par = NULL;
    }
 
    nmo = fint1_init (nw, nt, mute);
    
    for (ix = 0; ix < nx; ix++) {
    sf_warning("CMP %d of %d;",ix+1,nx);
 
    sf_floatread (vel,nt,velocity); 
    if (NULL != het) sf_floatread(par,nt,het);
    if (NULL != offset && noff != nh) sf_floatread (off,nh,offset);
    if (NULL != msk && nmask != nh) sf_intread (mask,nh,msk);
 
    for (ih = 0; ih < nh; ih++) {
        sf_floatread (trace,nt,cmp);
        
        /* skip dead traces */
        if (NULL != msk && 0==mask[ih]) {
        sf_floatwrite (trace,nt,nmod);
        continue;
        }
        
        fint1_set(nmo,trace);
 
        h = off[ih];
        if (NULL == offset) h += (dh/CDPtype)*(ix%CDPtype); 
        if (half) h *= 2;
        h = h*h - h0*h0;
        
        stretch(nmo,nmofunc,nt,dt,t0,nt,dt,t0,trace,str);
        sf_floatwrite (trace,nt,nmod);
    }
    }
    sf_warning(".");
 
    exit (0);
}
 
/*  $Id$     */
как мне изменить этот файл, чтобы velocity задавался в файле в котором я хочу? (см прикрепленные файлы, pick.rsf- как я хочу, picks.rsf-как есть на самом деле)
0
Вложения
Тип файла: rar demo.rar (5.7 Кб, 2 просмотров)
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru