MKOCTFILE ?= mkoctfile
FFC ?= ffc

OCTFILES= Mesh.oct \
             DirichletBC.oct \
             Expression.oct \
             Function.oct \
             save.oct \
             assemble.oct \
             assemble_system.oct \
             plot_func.oct \
             plot_mesh.oct \
             SubSpace.oct \
             feval.oct 


LIBS += -ldolfin

all: $(OCTFILES)

fem_init_env.oct: fem_init_env.o
	$(MKOCTFILE) $(CPPFLAGS) fem_init_env.o -o $@ $(LDFLAGS) $(LIBS)

fem_init_env.o:  fem_init_env.cc mesh.h functionspace.h boundarycondition.h function.h coefficient.h form.h
	$(MKOCTFILE) $(CPPFLAGS) -c fem_init_env.cc $(LDFLAGS) -o $@ -I.

Mesh.oct: Mesh.o
	$(MKOCTFILE) $(CPPFLAGS) Mesh.o -o $@ $(LDFLAGS) $(LIBS)

Mesh.o: Mesh.cc mesh.h
	$(MKOCTFILE) $(CPPFLAGS) -c Mesh.cc $(LDFLAGS) -o $@ -I.

fem_get_mesh.oct: fem_get_mesh.o
	$(MKOCTFILE) $(CPPFLAGS) fem_get_mesh.o -o $@ $(LDFLAGS) $(LIBS)

fem_get_mesh.o: fem_get_mesh.cc
	$(MKOCTFILE) $(CPPFLAGS) -c fem_get_mesh.cc $(LDFLAGS) -o $@ -I.

DirichletBC.oct: DirichletBC.o
	$(MKOCTFILE) $(CPPFLAGS) DirichletBC.o -o $@ $(LDFLAGS) $(LIBS)

DirichletBC.o: DirichletBC.cc functionspace.h boundarycondition.h expression.h
	$(MKOCTFILE) $(CPPFLAGS) -c DirichletBC.cc $(LDFLAGS) -o $@ -I.

Expression.oct: Expression.o
	$(MKOCTFILE) $(CPPFLAGS) Expression.o -o $@ $(LDFLAGS) $(LIBS)

Expression.o: Expression.cc expression.h coefficient.h
	$(MKOCTFILE) $(CPPFLAGS) -c Expression.cc $(LDFLAGS) -o $@ -I.

Function.oct: Function.o
	$(MKOCTFILE) $(CPPFLAGS) Function.o -o $@ $(LDFLAGS) $(LIBS)

Function.o: Function.cc function.h
	$(MKOCTFILE) $(CPPFLAGS) -c Function.cc $(LDFLAGS) -o $@ -I.

assemble.oct: assemble.o
	$(MKOCTFILE) $(CPPFLAGS) assemble.o -o $@ $(LDFLAGS) $(LIBS)

assemble.o: assemble.cc form.h boundarycondition.h
	$(MKOCTFILE) $(CPPFLAGS) -c assemble.cc $(LDFLAGS) -o $@ -I.

assemble_system.oct: assemble_system.o
	$(MKOCTFILE) $(CPPFLAGS) assemble_system.o -o $@ $(LDFLAGS) $(LIBS)

assemble_system.o: assemble_system.cc form.h boundarycondition.h
	$(MKOCTFILE) $(CPPFLAGS) -c assemble_system.cc $(LDFLAGS) -o $@ -I.

save.oct: save.o mkfunction
	$(MKOCTFILE) $(CPPFLAGS) save.o -o ./@function/save.oct $(LDFLAGS) $(LIBS)

save.o: save.cc 
	$(MKOCTFILE) $(CPPFLAGS) -c save.cc $(LDFLAGS) -o $@ -I.

mkfunction:
	 mkdir -p @function

mkmesh:
	 mkdir -p @mesh

plot_mesh.oct: plot_mesh.o mkmesh
	$(MKOCTFILE) $(CPPFLAGS) plot_mesh.o -o ./@mesh/plot.oct $(LDFLAGS) $(LIBS)

plot_mesh.o: plot_mesh.cc Plot_2d.h mesh.h Plot_3d.h
	$(MKOCTFILE) $(CPPFLAGS) -c plot_mesh.cc $(LDFLAGS) -o $@ -I.

Plot_2d.h: Plot_2d.ufl
	$(FFC) -l dolfin Plot_2d.ufl

Plot_3d.h: Plot_3d.ufl
	$(FFC) -l dolfin Plot_3d.ufl

plot_func.oct: plot_func.o mkfunction
	$(MKOCTFILE) $(CPPFLAGS) plot_func.o -o ./@function/plot.oct $(LDFLAGS) $(LIBS)

plot_func.o: plot_func.cc Plot_2d.h mesh.h Plot_3d.h
	$(MKOCTFILE) $(CPPFLAGS) -c plot_func.cc $(LDFLAGS) -o $@ -I.

SubSpace.oct: SubSpace.cc functionspace.h
	$(MKOCTFILE) $(CPPFLAGS) -I. SubSpace.cc $(LDFLAGS) $(LIBS)

feval.oct: feval.o mkfunction
	$(MKOCTFILE) $(CPPFLAGS) feval.o -o ./@function/feval.oct $(LDFLAGS) $(LIBS)

feval.o: feval.cc function.h
	$(MKOCTFILE) $(CPPFLAGS) -c feval.cc $(LDFLAGS) -o $@ -I.

clean:
	-rm -f *.o core octave-core *.oct *~ *.xml

cleanall:
	-rm -f *.o core octave-core *.oct *~ *.xml *.status *.log \
          octave-workspace configure *.pvd *.vtu
	-rm -r autom4te.cache
	-rm -f Makefile