/*
sample1.c

make EEM-STF data
*/

#include "stf_datalib.h"

int main(void)
{
	// (1) initialize

	stf_init();

	// (2) title

	stf_title("sample1");

	// (3) mesh

	stf_xsection(2, -5.0, +5.0);
	stf_xdivision(1, 50);

	stf_ysection(2, -5.0, +5.0);
	stf_ydivision(1, 50);

	stf_zsection(2, -5.0, +5.0);
	stf_zdivision(1, 50);

	// (4) material

	stf_epsr(2.0);

	// (5) volt

	stf_volt(-1);
	stf_volt(+1);

	// (6) geometry

	stf_unit6(1, 1, 1, -2.0, -2.0, -1.0, +2.0, +2.0, +1.0);		// dielectric
	stf_unit6(2, 1, 1, -2.0, -2.0, -1.0, +2.0, +2.0, -1.0);		// - electrode
	stf_unit6(2, 2, 1, -2.0, -2.0, +1.0, +2.0, +2.0, +1.0);		// + electrode

	// (7) charge

	// (8) solver

	stf_solver(1, 1e-5, 1000, 20, 0.9);

	// (9) output

	stf_outdat("sample1.stf");

	return 0;
}