The Computer Language
24.11 Benchmarks Game

k-nucleotide Java OpenJ9 #8 program

source code

/* The Computer Language Benchmarks Game
   https://salsa.debian.org/benchmarksgame-team/benchmarksgame/
   
   Naive transliteration from bearophile's program 
   contributed by Isaac Gouy
*/

import java.io.*;
import java.util.*;
import java.util.stream.*;

public final class knucleotide {

   static ArrayList<String> seqLines() throws IOException 
   {
      final var in = new BufferedReader(new InputStreamReader(System.in));
      String line;
      while ((line = in.readLine()) != null) {
         if (line.startsWith(">THREE")) break;
      }
      final var lines = new ArrayList<String>();
      while ((line = in.readLine()) != null) {
         if (line.startsWith(">")) break;
         lines.add(line);
      }    
      return lines;      
    }         

   static HashMap<String,Integer> baseCounts(int bases, String seq)
   {
      var counts = new HashMap<String,Integer>();  
      final int size = seq.length() + 1 - bases;     
      for (int i = 0; i < size; i++) {
         var nucleo = seq.substring(i,i+bases);
         Integer v;
         if ((v = counts.get(nucleo)) != null) {
            counts.put(nucleo, v+1);
         } else {
            counts.put(nucleo, 1);
         }                           
      }      
      return counts;
   }
   
   static List<String> sortedFreq(int bases, String seq)
   {
      final int size = seq.length() + 1 - bases;
      return baseCounts(bases, seq).entrySet()     
         .stream()
         .sorted(Map.Entry.comparingByValue(Comparator.reverseOrder()))
         .map(e -> String.format("%s %.3f", 
               e.getKey(), 100.0 * e.getValue() / size))
         .toList();
   }   
   
   static int specificCount(String code, String seq)
   {
      return baseCounts(code.length(),seq).getOrDefault(code,0);   
   }

   public static void main(String[] args) throws Exception 
   {
      final var seq = seqLines().stream()
         .map(s -> s.toUpperCase())
         .collect(Collectors.joining(""));
      
      for (int i : Arrays.asList(1,2)) {  
         for (String s : sortedFreq(i,seq)) {      
            System.out.println(s); 
         }
         System.out.println();          
      }
      
      for (String code : Arrays.asList("GGT", "GGTA", "GGTATT",
            "GGTATTTTAATT", "GGTATTTTAATTTATAGT")) {       
         System.out.println(specificCount(code, seq) + "\t" + code);                        
      }
   }
}
    

notes, command-line, and program output

NOTES:
64-bit Ubuntu quad core
openjdk 23 2024-09-17
IBM Semeru Runtimebr>
Eclipse OpenJ9
(build openj9-0.47.0) 
(JIT enabled, AOT enabled)


 Fri, 18 Oct 2024 22:46:37 GMT

MAKE:
mv knucleotide.openj9-8.openj9 knucleotide.java
/opt/src/ibm-semeru-23/bin/javac -d . -cp .:/opt/src/java-libs/fastutil-8.3.1.jar  knucleotide.java

3.49s to complete and log all make actions

COMMAND LINE:
 /opt/src/ibm-semeru-23/bin/java -Xshareclasses -XX:SharedCacheHardLimit=200m -Xscmx60m -Xtune:virtualized  -cp .:/opt/src/java-libs/fastutil-8.3.1.jar knucleotide 0 < knucleotide-input25000000.txt

PROGRAM OUTPUT:
A 30.295
T 30.151
C 19.800
G 19.754

AA 9.177
TA 9.132
AT 9.131
TT 9.091
CA 6.002
AC 6.001
AG 5.987
GA 5.984
CT 5.971
TC 5.971
GT 5.957
TG 5.956
CC 3.917
GC 3.911
CG 3.909
GG 3.902

1471758	GGT
446535	GGTA
47336	GGTATT
893	GGTATTTTAATT
893	GGTATTTTAATTTATAGT