2005-01-06  julien pommier  <pommier@gmm.insa-tlse.fr>

	* minor changes for MinGW (getfem-1.7.1)

2005-01-03
	* Getfem-1.7 released.

2004-12-10  Yves Renard  <Yves.Renard@insa-toulouse.fr>

	* src/getfem_mat_elem_type.h adding the unit normal to a boundary
	available for elementary computations.

2004-12-08  julien pommier  <pommier@gmm.insa-tlse.fr>

	* superlu/BLAS.c configure.in: added a f2c-ed BLAS, used when no
	other suitable blas is found.

2004-12-02  julien pommier  <pommier@gmm.insa-tlse.fr>

	* src/bgeot_geotrans_inv.C (bgeot): Made inversion geometric
	transformations more robust (linear search in the Newton
	algorithm).

2005-01-05  Yves Renard  <Yves.Renard@insa-toulouse.fr>

	* src/getfem_linearized_plate.h linearized plate brick

2004-11-25  Yves Renard  <Yves.Renard@insa-toulouse.fr>

	* src/getfem_linearized_plate.h linearized plate brick

2004-11-03  Yves Renard  <Yves.Renard@insa-toulouse.fr>

	* src/getfem_fem_sum.h introduction of a fem which is the direct
	sum of two fems..

2004-11-03  Yves Renard  <Yves.Renard@insa-toulouse.fr>

	* src/getfem_fem_virtual_link.h replaced by
	getfem_interpolated_fem.h redefinition of the interpolation of a
	fem with real element fems.

2004-10-19  Yves Renard  <Yves.Renard@insa-toulouse.fr>

	* src/gmm_precond_ilutp.h: New preconditionner, incomplete LU with
	  fill-in, threshold and column pivoting.

2004-10-15  Yves Renard  <Yves.Renard@insa-toulouse.fr>

	* src/gmm_vector.h: bug fixed and optimization of the copy toward
	a rsvector<T> : the unsorted indices were not taken into account.

2004-09-15  Yves Renard  <Yves.Renard@insa-toulouse.fr>

	* src/gmm_solver_bfgs.h: A BFGS and a DFP algorithm for
	unconstrained optimization.

2004-09-15  Julien Pommier  <pommier@insa-toulouse.fr>

	* src/gmm_blas.h: 
	gmm now handles (valid) operations mixing complex and scalars
	(e.g. multiplication of a complex vector with a real matrix)

2004-08-27  julien pommier  <pommier@gmm.insa-tlse.fr>

	* src/getfem_fem.h, src/getfem_interpolation.h: added hability to
	interpolate complex-valued fields.

2004-08-09  julien pommier  <pommier@gmm.insa-tlse.fr>

	* src/getfem_assembling_tensors.C, src/bgeot_sparse_tensors.C,
	src/getfem_mat_elem.C: Added the ability to perform reductions
	during the integration of elementary tensors, and daxpy is used
	when possible: huge speedup for the assembly of nonlinear tangent
	matrices.

2004-08-03  julien pommier  <pommier@gmm.insa-tlse.fr>
	* superlu: SuperLU is included into getfem, which does now depend
	on BLAS.

2004-08-02  Yves Renard  <Yves.Renard@gmm.insa-tlse.fr>
	* src/gmm_dense_qr.h: Bug fixed on the computation of Schur vectors.
	the reductible 2x2 blocks has to be reduced to obtain the Schur
	vectors (this is not necessary just to obtain eigenvalues).
	(Bug report from Kasper Peeters <kasper.peeters@aei.mpg.de>).

2004-07-20  julien pommier  <pommier@gmm.insa-tlse.fr>
	* src/getfem_fem.h etc.: interpolation_grad now stores the
	gradients in the rows of the matrix, no more in the columns.
	
2004-07-05  julien pommier  <pommier@gmm.insa-tlse.fr>
	* src/getfem_export.C: added VTK and OpenDX export of slices,
	getfem_mesh and mesh_fem objects.

2004-06-25 
	* src/gmm_real_part.h : added the functions gmm::real_part and
	gmm::imag_part.

2004-06-18  Yves Renard  <Yves.Renard@gmm.insa-tlse.fr>
	* src/getfem_fem.h : Now a dof has to be on the boundary if
	his corresponding shape function is non-zero on that boundary
	(a P0 node is on each boundary).

2004-06-18  Yves Renard  <Yves.Renard@gmm.insa-tlse.fr>
	* src/getfem_context.h : New structure to deal with interdependency
	of objects.

2004-06-18  Yves Renard  <Yves.Renard@gmm.insa-tlse.fr>
	* src/getfem_modeling.h : New structures defining elementary
	bricks to build complete models: linearized elasticity, boundary
	conditions ...

2004-06-17  julien pommier  <pommier@gmm.insa-tlse.fr> 
        * src/gmm_matrix.h: fixed the swap() method for csc_matrix and 
        csr_matrix. 

2004-06-04  Yves Renard  <Yves.Renard@gmm.insa-tlse.fr>
	* src/gmm_sub_index.h : the sturcture unsorted_sub_index is
	defined to represent a non-sorted sub-vector or sub-matrix, for
	instance a permutation of a vector.
	All linalg_traits have been modified, a new component called
	index_sorted has been added. The algorithms in gmm_blas.h have been
	modified to take into account the fact that a sparse index can be
	unsorted.

