remove set_acceleration from all interfaces

src/amuse/community/adaptb/interface.cc

@@ -443,5 +443,4 @@ int get_kinetic_energy(double* ek){return -2;}
 int get_potential_energy(double* ep){return -2;}
 int get_center_of_mass_velocity(double* vx, double* vy, double* vz){return -2;}
 int get_acceleration(int id, double* ax, double* ay, double* az){return -2;}
-int set_acceleration(int id, double ax, double ay, double az){return -2;}
 

src/amuse/community/bhtree/interface.cc

@@ -564,12 +564,6 @@ int get_number_of_particles(int *number_of_particles)
   return 0; 
 }
 
-int set_acceleration(int index_of_the_particle, double ax, double ay, double az)
-{
-  return -2;
-}
-
-
 int get_center_of_mass_position(double *x, double *y, double *z)
 {
   nbody_particle *np = bhtcs.get_particle_pointer();

src/amuse/community/bonsai/interface.py

@@ -143,31 +143,6 @@ def set_velocity():
         """
         return function        
 
-    @legacy_function
-    def set_acceleration():
-        """
-        Set the velocity vector of a particle.
-        """
-        function = LegacyFunctionSpecification()
-        function.must_handle_array = True
-        function.addParameter('index_of_the_particle', dtype='int32', direction=function.IN,
-            description = "Index of the particle to get the state from. This index must have been returned by an earlier call to :meth:`new_particle`")
-        function.addParameter('ax', dtype='float64', direction=function.IN, description = "The current x component of the velocity vector of the particle")
-        function.addParameter('ay', dtype='float64', direction=function.IN, description = "The current y component of the velocity vector of the particle")
-        function.addParameter('az', dtype='float64', direction=function.IN, description = "The current z component of the velocity vector of the particle")
-        function.addParameter('length', 'int32', function.LENGTH)
-        function.result_type = 'int32'
-        function.result_doc = """
-        0 - OK
-            current value was retrieved
-        -1 - ERROR
-            particle could not be found
-        -2 - ERROR
-            not yet implemented
-        """
-        return function        
-
-
     @legacy_function
     def get_theta_for_tree():
         """

src/amuse/community/bonsai/libraryInterface.cpp

@@ -550,24 +550,6 @@ int get_potential(int id, double * potential){
   return 0;
 }
 
-int set_acceleration(int *id, double *ax, double *ay, double *az, int length){
-  getCurrentStateToHost();
-
-  for(int i=0; i < length; i++)
-  {
-    int index_of_the_particle = getIdxFromId(id[i]);  
-    if(index_of_the_particle < 0)    return -3;
-
-    bonsai->localTree.bodies_acc1[index_of_the_particle].x = ax[i];
-    bonsai->localTree.bodies_acc1[index_of_the_particle].y = ay[i];
-    bonsai->localTree.bodies_acc1[index_of_the_particle].z = az[i];
-  }
-
-  bonsai->localTree.bodies_acc1.h2d();
-
-  return 0;
-}
-
 int get_acceleration(int id, double * ax, double * ay, double * az){
   getCurrentStateToHost();
   int index_of_the_particle = getIdxFromId(id);  

src/amuse/community/bonsai2/interface.py

@@ -136,31 +136,6 @@ def set_velocity():
         """
         return function        
 
-    @legacy_function
-    def set_acceleration():
-        """
-        Set the velocity vector of a particle.
-        """
-        function = LegacyFunctionSpecification()
-        function.must_handle_array = True
-        function.addParameter('index_of_the_particle', dtype='int32', direction=function.IN,
-            description = "Index of the particle to get the state from. This index must have been returned by an earlier call to :meth:`new_particle`")
-        function.addParameter('ax', dtype='float64', direction=function.IN, description = "The current x component of the velocity vector of the particle")
-        function.addParameter('ay', dtype='float64', direction=function.IN, description = "The current y component of the velocity vector of the particle")
-        function.addParameter('az', dtype='float64', direction=function.IN, description = "The current z component of the velocity vector of the particle")
-        function.addParameter('length', 'int32', function.LENGTH)
-        function.result_type = 'int32'
-        function.result_doc = """
-        0 - OK
-            current value was retrieved
-        -1 - ERROR
-            particle could not be found
-        -2 - ERROR
-            not yet implemented
-        """
-        return function        
-
-
     @legacy_function
     def get_theta_for_tree():
         """

src/amuse/community/bonsai2/libraryInterface.cpp

@@ -592,24 +592,6 @@ int get_potential(int id, double * potential){
   return 0;
 }
 
-int set_acceleration(int *id, double *ax, double *ay, double *az, int length){
-  getCurrentStateToHost();
-
-  for(int i=0; i < length; i++)
-  {
-    int index_of_the_particle = getIdxFromId(id[i]);  
-    if(index_of_the_particle < 0)    return -3;
-
-    bonsai->localTree.bodies_acc1[index_of_the_particle].x = ax[i];
-    bonsai->localTree.bodies_acc1[index_of_the_particle].y = ay[i];
-    bonsai->localTree.bodies_acc1[index_of_the_particle].z = az[i];
-  }
-
-  bonsai->localTree.bodies_acc1.h2d();
-
-  return 0;
-}
-
 int get_acceleration(int id, double * ax, double * ay, double * az){
   getCurrentStateToHost();
   int index_of_the_particle = getIdxFromId(id);  

src/amuse/community/brutus/interface.cc

@@ -675,5 +675,3 @@ int get_center_of_mass_velocity(double* vx, double* vy, double* vz){
 }
 
 int get_acceleration(int id, double* ax, double* ay, double* az){return -2;}
-int set_acceleration(int id, double ax, double ay, double az){return -2;}
-

src/amuse/community/etics/interface.cu

@@ -282,10 +282,6 @@ int recommit_particles() {
     return -2;
 }
 
-int set_acceleration(int index_of_the_particle, double ax, double ay, double az) {
-    return -2;
-}
-
 int get_center_of_mass_position(double *x, double *y, double *z) {
 //     vec3 CenterMass = thrust::transform_reduce(PPPPP.begin(), PPPPP.end(), CenterMassFunctor(mass), vec3(0,0,0), thrust::plus<vec3>());
 //     CenterMass = CenterMass * (1.0/N); //ugly should divide by the total mass
@@ -375,4 +371,4 @@ int update_force_potential_arrays(double tttt) {
 #warning time shouldnt be a parameter to this one
 //     ForceSCF(N, PTR(PotPotPot), PTR(PPPPP), PTR(F0xxxxxx));
     return 0;
-}
+}

src/amuse/community/fi/interface.py

@@ -36,7 +36,6 @@ class FiInterface(
     get_total_mass=None
     get_center_of_mass_position=None
     get_center_of_mass_velocity=None
-    set_acceleration=None
     get_acceleration=None
 
     use_modules=['StoppingConditions','AmuseInterface']

src/amuse/community/gadget2/interface.cc

@@ -2113,10 +2113,6 @@ int get_acceleration(int *index, double * ax, double * ay, double * az, int leng
     return result;
 }
 
-int set_acceleration(int index, double ax, double ay, double az){
-    return -2;
-}
-
 int get_internal_energy(int *index, double *internal_energy, int length){
     int errors = 0;
     double *buffer = new double[length];

src/amuse/community/hacs64/interface.cc

@@ -592,20 +592,6 @@ int get_velocity(int index_of_the_particle,
 
 /****************/
 
-int set_acceleration(int index_of_the_particle,
-    double ax, double ay, double az)
-{
-    try
-    {
-        assert(nbody_ptr != NULL);
-        assert(nbody_ptr->is_sane());
-        return -1;
-    }
-    catch(assert_failed ex)
-    {
-        return handle_assert_failed(ex);
-    }    
-}
 int get_acceleration(int index_of_the_particle,
     double * ax, double * ay, double * az)
 {

src/amuse/community/hermite/interface.cc

@@ -1287,24 +1287,6 @@ int get_acceleration(int id, double *ax, double *ay, double *az)
     }
 }
 
-int set_acceleration(int id, double ax, double ay, double az)
-{
-  unsigned int i = find(ident.begin(), ident.end(), id) - ident.begin();
-  if (i < ident.size())
-    {
-      acc[i][0] = ax;
-      acc[i][1] = ay;
-      acc[i][2] = az;
-      return 0;
-    }
-  else
-    {
-      return -2;
-    }
-
-  return -3;
-}
-
 int evolve_not_on_root() {
 
 #ifndef NOMPI

src/amuse/community/higpus/interface.cc

@@ -593,15 +593,6 @@ int recommit_particles(){
    return 0;
 }
 
-int set_acceleration(int index_of_the_particle, double ax, double ay, double az){
-	index_of_the_particle = 0;
-	ax = 0;
-	ay = 0;
-	az = 0;
-   if (rank_higpus == 0) cout<<"accelerations are stored only on GPU, you can't set it"<<endl;
-   return 0;
-}
-
 int get_center_of_mass_position(double * x, double * y, double * z){
    unsigned int ppG = (unsigned int) ceil ( (double) M / size_higpus);
 	double *sum = new double[4];

src/amuse/community/huayno/interface.c

@@ -574,11 +574,6 @@ int get_potential(int id, double *pot)
   return 0;
 }
 
-int set_acceleration(int id, double ax, double ay, double az)
-{
-  return -2;
-}
-
 int get_acceleration(int id, double *ax, double *ay, double *az)
 {
   return -2;

src/amuse/community/interface/gd.py

@@ -506,46 +506,6 @@ def get_acceleration():
         """
         return function
 
-    @legacy_function
-    def set_acceleration():
-        """
-        Update the acceleration of a particle.
-        *Defined for symetry with the get_acceleration function.*
-        *Should be removed if physaccily unsound*
-        *Maybe moved to snapshot support functionality*
-        """
-        function = LegacyFunctionSpecification()
-        function.addParameter(
-            'index_of_the_particle', dtype='int32', direction=function.IN,
-            description=(
-                "Index of the particle for which the state is to be updated. "
-                "This index must have been returned by an earlier call to "
-                ":meth:`new_particle`"
-            )
-        )
-        function.addParameter(
-            'ax', dtype='float64', direction=function.IN,
-            description="The new acceleration vector of the particle")
-        function.addParameter(
-            'ay', dtype='float64', direction=function.IN,
-            description="The new acceleration vector of the particle")
-        function.addParameter(
-            'az', dtype='float64', direction=function.IN,
-            description="The new acceleration vector of the particle")
-        function.result_type = 'int32'
-        function.can_handle_array = True
-        function.result_doc = """
-        0 - OK
-            particle was found in the model and the information was set
-        -1 - ERROR
-            particle could not be found
-        -2 - ERROR
-            code does not support updating of a particle
-        -3 - ERROR
-            not yet implemented
-        """
-        return function
-
     @legacy_function
     def get_potential():
         """

src/amuse/community/kepler_orbiters/interface.cc

@@ -69,10 +69,6 @@ int cleanup_code()
     return 0;
 }
 
-int set_acceleration(int particle_identifier, double ax, double ay, double az) {
-    return -2; // Not implemented
-}
-
 int get_acceleration(int particle_identifier, double *ax, double *ay, double *az) {
     return -2; // Not implemented
 }

src/amuse/community/mi6/interface.cc

@@ -849,10 +849,6 @@ int get_state(int particle_identifier, double *mass,
     return -3; // Not found!
 }
 
-int set_acceleration(int particle_identifier, double ax, double ay, double az) {
-    return -2; // Not implemented
-}
-
 int get_acceleration(int particle_identifier, double *ax, double *ay, double *az) {
     int index;
     if (found_particle(particle_identifier, &index)){

src/amuse/community/mikkola/interface.f90

@@ -596,14 +596,6 @@ FUNCTION get_number_of_particles(number_of_particles)
   get_number_of_particles=number_of_particles_allocated
 END FUNCTION
 
-FUNCTION set_acceleration(index_of_the_particle, ax, ay, az)
-  IMPLICIT NONE
-  INTEGER :: set_acceleration
-  INTEGER :: index_of_the_particle
-  DOUBLE PRECISION :: ax, ay, az
-  set_acceleration=0
-END FUNCTION
-
 FUNCTION get_indices_of_colliding_particles(index_of_particle1,  &
     index_of_particle2)
   IMPLICIT NONE

src/amuse/community/nbody6xx/interface.f90

@@ -301,13 +301,6 @@ FUNCTION get_acceleration(index_of_the_particle, ax_amuse,
         get_acceleration = 0
       END FUNCTION
 
-      FUNCTION set_acceleration(index_of_the_particle, ax, ay, az)
-        INTEGER :: index_of_the_particle
-        DOUBLE PRECISION :: ax, ay, az
-        INTEGER :: set_acceleration
-        set_acceleration = -2
-      END FUNCTION
-
       FUNCTION get_potential(index_of_the_particle, potential)
         INCLUDE 'src/common6.h'
         INTEGER :: index_of_the_particle

src/amuse/community/octgrav/interface.cc

@@ -704,11 +704,6 @@ int set_eps2(double epsilon_squared)
   return 0;
 }
 
-int set_acceleration(int index_of_the_particle, double ax, double ay, double az)
-{
-  return -2;
-}
-
 int synchronize_model()
 {
   return 0;

src/amuse/community/petar/interface.cc

@@ -533,24 +533,6 @@ extern "C" {
         return 0;
     }
 
-    int set_acceleration(int index_of_the_particle, double ax, double ay, double az) {
-        reconstruct_particle_list();
-        int index = ptr->getParticleAdrFromID(index_of_the_particle);
-        if (index>=0) {
-            FPSoft* p = &(ptr->system_soft[index]);
-            p->acc.x = ax; 
-            p->acc.y = ay; 
-            p->acc.z = az; 
-        }    
-#ifdef PARTICLE_SIMULATOR_MPI_PARALLEL
-        int check = PS::Comm::getMaxValue(index);
-        if (check==-1) return -1;
-#else
-        else return -1;
-#endif
-        return 0;
-    }
-
     int get_potential(int index_of_the_particle, double * potential) {
         reconstruct_particle_list();
         int index = ptr->getParticleAdrFromID(index_of_the_particle);

src/amuse/community/petar/interface.h

@@ -41,8 +41,6 @@ int set_velocity(int index_of_the_particle, double vx, double vy, double vz);
 
 int get_acceleration(int index_of_the_particle, double * ax, double * ay, double * az);
 
-int set_acceleration(int index_of_the_particle, double ax, double ay, double az);
-
 int get_potential(int index_of_the_particle, double * potential);
 
 int evolve_model(double time);