 News |
FAQ
|
| Search
|
Home
|
| Getting Started | Documentation | Glish | Learn More | Programming | Contact Us |
|
| Version 1.9 Build 1556 |
|
40s to
10mins, depending on whether the compiler
(for unknown reasons) decided to re-instantiate everything or not. Sun is
aware of this bug and a recent patch has alleviated this problem somewhat.
In g++ the template instantiation is done by hand (with assistance from the aips++ system) so this problem cannot occur. On a line by line basis the g++ compiler is slower than the Sun native compiler. In practice (my experience, and in the time taken for system rebuilds) it appears to be faster, because it does not do excessive recompilations.
Because of these problems with the standard Sun compiler and to aid porting of AIPS++ to other architectures, g++ has become the standard C++ compiler for AIPS++.
This note describes details on compiling AIPS++ code that are specific to using g++ and the ``do it yourself'' template instantiation mechanism that is used to specify which templates are created.
Code compiled and optimised under g++ is in my experience about 30% slower than code compiled and optimised under the Sun native compiler. Currently we are investigating using the using the ``do it yourself'' template instantiation mechanism in conjunction with the latest release of the Sun compiler and other compilers.
If not you can switch to to the g++ compiler using:
aipsinit gnuThis assumes that you have already run aipsinit to get the normal AIPS++ environment variables set.
Ensure that the following directories exist (otherwise create them)
${MYAIPS2ROOT}/sun4sol_gnu
${MYAIPS2ROOT}/sun4sol_gnu/aux
${MYAIPS2ROOT}/sun4sol_gnu/bin
${MYAIPS2ROOT}/sun4sol_gnu/lib
${MYAIPS2ROOT}/sun4sol_gnu/bindbg
${MYAIPS2ROOT}/sun4sol_gnu/libdbg
where $MYAIPS2ROOT is the root of your workspace.
(eg. $HOME/aips++).
At Socorro where the code is on one partition and the executables are on another (which is not backed up), you should substitute for $MYAIPS2ROOT the path /tarzan2/$USERNAME, and create a symbolic link pointing from /tarzan/$USERNAME/sun4sol_gnu to /tarzan2/$USERNAME/sun4sol_gnu
Debugging should be done using gdb (use ddd or xxgdb for a graphical interface), rather than the normal system debuggers (ie. dbx and friends). Note that the commands differ from the standard ones (e.g. no 'stop' but rather 'break'). Help can be obtained by typing 'help' or 'help subject', while information is available with the 'info' command (see 'help info').
If gmake dumps core it may be because you are using a very old version that exists in /opt/local1/bin (Socorro only).
This problem is particularly endemic in the Arrays module and for these classes a particular solution has been devised. The Array class (and hence derived classes) have an upcasting member function called arrayCast() (or ac() for short).
For example in the aips/Arrays/ArrayMath.h class there is a function
T max(const Array<T> &a)which will not compile when called with a vector as an argument. ie.
Vector<Float> vec(4); cout << max(vec) << endl; // error flagged on this line!but will compile if rewritten as:
cout << max(vec.arrayCast()) << endl;
This solution gets particularly useful when doing arithmetic on Vectors because g++ will complain about:
Vector<Float> v1(4), v2(4); v1 += v2;And this should be changed to
v1.arrayCast() += v2.arrayCast(); // make the cast explicit or v1.ac() += v2.ac(); // emphasize the arithmetic not the cast
This problem has also been recently seen in the Lattice classes and a corresponding set of member functions called latticeCast(), or lc() for short, have been written.
On rare occasions you may need to selectively do an upcast when a templated argument can be either a builtin data type (Int/Float etc.) or an Array (Vector/Matrix etc.). For these occasions the functions at_c() & at_cc() exist. These global functions only upcast Arrays and leave the builtin data types alone. They are defined in the Array, Vector, Matrix & Cube classes7.1, and in aips.h for the builtin data types (and in String.h also).
The difference between the two functions is that at_c() is takes a non-constant input arguement while at_cc() takes a constant input arguement. Two function names are used rather than overloading just one because of what we think is a bug in the g++ compiler.
An example of this function in use in is the Measures module where a Quantum has standard templates of Double and Vector<Double>. Then the floor function is defined in QMath.cc as:
template <class Qtype>
Quantum<Qtype> floor(const Quantum<Qtype> &left) {
Qtype tmp = left.getValue();
Qtype ret;
at_c(ret) = floor(at_c(tmp));
return (Quantum<Qtype>(ret,left));
}
If you see this problem in other class hierarchies, you can do an explicit type cast but are instead urged to create an upcasting member function in the base class. The advantage of using member functions, rather than explicit casts, is that if g++ is ever able to do upcasts in the future it will then be easier to search for the upcasting functions and remove them. Please tell me if you do this, so I can add it to the documentation.
It is expected that this bug will disappear in g++ version 2.8.0 (private communication between Brian and one of the g++ developers).
typedef Vector<Int> gpp_vector_int;near the top of the .cc file. the type gpp_vector_int is arbitrary but this seems like a good convention (Brian's convention). This appears to occur when a specific templated types used in the context of a general template definition. eg.
template <class T> const Lattice<Bool> & PagedImage<T>::mask() constwill trigger this bug because requiring a typedef for Lattice<Bool>
I enclose this change and others that are specific to the gnu compiler in ifdefs's to avoid getting other compilers off side. eg.
#ifdef __GNUG__
typedef Vector<Int> gpp_vector_int;
#endif
It is not necessary to do this when using the array upcasting functions
(arrayCast() etc.) described in section 7.6.
The error messages that indicate you are being affected by this bug can be quite obscure and some are listed in section 7.11
The format of the template entries is described in detail in Chapter 9 of the AIPS++ system manual (the mkinst section). I will only outline the basics here.
The templates file should contain entries like:
1010 trial/Arrays/Convolver.cc
template class Convolver<Float>
1000 aips/Arrays/ArrayMath.cc aips/Arrays/Array.h
template Double min(Array<Double> const &)
The entries are:
Templates that are used by the test programs should not be put in the system file (unless there is an overriding reason), but stored in separate templates files that exist in the test directories. This keeps the aips++ libraries at a manageable size and avoids excessive link times.
If your templates file is getting long and unmanageable there is a utility called reident to sort/uniq, re-sequence and ``canonicalise''the templates file.
In the default mode7.2 reident will do the following:
#ifdef and #else and #endif are
properly sequenced.
Its output is always to the file called templates (and input and output can be the same file) so I usually just type
reidentas the default input file is also called templates and I have already appended my changes tothis file.
reident can also merge multiple template files together so that
reident templates templates.additionswill merge both input files into the output file.
The shell comment character # can be used to comment entries in
your templates file and the preprocessor directives #if,
#if defined(symbol), #if defined(symbol)!, #ifdef,
#ifndef, #else, #endif can also be used.
Each entry in the templates file gets compiled into an object file
called myClass_Ident.o where the myClass is the name of the first .cc
file in the template line and Ident is the indentification number.
Normally this object file contains only one instantiation. But often
templates instaitions are required in groups. For example to use a
CountedPtr<SkyCompRep> requires you instantiate four other
classes, namely: CountedConstPtr<SkyCompRep>,
PtrRep<SkyCompRep>, SimpleCountedConstPtr<SkyCompRep> &
SimpleCountedPtr<SkyCompRep>. In this case it is better to
group all five instantiations in one object file. This is done by
adding multiple template lines as shown below.
1000 aips/Utilities/CountedPtr.cc trial/ComponentModels/SkyCompRep.h
template class CountedConstPtr<SkyCompRep>
template class CountedPtr<SkyCompRep>
template class PtrRep<SkyCompRep>
template class SimpleCountedConstPtr<SkyCompRep>
template class SimpleCountedPtr<SkyCompRep>
gmake test.cc >&! g++.out /usr/local/gnu/bin/c++filt < g++.outor using sh
gmake test.cc 2>&1 | /usr/local/gnu/bin/c++filt | tee missing-templates
The c++filt command will translate the following mangled name
getArray__Ct11MaskedArray1Z14PointComponento
MaskedArray<PointComponent>::getArray(void)and because the getArray function is a member function of the MaskedArray class we need to instantiate the whole class (MaskedArray<PointComponent>)
When code is checked in, the templates entries should be added to the relevant system/test template instantiation file:
.../code/{aips,trial}/implement/\_ReposFiller/templates
if the templates are associated with an application or
.../code/{aips,trial}/implement/yourModule/test/templates
if the templates are associated with a test program. 7.3
Be very careful when adding your templates to the system file as any duplication of the indent. number between your added templates and the already defined ones will cause the next system rebuild to fail (a sure way to become popular). Mark Calabretta suggests:
``When you add idents to the system templates file you should always run 'reident' before checking it back in. A good strategy is to give your new entries a 0000 ident.''
A question that is often asked is ``when do I put a template in the aips package and when do I put it in the trial package?''. If any template required for an application uses code (.cc or .h files) only from the aips package it should go in the aips/implement/_Reposfiller/template file. If it uses anything from the trial package it should go in the trial/implement/_Reposfiller/template template file. This avoids the problem of having undefined symbols when some purely aips templates are defined in the trial library (at the expense of having a larger than necessary aips library).
Pure aips templates that are needed for test programs in the aips package should still be in aips/implement/yourModule/test/templates. The aips test programs should not rely on anything in trial, they are not linked against the trial libraries, and should not use any templates from the trial package.
Any templates (pure aips or involving trial) that are needed for test programs in the trial package should still be in the trial/implement/yourModule/test/templates file. This may mean that you may have to use the EXTRA_PGMRLIBS flag when compiling these test programs, although we will attempt to automate this.
You should also ensure that the templates you are adding to the system are not defined alsewhere. Duplicate templates lead to link time errors on many systems (eg. Dec Alpha). The policy for checking for duplicates is:
Brian Glendenning has a tool that can automatically check for template duplicates. It is checked into the system and information on how to use it can be obtained by typing duplicates -help
The upcasting bug has also been known to cause an exception to be thrown. This occured in the following piece of code:
template <class Qtype>
ostream& operator<< (ostream &os, const Quantum<Qtype> &ku) {
os << ku.qVal;
return os;
}
where ku.qVal could be a Double or a Vector<Double> (but
the operator<< function is defined for Arrays not Vectors). The
compiler did not detect that there was an upcasting problem and examination
with a debugger revealed it was getting into an infinite constructor loop!
Further investigation revealed that the compiler did not find a
function match using an upcast but was able to find another unintended
match, that was further resolved to the original function!