demo for polarization ordering

extensions/demo/other/init/polarizationField.dat.2d

+dimension of world: 2
+
+% =================== MESH ================================
+mesh_name: mesh
+mesh->macro file name:                   ./macro/macro.square.2d
+mesh->global refinements:                10
+mesh->check:                             0
+
+
+% ============== USER-PARAMETER ==========================
+p->alpha2: 1
+p->alpha4: 1
+p->epsilon: 0.05
+p->K:      (${p->epsilon}) * 1
+
+% =========== OUTPUT ==============================================
+p->space->output->filename: ./output/p_
+p->vectorField->output->filename: ./output/vec_
+
+% ==================== TIMESTEPS ===============================
+adapt->timestep:      0.00001
+adapt->max timestep:  1e+1
+adapt->min timestep:  1e-6
+adapt->start time:    0.0
+adapt->end time:      10000.0
+
+p->space->rosenbrock->method: ros3Pw
+
+p->space->rosenbrock->error weights: [1,1,0,0]
+
+adapt->fix first timesteps: 10
+adapt->rosenbrock->timestep study: 0
+adapt->rosenbrock->timestep study steps: 0
+
+adapt[0]->time tolerance: 1.e-2
+
+% ============= PROBLEM-SPACES ==================================
+p->space->components:           4
+p->space->polynomial degree[0]: 2
+p->space->polynomial degree[1]: 2
+p->space->polynomial degree[2]: 2
+p->space->polynomial degree[3]: 2
+p->space->dim:                  2
+p->space->mesh:                 mesh
+
+
+% ================== SOLVER ======================================
+p->space->solver->backend: mtl
+p->space->solver: umfpack
+p->space->solver->use old initial guess: 1
+p->space->solver->tolerance:          1.e-8
+p->space->solver->info:               10
+
+% =================== OUTPUT =========================================
+p->space->output->ParaView animation:        1
+p->space->output->ParaView format:           1
+p->space->output->write every i-th timestep: 10
+p->space->output->append index:              1
+p->space->output->index length:              9
+p->space->output->index decimals:            7
+
+p->vectorField->output->ParaView animation:        1
+p->vectorField->output->ParaView vector format:    1
+p->vectorField->output->write as 3d vector:        1
+p->vectorField->output->write every i-th timestep: 10
+p->vectorField->output->append index:              1
+p->vectorField->output->index length:              9
+p->vectorField->output->index decimals:            7