aboutsummaryrefslogtreecommitdiff
path: root/dist/des.patch
diff options
context:
space:
mode:
Diffstat (limited to 'dist/des.patch')
-rw-r--r--dist/des.patch109
1 files changed, 62 insertions, 47 deletions
diff --git a/dist/des.patch b/dist/des.patch
index 2434daa..e967ee1 100644
--- a/dist/des.patch
+++ b/dist/des.patch
@@ -4,16 +4,25 @@ RCS file: /cvsroot/Lorene/C++/Source/Bin_hor/binhor_hh.C,v
retrieving revision 1.3
diff -u -r1.3 binhor_hh.C
--- Lorene/C++/Source/Bin_hor/binhor_hh.C 9 Jan 2008 14:28:58 -0000 1.3
-+++ Lorene/C++/Source/Bin_hor/binhor_hh.C 30 Nov 2009 03:47:22 -0000
-@@ -892,7 +892,7 @@
- // des_profile(hh1(i,j), 0., 20., M_PI/2, M_PI) ;
- //des_profile(hh1(i,j), 0., 20., M_PI/2, 0) ;
++++ Lorene/C++/Source/Bin_hor/binhor_hh.C 27 Mar 2014 13:38:39 -0000
+@@ -885,7 +885,7 @@
+ cout << hole1.mp.phi << endl ;
+ cout << hole1.mp.tet << endl ;
+
+-
++ /*
+ //des_meridian(hh1, 0., 20., "hh1 cart", 20) ;
+ for (int i=1 ; i<= 3 ; i++)
+ for (int j=i ; j<= 3 ; j++){
+@@ -894,7 +894,7 @@
//des_profile(hh1(i,j), 0., 20., 0, M_PI) ;
-- des_coupe_z (hh1(i,j), 0., 5) ;
-+ //des_coupe_z (hh1(i,j), 0., 5) ;
+ des_coupe_z (hh1(i,j), 0., 5) ;
}
-
+-
++ */
hh1.change_triad(hole1.mp.get_bvect_spher()) ;
+ hh2.change_triad(hole2.mp.get_bvect_spher()) ;
+
@@ -908,7 +908,7 @@
hole2.tgam = tgam_2 ;
@@ -29,41 +38,45 @@ RCS file: /cvsroot/Lorene/C++/Source/Isol_hor/init_data.C,v
retrieving revision 1.28
diff -u -r1.28 init_data.C
--- Lorene/C++/Source/Isol_hor/init_data.C 19 Aug 2008 06:42:00 -0000 1.28
-+++ Lorene/C++/Source/Isol_hor/init_data.C 30 Nov 2009 03:47:22 -0000
-@@ -1262,10 +1262,10 @@
++++ Lorene/C++/Source/Isol_hor/init_data.C 27 Mar 2014 13:38:39 -0000
+@@ -1261,13 +1261,13 @@
+
for (int mer=0; mer<niter; mer++) {
-
-- des_meridian(psi_j, 1, 10., "psi", 0) ;
-- des_meridian(nn_j, 1, 10., "nn", 1) ;
-- des_meridian(theta_j, 1, 10., "Theta", 2) ;
-- des_meridian(chi_j, 1, 10., "chi", 3) ;
-+ //des_meridian(psi_j, 1, 10., "psi", 0) ;
-+ //des_meridian(nn_j, 1, 10., "nn", 1) ;
-+ //des_meridian(theta_j, 1, 10., "Theta", 2) ;
-+ //des_meridian(chi_j, 1, 10., "chi", 3) ;
+-
++ /*
+ des_meridian(psi_j, 1, 10., "psi", 0) ;
+ des_meridian(nn_j, 1, 10., "nn", 1) ;
+ des_meridian(theta_j, 1, 10., "Theta", 2) ;
+ des_meridian(chi_j, 1, 10., "chi", 3) ;
arrete() ;
+-
++ */
-
-@@ -1409,8 +1409,8 @@
+ //========
+ // Sources
+@@ -1408,11 +1408,11 @@
+
// Iteration in theta/chi
for (int i=0 ; i<niter ; i++) {
-
-- des_meridian(theta_i, 1, 10., "Theta", 2) ;
-- des_meridian(chi_i, 1, 10., "chi", 3) ;
-+ //des_meridian(theta_i, 1, 10., "Theta", 2) ;
-+ //des_meridian(chi_i, 1, 10., "chi", 3) ;
+-
++ /*
+ des_meridian(theta_i, 1, 10., "Theta", 2) ;
+ des_meridian(chi_i, 1, 10., "chi", 3) ;
arrete() ;
+-
++ */
+ //Sources
Index: Lorene/C++/Source/Time_slice/tslice_conf_init.C
===================================================================
RCS file: /cvsroot/Lorene/C++/Source/Time_slice/tslice_conf_init.C,v
-retrieving revision 1.10
-diff -u -r1.10 tslice_conf_init.C
---- Lorene/C++/Source/Time_slice/tslice_conf_init.C 4 Dec 2008 18:22:49 -0000 1.10
-+++ Lorene/C++/Source/Time_slice/tslice_conf_init.C 30 Nov 2009 03:47:23 -0000
-@@ -223,7 +223,7 @@
+retrieving revision 1.11
+diff -u -r1.11 tslice_conf_init.C
+--- Lorene/C++/Source/Time_slice/tslice_conf_init.C 20 Oct 2010 07:58:09 -0000 1.11
++++ Lorene/C++/Source/Time_slice/tslice_conf_init.C 27 Mar 2014 13:38:39 -0000
+@@ -225,7 +225,7 @@
maxabs(psi_jp1.laplacian() - source_psi,
"Absolute error in the resolution of the equation for Psi") ;
@@ -72,7 +85,7 @@ diff -u -r1.10 tslice_conf_init.C
// Resolution of the Poisson equation for the lapse
// ------------------------------------------------
-@@ -236,7 +236,7 @@
+@@ -238,21 +238,21 @@
maxabs(nn_jp1.laplacian() - source_nn,
"Absolute error in the resolution of the equation for N") ;
@@ -81,10 +94,10 @@ diff -u -r1.10 tslice_conf_init.C
// Resolution of the vector Poisson equation for the shift
//---------------------------------------------------------
-@@ -244,12 +244,14 @@
+
Vector beta_jp1 = source_beta.poisson(0.3333333333333333, ff,
method_poisson_vect) ;
-
+-
+ /*
des_meridian(beta_jp1(1), 0., ray_des, "\\gb\\ur\\d", ngraph0+2,
graph_device) ;
@@ -92,30 +105,32 @@ diff -u -r1.10 tslice_conf_init.C
graph_device) ;
des_meridian(beta_jp1(3), 0., ray_des, "\\gb\\u\\gf\\d", ngraph0+4,
graph_device) ;
-+ */
-
+-
++ */
// Test:
Vector test_beta = (beta_jp1.derive_con(ff)).divergence(ff)
+ + 0.3333333333333333 * (beta_jp1.divergence(ff)).derive_con(ff) ;
Index: Lorene/C++/Source/Time_slice/tslice_dirac_max_evolve.C
===================================================================
RCS file: /cvsroot/Lorene/C++/Source/Time_slice/tslice_dirac_max_evolve.C,v
-retrieving revision 1.16
-diff -u -r1.16 tslice_dirac_max_evolve.C
---- Lorene/C++/Source/Time_slice/tslice_dirac_max_evolve.C 4 Dec 2008 18:22:49 -0000 1.16
-+++ Lorene/C++/Source/Time_slice/tslice_dirac_max_evolve.C 30 Nov 2009 03:47:23 -0000
-@@ -673,6 +673,7 @@
+retrieving revision 1.20
+diff -u -r1.20 tslice_dirac_max_evolve.C
+--- Lorene/C++/Source/Time_slice/tslice_dirac_max_evolve.C 24 Jan 2013 12:55:18 -0000 1.20
++++ Lorene/C++/Source/Time_slice/tslice_dirac_max_evolve.C 27 Mar 2014 13:38:39 -0000
+@@ -688,7 +688,7 @@
+ beta_new = relax * beta_new + (1.-relax) * beta() ;
beta_evol.update(beta_new, jtime, ttime) ;
}
-
-+ /*
- des_meridian(beta()(1), 0., ray_des, "\\gb\\ur\\d", ngraph0+6,
+-
++ /*
+ des_meridian(vec_X()(1), 0., ray_des, "\\gb\\ur\\d", ngraph0+6,
graph_device) ;
- des_meridian(beta()(2), 0., ray_des, "\\gb\\u\\gh\\d", ngraph0+7,
-@@ -695,6 +696,7 @@
+ des_meridian(vec_X()(2), 0., ray_des, "\\gb\\u\\gh\\d", ngraph0+7,
+@@ -713,6 +713,7 @@
graph_device) ;
arrete(nopause) ;
-+ */
-
++ */
}
+ par_A.clean_all() ;
generated by cgit v1.2.3 (git 2.39.1) at 2024-06-10 12:56:10 +0000