This repository has been archived by the owner on Jul 30, 2018. It is now read-only.
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathINNOVATION
1437 lines (1333 loc) · 66.5 KB
/
INNOVATION
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
- add spline interface of gsl library in gsl plugin see
gsl.edp examples.
version 3.47 (1/june/2016)
- tools to build new FiniteElement2d and FiniteElement3d from quadrature formula
FiniteElement2d EFQF(qf23pT);
fespace Zh(Th2d,EFQF);
FiniteElement3d EFQF3(qfVp14);
fespace Zh(Th3d,EFQF3);
- add solveurs FETI and BDD in HPDDM
- add savegmsh in gmsh plug in to write file for gmsh (thank to Lo Sala [email protected])
see ttestio.edp for the examples
version 3.46
version 3.45-1
- add new finite element to store value a quadradure point.
in plugin Element_QF (see exemples++load/Element_QF.edp)
in 2d FEQF1,FEQF2,FEQF5,FEQF7,FEQF9 corresponding resp. to quadrature
qf1pT,qf2pT,qf5pT,qf7pT,qf9pT
in 3d FEQF13d,FEQF23d,FEQF53d corresponding resp. to quadrature
qfV1 ,qfV2 , qfV5
and FEQF and FEQF3d for the default quadrature.
- ajout gluemesh of 3d mesh
- ajout operator b= A* c; // where is real matrix avec b and c complexe vector.
version 3.45 (10 march 2016)
- correct veru old bug in ffglut, some time we miss wait due to asyncronuos stuff
and we can hav e memory error for the same reason .
- try use lot of preinstall software (hard job for ubuntu and debian)
- pass to tetgen version 1.5.1-beta (the last one)
- add Th.hmin and Th.hmax methode to Mesh type
version 3.44-1
- add new hpddm example for Stokes probleme (Thank to P. Jolivet )
- put the new convect 3d as the default, (newconvect do nothing now). 25/02/16
- try to add new version do 3d convect ( add newconvect=1;) in test ...
version 3.44 (22 feb. 2016)
- add tools to extract border after movemesh in Curvature plugin)
see curvature.edp, can be useful to remesh after movemesh.
- restruct the hpddm and petsc plugin ( thanks to P. Jolivet)
add lots hpddm and petsc examples in examples++-hpddm
version 3.43-33
- remove some spurious print when verbosity = 0
- add initialisation of array which string parameter like in
string[string] b= ["-1",1,"13","qsdqdq"];
int[string] a=["2",1,"12",4];
- add setw iomap stuff like in:
cout << setw(10) << 1 << " " << 3;
version 3.34-2
- correct examples++-load/ff-get-dep.awk in case of choose of lib in syntaxe of library dependes in
WHERE-LIBRARY stuff
like in :
//ff-c++-LIBRARY-dep: cxx11 [petsc|mumps parmetis ptscotch scotch] scalapack blas [mkl] mpifc fc mpi pthread
version 3.43-1
- correct PB with hpddm (pass to git version d2a99d89ee37777d2423da218b29e2080950da1e Jan 20 22:03:12 2016)
- correct PB of FLIBS
- pass to evrsion 3.43-1 3 feb 2015
version 3.43 (version 3.43 (3626:8e3d8953975d), 2 feb 2016 16;00 CET))
- change RT03d finite element interpolation (put 1 point integration by face)
- add renum=false in all trunc fonction of DDM examples
- warning change the defaut no renumbering to be sure trunc commute
the two mesh must be the same truc(trunc(Th,I),J) == truc(trunc(Th,J),I)
remove this version (version 3.43 (3623:882bfbae1bf6, feb 2016 15:30 CET))
- change the trunc function in 2d to keep the ordre of old vertices to remove signe probleme
in some restriction operator of RT finite element for example
- correct of the formal operator trace (bug was comute a00 + trace(A), thank to O. Pironneau )
- put hppdm source include in the distrubition to have more safe compile process.
- update a little the documentation (add new finite element)
- add new finite elemnt just for interpolate partition of unity in DDM / HPDDM method
Edge23ds0, Edge13ds0,Edge03d in plugin Element_Mixte3d
- add new finite element P1bl, P1bl3d in plugin Element_P1bl
- add new finite element Edge23d, P1bl,in plugin Element_Mixte3d
- change EPSD in eigenv because touble in some very rare case (see G. Vergez)
- correct problem for vectoral solution on implicity declaration of product of simple EF like P1b
- add tools to compute curvature of border in Cuvature plugin
- all implicity loop for lot of array : real[ini], real[int,int], real[string] ,matrix , string[string]
the syntaxe is
if B is like a Array :
for [i,bi: A] bi=f(i);
if A is like a matrix :
for [i,j,aij: A] aij=g(i,j);
see forall.edp example file for more exemple
- change the chinise documentation (thank to Helin GONG)
version 3.42
- correct some Makefile
- correct optimisation of convect in case on windows bug)
- add tool do build quadrature formule when we split triangle in 3 triangle for HCT element.
- add new quadrature formula 3d exact form degre 6 to 14 in qf11to25.cpp
the name of qfV quadrature are qfVp1 .. qfVp14
- correct optimisation in convect for navier stokes problem..
for this kind of expression
+ int2d(Th) ( -alpha*convect([up1,up2],-dt,up1)*v1 -alpha*convect([up1,up2],-dt,up2)*v2 )
in this case now only on characteristics is computed
- correct trouble in configure this /usr/bin/machine
- map&key to know if a key exist in a map int[string], ....
version 3.41
- add tools to get mesh of fespace Vh in 2d and 3d
mesh th2= Vh2d.Th ; mesh3 th3= Vh3d.Th ;
not the type of the mesh is a const pointeur => lot of cosmetic change
version 3.40-2
- try to resolve problem with petsc metis and freefem++ metis
now freefem++ use petsc lib if petsc is find for
parmetis, metis,ptscotch,scotch,fftw3_mpi, mumps,scalapack
- correct hpddm examples
- rename plugin schwarz in hpddm
- correct problem on compilation do scotch
version 3.40-1
- move of all hpddm code un examples examples++-hpddm ( Own by P. Jolivet and F. Hecht)
version 3.40 (14 sep. 2015)
- add new C1 Finite element HCT see examples++-load/bilabHCT.edp
- change permission in install example++ for read and write under windows
- try to build a windows 64 version under msys64 system ...
- Add use ArrowSize= of size
version 3.39
- remove compile flags -fPIC under windows
- pas to version MUMPS-5.0.1
- correct some problem of memory leak, now the decrement method in RefCount call delete
and all decrement are change in Add2StackOfPtr2FreeRC(stack,m); (more safe).
so the previous plugin become incompatible.
version 3.38-2
- add hppdm explain in the doc.
- add function cube to simplify the construction of mesh of cube like square in 2d
see cube.edp example.
- add function ClosePoint in ClosePoint plugin
- add += example of varf
- correct mortar example: replace P0Egde in P0 => better result
- correct map func if dfft plugin
- add splitComm for hpddm
version 3.38-1
- add finite element Edge13d (Finite Element of degrees 1) Thanks to [email protected]
and exemples of wave guide in waveguide.edp
- correct renumbering function in case of rhs
- correct new type of optimisation in integral, (version 3.38 25june 2015)
optimized=2 => do optimisation without check
can be useful in case of random problem.
- correct convect operator in 3d (remove randow part,
to remove problem of optimisation in varf.
28/7/2015 (ALH)
- fixed behavior of "verbosity" to use value from command line ("-v" option)
- synchronized lg.ypp with lg.tab.cpp
version 3.38
- add parameter in renumbering function.
version 3.37-1 (22 mai 2015)
- clean examples-mpi, remove of all usage to RemoveDOF, bb2d bb3d, findDiff (see version 3.36)
- update schwarz.cpp for new hpddm version
- add read of real[int,int], int[int,int], complex[int,int] ..
version 30 april 2015 (3319:8657e0526391)
- correct freefem++ launch in windows (remove wait when launch vai freefem++)
add wait in case of launch through launchff++.exe
- in DG linear form with jump or mean in test fonction was wrong
like varf a(u,v)=intalledges(Th)(jump(v)*u) ;
- correct problem of try/catch in freefem++ func
version 3.36-2
- correct problem some null ref with compile Apple LLVM version 6.1.0 (clang-602.0.49) (based on LLVM 3.6.0svn) xcode 6.3
some warning exist but all test pass.
- correct problem of cleaning tmp memory for small array
- correct configure and makefile for windows compilation process
version 3.36
- add tools to remove (now build in)
utility.dylib removeDOF.dylib symmetrizeCSR.dylib
bb2d bb3d -> boundingbox
fast -> ltab=lround(dtab)
Unique to get all unique value of a array (findDiff)
- after version 3.33 the some compilation flags a lose, correct configure.ac
Please do not use version from [3.33 .. 3.35] included.
- add disable-gmm configure flags
version 3.35 dist (12/3/15) (rev 3246:664a6473d705) (warning slow version)
- optional lib seach [toto|tyty] in WHERE-LIBRARY seach lib
- change the metis to scotch-metis interface in hips and parms ..
- change the current version of scotch metis mumps parmetis
- remark: .dll , .so, .dylib are incompatible with previous version.
- ff-run.in under windows (use src/bin-win32/FreeFem++-mpi
and not src/mpi/FreeFem++-mpi for dll problem).
- correction init plugin missing static in int DoLoadInit() ;
version 3.34-2 (20/2/2015) (warning slow version)
- do correction do compilation for mpi version under window
big change in ffapi.cpp rules ( use now call back) end init in mymain.
schwarz plugin run under windows without petsc lob.
- add some function to then gsl interface (airy, biry, ...)
version 3.34-1
- correct lot of mistake for simple compilation of
hpddm interface ..
- add no mandatory lib for petsc
write the WHERE-LIBRARY search lib in awk (more simple )
version 3.34 (warning slow version)
- configure : 6/02/ 2015
correct mpi for sgi uv computer
find gsl lib
add find petsc (in progress)
add missing file in distribution
- correct for compilation with g++-4.9.1 -std=c++11 ( without downlaod)
BUG in configure => slow version after until version 3.36 , 25 mars 2015 ... (dur Dur)
no optimisation flags in compilation
- add hd5 interface (13/01/2015)
Thank to Mathieu Cloirec CINES - http://www.cines.fr
voir example iohd5-beam-2d.edp iohd5-beam-3d.edp
- add find of libgsl in configure script
- correct pb of memory leak in case
matrix A = ...; in loop ( this occur in lot of case, tanks to P. Jolivet )
correct small memory leak in use of routine due to debugstack.
- correct bug in periodic condition in case common dof with periodic.
- correct big bug in memory management of sparse matrix
4/3/2015 (ALH)
- Edited for compatibility with Emscripten
- Added calls for Javascript graphical output along with postscript output
15/10/2014 (ALH)
- Enabled monothreaded downloads in download/getall for platforms that do not support threads (eg Mingw)
- Corrected --disable-pdf configuration option (useful for FFCS build)
- Removed obsolete cvs history directives
- [[file:download/getall]] now checks downloaded packages against known MD5 values
version 3.32
- correct of problem of plugin and mpi,
build all dynamics lib with and without mpi,
the mpi version is install dir lib/mpi
- correct of plugin MUMPS.cpp for complex value.
- add vectorial operator a/v and v/a where a is scalar and v vector like real[int], ...
version 3.31-3 (rev 3053:4e164226411d 12 Aug 2014)
- correct the problem of size of arrow in 2d plot
version 3.31-2 (rev 3052, 11 July 2014)
- correct stop test function in LinearGC (for [email protected])
build tar.gz distribution (rev 3050)
build version MacOs 3.31-1
- correct bug put in DG formulation (rev 3044)
jump, mean , was wrong from Sun Jun 29 22:39:20 2014 +0200 rev 3028
version 3.31-1 (rev 3042, 10 july 2014)
- function to put your stop test in LinearGC and NLGC
the prototype is
func bool stop(int iter,real[int] u,real[int] g) ;
{ return g.linfty < 1e-5 || iter > 15;}
LinearCG(DJ,u,eps=1.e-15,nbiter=20,precon=matId,verbosity=50,stop=stop);
- add functionnal interface to arpack (Eigen Value)
func real[int] FA(real[int] & u) { real[int] Au=A^-1*u;return Au;}
func real[int] FB(real[int] & u) { real[int] Au=B*u;return Au;}
int k=EigenValue(n,FA,FB,sym=true,sigma=sigma,value=ev,vector=eV,tol=1e-10,maxit=0,ncv=0);
see examples++-eigen/LapEigenValueFunc.edp for a true example
version 3.31 (rev 3037, 1 july 2014)
- re-add tan function for complex number
- correct a big mistake in LinearGMRES , the result are completely wrong,
correct also the algo.edp
- add sqr functon of O. Pironneau
- correct update of mercurial depot (rev 3034, 1 july 2014)
- correct misstake in examples++-3d/MeshSurface.idp
about computation metric in function build surface Sphere mesh
and add a function to build Ellipsoide surface mesh
func mesh3 Ellipsoide (real RX,real RY, real RZ,real h,int L,int orientation)
- correct a bug the DG with periodic boundary condition with only one
layer of element.
- add pluging "bfstream" to write and read in binary file (long, double, complex<double> and array)
see bfstream.edp for an example.
version 3.30-1 may/2014 ( hg rev: 3017)
- add levelset integral on 3d case ( on levelset and under level set)
- correct problem with Ipopt / lapack configure ...
- add BEC plugin of Bose-Enstein Optimisation
- standardisation movemesh3 -> movemesh ( same parameter of 2d version )
- correct jump in basic integrale to be compatible with varf definition
jump is not zero on boundary, it is - boundary value ..
fespace Ph(Th,P0);u= 1;int1d(Th)( jump(u) ) = - length border .
30/6/14 (ALH)
- [[file:download/getall]]
- now checks downloaded packages against known MD5 values
- new -h inline help option
version 3.30. ( for windows hg rev : 3013)
- add binary ios:mode constant, to open file in binary mode under window
to solve pb of seekg under windows
- add multi border april 23 2014 , (hg rev : 3004)
syntaxe example:
// i is the index of the the multi index
// the number of sub border is given by teh size if the array to set the number of seg on the border ..
// so for multi border the number of seg of the border , must be in int[int] array
real[int] RC=[ 0.1, 0.05, 0.05, 0.1],
XC= [0.2,0.8,0.2,0.8],
YC= [0.2,0.8,0.8,0.2];
int[int] NC=[-10,-11,-12,13];
border bb(t=0,1;i)
{
// cout << " i = " << i << endl;
int ii = (i+1)%4; real t1 = 1-t;
x = xx[i]*t1 + xx[ii]*t;
y = yy[i]*t1 + yy[ii]*t;
label = 0; ;
}
border cc(t=0,2*pi;i)
{
x = RC[i]*cos(t)+XC[i];
y = RC[i]*sin(t)+YC[i];
label = i+1;
}
int[int] nn=[4,4,5,7];
plot(bb(nn),cc(NC),wait=1);
mesh th= buildmesh(bb(nn)+cc(NC)) ;
plot(th,wait=1);
- add ltime() (rev 2982) function returns the value of time in seconds
since 0 hours, 0 minutes, 0 seconds, January 1, 1970, (int)
- add new macro tool like in C (rev 2980)
FILE,LINE,Stringification() to get line number and edp filename,
add quote to a parameter after macro generation …
Like in
cout << "in " << FILE << " line " << LINE << " -- '"
<< Stringification( "zzz" aa () {} + /* */ bb cc) << "'" << endl;
add new int2d on levelset in 3d (int test)
add basic func mesh3 Cube(int nx,int ny,int nz) in cube.idp file.
version 3.29 (hg rev 2973)
- add int storagetotal(); and int storageused(); function
to get static of malloc ( memory storage usage)
- correct problem of region evalution in jump and mean function
version 3.28 ( merge with freefem++-cs tool)
- add download/getall perl script to download all related soft
- add int1d on isoline for matrix ...
version 3.27
- correct bug in display of P1b finite element in 3d
error in SplitMesh<R3> function (Thank to O. Pironneau)
- add AddLayers(Th,suppi[],sizeoverlaps,unssd[]);
- add tool to trunc to get element numbering for Thn to Tho
int[int] In2o(1),Io2n(1);
Tho = trunc(Tho, x<0, new2old=In2o, old2new=Io2n);
- add restrict function for get dof numbering old to new
fespace Vnh(Thn,Pk), Voh(Tho,Pk);
int[int] n2o=restrict(Vnh,Voh,n2ok);
- correct mistake in gsl interface random number (9/1/14)
missing all random distribution ..
- add interface with gsl random number generation in test ..
- correct pb of compilation under window (mingw32/msys)
6/2/14 (ALH)
- Creating documentation hyperlinks (in emacs orgtime and in Doxygen format)
- New script build/orgindex and file index.org list all hyperlinks in the FF source in Emacs org-mode format
- Set all hyperlinks to relative paths to make sure that they work from any location
- Separate download script for all third-party software [[file:download/getall]] (request from FH)
- created common makefile goals in [[file:download/common.mak]] for all downloaded packages (request from FH)
- changed [[file:download/arpack/Makefile.am]] to make use of [[file:download/common.mak]]
- done for [[file:download/blas/Makefile.am]] as well
- added option --enable-mkl-mlt in [[file:configure.ac::enable_mkl_mlt]] to allow FF to be linked with the
multithreaded MKL when an external library requires it (request from Atsushi Suzuki)
- cleaned up old comments from the documentation to enable automatic convertion into other formats (eg wiki format)
- many small corrections for Windows compilation
- added "NbColors" in plot() parameters
version 3.26-3 09/12/2013
correct problem with openblas need pthread lib
01/12/13
install dowload lib
update lapack interface in examples.
add WHERE_LIBRARY-file in install
change the order of seach lib with WHERE_LIBRARY-download
WHERE_LIBRARY-config WHERE_LIBRARY
remove all relative in WHERE_LIBRARY-config
correct of superlu compile problem
13/11/203 Pass to version 3.26 ...
- correct pb of ambiguity with new complex lib
31/10/2013 compilation on MacOS 10.9 with compile of xcode 5.0.1
- correct the pipe.cpp to compile with clang-500.2.79
- remove of warning message genrate by clang++ -std=c++11 (version clang-500.2.79)
- correct also compilation problem
- correct error in case of block matrix with block of 0 row or 0 column.
21/10/2013 (ALH)
- umfpack configuration cleanup (request from Fred)
- blas configuration cleanup (request from Cico)
- configuration option --enable-generic-blas has no effect anymore, removed
- Atlas compilation directives have no effect anymore, removed
- Moved OpenBLAS compilation from FFCS source tree to FF source tree
- Added configure option --disable-system-fftw to skip any default FFTW library installed on the system
- Changed location of lapack WHERE-library setup in [[file:configure.ac::WHERE_lapack]] to insure that it is always
defined (request from Fred)
- changed FFCS Fortran MPI configuration for MPICH2
- split DOC makefiles into Makefile.am and figs.mak to suppress automake warnings about portability
- moved build/download script from FFCS to FF before reorganizing FF downloads (request from FH)
10/9/2013 (ALH)
- Corrected pastix compilation for FFCS on MacOS 10.6
- Removed dependency from lg.ypp to lg.tab.?pp in [[file:src/lglib/Makefile.am::lg.tab.?pp]] to avoid automatically
generated conflicts
05/09/2013 put version .tar.gz on the web
- correct makefile wget in pastis and superludist
- correct compile of load example : paradiso, gsl automaticaly
- add AutoGeneratedFile.tar.gz a file contening all
file build by autoreconf in case of non automake tools (1.13)
- add missing file in dist or in mercurial distribution
- correct download of scotch un curl
- correct problem in a*[b,c, ... ]' in case of complex value
5/9/2013 (ALH)
- force Umfpack build to run sequentially (parallel make crashes with "pipe from processes is a directory"?)
- pARMS download URL changed to http://www-users.cs.umn.edu/~saad/software/pARMS/pARMS_2.2.php (thanks Fred)
- corrected download/mumps parallel compilation (WHERE was not built properly)
- Removed all Mercurial-tracked files from .hgignore
- added options --with-[package]-include= and --with-[package]-ldflags= to avoid downloading existing packages
- Parallelized MUMPS compilation in download/mumps (and mumps-seq)
- Applied a patch from Fred for compiling SuiteSparse on Slackware64-14.0
- Added configuration option --enable-hypre (disabled by default, contrary to other tools)
- Deactivate FFTW download when a local version is found (request from Helmut Jarausch)
10/7/2013 (version 3.25)
- remove of Makefile.in configure for the hg distrubion
use : autoreconf -i # too build Makefile.in
before to configure
need automake version 1.13 ...
- merge FFCS (ALH) version and ff++ version (FH) of freefem++ programs
- remove all automake file form the hg data base
- add new parameter to ffglut for demo of freefem++
ffglut [-nv|-v|-vv|-vvv] [-wait 0.5] [-g 512x300+10+10] [-t title] [file]
all number can be change the wait is in second 0 is default value nowait bewteew plot
-g 512x300+10+10 is the geometry of the graphic window
-t the title of the windows
26/6/2013 (ALH)
- created a build/ subdirectory for build tools
- enabled parallel make ("make -j")
- created a separate file (acmpi.m4) for complex MPI configuration options
- added configuration option --enable-ffcs to make the FF source compatible with FFCS
- backported all current FreeFem++ patches for FreeFem++-cs into the FreeFem++
- started main doxygen page mainpage.dox
9/06/2013
- correct extract function of mesh Lo Sala <[email protected]>
- correct int2d on levelset see example intlevelset.edp
- correct automake TESTING part (in progress)
- correct typo on the doc with .*= ./= operator
- correct bug in RT0 3d , code in the construction of the DOF.
the bug is all dof on border face of same element have same dof number.
thank to Laurent Bernard <[email protected]>
version 3.23
- do cleanning in version remove x11, glx, std : freefem++
clean compile option for clang and clang++ compiler
- add flags to remove internal boundary in 2d,3d in function change
rmInternalEdges=1
- correct glumesh in case of no mesh in 2d
version 3.22
- add multi windows graphics ; WindowIndex=0 in plot function
add new event in graphic windows
* to set/unset default graphics stat to previous plot
- add getenv, setenv , unsetenv function in shell plugins for gestion of environnemnt variable for
openmp.
- correct pb un trunc for 3d mesh with too flat element (sliver) , and cleaning code .
- correct bug in gestion off outsite flag is 3d in case of brute force (searchMethod>0)
version 3.21-1
- correct bug a=b' ; of full matrix
- bug in assert on 3d mesh of huge dimension with periodic condition.
- correct build Delaunay 3d mesh of set of point.
see examples++-load/convexehull3d.edp
version 3.21 frev 2013
- n, resize in array of finite element function ...
see edp array.edp file
- correct pb of compilation of Ipopt with clang++ -std=c++11
- correct une NSCahouetChabart.epd examples , and correct un doc.
this is a complete rewritting. FH. 12/02/2013.
- correct in change function the "flabel=" parameter in 2d and 3d
version 3.20-3
- add master= in solver of MUMPS interface to set the master mpirank
<0 => distributed matrix.
- correct problem in label generation in freeyams (18/01/2013)
all label border was set on 1 before.
version 3.20-2
- add MUMPS parallel version (in test)
- add paradiso seq solver ..
version 3.20
- correct isoline plug in case of sadlle point
- change the compilation tools under windows gcc 4/7 + freeglut / ...)
- change the compile tools on mac pass to clang++
- compile ok of c++11 compiler
- add formal tools on array [] or matrix [[],[],] for elastic problem.
let A a formal array
if A is matrix 2x2 or 3x3 : trace(A) , det(A) , Cofactor(A) ;
A:A = sum_ij A_ij * A_ij
2*A, A*2 // muliplication by a scalar
-add integration on levelset line (in test)
version 3.19-2
- correct pb of C in/output in pluging (in test)
- correct bugs mshmet pluging in case of double Eigen value
- correct typo problem (string size) when a change the default
window size in postscript in version 3.19-1 to have more
precise postscript plot
version 3.19-1
- add tool to create Quadrature formular 1d,2d,3d with
plugind: load "qf11to25"
real[int,int] tab1(np,2),tab2(np,3),tab3(np,4);
QF1 qfe1(norder,tab1);// 1d
QF2 qfe1(norder,tab2);// 2D
QF3 qfe1(norder,tab3);// 3D
where tab(0,i) = weigth , tab(j,i) = j coord, norder order of the
quadrature
see examples++-load/LaplaceP4.edp for example
- correct download auto compile of
mmg3d, mshmet, scotch, ...
version 3.19 (20 april 2012)
- scotch partitionner interface see scotch.edp in examples++-load
- add isNaN(x), isInf(x), IsNormal(x) function to check floating point
of real value x, see ISO C99.
- add function : NaN() and NaN("") to build NaN real number (double in C) .
- correct error in macro with operator ./= and .*=
- add Ipopt doc (thanks to Sylvain Auliac)
- add Ipopt interface (thanks to Sylvain Auliac)
- correct 3d trunc bug in case of internal boundary
thank to Yoshihiro Tomita <[email protected]>.
- add new type of array , array of array
see taboftab:edp in examples++-tuturial
real[int] a;
real[int,int][int] m(10);
real[int][int] v(10);
not well tested.
version 3.18-2
- add plugins with sequentiel mumps without mpi
- add conversion of re and im part of complex sparse matrix
A = C.im;
A = C.re;
version 3.18-1
- correct Typo error in example
- add generation of error in case of periodic boundary condition
non scalar problem.
- add tools for adaptation of P2 and P3 finite elements with metrics
see APk-AdaptEpsDeltaPk.edp APk-FreeFemQA.edp APk-MetricPk.edp
APk-ExplicitPkTest.edp APk-LaplaceDirac.edp
- New example in in Chapter 3 Navier Stokes Newton NSNewton.edp
- add cod to build matrix and vector form varf in 3d in case of
different meshes (in test )
- correct NSprojection.edp chap3 example (PB in out flow BC.)
- correct compile of mmg3d v4 plugins (small change un distrib archive)
v 3.18 (11/01/2012)
- rewrite the isoline-P1 plugins (new name isoline.{dll,so,dylib} )
see example Leman-mesh.edp and isoline.edp in
examples++-load directory
- correct bug in cas of resize of array with 2 index (full matrix)
- correct assert in MPI in gather and scatter
- correct bug in case of return in loop for or while.
- correct a=int2d(Th,l)(1) in case simple expression when l is a array.
- build a pkg under macos for distribution .
v 3.17 (17/11/2011)
- correct PB of pugins: MUMPS, parmetis, metis, mmg34, mshmet
- the new load interface (for more safe IO)
- build mpi for windows with msmpi ( begin)
memory error after end ..
- remove add by default $LIB_MPI in ff-c++
- change the way add thing in freefem++ in case of dynamics load
to be compatible plugin with freefem++cs ( solve pb with
cin,cout, cerr in dll may be..).
LOADINIT(Init); // where init is a class
or do
addingInitFunct FFinit(100,ffinit,"MUMPS_FreeFem"); // where ffinit the init function ...
v 3.16-1
- cmaes interface in scalar and MPI case (thank to S. Auliac)
see doc and examples : cmaes-mpi-VarIneq.edp , cmaes--VarIneq.edp
- add NLopt interface (thank to S. Auliac)
v 3.16
technical stuff:
- correct auto load of mumps
- add tool to do automatic load or static load (for testing)
see MUMPS_*cpp example
v 3.15
- correct all examples++ load and 3d rebuilde the all.edp
and test this examples.
- correct version old bug when full matrix (arry)
A=A;' is now correct (set and initialization )
and add A+=A'; A-=A';
- reput metis 4.0 form netlib ... (pb of compatibilite with other // soft hips)
v 3.14-1
- change interface with metis 5.0.1
- Complete writing of mmg3d interface with version 4 of mmg3d
le plugin is "mmg3d-v4.0", the parametre are the same have
command line mmg3d , execpt the file name are not given.
- remove old bug 3D in interpolation P1 operator
correct bug 29/08/2011 (thanks to [email protected])
remove wrong bulid of KHat (memory out of bound)
v 3.14
- correct in configure remove the default -O2 -g autoconf value
and add -g in case of --enable-debug
- a very old error in the on() functionnal
the bug is with a vectorial finite element like RT0, ..
the name of vectorial must be in lexicographic order
so u1,u2 is ok but v,u is wrong
is correct in version 3.14 the 24/08/2011
- correct trap in check convect-apt.edp example
- add the existance of patch in configure
v 3.13-3 ( 30 june 2011 Seville)
- correct the Hips interface (not to bad , in test)
load "hips_FreeFem"
int[int] iparm(1);real[int] dparm(1);
HipsDefaults(iparm,dparm);) // set def option
( limit to 100 Hips active linear system ).
- add interface with MUMPS_4.10.0 version (with automatic download )
- add P1dc3d finite element in pluging "Element_P1dc1"
- correct bug in gibbs renumbering in 1 case veru sample mesh)
- correct bug in interpolation operator with periodic BC
mesh Th1=square(2,1), Th2=square(2,2);
fespace Vh2(Th2, P1),Ph2(Th&,P1,periodic=[[1,x],[3,x]]);// une couche
matrix Jh=interpolate(Vh2,Ph2,op=0,inside=0);
Ph2 w=1.;Vh2 wi; wi[]= Jh*w; // wi must be 1, now this ok.
v 3.13-1
- correct compilation problem on fedora 13 WITH MPI
v 3.13 (25 may 2011)
- update the finite element list in documentation
- add (load "Element-Mixte") NEW FINITE ELEMENT 2D
TDNSS1 sym matrix 2x2 conforme in $\{H(div div) / div(div s)) \in H^{-1} \}$
RT1 and BDM1 conforme in H(div) Raviert Thomas of Degre 1 and Bezzi, Douglas, Marini
RT1 and BDM1ortho conforme in H(curl) 2 Nedelec Finite Element
v 3.12-3
- new finite element in 2d for elasticity in test. TD-NSS0 (see these of Astrid Sabine Sinwel)
A New Family of Mixed Finite Elements for Elasticity
- add matrice matrix multiply in lapack interface
- add tool to change the region and label number
change(Th,fregion= integer function to set the region number )
- nuTriangle now given the tet number in 3D.
- add mpireduce of matrix to build parallel matrix
v 3.12-1 april/2011
see remove examples++-mpi/chaleur3D-superludist.edp
- correct the precond of gmres algo in A^-1 : (29 mars 2011) to build a correct Navier-Stokes
- add cast TypeOfSolver on int to parametrize the choice of linear solver.
- correct intersection of given metrix in adaptmesh (bug introduc in 07/10 verion 3.9
in case of convect-apt example
v 3.12-1 10 10 fevr 2011
- add VTK write for paraview, examples++-load/VTK_writer.cpp examples++-load/VTK_writer_3d.cpp (see source for moe details)
- put comment in the documentation about negative tgv
- correct pb blacs mkl (unde with mumps)
- correct mpiReduce for complex and real data.
- add mpiAllReduce(umax,dmaxg,comm,mpiMAX); where for real umax,dmaxg;
- add tag verson 3.12-win32
- add inferface for mpi in win32 architecture (form 32 version)
v 3.12 17 jan 2011
- correct probleme of comm world on SuperLuDist (complex version)
- correct medit Makefile.am in case of no compilation of medit ..
- correct a lot of MPI parallel solver example
a bug pastix interface ???
MUMPS, superludist, hypre , hips works
- correct link problem in hips and hypre under linux
- Add thresholdings.cpp thresholdings.edp in examples++-load to remove to small coef
in a matrix .
- Add lots of DDM Schwarz GMRES precondionned with a coarse grid solver por overlapping :
DDM-Schwarz-Lame-2d.edp DDM-Schwarz-Lap-2dd.edp DDM-Schwarz-Stokes-2d.edp
DDM-funcs-v2.idp DDM-Schwarz-Lame-3d.edp DDM-Schwarz-Lap-3d.edp
DDM-Schwarz-macro.idp
I will add the explanation in the doc in the future;
- Add a true exemple to buld mesh form a image (lg.pgm) , Leman-mesh.edp
- Add New syntaxe in macro generation
NewMacro a(i)
EndMacro
with // comment and macro definition, first trik
- Add interface with special function of gls http://www.gnu.org/software/gsl/
the full list is in the example examples++-loal/gsl.edp
v 3.11-1 25 dec 2010
- Add coarse preconditioner in MPIGMRES[23]d.edp (to be optimal)
now the number of iteration is close to 10.
- Show ff++ line number in case of assertion in RNM class.
- add Coarse Preconditionner for MPIGMRES[23]d.edp (Good)
- solve bug For MPIGMRES2D.edp due to interpolation bug in rare case
by Add brute force for seach of point in 2d like in 3d in the search in find ouside
set global variable : searchMethod=1; // more safe search algo (can be expensive in case of lot of outside point)
- add hack of ILU precond, if tgv is < 0 then
we remove all the line and put 1 on the diag term ..
v 3.11 9 dec 2010
- update the documentation of 3d adaption process
and / correct freeyams,mshmet, mmg3d interface and associaded example
- configure is compatible with MKL lib (on gnome)
- add quoting argument in macro argument with { } for mpi plot ..
v 3.10-2
- try to compile with MKL libs.
- correct MPIGMRES[23]d.edp example
add doc on this example.
v 3.10-2
- add operator to inverse permutation to set or initial int[int] array
I=J^-1;
if J is a permuation of 0:n-1 then we have : I[j[i]]=i
- correct comment problem of periodic boundary condition in 3D
see examples++-3d/periodic-3d.edp
- correct configure to scotch compilation (phtread)
v 3.10-1
- at convection function form formal array to int,real,complex array
(resp. toZarray, toRarray, toCarray)
- correct ffrandom.cpp to read /dev/random to get a true random seed
srandomdev
- add ff-mpirun script to simplify the launch of FreeFem++-mpi version
- correct MPI for MPI icc on gnome
- correct pb of computation of area, lenbord in mesh type in some case .
- correct WHERE_LIBRARY-config for blas find in configure
- correct atof
v 3.10
- add true Domain Decomposition example in 2d and 3d.
see MPIGMRES2D.edp MPIGMRES3D.edp in MPI examples
- add mpi Isend/Irevd for complex struct like matrix, meshes
for send/recv Huge objet
- add named parameter verbosity= in add GMRES and GC function
for show algo evolution .
- add MPIGC dynamic LIB for // GC and GMRES , add fully // scharwz RAS
in MPIGMRES2d.edp
- correct problem of a=A^-1*b when a, or b is not consecutive array
ie. expression like M(2,:)
- correct problem of Makefile in download clean,install,WHERE target
v 3.10 ( 5 oct 2010)
- add install of missing MPI dynamics lib.
add argument "-cd" to FreeFem++ command to change current directory to edp script directory
- clean configure.ac
- do correction for g++-4.6 compiler (lambda expressions , and trap ..)
- add MPICG for Parallel Conjugate Gradient for full split matrix
see MPICG
- correct Makefiles in download ( add WHERE interface)
- remove wait option in medit because the code is wrong , always waiting now
- correct bug renumbering of matrix : B= A(I,J),
the last term N,M was force to zero =+ -> += (line 1885 of file lgmat.cpp)
where N= I.max, M=J.max
v 3.9-3
- correct compilation of gmm, mumps on linux
- add parameter -ne in FreeFem++ commands to remove edp script print
- correct mistake in --enable-m64 or --enable-m32 (suppress the configure warning)
- correct of ff-c++ script if whith space in path (for windows)
- add compilation of gmm library
v 3.9-2
- correct compilation of mshmet dynamic library
- correct pb of compile of superludist ( add CNOFLAGS no optimize CFLAGS)
- correction of lib on win32 for freeyams and mmg3d
- correction download/Makefile.am to be sure than bin is a directory.
- add fftw of win32
v 3.9-1
- correction mmg3d interface (J Morice)
- correct of mmg3d and freeyams under wind32 (ld problem)
v 3.9-1 August 2009 (For FreeFem++ days )
- correct configure (find lapack lib change $ll_lapack_libs in $ll_lapack_lib)
- correct mistake in mpi
add: gatherv, allgatherv, scatterv, alltoallv
and complex data type in: allgather, gather, scatter, alltoall, ..
correction bug : gather, scatter
correct essai.edp mpi example
- correct problem of compatibility of dynamic lib with and without MPI
change ff-cc++ to add all MPI libs if MPI version exist.
- correct default region number in square build mesh function now 0 , (2 between v3.8 -- 3.9)
- change in all example reffacexx= in labelxx= ..., etc to be correct with v 3.8 change.
- correct a mistake of type region= parameter in tetgen all functions.
- correct the mpi configure search tool
miss when no full path are given ( only change in configure.ac)
v 3.9 July 2009.
- add lots of automatic compilation of download soffware,
tetgen superlu fftw metis yams mshmet
blacs parmetis scalapack scotch superludist MUMPS pastix hypre hips
For the link with mmg3d software put the tar.gz archive in .../dowload/pgk directory.
add interface with freeyams, mmg3d, mshmet (3d mesh adatation) software,
add automatique compilation of // solver with
flags -auto in ff-c++ commands
v 3.8-2
- add cast operator from SubString to String
- correct the SubString tools
v 3.8
- add in change mesh2 tool to make a renumbering of vertex for periodic 3d mesh.
Th= change(Th,renumv=old2new);
- correct ' operator do alway in complex case the conj and trans (Hermitian stuff)
in formal array with [ ]'
in linear comb of complex matrix
now A + (-1)*A' is zero is A is hermitian .
in varf term (a*dx(u)*dy(v) )' <=> a' dx(u)*dy(v)
now you can [a*dx(v),dy(v)] ' *[a*dx(u),dy(u)] in complex case..
- plot of complex field and 3d vector ..
- upgrade the documentation in 3d mesh example
- change configure.ac for mpifc, mpif77, mpicc, ...
- correct include seach path
- uniformize named parameter in change, movemesh, in load "msh" , glumesh, ...
add synonyme:
refface= -> label=
reftet= -> region=
in buildlayer:
labebup=
labebdown=
labelmib=
- writing schwarz-nm-3d.edp examples
- add array of 3d mesh
- add word volume to get the volume of the current tet element
- correct metis.cpp in 3d case
- correct in fflapack.cpp (example++-load, computation of eigenvalue of complex matrix)
-llapack search in configure.ac
and unify fflapack.cpp lapack.cpp are the same file.
- add seekp, tellp method on ostream type
seekg, teeg method on istream type
see examples++-tutorial/readmesh.edp examples
v 3.7-1 13 jan 2010..
- EqPoisson.edp (Solve Poison equation in 3d fish)
- add possibility to put array (int[int] ) to set a set of label/region in
in integral (int1d,int2d,int3d) or "on" key word
example : int[int] l=[1,2,3] , ... on(l,u=1);
- add tool to compile under mingw/msys under windows 32 (freefem++-nw works,
- add bug (random trap) in bamg in sub domain computation in case of
internal edge in same region (fist time in 15 years).
- add bug in varf interpretation in complex case (miss some convertion)
- add build interpolation matrix in 3d (see examples++-3d/mat_interpole.edp)
- correct CFLAGS (add -fFIP) in superlu , umfpack under 64 architecture.
v 3.7 8 december 2009
add 3d beam examples examples++-3d/beam.edp
correct install problem due to metis directory
add dynamic load interface with newuoa fortran optimizer without derivative
see ffnewuoa.edp example ins examples++-load
correct problem of free of mesh in case of gluing meshes
Th= Tha+Thb; // now just the ref counter on Th is decrease
// and before Th is delete
add .im, .real method on complex [int,int] matrix array
add missing matrix<complex>[int] type in grammar.
correct mistake in generation of file WHERE_LIBRARY-config
add initialization of scalar variable at zero
add warning in case of use of variable with hide freefem++ name
add missing file in download version
v 3.6.1
correct font problem in documentation
lot of mistake in examples++-load
add tool to compile all examples
add PICHON stuff for brute force search point
add umpack long interface to by pass the 32bit limits
see UMFPACK64.cpp and LapUmfpack64.edp in examples++-load.
add some interface with lapack (inverse of full matrix,
eigenvalue of full matrix)
correct ::Draw undef ref. ubuntu 10. (gcc 4.4)
correct typo error image tools
see UMPACK64 (code and example in examples++-load)
v 3.6 (4 nov 2009)
change the version of tetgen to 1.4.3 download
correct in ffglut max value of arrow
add examples of function with stack parameter in
examples++-load/funcTemplate.cpp load facility
and
examples++-load/funcTemplate.edp
add tools to read .pcm files for optic-flow computation (pcm2rnm.cpp)
see http://www.cs.otago.ac.nz/research/vision/Research/OpticalFlow/pcm.html
correct bug in GMRES without preconditionner assert error.
correct a rare and ramdomly bug in the interpolation 2d process when the
the element 0 touch the boundary
correct on versy old bug in bamg in case of internal boundary edge with
two extremity in true boundary (symptoms FillHoleInMesh trap).
correct operator subsript a:b:c operator (in some case the last value was miss)
correct the ?: operator in case of constant array
- add -enable-m32 flag to build 32bits freefem++ on 64 architecture computer
- correct eigenvalue.cpp in case of complex finite element
v 3.5-2 (28 sept. 2009)
- correct install problem under windows (metis.dll do not compile)
- correct MacOS install missing /usr/local/lib/ff++/3.xx-yy link
v 3.5-1 (24 sept. 2009) not to bad ...
- add interface to metis
- correct display of 3d mesh (with light today)
- add trap in case of mesh with negative volume
- correct bug in case 3d finite element with constant number of df / node
example fespace Vh(Th2,[P1,P1]); or fespace Wh(Th3,[P0,P0,P0 ])
- correct problem with tetgen under windows and ubuntu
(remove optimization during compilation of the library)
- correct setting array of vectorial finite element Wh
Wh[int] [u,v](10); // array of 10 value
u[1] or v[1] 2 componante function for indice 1.
[u[2],v[2]] = [1,2]; // set the function for indices 2.
- correct in matrix size of matrix B=A(I,J) to I.n x I.m
where I, J are array of int
- add resize of complex sparse matrix
- Add inferface with metis in examples++-load (metis.cpp, metis.edp)
v 3.5 (27 aug. 2009)
Complet writing to MPI interface, tested with Open MPI 1.2.3 and mpich 2.1
see chapter 10, page 229, of the freefem++doc.pdf and
see example++-mpi/essai.edp.
correct
operator a(:)= 1:4; // before change a copy of a and so do nothing
string[string] s;
s[i]; // before trap (two delete)
v 3.4-2
add boundary mesh inquire in 2d and 3d
Th.nbe ; // return the number of boundary element
Th.be(k); // return the boundary element k $\in \{0,...,Th.nbe-1\}$
Th.be(k)[l]; // return the vertices l $\in \{0,1\}$ of boundary elmt k
Th.be(k).Element ; // return the triangle containing the boundary elmt k
Th.be(k).whoinElement ; // return the edge/face number of element containing the boundary elmt k
Th[k].adj(e) ; // return adjacent element to k by edge/face e, and change
// the value of e to the corresponding edge in the adjacent element
Th[k] == Th[k].adj(e) // non adjacent element
Th[k] != Th[k].adj(e) // true adjacent element
do small change in msh3.cpp file
v 3.4-1
rewrite of the mpi interface in vue of use group and communicator
now we only use mpi.h and not mpi++.h
add plot of array of mesh and array of finite element functions
add new command in ffglut graphic interface
n,N,i,I to change the number of iso-value and to set the min,max value to the
graphic
add left button motion to zoom , and + alt to translate
V 3.4
add in parallel version mpi send, receive and broadcast of 3d mesh and matrix ( sparse matrix).
build a 3d mpi schwarz example (see example++-mpi/schwarz-3.edp)
make a wrapper (launchff++.exe) to launch freefem++ under windows/OS in the edp file directory.
correct in interpolation matrix ( unset variable in case of same mesh
=> missing some term in the matrix)
V 3.3-3
add missing file in distribution examples++-load/DxWriter.cpp
V 3.3-2 8 June 2009
add 2 new finite element:
Edge03d : the Nedelec 3d Edge finite elements.
RT03d : the Raviart-Thomas 3d finite elements
correct the size of anytype for 64 bits architecture.
install the new .idp file (include of edp file) in
directory idp
V 3.3-1 1 June 2009
to correct the problem of the destruction of parameter in case:
func real[int] a2(real[int] a) {a[0]=2; return a;}
so now the parameter cast to KN_<double> so the array is not
duplicate like in C++, and the value of the parameter is change.
- simplify the internal (C++) definition of array now of 2 type
KN<> * of variable (left expression)
and KN_<> of right expression
- correct UnRef<A,B> where B ~ A*
v 3.3 29 may 2009
- correct some problem in return in function :
the following function prototype now work without memory fault.
func real[int] a1(int n) {real[int] a(n);return a;}
func real[int] a2(real[int] &a) { return a;}
real[int] G=[7,8,9];
func real[int] a3(real[int] &a) { return G;}
func mesh carre(int n) { mesh th=square(n,n); return th;}
(lot of change, warning: now KN<K> is not a freefem++ internal type).
warning :
- in int3d flag qfV= a 3d quadrature formula
the 3d quadrature formula are
qfV1 , qfV1lump, qfV2 , qfV5 (where the number is the order of the quadrature)
- correct sign of the normal in 3d (now the normal is exterior)
-------------------------------------------------------------
- correct in adapted mesh the flag splitin2= (no tested)
- correct the convect operator in 3d
- add EigenValue tool for 3d case May 15, 2009.
v 3.2
- add tools to set array like in matlab/scilab
real[int] a(1:2:10), def a to 1,3,.. 10.
int[int] I(0:n-1) , set array of size n to 0, .., n-1
real[int] b(0.1:0.5:9.99) , set 0.1, 0.6, 1.1,..., 9.6
add : b.n == (9.99-0.1)/0.5 +1
- sort(a,p); // sort a is and array and p in parallel
where is a integer array same size of a.
- add interface with MUMPS parallel solver (in progress)
http://graal.ens-lyon.fr/MUMPS/
- rename movemesh2D3Dsurf in movemesh23 and change named
parameter normal= in orientation=
- add periodic boundary in condition like in 2d
example for a cube with label face numbering
1 : ( x == xmin) 2 : ( x == xmax)
3 : ( y == ymin) 4 : ( y == ymax)
5 : ( z == zmin) 6 : ( z == zmax)
fespace Vh(Th,P2,periodic=[[1,y,z],[2,y,z],[3,x,z],[4,x,z],[5,x,y],[6,x,y]]);
- automatic: Finite element name conversion beetween 2d and 3d
old : fespace(Th3,P13d);
new : fespace(Th3,P1);
- add complete rewrite of the way to set parameter of the sparse solver
for parallel solver ( no progress)
- add --enable-m64 for 64 bit compilation flags
- correct f2c of take the -m64 flags
- build the universal f2c library with 4 architecture ppc i386 ppc64 x86_64