# a bit of history, code developed before pde_spike
#                   just the C code, also Fortran, Ada, Java, ...


# generate spike.dat from sphere
sphere_div.out: sphere_div.c datread.c datread.h
	gcc -o sphere_div sphere_div.c datread.c -lm
	./sphere_div spike.dat 10 > sphere_div.out


# convert spike.dat to spike.inp UCD format
make_spike_inp_c.out: make_spike_inp.c datread.h datread.c
	gcc -o make_spike_inp make_spike_inp.c datread.c -lm
	./make_spike_inp spike.dat > make_spike_inp_c.out


# test  pde_read_ucd 
test_pde_read_ucd_c.out: test_pde_read_ucd.c pde_read_ucd.h pde_read_ucd.c \
                         point_in_poly.c point_in_poly.h simeq.c simeq.h \
                         spike.inp
	gcc -o test_pde_read_ucd test_pde_read_ucd.c pde_read_ucd.c \
                                 point_in_poly.c simeq.c -lm
	./test_pde_read_ucd  spike.inp > test_pde_read_ucd_c.out


# test  pde_inside3_ucd
test_pde_inside3_ucd_c.out: test_pde_inside3_ucd.c pde_read_ucd.h \
                            pde_read_ucd.c simeq.h simeq.c \
                            point_in_poly.c point_in_poly.h spike.inp
	gcc -g3 -o test_pde_inside3_ucd test_pde_inside3_ucd.c \
                            pde_read_ucd.c simeq.c point_in_poly.c -lm
	./test_pde_inside3_ucd spike.inp  > test_pde_inside3_ucd_c.out


# test  lsfit
test_lsfit_c.out: test_lsfit.c lsfit.c lsfit.h
	gcc -o test_lsfit test_lsfit.c lsfit.c -lm
	./test_lsfit > test_lsfit_c.out


# test  nuderiv3dg
test_nuderiv3dg_c.out: test_nuderiv3dg.c nuderiv3dg.c nuderiv3dg.h
        gcc -g3 -o test_nuderiv3dg test_nuderiv3dg.c nuderiv3dg.c  -lm
        test_nuderiv3dg > test_nuderiv3dg_c.out


# run  pde_spike  solving for each unknown independently
pde_spike_c.out: pde_spike.c pde_read_ucd.h pde_read_ucd.c \
                 lsfit.h lsfit.c simeq.h simeq.c \
                 point_in_poly.c point_in_poly.h \
                 nuderiv3dg.h nuderiv3dg.c spike.inp
	gcc -g3 -o pde_spike pde_spike.c pde_read_ucd.c lsfit.c \
                          simeq.c nuderiv3dg.c point_in_poly.c -lm
	./pde_spike > pde_spike_c.out