2004-06-01  Yves Renard  <Yves.Renard@gmm.insa-tlse.fr>
	* src/getfem_export.h : the interpolation function has been moved
	in the new file src/getfem_interpolation.h. An option for the
	extrapolation of points just outside the boundary has bee added.
	the function interpolation_solution is now just called interpolation.

2004-05-28  julien pommier  <pommier@gmm.insa-tlse.fr>

	* src/getfem_assembling.h (getfem): minor speedup for the assembly
	of linear elasticity.
	
	* src/bgeot_small_vector.h : small_vectors are now initialized
	with zeros.

	* src/dal_singleton.h, *.C: added ability to cleanup allocated
	static data on exit -- mostly for valgrind --leak-check=yes.
	
2004-05-14  Yves Renard  <Yves.Renard@gmm.insa-tlse.fr>

	* src/gmm_solvers_Schwarz_additive.h: A sequential version of
	the domain decomposition / precnditionner is finalized. Gmres,
	Bicgstag or CG can be used as local or global solvers. Superlu
	can be used as local solver.
	
2004-05-14  julien pommier  <pommier@gmm.insa-tlse.fr>

	* src/gmm_solver_gmres.h: fixed an incorrect computation of the
	residu in gmres in cases where a preconditioner was used.

	* src/gmm_vector.h: fixed a bug in the constructor of slvector
	(thanks to li zhipeng).
	
2004-04-02  julien pommier  <pommier@gmm.insa-tlse.fr>

	* tests/elastostatic.C tests/laplacian.C: now take into account
	the GENERIC_DIRICHLET parameter

2004-03-31  Yves Renard  <Yves.Renard@gmm.insa-tlse.fr>

	* src/getfem_fem.h and src/getfem_import.h: added a method giving
	an interpolation matrix, useful for repeated interpolations of one
	field onto another.

2004-03-31  Yves Renard  <Yves.Renard@gmm.insa-tlse.fr>

	* src/gmm_superlu_interface.h: The superlu interface is now
	uptodate with the lastest release of superlu (v.3.0).  A new
	structure was introduced for storing an LU decomposition and doing
	tri_solves later.

2004-03-26  Yves Renard  <Yves.Renard@gmm.insa-tlse.fr>

	* src/gmm_precond_ildltt.h: fix bug which slowed the ildltt from
	a factor of about 50.

2004-03-11  julien pommier  <pommier@gmm.insa-tlse.fr>

	* src/getfem_fem.C (getfem): added a classical_discontinuous_fem
	similar to the classical_fem already available.

	* src/getfem_integration.C (getfem): added the classical_approx_im
	function, which returns an approx. integration method for a given
	{geometric_trans,degree}. Renamed exact_classical_im to
	classical_exact_im.
	
2004-03-09  julien pommier  <pommier@gmm.insa-tlse.fr>

	* src/bgeot_geometric_trans.C (bgeot): added compute_local_basis
	which returns the normal of a face + tangents

2004-02-26  julien pommier  <pommier@gmm.insa-tlse.fr>

	* src/getfem_fem.C (getfem): correction the wrong quadrature rules
	for IM_GAUSSLOBATTO1D. Addition of a new fem:
	FEM_PK_GAUSSLOBATTO1D. Integration points are now sorted.

2004-02-09  julien pommier  <pommier@gmm.insa-tlse.fr>

	* src/getfem_fem.[Ch] src/getfem_precomp.[Ch],
	src/bgeot_precomp.[Ch], etc: large modifications of the fem
	interpolation interface. Now uses fem_interpolation_context (and
	geotrans_interpolation_context). geotrans_precomp is now part of
	namespace bgeot.

2004-02-04  julien pommier  <pommier@gmm.insa-tlse.fr>

	* src/getfem_precomp.[Ch] (getfem): cleaned up getfem_precomp,
	replaced the ugly fem_precomp_not_stored and
	geotrans_precomp_not_stored with the much better fem_precomp_pool
	and geotrans_precomp_pool.

	* src/getfem_Xfem.C (getfem): small perf improvement (better use
	of pgeotrans_precomp).

2004-02-02  julien pommier  <pommier@gmm.insa-tlse.fr>

	* src/getfem_mesh_slice.[Ch] src/getfem_mesh_slicers.[Ch]: the
	mesh_slices have been splitted in two parts: 

	- the mesh_slicer object, which takes a mesh, simplexifies it and
	sends thes simplexes to a list of "slicer_action" objects

	- the stored_mesh_slice object, which contains all simplexes that
	have been constructed by the mesh_slicer. This recording of
	simplexes is done via the slicer_build_stored_mesh_slice (derived
	from slicer_action) class.

	src/getfem_export.h: interpolation to a set of sparse points can 
	now interpolate the gradient of a field.
	
2004-01-29  julien pommier  <pommier@gmm.insa-tlse.fr>

	* src/getfem_export.h: removed some obsolete interpolation functions. 

2004-01-27  Yves Renard <renard@gmm.insa-tlse.fr>
	
	* tests/elastostatic.C: simplification, copied from laplacian.C 
	(elastostatique -> elastostatic)
	
2004-01-27  julien pommier  <pommier@gmm.insa-tlse.fr>

	* tests/laplacian.C: huge simplification of the sources..

	* removed obsolete getfem_norm.h

	* src/getfem_mesh_slice.C: fixed the problem (missing simplexes)
	with 'union' operation on slice. Moreover, the slice boundary
	flagging is now better than ever. Added a bounding_box member
	function in geometric_trans.

2004-01-23
	* Getfem-1.6 released. The ChangeLog is started!
