Loading src/include/Commons.h +39 −2 Original line number Diff line number Diff line Loading @@ -78,54 +78,91 @@ public: ParticleDescriptor *c1; //! \brief Vector of geometric asymmetry factors. double *gaps; //! \brief Components of extinction contribution to radiation torque on a single sphere along k. double **tqse; //! \brief Components of polarized extinction contribution to radiation torque on a single sphere along k. dcomplex **tqspe; //! \brief Components of scattering contribution to radiation torque on a single sphere along k. double **tqss; //! \brief Components of polarized scattering contribution to radiation torque on a single sphere along k. dcomplex **tqsps; double ****zpv; //! \brief Mean geometric asymmetry parameters. double **gapm; //! \brief Mean polarized geometric asymmetry parameters. dcomplex **gappm; double *argi; double *args; //! \brief Geometric asymmetry parameters. double **gap; //! \brief Polarized geometric asymmetry parameters. dcomplex **gapp; //! \brief Components of extinction contribution to radiation torque on the cluster along k. double **tqce; //! \brief Components of polarized extinction contribution to radiation torque on the cluster along k. dcomplex **tqcpe; //! \brief Components of scattering contribution to radiation torque on the cluster along k. double **tqcs; //! \brief Components of polarized scattering contribution to radiation torque on the cluster along k. dcomplex **tqcps; double *duk; //! \brief Cluster extinction cross-section components. double **cextlr; //! \brief Cluster extinction cross-section components. double **cext; //! \brief Cluster Mueller Transformation Matrix components. double **cmullr; //! \brief Cluster Mueller Transformation Matrix components. double **cmul; //! \brief Geometric asymmetry parameter components. double *gapv; //! \brief Radiation extinction torque components. double *tqev; //! \brief Radiation scattering torque components. double *tqsv; //! \brief Incident unitary vector components. double *u; //! \brief Scattered unitary vector components. double *us; //! \brief Normal unitary vector components. double *un; //! \brief Normal scattered unitary vector components. double *uns; //! \brief Incident unitary vector components on polarization plane. double *up; //! \brief Scattered unitary vector components on polarization plane. double *ups; //! \brief Mean unitary vector components normal to polarization plane. double *unmp; //! \brief Mean scattered unitary vector components normal to polarization plane. double *unsmp; //! \brief Mean incident unitary vector components on polarization plane. double *upmp; //! \brief Mean scattered unitary vector components on polarization plane. double *upsmp; //! \brief Scattering angle. double scan; //! \brief Control parameter on incidence direction referred to meridional plane. double cfmp; //! \brief Control parameter on scattering direction referred to meridional plane. double sfmp; //! \brief Control parameter on incidence direction referred to scattering plane. double cfsp; //! \brief Control parameter on scattering direction referred to scattering plane. double sfsp; double qsfi; //! \brief SQSFI = XI^-2 double sqsfi; //! \brief Vectorized scattering coefficient matrix. dcomplex *am_vector; //! \brief Scattering coefficient matrix. dcomplex **am; //! \brief Argument of harmonic functions. QUESTION: correct? dcomplex arg; //! \brief Vacuum magnitude of wave vector. double vk; //! \brief Wave number. double wn; //! \brief Normalization scale. QUESTION: correct? double xip; //! \brief Number of scales (wavelengths) to be computed. int number_of_scales; Loading Loading @@ -389,7 +426,7 @@ public: dcomplex **gis; //! \brief TBD. dcomplex **gls; //! \brief TBD. //! \brief Mean scattering amplitude components. dcomplex **sam; // >>> END OF SECTION NEEDED BY CLUSTER <<< // Loading src/libnptm/Commons.cpp +0 −4 Original line number Diff line number Diff line Loading @@ -143,7 +143,6 @@ ClusterIterationData::ClusterIterationData(GeometryConfiguration *gconf, Scatter sfmp = 0.0; cfsp = 0.0; sfsp = 0.0; qsfi = 0.0; // End of suspect initializations wn = sconf->wp / 3.0e8; xip = sconf->xip; Loading Loading @@ -298,7 +297,6 @@ ClusterIterationData::ClusterIterationData(const ClusterIterationData& rhs) { sfmp = rhs.sfmp; cfsp = rhs.cfsp; sfsp = rhs.sfsp; qsfi = rhs.qsfi; // End of suspect initializations wn = rhs.wn; xip = rhs.xip; Loading Loading @@ -431,7 +429,6 @@ ClusterIterationData::ClusterIterationData(const mixMPI *mpidata, const int devi MPI_Bcast(&sfmp, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&cfsp, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&sfsp, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&qsfi, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&wn, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&xip, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&sqsfi, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); Loading Loading @@ -512,7 +509,6 @@ void ClusterIterationData::mpibcast(const mixMPI *mpidata) { MPI_Bcast(&sfmp, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&cfsp, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&sfsp, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&qsfi, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&wn, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&xip, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&sqsfi, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); Loading Loading
src/include/Commons.h +39 −2 Original line number Diff line number Diff line Loading @@ -78,54 +78,91 @@ public: ParticleDescriptor *c1; //! \brief Vector of geometric asymmetry factors. double *gaps; //! \brief Components of extinction contribution to radiation torque on a single sphere along k. double **tqse; //! \brief Components of polarized extinction contribution to radiation torque on a single sphere along k. dcomplex **tqspe; //! \brief Components of scattering contribution to radiation torque on a single sphere along k. double **tqss; //! \brief Components of polarized scattering contribution to radiation torque on a single sphere along k. dcomplex **tqsps; double ****zpv; //! \brief Mean geometric asymmetry parameters. double **gapm; //! \brief Mean polarized geometric asymmetry parameters. dcomplex **gappm; double *argi; double *args; //! \brief Geometric asymmetry parameters. double **gap; //! \brief Polarized geometric asymmetry parameters. dcomplex **gapp; //! \brief Components of extinction contribution to radiation torque on the cluster along k. double **tqce; //! \brief Components of polarized extinction contribution to radiation torque on the cluster along k. dcomplex **tqcpe; //! \brief Components of scattering contribution to radiation torque on the cluster along k. double **tqcs; //! \brief Components of polarized scattering contribution to radiation torque on the cluster along k. dcomplex **tqcps; double *duk; //! \brief Cluster extinction cross-section components. double **cextlr; //! \brief Cluster extinction cross-section components. double **cext; //! \brief Cluster Mueller Transformation Matrix components. double **cmullr; //! \brief Cluster Mueller Transformation Matrix components. double **cmul; //! \brief Geometric asymmetry parameter components. double *gapv; //! \brief Radiation extinction torque components. double *tqev; //! \brief Radiation scattering torque components. double *tqsv; //! \brief Incident unitary vector components. double *u; //! \brief Scattered unitary vector components. double *us; //! \brief Normal unitary vector components. double *un; //! \brief Normal scattered unitary vector components. double *uns; //! \brief Incident unitary vector components on polarization plane. double *up; //! \brief Scattered unitary vector components on polarization plane. double *ups; //! \brief Mean unitary vector components normal to polarization plane. double *unmp; //! \brief Mean scattered unitary vector components normal to polarization plane. double *unsmp; //! \brief Mean incident unitary vector components on polarization plane. double *upmp; //! \brief Mean scattered unitary vector components on polarization plane. double *upsmp; //! \brief Scattering angle. double scan; //! \brief Control parameter on incidence direction referred to meridional plane. double cfmp; //! \brief Control parameter on scattering direction referred to meridional plane. double sfmp; //! \brief Control parameter on incidence direction referred to scattering plane. double cfsp; //! \brief Control parameter on scattering direction referred to scattering plane. double sfsp; double qsfi; //! \brief SQSFI = XI^-2 double sqsfi; //! \brief Vectorized scattering coefficient matrix. dcomplex *am_vector; //! \brief Scattering coefficient matrix. dcomplex **am; //! \brief Argument of harmonic functions. QUESTION: correct? dcomplex arg; //! \brief Vacuum magnitude of wave vector. double vk; //! \brief Wave number. double wn; //! \brief Normalization scale. QUESTION: correct? double xip; //! \brief Number of scales (wavelengths) to be computed. int number_of_scales; Loading Loading @@ -389,7 +426,7 @@ public: dcomplex **gis; //! \brief TBD. dcomplex **gls; //! \brief TBD. //! \brief Mean scattering amplitude components. dcomplex **sam; // >>> END OF SECTION NEEDED BY CLUSTER <<< // Loading
src/libnptm/Commons.cpp +0 −4 Original line number Diff line number Diff line Loading @@ -143,7 +143,6 @@ ClusterIterationData::ClusterIterationData(GeometryConfiguration *gconf, Scatter sfmp = 0.0; cfsp = 0.0; sfsp = 0.0; qsfi = 0.0; // End of suspect initializations wn = sconf->wp / 3.0e8; xip = sconf->xip; Loading Loading @@ -298,7 +297,6 @@ ClusterIterationData::ClusterIterationData(const ClusterIterationData& rhs) { sfmp = rhs.sfmp; cfsp = rhs.cfsp; sfsp = rhs.sfsp; qsfi = rhs.qsfi; // End of suspect initializations wn = rhs.wn; xip = rhs.xip; Loading Loading @@ -431,7 +429,6 @@ ClusterIterationData::ClusterIterationData(const mixMPI *mpidata, const int devi MPI_Bcast(&sfmp, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&cfsp, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&sfsp, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&qsfi, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&wn, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&xip, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&sqsfi, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); Loading Loading @@ -512,7 +509,6 @@ void ClusterIterationData::mpibcast(const mixMPI *mpidata) { MPI_Bcast(&sfmp, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&cfsp, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&sfsp, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&qsfi, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&wn, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&xip, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); MPI_Bcast(&sqsfi, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD); Loading