-
Notifications
You must be signed in to change notification settings - Fork 26
/
ChangeLog
14103 lines (9872 loc) · 460 KB
/
ChangeLog
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
2014-09-05 11:05 Rajesh Kommu
* multiple files changed
Removed Python sections of CitcomS
Wrote Python to standard parameter file conversion tool (Py2C)
Converted Python cookbooks parameter files to original style parameter files
Added output of parameters to a pidXXXXXX style file
Updated the manual to correspond to the new changes
2012-01-24 21:54 becker
* lib/Ggrd_handling.c, lib/Instructions.c, lib/global_defs.h: Fixed
minor ggrd assignment bug.
2012-01-24 18:44 becker
* lib/Ggrd_handling.c, lib/Instructions.c: Debugging GGrd
2012-01-24 18:03 becker
* lib/Ggrd_handling.c, lib/Instructions.c, lib/global_defs.h:
Updating to latest version for sync.
2012-01-19 19:24 tan2
* configure.ac: Stop configure if pyre is enabled and python 2.7 is
used. See Issue622.
2012-01-17 04:57 tan2
* INSTALL: updated installation instruction
2012-01-17 02:41 tan2
* configure.ac: When python 2.7 is used, ask users to download
pythia manually. See issue622 and issue606.
2012-01-14 08:00 becker
* lib/Advection_diffusion.c, lib/Composition_related.c,
lib/Drive_solvers.c, lib/Ggrd_handling.c, lib/Output_gzdir.c,
lib/Pan_problem_misc_functions.c, lib/Topo_gravity.c,
lib/Viscosity_structures.c: Made sure that compositional
viscosity without compositional buoyancy can be assigned without
setting ibuoy=1 and the buoyancy number to zero.
2011-12-20 18:04 becker
* lib/Ggrd_handling.c, lib/Instructions.c, lib/Topo_gravity.c,
lib/global_defs.h: Merged with recent check in, minor changes.
2011-11-14 05:06 tan2
* lib/Regional_boundary_conditions.c: Fixed a bug in reflective TBC
where y-index should be used, not z-index. The bug is harmless
though.
I checked the code history. This bug exists since CitcomS 1.0,
but not in CitcomCU nor CitcomT.
2011-10-14 18:27 tan2
* lib/Output_h5.c: fixed a bug that h5close_field() didn't have a
return statement, reported by Nicolas Flament in issue611.
2011-09-18 03:54 tan2
* lib/Initial_temperature.c: fixed a bug that when the perturblayer
is on the internal boundary, the temperature perturbation wasn't
added to the lower processor.
2011-07-31 20:38 becker
* lib/Ggrd_handling.c, lib/Initial_temperature.c,
lib/Instructions.c, lib/Stokes_flow_Incomp.c, lib/global_defs.h:
Updating to newest version, don't think that there should be any
significant changes on my side.
2011-07-12 16:27 tan2
* CitcomS/Components/Output.py, lib/Instructions.c, lib/Output.c,
lib/Output_vtk.c, lib/global_defs.h, module/setProperties.c:
enhancing VTK output, mostly contributed by Rene Gassmoeller at
GFZ.
- A new input parameter 'vtk_format', which can be either 'ascii'
(default) or 'binary'. When 'binary', gzip and base64 encoding
are used. - The node ordering is changed to CitcomS native
ordering. The consequence is the axis are rotated as z in CitcomS
is mapped to X in vtk, and x to Y, y to Z. - More data can be
outputted in vtk format - Using .pvts file for regional model;
.vtm and .visit file for global model (.vtm format provided by
Tobias Hoeink at Rice, .visit format provided by Kat at
visitusers.org/forum)
2011-03-21 00:13 becker
* lib/Topo_gravity.c: Fixed typo for anisotropic stress
computation.
2011-02-22 04:13 becker
* lib/Drive_solvers.c, lib/Instructions.c, lib/global_defs.h: Added
a flag force_iteration which if on (off by default), will enforce
an iterative solution to the velocity fields even if linear
rheologies were chosen. This is useful for debugging/robustness
tests.
2011-02-18 19:45 becker
* lib, lib/Anisotropic_viscosity.c, lib/Output_gzdir.c,
lib/anisotropic_viscosity.h: Modified storage of anisotropic
viscosity files. Those are now
2011-02-14 05:45 becker
* lib/Output_vtk.c: Improved vtk output, still for Dan to test.
2011-02-12 05:30 becker
* lib/Anisotropic_viscosity.c, lib/anisotropic_viscosity.h: Aedded
the actual files because I didn't know how to link from CitcomCU.
Those are didentical
2011-02-12 05:20 becker
* lib/Anisotropic_viscosity.c, lib/anisotropic_viscosity.h:
2011-02-08 04:06 becker
* lib/Citcom_init.c, lib/Tracer_setup.c: Removed the srand() call
that seemed like a good idea (rev. 17175) but screwed up
tracer-based runs with large number of processors (12x3x3x2).
Weird, but oh well. Someone advise why this didn't work?
2011-02-07 19:00 becker
* lib/Ggrd_handling.c: Added one missing line for anisotropic
viscosity init.
2011-02-07 00:32 becker
* lib/Ggrd_handling.c: Reversed, temporarily, the use of
nearneighbor because it makes it easier to compare model
computations with previous versions.
2011-02-06 23:57 becker
* lib/Full_boundary_conditions.c: Fixed typo that was introduced in
17555 and made temperature field assignment apparently screwed
up.
2011-02-05 23:45 becker
* lib/Instructions.c: Backward compatibility of neglect of both
divergence and pressure for convergence checks.
2011-02-04 05:03 becker
* lib/Viscosity_structures.c: Temporarily moved velocity gradient
computation to ensure compilation.
2011-02-03 21:06 becker
* lib/Ggrd_handling.c, lib/Parsing.c, lib/Viscosity_structures.c,
lib/global_defs.h, lib/prototypes.h: Uptdated prototypes file.
Fixed a few compiler complaints.
Made sure that atemp is always set in Viscosity_structures,
making mat_control > 0 feasible, not just 0 or 1
2011-02-03 17:11 becker
* lib/Ggrd_handling.c, lib/Pan_problem_misc_functions.c,
lib/Viscosity_structures.c, lib/ggrd_handling.h,
lib/prototypes.h: Forgot some functions, not sure why the linker
didn't complain previously.
2011-02-02 01:09 becker
* lib/Instructions.c, lib/Topo_gravity.c,
lib/Viscosity_structures.c: Futher synchronization with CitcomCU
anisotropic viscosity.
2011-02-02 00:30 becker
* lib/Anisotropic_viscosity.c, lib/Element_calculations.c,
lib/Topo_gravity.c, lib/anisotropic_viscosity.h,
lib/global_defs.h, lib/viscosity_descriptions.h: Trial
synchronization between CitcomCU and CitcomS anisotropic
viscosity.
2011-01-26 00:14 becker
* lib/Full_parallel_related.c, lib/General_matrix_functions.c,
lib/Ggrd_handling.c: Added parentheses for logic statements as
per compiler suggestion
2011-01-19 21:00 becker
* lib/Advection_diffusion.c: Checking in Shijie Zhong's revised
Advection_diffusion routine.
2010-12-29 20:47 becker
* lib/Ggrd_handling.c: Added support for non-interpolated grd
input.
2010-12-22 21:53 becker
* lib/Anisotropic_viscosity.c, lib/BC_util.c,
lib/Full_boundary_conditions.c, lib/Ggrd_handling.c,
lib/Regional_boundary_conditions.c, lib/Solver_multigrid.c,
lib/prototypes.h: Minor fixes for compile, still hunting for
convergence problem.
2010-12-08 21:06 becker
* lib/BC_util.c, lib/Full_boundary_conditions.c,
lib/Regional_boundary_conditions.c, lib/prototypes.h: Reverted to
separate horizontal_bc routines for regional and full. the former
has a check for
ROW==E->lmesh.NOZ[level] &&
E->parallel.me_loc[3]==E->parallel.nprocz-1
while the full version has
ROW==E->mesh.NOZ[level]) &&
(E->parallel.me_loc[3]==E->parallel.nprocz-1) )
2010-12-08 20:14 becker
* lib/BC_util.c: Reverted the rowl adjustment since it didn't fix
the problem reported by Dan.
2010-12-08 18:43 becker
* lib/global_defs.h: Added definition of TRUE (1) and FALSE (0) in
case undefined (was defined for GMT/ggrd compile before).
2010-12-08 00:54 becker
* lib/BC_util.c: Addded a test for top or bottom nodal row back
into horizontal_bc (earlier moved from
Regional_boundary_conditions and Full_boundary_conditions to
BC_util.c) to see if this helps fix a bug reported by Dan that I
could not reproduce in the non-pyre version.
2010-10-12 06:36 becker
* lib/Anisotropic_viscosity.c, lib/Drive_solvers.c,
lib/Ggrd_handling.c, lib/Instructions.c: More accurate
computation of element coordinates in anisotropic part.
2010-09-24 03:35 becker
* lib/Anisotropic_viscosity.c, lib/Element_calculations.c,
lib/Instructions.c: Augmented G matrix with anisotropic
computation, may or may not be a good idea.
2010-09-15 10:40 becker
* lib/General_matrix_functions.c, lib/prototypes.h: Missed one
routine.
2010-09-15 10:36 becker
* lib/Anisotropic_viscosity.c, lib/BC_util.c, lib/Convection.c,
lib/Element_calculations.c, lib/General_matrix_functions.c,
lib/Instructions.c, lib/Obsolete.c, lib/Output.c,
lib/Size_does_matter.c, lib/Solver_multigrid.c,
lib/Topo_gravity.c, lib/anisotropic_viscosity.h,
lib/global_defs.h, lib/prototypes.h: Suggested use of general
prototypes
Retired N>3 matrix determinant functions (because they caused
compile problems)
2010-09-14 01:08 becker
* lib/Viscosity_structures.c: Read in sdepv_misfit regardless of
rheology.
2010-09-13 14:07 becker
* lib/Anisotropic_viscosity.c, lib/Drive_solvers.c,
lib/Element_calculations.c, lib/Ggrd_handling.c,
lib/Topo_gravity.c, lib/Viscosity_structures.c,
lib/anisotropic_viscosity.h, lib/global_defs.h,
lib/viscosity_descriptions.h: Added option to start with
isotropic viscosity during first anisotropic iteration.
2010-09-13 12:01 becker
* lib/Anisotropic_viscosity.c, lib/Element_calculations.c,
lib/Ggrd_handling.c, lib/Topo_gravity.c,
lib/anisotropic_viscosity.h: Forgot part of the rotation matrix,
minor fixes else, still experimental.
2010-09-13 07:13 becker
* configure.ac, lib/Anisotropic_viscosity.c,
lib/Element_calculations.c, lib/Ggrd_handling.c,
lib/Instructions.c, lib/Output.c, lib/Output_gzdir.c,
lib/Solver_multigrid.c, lib/Topo_gravity.c,
lib/Viscosity_structures.c, lib/anisotropic_viscosity.h,
lib/global_defs.h, lib/viscosity_descriptions.h: Made sure
current version compiled nicely without additional flags.
Anisotropic viscosity is set up more flexibly.
2010-09-07 07:05 becker
* lib/Ggrd_handling.c, lib/Instructions.c, lib/Output_gzdir.c,
lib/Viscosity_structures.c, lib/global_defs.h: Made sure that
viscosity prefactors get assigned for material dependent
viscosity only (TDEPV = off)
2010-09-07 02:22 becker
* lib/Citcom_init.c, lib/Tracer_setup.c: Added random number
generator initialization to ensure exact reproducibility of all
runs, including those with tracers.
Earlier reported solver issues might be related to our cluster,
at least all revisions seems to give the same answer now.
2010-09-07 01:40 becker
* lib/Anisotropic_viscosity.c, lib/Ggrd_handling.c,
lib/Nodal_mesh.c, lib/Viscosity_structures.c,
lib/viscosity_descriptions.h: Checking in current version because
of need for debugging.
Revision 17172, compared to (presumably) 17153 give strange
multigrid convergence behavior. Checking why that is now.
2010-09-06 05:20 becker
* configure.ac, lib/Anisotropic_viscosity.c, lib/Drive_solvers.c,
lib/Element_calculations.c, lib/Ggrd_handling.c,
lib/Instructions.c, lib/Makefile.am, lib/Nodal_mesh.c,
lib/Output.c, lib/Output_gzdir.c,
lib/Pan_problem_misc_functions.c, lib/Solver_multigrid.c,
lib/Topo_gravity.c, lib/Tracer_setup.c,
lib/Viscosity_structures.c, lib/anisotropic_viscosity.h,
lib/ggrd_handling.h, lib/global_defs.h,
lib/viscosity_descriptions.h: Trial implementation of anisotropic
viscosity.
2010-08-30 20:41 tan2
* CitcomS/Coupler/Coupler.py, module/Exchanger/BaseSVTInlet.cc,
module/Exchanger/BaseSVTInlet.h: Added support of multigrid
solver in Exchanger. Lifted the solver assertion in Coupler.
Cookbook9 can be converted to multigrid solver once appropriate
parameters (*solver.mesher.levels=2 and
*solver.vsolver.Solver=multigrid) are added. However, the
embedded solver will converge incorrectly. Setting
esolver.vsolver.accuracy=6e-3 seems to fix the problem.
2010-05-25 02:47 becker
* lib/BC_util.c: And a typo fix.
2010-05-24 18:00 becker
* lib/BC_util.c, lib/Ggrd_handling.c, lib/Initial_temperature.c,
lib/Instructions.c, lib/ggrd_handling.h, lib/global_defs.h:
toplayerbc > 0 will now assign BCs for all noes with r >=
toplayerbc_r, which defaults to 0.984303876942 i.e. 100 km depth.
(Still experimental.)
2010-03-12 00:41 becker
* lib/BC_util.c, lib/Ggrd_handling.c: Improvements for internal
BCs.
2010-03-10 19:55 becker
* lib/BC_util.c, lib/Ggrd_handling.c, lib/Topo_gravity.c: This
logic should make more sense for internal BCs, but still testing.
2010-03-10 19:38 becker
* lib/BC_util.c, lib/Full_boundary_conditions.c,
lib/Ggrd_handling.c, lib/Instructions.c,
lib/Regional_boundary_conditions.c, lib/Topo_gravity.c,
lib/global_defs.h: Experimental implementation of a single node
layer for internal boundary condition assignment, assigned if
toplayerbc < 0, to nodes at layer noz+toplayerbc, but only for
the top level multigrid, for now.
2010-02-26 21:04 walter
* lib/Topo_gravity.c: Fix a misleading comment
2010-02-02 23:41 walter
* Makefile.am: Fix an error when running autoreconf -i
2010-02-01 19:43 becker
* lib/Viscosity_structures.c, lib/viscosity_descriptions.h: Changed
internal variable name for use_se_smoothing, as per Eh's
suggestion.
2010-01-29 14:48 becker
* lib/Viscosity_structures.c, lib/viscosity_descriptions.h: Added
old CitcomS element -> nodes -> element viscosity smoothing with
new use_ne_visc_smooth parameter, default is off.
2010-01-21 21:23 tan2
* CitcomS/Components/Stokes_solver/Incompressible.py,
module/setProperties.c: Added pyre input for the new parameters
added in r16138.
2010-01-21 21:22 tan2
* module/Exchanger/CitcomInterpolator.cc: Replaced fixed accuracy
setting. Bug reported by Rob Moucha.
2010-01-21 20:11 tan2
* visual/OpenDXMacro/CitcomSImportOpt.net, visual/visFull.cfg,
visual/visFull.net.in: Added opendx macro to visualize opt (e.g.
composition) fields. This fixed Issue99.
2010-01-17 01:19 danb
* lib/Composition_related.c, lib/Tracer_setup.c, lib/tracer_defs.h:
Preliminary support for the absolute tracer method. Only tested
with one flavor, reading in from a file, and the initial volume
of the tracers is HARD-CODED! In addition, the code automaticaly
truncates the composition to one, although this should be a
switch for the user to specify.
2010-01-16 02:47 becker
* lib/Instructions.c, lib/Nodal_mesh.c, lib/Stokes_flow_Incomp.c,
lib/Viscosity_structures.c, lib/global_defs.h: Implementation of
two, perhaps entirely unnecessary features, in search of why
certain prescribed velocity models do not converge like they
should.
- allow for each inner loop solution to have a higher accuracy
than the outer, Uzawa loop. This is set by new parameter
"inner_accuracy_scale", which pre-multiplies the inner accuracy
(imp) setting (defaults to unity)
- remove any rigid body rotation at each Uzawa iteration step.
Parmaeter "inner_remove_rigid_rotation", by default off. For
this, I also added an experimental routine assign_v_to_vector()
in order to reasign the V velocities to the U solution vector
after removing the net rotations from V (how is supposed to work
for the regular operational mode?)
2009-12-09 02:25 tan2
* lib/BC_util.c: Replaced TRUE and FALSE. Added function
prototypes.
2009-12-08 21:44 tan2
* CitcomS/Components/Stokes_solver/Incompressible.py,
lib/Instructions.c, lib/Stokes_flow_Incomp.c, lib/global_defs.h,
module/setProperties.c: Added two boolean parameters for
controlling convergence criterion:
check_continuity_convergence (default=on): include div/v in the
criterion. check_pressure_convergence (default=on): include dp/p
in the criterion.
Retired the parameter: only_check_vel_convergence.
2009-12-08 21:43 tan2
* lib/Drive_solvers.c: removed unused variables in Drive_solvers.c
2009-12-08 21:43 tan2
* lib/Drive_solvers.c, lib/Element_calculations.c,
lib/Stokes_flow_Incomp.c, module/bindings.c,
module/stokes_solver.c, module/stokes_solver.h: Merged
solve_constrained_flow_iterative_pseudo_surf() with void
solve_constrained_flow_iterative(); assemble_forces_pseudo_surf()
wi
Removed used python bindings in module/stokes_solver.c
2009-12-08 17:19 becker
* lib/BC_util.c, lib/Ggrd_handling.c, lib/Topo_gravity.c: Truly
apply internal boundary conditions on horizontal surfaces for
E->mesh.toplayerbc > 0 (default, of course, 0)
2009-12-04 23:16 becker
* lib/BC_util.c: Modified handling of internal BCs, as per Eh's
suggestion, now handled in a separate function.
2009-12-04 21:57 tan2
* README: Removed metioning the version number in README. Added the
url of the manual.
2009-12-04 17:44 becker
* lib/BC_util.c, lib/Full_boundary_conditions.c,
lib/Full_read_input_from_files.c, lib/Ggrd_handling.c,
lib/Regional_boundary_conditions.c,
lib/Regional_read_input_from_files.c, lib/ggrd_handling.h:
Suggested unified treatment of horizontal_bc (moved to BC_util)
for review by Eh
Improved handling of grd-read velocity and traction boundary
conditions.
2009-12-04 02:33 becker
* lib/Full_boundary_conditions.c, lib/Ggrd_handling.c,
lib/Regional_boundary_conditions.c: Temporary fix for regular
operations, for Eh's review.
Ggrd assignment still needs tending to.
2009-12-03 23:51 becker
* lib/BC_util.c, lib/Full_boundary_conditions.c,
lib/Full_read_input_from_files.c, lib/Ggrd_handling.c,
lib/Instructions.c, lib/Regional_boundary_conditions.c,
lib/Regional_read_input_from_files.c, lib/ggrd_handling.h,
lib/global_defs.h: Experimental implementation of internal stress
or velocity boundary conditions. As of right now, toplayerbc > 0
(default = 0) will assign the surface boundary condition to all
nodes within the first toplayerbc layers. If BCs are assigned
from file, this will only work for the grd method, not the
regular file I/O.
2009-11-21 00:15 becker
* lib/Instructions.c, lib/Stokes_flow_Incomp.c, lib/global_defs.h:
Exit criterion for incompressibility iteration allowed for
velocity and pressure convergence at finite compressibility >
accuracy. Changed criterion to loop while (convergence < 2) OR
(compressibility < imp)
Please review
Also modified the sloppy iteration branch, where pressure isn't
checked.
2009-11-20 21:58 tan2
* lib/Convection.c, lib/Parsing.c: Clean up #include
2009-11-20 21:58 tan2
* lib/Material_properties.c: Fixed two minor errors
2009-11-19 23:11 tan2
* lib/Full_version_dependent.c, lib/Global_operations.c,
lib/Instructions.c, lib/Material_properties.c,
lib/Regional_version_dependent.c, lib/global_defs.h: Created a
new array E->sphere.gr for the radii of all layers. This makes
gather_r() obsolete.
2009-11-19 23:09 tan2
* lib/Viscosity_structures.c, lib/viscosity_descriptions.h:
refactor viscosity definition
2009-11-19 23:08 tan2
* visual/batchcombine.py: 'localhost' can abbr. as '-'
2009-10-14 18:26 becker
* lib/Instructions.c: Changed defaults for z_410, z_lmantle, z_lith
to more Earth like parameters.
Does this break any tests?
2009-10-13 16:37 becker
* lib/Ggrd_handling.c: Fixed typo.
2009-10-08 18:01 tan2
* CitcomS/Components/Visc.py: z_layer needs 4 elements
2009-10-07 04:39 danb
* CitcomS/Components/Stokes_solver/Incompressible.py: added
only_check_vel_convergence to pyre version.
2009-10-06 18:48 becker
* lib/Full_read_input_from_files.c, lib/Ggrd_handling.c,
lib/Instructions.c, lib/Regional_read_input_from_files.c:
Modified handling of ggrd_mat_control, > 0 will select all layers
above, < 0 will select individual layer == - ggrd_mat_control
2009-10-05 19:06 tan2
* CitcomS/Components/Sphere/Sphere.py, CitcomS/Components/Visc.py,
module/setProperties.c: Add pyre binding of z_layer,
r_grid_layers, rr, and nr.
replace MAX_MAT with CITCOM_MAX_VISC_LAYER.
2009-10-05 19:05 tan2
* lib/Instructions.c: Refactoring z_layer input validation.
2009-10-05 19:05 tan2
* lib/Full_tracer_advection.c, lib/Tracer_setup.c: Dealing with
tracers that belong to the current processor but cannot find its
element.
When itracer_warning is off, the code will ignore the lost
tracers. When itracer_warning is on (default), the code will
abort.
2009-10-02 23:16 becker
* lib/Full_boundary_conditions.c, lib/Ggrd_handling.c,
lib/Instructions.c, lib/Regional_boundary_conditions.c:
Experimental traction surface BC from grds.
2009-09-28 20:08 tan2
* CitcomS/Components/Visc.py, lib/Construct_arrays.c,
lib/Pan_problem_misc_functions.c, lib/Viscosity_structures.c,
lib/viscosity_descriptions.h, module/setProperties.c:
visc_layer_control: allowing each radial layer has its own
viscosity parameters
When "visc_layer_control" is true, the file "visc_layer_file" is
read in. The format of visc_layer_file:
param_nameA Aval_1 Aval_2 ... Aval_nel
param_nameB Bval_1 Bval_2 ... Bval_nel
. . .
The first line is assigned to the basal layer (CMB), and the last
line to the surface layer. Note that each parameter has nodez-1
values.
2009-09-28 20:05 tan2
* lib/viscosity_descriptions.h: clean up unused viscosity variables
2009-09-28 20:05 tan2
* lib/Element_calculations.c: Fixed a bug in pseudo free surface,
reported by Robert Moucha.
2009-09-28 20:04 tan2
* lib/Topo_gravity.c: remove duplicating prototype declaration
2009-09-03 17:24 tan2
* CitcomS/Components/Output.py, CitcomS/Solver/Solver.py,
module/bindings.c, module/outputs.c, module/outputs.h,
module/setProperties.c: Added pyre binding for write_q_files
2009-09-03 17:24 tan2
* CitcomS/Solver/CoupledSolver.py: Added a warning message about
tracer + coupling
2009-08-25 19:54 becker
* lib/Stokes_flow_Incomp.c: Forgot one valid check for Uzawa
iterations.
2009-08-25 17:31 becker
* lib/Stokes_flow_Incomp.c: Ensure that Uzawa solver does not
converge if low level, MG solver does not yield a valid solution.
2009-08-21 20:32 becker
* lib/Viscosity_structures.c, lib/viscosity_descriptions.h: test
implementation of steady state strain-rate weakending via PDEPV
psrw=on parameter.
2009-08-11 17:48 tan2
* lib/Material_properties.c: Add Murnaghan's integrated linear EoS
as reference_state=2
2009-08-11 17:47 tan2
* lib/Global_operations.c: Add gather_r() to gather all radial
coordinates
2009-08-11 17:35 tan2
* lib/Topo_gravity.c: Adding augmented Lagrangian component of
stiffness matrix in cbf topo
2009-08-04 01:45 leif
* lib/multigrid_kernel.cu: The function collect_terms() turns out
to be quite expensive (over 70% according to gprof). Perform it
once per level, caching the result.
2009-07-29 20:16 leif
* lib/multigrid_kernel.cu: Fixed the remaining bug in
n_assemble_del2_u(). Glancing at "global_defs.h", I had assumed
that VBX, VBY, and VBZ were defined in the natrual way:
#define VBX 0x2 #define VBY 0x4 #define VBZ 0x8
I.e., "0x1 << doff"... but today I looked at it more closely; Y &
Z are reversed!!!
#define VBX 0x2 #define VBZ 0x4 #define VBY 0x8
With this fix, the CUDA multigrid solver seems to work correctly.
2009-07-29 04:05 becker
* lib/Ggrd_handling.c: Ensure debugging output for ggrd_vtop mixed
BC if more than one CPU with depth is used.
Previous check in of precise_strain_rate=on mode seems to be OK,
stress, geoid, and dynamic topo output (for cbf topo off) have
been tested.
2009-07-28 01:00 leif
* lib/multigrid_kernel.cu: Experiments show that using a block size
of 1 leads to poor performance. Using the recommended minimum of
64, gauss_seidel_0 and gauss_seidel_1 are about 3x to 4x faster.
2009-07-27 23:44 leif
* lib/multigrid_kernel.cu: Fixed stupid reduce bug that affected
several routines. (This is what I get for coping & pasting
example code without thinking.) Now the only remaining broken
routine is n_assemble_del2_u().
2009-07-25 02:48 leif
* lib/multigrid_kernel.cu: Wrote two new versions of
gauss_seidel_2() -- neither of which is very parallel. Both
produce similar wrong answers.
2009-07-24 01:03 becker
* lib/Topo_gravity.c: This version of precise strain-rates may now
be worth testing out.
2009-07-23 18:40 becker
* lib/Topo_gravity.c, lib/Viscosity_structures.c: Corrected call to
get_ba in Topo_gravity for precise_strain_rate, still testing.
Rest should be cosmetic.
2009-07-23 03:50 leif
* lib/multigrid_kernel.cu: This junk is my effort to understand
gauss_seidel_2(). Because of the way E->temp[] is used, it is
quite difficult (impossible?) to parallelize. I even wrote code
to generate a "dot" graph of the expression tree. Of course, this
tree is too big to plot in its entirety. With neq == 2187, the
tree has 87479 nodes and is 355 levels deep; each level can have
anywhere from a handful to several hundred nodes, with over a
thousand nodes in each of the bottommost layers.
2009-07-22 20:17 becker
* lib/Topo_gravity.c: Tentative check in of stress computation
modification such that Eh can take a look of what we're working
on.
2009-07-17 21:22 leif
* lib/multigrid_kernel.cu: Created drop-in, host-side replacements
for each of the multigrid kernel functions. After swapping them
in/out, it seems n_assemble_del2_u(), gauss_seidel_2(), and
gauss_seidel_3() are all buggy :-( Both gauss_seidel_0() and
gauss_seidel_1() appear to be OK... but that isn't worth writing
home about, because they are both trivial.
2009-07-16 03:13 leif
* lib/multigrid_kernel.cu: New scaffolding with annotations.
2009-07-16 00:49 leif
* lib/multigrid_kernel.cu: Bug fixes... still doesn't work right.
2009-07-15 21:54 leif
* lib/multigrid_kernel.cu: Wrote code to allocate & initialize
device memory, and copy the result from the device.
2009-07-15 01:58 leif
* lib/multigrid_kernel.cu: Parallelized the last fragment of the
gauss_seidel() main loop as gauss_seidel_3(). It turns out to be
very similar to the second half of n_assemble_del2_u(). Tore down
scaffolding.
2009-07-15 00:14 leif
* lib/multigrid_kernel.cu: Parallel n_assemble_del2_u() takes
shape.
2009-07-10 22:15 leif
* lib/multigrid_kernel.cu: Pardon my dust...
tally_n_assemble_del2_u() is instrumention and printfs I wrote in
order to understand n_assemble_del2_u(). As with
e_assemble_del2_u(), it appears I need one block for each Au[i]
element, with a reduction in shared memory.
2009-07-10 21:08 danb
* lib/Process_buoyancy.c: include density (rho) scaling in
convective heat flux (necessary for compressible flow)
2009-07-09 20:40 danb
* lib/Regional_read_input_from_files.c: Second material element
file only read in and processed when (pos_age) criteria met (i.e.
consistent with vel b.c. approach).
2009-07-09 17:27 becker
* lib/Instructions.c: Only one processor complains about removing
rigid rotations while prescribing plate motions.
2009-07-06 22:13 leif
* lib/General_matrix_functions.c, lib/cgrad_kernel.cu,
lib/multigrid_kernel.cu: Make sure "-arch sm_13" is given -- for
double-precision floating-point. Added missing #ifndef USE_CUDA.
2009-07-02 18:41 leif
* lib/multigrid_kernel.cu: Sketched data-parallel version of
gauss_seidel(). I haven't figured out how to parallelize the
following:
Au[C[i]] += ... Ad[C[j]] += ...
2009-07-01 19:09 leif
* lib/General_matrix_functions.c, lib/cgrad_kernel.cu,
lib/multigrid_kernel.cu: Fixed a pair of mistakes I made: under
cgrad, don't replace solve_del2_u(); under multigrid, 'level' is
not a constant.
2009-06-30 22:33 leif
* lib/multigrid_kernel.cu: The 'temp1' array is unnecessary since
nproc == 1.
2009-06-29 22:10 leif
* lib/multigrid_kernel.cu: For my own sanity, simplified code
assuming the following:
E->sphere.caps_per_proc == 1 E->parallel.nproc == 1
E->mesh.levmax == 0 E->mesh.nsd == 3
2009-06-29 20:18 leif
* lib/General_matrix_functions.c, lib/Makefile.am,
lib/multigrid_kernel.cu: Proposed CUDA kernel for multigrid
solver.
2009-06-25 16:51 tan2
* examples/Cookbook10/cmt-solution.cookbook10,
examples/Cookbook10/stations.cookbook10: forgot to add these
files
2009-06-15 18:21 tan2
* ChangeLog: Update ChangeLog to r15180
2009-06-10 20:42 tan2
* configure.ac: Added COND_GGRD flags; removed redundant LIBS
2009-06-10 20:41 tan2
* configure.ac: Checking malloc.h earlier to initialize even when
not using pyre
2009-06-10 17:49 tan2
* configure.ac: checking netcdf header and library
2009-06-10 04:11 tan2
* configure.ac: Append GMTHOME/lib to LD path
2009-06-09 23:47 tan2
* configure.ac: Adding ggrd compilation support
2009-06-08 21:03 tan2
* NEWS: One more feature and mention where to find the features in
the cookbooks.
2009-06-08 18:11 tan2
* CitcomS/Coupler/EmbeddedCoupler.py: Added a safe guard
2009-06-04 21:04 tan2
* examples/Cookbook9/cookbook9.cfg: Reordering parameter
2009-06-03 21:08 tan2
* module/Exchanger/BoundarySVTInlet.cc,
module/Exchanger/SVTInlet.cc, module/Exchanger/SVTInlet.h: Update
SVTInlet to inherit from BaseSVTInlet
2009-06-03 21:07 tan2
* module/Exchanger/Makefile.am: adding files to makefile
2009-06-02 22:58 tan2
* tests/coupled.cfg, tests/exchange.py: update to current coupler
2009-06-02 22:58 tan2
* examples/Cookbook9/cookbook9.cfg: reducing the accuracy of
solvers, otherwise, they will converge to the wrong result
2009-06-02 22:56 tan2
* CitcomS/Coupler/ContainingCoupler.py, CitcomS/Coupler/Coupler.py,
CitcomS/Coupler/EmbeddedCoupler.py, CitcomS/Coupler/Inlet.py,
CitcomS/Coupler/Outlet.py, CitcomS/Solver/CoupledSolver.py,
module/Exchanger/AreaWeightedNormal.cc,
module/Exchanger/AreaWeightedNormal.h,
module/Exchanger/BaseSVTInlet.cc,
module/Exchanger/BaseSVTInlet.h,
module/Exchanger/BoundarySVTInlet.cc,
module/Exchanger/BoundarySVTInlet.h,
module/Exchanger/BoundaryVTInlet.cc,
module/Exchanger/CitcomInterpolator.cc,
module/Exchanger/Makefile.am, module/Exchanger/PInlet.cc,
module/Exchanger/PInlet.h, module/Exchanger/PInterior.cc,
module/Exchanger/PInterior.h, module/Exchanger/POutlet.cc,
module/Exchanger/POutlet.h, module/Exchanger/bindings.cc,
module/Exchanger/exchangers.cc, module/Exchanger/exchangers.h,
module/Exchanger/inlets_outlets.cc,
module/Exchanger/inlets_outlets.h: Added parameters
'amending_outflow' and 'exchange_pressure' to help the
convergence of esolver
When 'amending_outflow' is set to true, the imposed velocity BC
will be amended slightly to be divergence-free. (The divergence,
e.g. outflow, is caused by the combination of solver inaccuracy
and interpolation inaccuracy.)
When 'exchange_pressure' is set to true, the initial pressure (at
element level) of the embedded solver is taken from the pressure
solution of the containing solver.
2009-06-02 22:09 tan2
* module/Exchanger/Boundary.cc, module/Exchanger/Boundary.h:
Boundary for global mesh, not sure whether it works with
interpolator or not
2009-06-02 21:58 tan2
* lib/Instructions.c: init'd E->output.seismic and
E->output.coord_bin
2009-06-02 21:58 tan2
* examples/Cookbook10/cookbook10.cfg: refining cookbook10