#include "def21060.h"
#define N 167
#define c 1.0
#define c_star 1.0
#define uc 1.0 //1/c

.SECTION/DM seg_dmda;
//dmdata
.VAR input[N] = "test.txt";
.VAR mas1[N];

.SECTION/PM seg_pmda;
.VAR out[N];
.VAR nabla = 0.0;
//pmdata


.SECTION/PM seg_int_code;
	nop;
	jump start;
	nop;
	nop;


.SECTION/PM seg_pmco;
//pmcode
start:
	I0 = input;
	I1 = mas1;
	I8 = out;
	i10 = nabla;
	f0 = 0.0; //Yi
	f1 = 0.0;//yi
	f2 = 0.0; //zi
	f3 = 0.0; //dzi
	f4 = 0.0; //zi1
	F12 = 0; //F
	f15 = 0.0; //äëÿóñëîâèÿ!!!!
	f13 = 0.0;//äåëüòà
	f5 = 0.0;


	LCNTR = N, DO end UNTIL LCE;

	f1 = dm(i0, 1);
	f2 = f0 - f1; //âû÷èñëåíèå zi
	f3 = f2 - f4; //dzi = zi -zi+1
	f4 = f2;

	f6 = f3 * f3; //dz*dz
	f10 = 1.5;
	f7 = f10 * f3;//1.5dz
	f10 = uc;
	f6 = f6 * f10;//dz*dz/c
	f10 = 0.125;
	f11 = c;
	f10 = f10 * f11;
	f6 = f6 - f10;//dz*dz/c-0.125c
	f12 = f2 + f7; //zi+1.5dz
	f12 = f12 + f6;//tot samii fi

	f14 = 1.0;
	comp(f3, f15); //znak esli bolshe to 1
	if GT jump xxx;
	f14 = -1.0;
xxx: f8 = -1.0;
    f14 = f14*f8;
	f12 = f12 * f14; //f ended


	f13 = 1.0;
	comp(f12, f15); //delta -sign
	if GT jump uuu;
	f13 = -1.0;

uuu: f8 = -1.0;
	f13 = f8*f13;


zzz:
	DM(i1, 1) = f13;
	f10 = c_star;//áóäåò 2íàáëà
	f10 = f10 * f13;//â ô10 2íàáëà èëè íàáëà 2 õç êîðî÷

	f5 = f5 + f10; //íàáëà íàáëà íà ïðåäûäóùåì + 2 íàáëà



end: f0 = f0 + f5;//Yi count

	nop;nop;F1 = 0; //íàáëà 2
	F2 = 0; //íàáëà
	f0 = 0; //Y
	i4 = mas1;
	LCNTR = N, DO aaa UNTIL LCE;

		f11 = dm(i4, 1); //deltas
		f10 = c_star;
		f1 = f10 * f11;
		f2 = f2 + f1;
		f0 = f0 + f2;

aaa:	PM(i8, 1) = f0;